Hyaluronic acid-based hydrogels with tunable mechanics improved structural and contractile properties of cells

Extracellular matrix (ECM) regulates cellular responses through mechanotransduction. The standard approach of in vitro culturing on plastic surfaces overlooks this phenomenon, so there is a need for biocompatible materials that exhibit adjustable mechanical and structural properties, promote cell adhesion and proliferation at low cost and for use in 2D or 3D cell cultures. This study presents a new tunable hydrogel system prepared from high-molecular hyaluronic acid (HA), Bovine serum albumin (BSA), and gelatin cross-linked using EDC/NHS. Hydrogels with Young's moduli (E) ranging from subunit to units of kilopascals were prepared by gradually increasing HA and BSA concentrations. Concentrated high-molecular HA network led to stiffer hydrogel with lower cluster size and swelling capacity. Medium and oxygen diffusion capability of all hydrogels showed they are suitable for 3D cell cultures. Mechanical and structural changes of mouse embryonic fibroblasts (MEFs) on hydrogels were compared with cells on standard cultivation surfaces. Experiments showed that hydrogels have suitable mechanical and cell adhesion capabilities, resulting in structural changes of actin filaments. Lastly, applying hydrogel for a more complex HL-1 cell line revealed improved mechanical and electrophysiological contractile properties.

Puncturing lipid membranes: onset of pore formation and the role of hydrogen bonding in the presence of flavonoids

Products of lipid peroxidation induce detrimental structural changes in cell membranes, such as the formation of water pores, which occur in the presence of lipids with partially oxidized chains. However, the influence of another class of products, dicarboxylic acids, is still unclear. These products have greater mobility in the lipid bilayer, which enables their aggregation and the formation of favorable sites for the appearance of pores. Therefore, dodecanedioic acid (DDA) was selected as a model product. Additionally, the influence of several structurally different flavonoids on DDA aggregation via formation of hydrogen bonds with carboxyl groups was investigated. The molecular dynamics of DDA in DOPC lipid bilayer revealed the formation of aggregates extending over the hydrophobic region of the bilayer and increasing its polarity. Consequently, water penetration and the appearance of water wires was observed, representing a new step in the mechanism of pore formation. Furthermore, DDA molecules were found to interact with lipid polar groups, causing them to be buried in the bilayer. The addition of flavonoids to the system disrupted aggregate formation, resulting in the displacement of DDA molecules from the center of the bilayer. The placement of DDA and flavonoids in the lipid bilayer was confirmed by small-angle X-ray scattering. Atomic force microscopy and electron paramagnetic resonance were used to characterize the structural properties. The presence of DDA increased bilayer roughness and decreased the ordering of lipid chains, confirming its detrimental effects on the membrane surface, while flavonoids were found to reduce or reverse these changes.

New insight into the biocompatibility/toxicity of graphene oxides and their reduced forms on Chlamydomonas reinhardtii

Graphene oxides (GOs) and their reduced forms are often discussed both positively and negatively due to the lack of information about their chemistry and structure. This study utilized GOs with two sheet sizes that were further reduced by two reducing agents (sodium borohydride and hydrazine) to obtain two different degrees of reduction. The synthesized nanomaterials were characterized using scanning electron microscopy (SEM), atomic force microscopy (AFM), X-ray photoelectron spectroscopy (XPS), elemental analysis (EA), Fourier transform infrared (FTIR) spectroscopy, and Raman spectroscopy (RA) to understand their chemistry and structure. The second focus of our investigation included in vitro testing of the biocompatibility/toxicity of these materials on a model organism, the freshwater microalga Chlamydomonas reinhardtii. The effects were studied on the basis of biological endpoints complemented by biomass investigation (FTIR spectroscopy, EA, and atomic absorption spectrometry (AAS)). The results showed that the biocompatibility/toxicity of GOs is dependent on their chemistry and structure and that it is impossible to generalize the toxicity of graphene-based nanomaterials.

Effect of titanium nanostructured surface on fibroblast behavior

As the consumption of implants increases, so do the requirements for individual types of implants, for example, improved biocompatibility or longevity. Therefore, the nano-modification of the titanium surface is often chosen. The aim was to characterize the modified surface with a focus on medical applications. The titanium surface was modified by the anodic oxidation method to form nanotubes. Subsequently, the material was characterized and analyzed for medical applications-surface morphology, surface wettability, chemical composition, and release of ions into biological fluids. A human gingival fibroblasts (HGFb) cell line was used in the viability study. A homogeneous layer of nanotubes of defined dimensions was formed on the titanium surface, ensuring the material's biocompatibility-the preparation conditions influence the resulting properties of the nanostructured surface. Nanostructured titanium exhibited more suitable characteristics (e.g., wettability, roughness, ion release) for biological applications than compared to pure titanium. It was possible to understand the behavior of the modified layer on the titanium surface and its effect on cell behavior. Another contribution of this work is the combination of material characterization (ion release) with the study of cytocompatibility (direct contact of cells with metals).

Atomic force microscopy and surface plasmon resonance for real-time single-cell monitoring of bacteriophage-mediated lysis of bacteria

The growing incidence of multidrug-resistant bacterial strains presents a major challenge in modern medicine. Antibiotic resistance is often exhibited by Staphylococcus aureus, which causes severe infections in human and animal hosts and leads to significant economic losses. Antimicrobial agents with enzymatic activity (enzybiotics) and phage therapy represent promising and effective alternatives to classic antibiotics. However, new tools are needed to study phage-bacteria interactions and bacterial lysis with high resolution and in real-time. Here, we introduce a method for studying the lysis of S. aureus at the single-cell level in real-time using atomic force microscopy (AFM) in liquid. We demonstrate the ability of the method to monitor the effect of the enzyme lysostaphin on S. aureus and the lytic action of the Podoviridae phage P68. AFM allowed the topographic and biomechanical properties of individual bacterial cells to be monitored at high resolution over the course of their lysis, under near-physiological conditions. Changes in the stiffness of S. aureus cells during lysis were studied by analyzing force-distance curves to determine Young's modulus. This allowed observing a progressive decline in cellular stiffness corresponding to the disintegration of the cell envelope. The AFM experiments were complemented by surface plasmon resonance (SPR) experiments that provided information on the kinetics of phage-bacterium binding and the subsequent lytic processes. This approach forms the foundation of an innovative framework for studying the lysis of individual bacteria that may facilitate the further development of phage therapy.

Myeloperoxidase mediated alteration of endothelial function is dependent on its cationic charge

Endothelial cell (EC) glycocalyx (GLX) comprise a multicomponent layer of proteoglycans and glycoproteins. Alteration of its integrity contributes to chronic vascular inflammation and leads to the development of cardiovascular diseases. Myeloperoxidase (MPO), a highly abundant enzyme released by polymorphonuclear neutrophils, binds to the GLX and deleteriously affects vascular EC functions. The focus of this study was to elucidate the mechanisms of MPO-mediated alteration of GLX molecules, and to unravel subsequent changes in endothelial integrity and function. MPO binding to GLX of human ECs and subsequent internalization was mediated by cell surface heparan sulfate chains. Moreover, interaction of MPO, which is carrying a cationic charge, with anionic glycosaminoglycans (GAGs) resulted in reduction of their relative charge. By means of micro-viscometry and atomic force microscopy, we disclosed that MPO can crosslink GAG chains. MPO-dependent modulation of GLX structure was further supported by alteration of wheat germ agglutinin staining. Increased expression of ICAM-1 documented endothelial cell activation by both catalytically active and also inactive MPO. Furthermore, MPO increased vascular permeability connected with reorganization of intracellular junctions, however, this was dependent on MPO's catalytic activity. Novel proteins interacting with MPO during transcytosis were identified by proteomic analysis. Altogether, these findings provide evidence that MPO through interaction with GAGs modulates overall charge of the GLX, causing modification of its structure and thus affecting EC function. Importantly, our results also suggest a number of proteins interacting with MPO that possess a variety of cellular localizations and functions.

Cell type specific adhesion to surfaces functionalised by amine plasma polymers

Our previously-obtained impressive results of highly increased C2C12 mouse myoblast adhesion to amine plasma polymers (PPs) motivated current detailed studies of cell resistance to trypsinization, cell proliferation, motility, and the rate of attachment carried out for fibroblasts (LF), keratinocytes (HaCaT), rat vascular smooth muscle cells (VSMC), and endothelial cells (HUVEC, HSVEC, and CPAE) on three different amine PPs. We demonstrated the striking difference in the resistance to trypsin treatment between endothelial and non-endothelial cells. The increased resistance observed for the non-endothelial cell types was accompanied by an increased rate of cellular attachment, even though spontaneous migration was comparable to the control, i.e., to the standard cultivation surface. As demonstrated on LF fibroblasts, the resistance to trypsin was similar in serum-supplemented and serum-free media, i.e., medium without cell adhesion-mediating proteins. The increased cell adhesion was also confirmed for LF cells by an independent technique, single-cell force spectroscopy. This method, as well as the cell attachment rate, proved the difference among the plasma polymers with different amounts of amine groups, but other investigated techniques could not reveal the differences in the cell behaviour on different amine PPs. Based on all the results, the increased resistance to trypsinization of C2C12, LF, HaCaT, and VSMC cells on amine PPs can be explained most probably by a non-specific cell adhesion such as electrostatic interaction between the cells and amine groups on the material surface, rather than by the receptor-mediated adhesion through serum-derived proteins adsorbed on the PPs.

The Structural Integrity of the Model Lipid Membrane during Induced Lipid Peroxidation: The Role of Flavonols in the Inhibition of Lipid Peroxidation

The structural integrity, elasticity, and fluidity of lipid membranes are critical for cellular activities such as communication between cells, exocytosis, and endocytosis. Unsaturated lipids, the main components of biological membranes, are particularly susceptible to the oxidative attack of reactive oxygen species. The peroxidation of unsaturated lipids, in our case 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), induces the structural reorganization of the membrane. We have employed a multi-technique approach to analyze typical properties of lipid bilayers, i.e., roughness, thickness, elasticity, and fluidity. We compared the alteration of the membrane properties upon initiated lipid peroxidation and examined the ability of flavonols, namely quercetin (QUE), myricetin (MCE), and myricitrin (MCI) at different molar fractions, to inhibit this change. Using Mass Spectrometry (MS) and Fourier Transform Infrared Spectroscopy (FTIR), we identified various carbonyl products and examined the extent of the reaction. From Atomic Force Microscopy (AFM), Force Spectroscopy (FS), Small Angle X-Ray Scattering (SAXS), and Electron Paramagnetic Resonance (EPR) experiments, we concluded that the membranes with inserted flavonols exhibit resistance against the structural changes induced by the oxidative attack, which is a finding with multiple biological implications. Our approach reveals the interplay between the flavonol molecular structure and the crucial membrane properties under oxidative attack and provides insight into the pathophysiology of cellular oxidative injury.

Application of the Enzymatic Electrochemical Biosensors for Monitoring Non-Competitive Inhibition of Enzyme Activity by Heavy Metals

The inhibition effect of the selected heavy metals (Ag⁺, Cd²⁺, Cu²⁺, and Hg²⁺) on glucose oxidase (GOx) enzyme from Aspergillus niger (EC 1.1.3.4.) was studied using a new amperometric biosensor with an electrochemical transducer based on a glassy carbon electrode (GCE) covered with a thin layer of multi-wall carbon nanotubes (MWCNTs) incorporated with ruthenium(IV) oxide as a redox mediator. Direct adsorption of multi-wall carbon nanotubes (MWCNTs) and subsequent covering with Nafion^® layer was used for immobilization of Gox. The analytical figures of merit of the developed glucose (Glc) biosensor are sufficient for determination of Glc in body fluids in clinical analysis. From all tested heavy metals, mercury(II) has the highest inhibition effect. However, it is necessary to remember that cadmium and silver ions also significantly inhibit the catalytic activity of Gox. Therefore, the development of Gox biosensors for selective indirect determination of each heavy metal still represents a challenge in the field of bioelectroanalysis. It can be concluded that amperometric biosensors, differing in the utilized enzyme, could find their application in the toxicity studies of various poisons.

Post-Translational Modifications and Diastolic Calcium Leak Associated to the Novel RyR2-D3638A Mutation Lead to CPVT in Patient-Specific hiPSC-Derived Cardiomyocytes

Background: Sarcoplasmic reticulum Ca²⁺ leak and post-translational modifications under stress have been implicated in catecholaminergic polymorphic ventricular tachycardia (CPVT), a highly lethal inherited arrhythmogenic disorder. Human induced pluripotent stem cells (hiPSCs) offer a unique opportunity for disease modeling. Objective: The aims were to obtain functional hiPSC-derived cardiomyocytes from a CPVT patient harboring a novel ryanodine receptor (RyR2) mutation and model the syndrome, drug responses and investigate the molecular mechanisms associated to the CPVT syndrome. Methods: Patient-specific cardiomyocytes were generated from a young athletic female diagnosed with CPVT. The contractile, intracellular Ca²⁺ handling and electrophysiological properties as well as the RyR2 macromolecular remodeling were studied. Results: Exercise stress electrocardiography revealed polymorphic ventricular tachycardia when treated with metoprolol and marked improvement with flecainide alone. We found abnormal stress-induced contractile and electrophysiological properties associated with sarcoplasmic reticulum Ca²⁺ leak in CPVT hiPSC-derived cardiomyocytes. We found inadequate response to metoprolol and a potent response of flecainide. Stabilizing RyR2 with a Rycal compound prevents those abnormalities specifically in CPVT hiPSC-derived cardiomyocytes. The RyR2-D3638A mutation is located in the conformational change inducing-central core domain and leads to RyR2 macromolecular remodeling including depletion of PP2A and Calstabin2. Conclusion: We identified a novel RyR2-D3638A mutation causing 3D conformational defects and aberrant biophysical properties associated to RyR2 macromolecular complex post-translational remodeling. The molecular remodeling is for the first time revealed using patient-specific hiPSC-derived cardiomyocytes which may explain the CPVT proband’s resistance. Our study promotes hiPSC-derived cardiomyocytes as a suitable model for disease modeling, testing new therapeutic compounds, personalized medicine and deciphering underlying molecular mechanisms.

Short communication: Reliability of single-step genomic BLUP breeding values by multi-trait test-day model analysis

The purpose of our study was to develop an approximation procedure to estimate reliabilities of single-step genomic BLUP breeding values in a test-day model for routine evaluation of milk yield in a dairy cattle population. Input data consisted of 20,220,047 first-, second-, and third-lactation test-day milk yield records of 1,126,102 Czech Holstein cows (each lactation being considered a separate trait), with 1,844,679 animals in the pedigree file and with genomic data from 2,236 bulls. Evaluation was according to a multi-lactation model. The procedure was based on the effective number of records per animal from milk recording as well as from genomic and pedigree relationships. Traits were analyzed individually, and genetic covariances among traits were subsequently taken into account. The use of genomic information increased average reliability in young bulls from 0.276 to 0.505, but increased reliability in proven bulls only from 0.828 to 0.855. The reliabilities of genomic breeding values in multi-trait evaluation for first, second and third lactations, respectively, averaged 0.652, 0.673, and 0.633 for young bulls and 0.907, 0.894, and 0.852 for proven bulls. For an index combining all 3 lactations, the average reliability of a single-step genomic BLUP prediction was 0.712 and 0.925 for younger and proven bulls, respectively. Increased reliability due to genotyping in the population of all genotyped and nongenotyped animals was very small (<0.01) because of the small proportion of genotyped animals in the population.

Contribution of domestic production records, Interbull estimated breeding values, and single nucleotide polymorphism genetic markers to the single-step genomic evaluation of milk production

Estimated breeding values (EBV) for first-lactation milk production of Holstein cattle in the Czech Republic were calculated using a conventional animal model and by single-step prediction of the genomic enhanced breeding value. Two overlapping data sets of milk production data were evaluated: (1) calving years 1991 to 2006, with 861,429 lactations and 1,918,901 animals in the pedigree and (2) calving years 1991 to 2010, with 1,097,319 lactations and 1,906,576 animals in the pedigree. Global Interbull (Uppsala, Sweden) deregressed proofs of 114,189 bulls were used in the analyses. Reliabilities of Interbull values were equivalent to an average of 8.53 effective records, which were used in a weighted analysis. A total of 1,341 bulls were genotyped using the Illumina BovineSNP50 BeadChip V2 (Illumina Inc., San Diego, CA). Among the genotyped bulls were 332 young bulls with no daughters in the first data set but more than 50 daughters (88.41, on average) with performance records in the second data set. For young bulls, correlations of EBV and genomic enhanced breeding value before and after progeny testing, corresponding average expected reliabilities, and effective daughter contributions (EDC) were calculated. The reliability of prediction pedigree EBV of young bulls was 0.41, corresponding to EDC=10.6. Including Interbull deregressed proofs improved the reliability of prediction by EDC=13.4 and including genotyping improved prediction reliability by EDC=6.2. Total average expected reliability of prediction reached 0.67, corresponding to EDC=30.2. The combination of domestic and Interbull sources for both genotyped and nongenotyped animals is valuable for improving the accuracy of genetic prediction in small populations of dairy cattle.

Selection of a suitable data set and model for the genetic evaluation of the linear description of conformation and type description in Old Kladruber horses

Abstract. The objective of this paper was to select a suitable data subset and statistical model for the estimation of genetic parameters for 36 traits of the linear type in 977 Old Kladruber horses. Two subsets were tested to identify a suitable subset for analysis. One subset included repeated evaluation of certain individuals, whereas the other did not. The most suitable subset included repeated evaluation (n=1 390). The selection of a suitable model was made from 4 candidate models. These models comprised a number of random effects (direct individual effect and animal permanent environmental effect of the animal) and a number of fixed effects (colour variant, stud, colour variant × stud interaction, sex, age at description, year of birth, year of description). The model was selected based on the Akaike information criterion (AIC, Akaike 1974), residual variance and heritability coefficient. The model that included colour variant, stud, colour variant × stud interaction, sex, age at description, and year of description as fixed effects and direct individual and animal permanent environment as random effects was the most suitable model for the estimation of genetic parameters and for the subsequent estimation of breeding values.

Genetic parameters for growth of young beef bulls

Abstract. The average daily gains of young bulls on test stations (ADGT) were analysed for the most frequent breeds of beef cattle in the Czech Republic using a multiple-trait animal model. Body weights at birth (W0), at 120 days of age (W120) and at weaning at 210 days (WW) were considered in this model as pre-weaning growth. The tested models comprised some of the random effects: direct genetic effect, maternal genetic effect, permanent animal environment effect, permanent maternal environment effect, and some of the fixed effects: dam’s age, sex, herd-year-season, linear and quadratic regression on age at the beginning of the test. For optimization of the models Akaike information criterion (AIC), residual variance and likelihood ratio test were used. Coefficients of direct and maternal heritability across breeds of about 0.25 for W0, about 0.17 for W120, about 0.17 for WW and about 0.29 for ADGT were estimated by all models. All criteria selected models including the permanent animal environment effect, which was the most important effect in the model.

Population structure of Czech cold-blooded breeds of horses

Abstract. In order to estimate effective population size, generation interval and the development of inbreeding coefficients (Fx) in three original breeds of cold-blooded horses kept in the Czech Republic: Silesian Noriker (SN), Noriker (N) and Czech-Moravian Belgian horse (CMB) all animals of the particular breeds born from 1990 to 2007 were analysed. The average values of generation interval between parents and their offspring were: 8.53 in SN, 8.88 in N and 8.56 in CMB. Average values of effective population size were estimated to be: 86.3 in SN, 162.3 in N and 104.4 in CMB. The average values of inbreeding coefficient were 3.13 % in SN stallions and 3.39 % in SN mares, in the N breed 1.76 % and 1.26 % and in the CMB breed 3.84 % and 3.26 % respectively. Overall averages of Fx were: 3.23 %, 1.51 % and 3.55 % for the breeds SN, N and CMB. The average value of inbreeding coefficient Fx increased by 1.22 % in SN, by 0.35 % in N and by 1.01 % in CMB, respectively. This may lead to a reduction in genetic variability. Reduction in genetic variability could be either controlled in cooperation with corresponding populations of cold-blooded breeds in other European countries or controlled by number of sires used in population

Estimation of animal × environment interaction in Czech beef cattle

Abstract. The objective was the estimation of animal × environment interaction for weaning weight (WWT) of the most frequent breeds of beef cattle in the Czech Republic. The environment was represented by 3 areas: Mountain (1), Foothills (2) and Lowland (3). A pairwise analysis was carried out between the areas: 1 vs 2 (n=5 149) and 2 vs 3 (n=4 971) for the estimation of (co)variance components and genetic correlations between the environments. A multiple-trait animal model was used where weaning weights in different areas were considered as different traits. The direct and maternal heritability estimates of WWT were for area (1): 0.12 (direct) and 0.05 (maternal); for (2): 0.26 and 0.11 (direct), 0.10 and 0.11 (maternal) and for (3): 0.26 (direct) 0.03 (maternal). Direct and maternal genetic correlations were between: 1 and 2 0.98 and 0.65 respectively, with standard deviations 0.03 and 0.37 respectively, and between 2 and 3: 0.82 and 0.50 respectively, with standard deviations 0.12 and 0.31 respectivly. The magnitude of the across-areas genetic correlation indicates that the animal by environment interaction was not biologically important and can be ignored in the evaluation of beef cattle in the Czech Republic.

Comparison of Breeding Values for Daily Gains of Bulls Estimated with Multi-Trait and Random Regression Models

Abstract. In this study, random regression models with Legendre polynomials of the 2nd, 3rd and 4th degree (RR2, RR3 and RR4) are compared with regard to the estimation of breeding values for the average daily gain of Czech Pied bulls (Simmental type). The data were prepared such that a multi-trait model (MTM) could be used as reference model. For each bull, 8 repeated records or fewer were available for the testing period from the 12th to the 420th day of life. For the modeling of the expected value structure, the fixed regression coefficients of the Legendre polynomials were subordinated hierarchically to the herd-year-season effects (HYS). For the comparison of the random regression models with the reference model, rank correlations between the estimated breeding values of various animal groups were calculated and a variety of top-lists were analyzed. In general, models RR3 and RR4 returned higher rank correlations with MTM in comparison to model RR2. Additionally, the number of common animals in the 1% and 10% top-lists showed that models RR3 and RR4 are to be preferred over RR2 when it comes to the estimation of breeding values.

Anpassung eines Fixed Regression Modells für die tägliche Zunahme von Fleckviehbullen mit Hilfe von Informationskriterien

Abstract. Title of the paper: Fitting a fixed regression model for daily gain of bulls using information criterion In this study the model choice is demonstrated exemplarily on data of 6405 Czech Simmental bulls using information criterion. Per bull up to 8 observations were available for the trait daily gain. Because the animals showed different age on control day, the expected gain curves were described in the population and within the herd*year*season-classes by second, third or fourth order Legendre polynomials of age. For optimization of the fixed effects and to choice the covariance structure of the repeated records the information criteria of Akaike (AIC), the Bayesian criteria (BIC) and the ICOMP-criteria, developed mainly from Bozdogan, were used. Within and over all covariance structures AIC selected generally the most complex model. On the other hand, BIC and ICOMP favoured a model with second order polynomials of age nested within the head*year*seasonclasses. All criterion selected models with nested second order polynomials within the herd*year*season-classes in comparison to models with non-nested polynomials of age.

Estimation of Genetic Parameters for Daily Gains of Bulls with Multi-Trait and Random Regression Models

Abstract. The average daily gains of 6,420 Czech Pied bulls (dual-purpose, Simmental type) from 7 breeding stations were analyzed using single-trait animal models, a multi-trait animal model and random regression models. The effects of station, year and season were taken into account by creating herd-year-season classes (HYS) with the season being defined as a 3-month class starting with December. Legendre polynomials of the 1st to the 4th degree were used to describe the daily gains within the HYS classes as well as to model bull-specific gain curves. The comparison of the h2-values estimated with single-trait models and those gained with a multi-trait model returned only insignificant differences. The comparison of genetic parameters based on the multi-trait model to those from different random regression models shows that polynomials of at least the 2nd degree are to be used for the genetic analysis of daily gains.

Selection of a suitable data set and model for the estimation of genetic parameters of the weaning weight in beef cattle

Abstract. The objective of this paper was to select a suitable data subset and statistical model for the estimation of genetic parameters for weaning weight of beef cattle in the Czech Republic. Nine subsets were tested for the selection of a suitable subset. The subsets differed from each other in the limit of sampling criteria. The most suitable subset satisfied these conditions: at least 5 individuals per each sire, 5 individuals per HYS (herd, year, season), 2 sires per HYS, and individuals per dams that have at least one half-sister and two offspring (n = 4 806). The selection of a suitable model was carried out from 10 models. These models comprised some of the random effects: direct genetic effect, maternal genetic effect, permanent maternal environment effect, HYS, sire × herd or sire × year interaction, and some of the fixed effects: dam’s age, sex (young bull, heifer × single, twin born), HYS, year, herd. The direct heritability (h2a) ranged from 0.06 to 0.17, of maternal heritability (h2m) from 0.03 to 0.06. The genetic correlations between the direct and maternal effect (ram) were in the range of –0.15 –0.42.