Search

Social interactions between cows play a fundamental role in the daily activities of dairy cattle. Real-time location systems provide on a continuous and automated basis information about the position of individual cows inside barns, offering a valuable opportunity to monitor dyadic social contacts. Understanding dyadic social interactions could be applied to enhance the stability of the social structure promoting animal welfare and to model disease transmission in dairy cattle. This study aimed to identify the impact of different cow characteristics on the likelihood of the formation and persistence of social contacts in dairy cattle. The individual position of the lactating cows was automatically collected once per second for 2 weeks, using an ultra-wideband system on a Swedish commercial farm consisting of nearly 200 dairy cows inside a free-stall barn. Social networks were constructed using the position data of 149 cows with available information on all characteristics during the study period. Social contacts were considered as a binary variable indicating whether a cow pair was within 2.5 m of each other for at least 10 min per day. The role of cow characteristics in social networks was studied by applying separable temporal exponential random graph models. Our results revealed that cows of the same parity interacted more consistently, as well as those born within 7 d of each other or are closely related by pedigree. The repeatability of the topological parameters indicated a consistent short-term stability of the individual animal roles within the social network structure. Additional research is required to elucidate the underlying mechanisms governing the long-term evolution of social contacts among dairy cattle and to investigate the relationship between these networks and the transmission of diseases in the dairy cattle population.

In modern freestall barns where large groups of cows are housed together, the behavior displayed by herd mates can influence the welfare and production of other individuals. Therefore, understanding social interactions in groups of dairy cows is important to enhance herd management and optimize the outcomes of both animal health and welfare in the future. Many factors can affect the number of social contacts in a group. This study aimed to identify which characteristics of a cow are associated with the number of contacts it has with other group members in 2 different functional areas (feeding and resting area) to increase our understanding of the social behavior of dairy cows. Inside 2 herds housed in freestall barns with around 200 lactating cows each, cow positions were recorded with an ultra-wideband real-time location system collecting all cows' positions every second over 2 wk. Using the positioning data of the cows, we quantified the number of contacts between them, assuming that cows spending time in proximity to one another (within a distance of 2.5 m for at least 10 min per day) were interacting socially. We documented in which barn areas these interactions occurred and used linear mixed models to investigate if lactation stage, parity, breed, pregnancy status, estrus, udder health, and claw health affect the number of contacts. We found variation in the number of contacts a cow had between individuals in both functional areas. Cows in later lactation had more contacts in the feeding area than cows in early lactation. Furthermore, in one herd, higher parity cows had fewer contacts in the feeding area than first parity cows, and in the other herd, cows in third parity or higher had more contacts in the resting area. This study indicates that cow characteristics such as parity and days in milk are associated with the number of contacts a cow has daily to its herd mates and provides useful information for further research on social interactions of dairy cows.

In modern freestall barns where large groups of cows are housed together, the behavior displayed by herd mates can influence the welfare and production of other individuals. Therefore, understanding social interactions in groups of dairy cows is important to enhance herd management and optimize the outcomes of both animal health and welfare in the future. Many factors can affect the number of social contacts in a group. This study aimed to identify which characteristics of a cow are associated with the number of contacts it has with other group members in 2 different functional areas (feeding and resting area) to increase our understanding of the social behavior of dairy cows. Inside 2 herds housed in freestall barns with around 200 lactating cows each, cow positions were recorded with an ultra-wideband real-time location system collecting all cows' positions every second over 2 wk. Using the positioning data of the cows, we quantified the number of contacts between them, assuming that cows spending time in proximity to one another (within a distance of 2.5 m for at least 10 min per day) were interacting socially. We documented in which barn areas these interactions occurred and used linear mixed models to investigate if lactation stage, parity, breed, pregnancy status, estrus, udder health, and claw health affect the number of contacts. We found variation in the number of contacts a cow had between individuals in both functional areas. Cows in later lactation had more contacts in the feeding area than cows in early lactation. Furthermore, in one herd, higher parity cows had fewer contacts in the feeding area than first parity cows, and in the other herd, cows in third parity or higher had more contacts in the resting area. This study indicates that cow characteristics such as parity and days in milk are associated with the number of contacts a cow has daily to its herd mates and provides useful information for further research on social interactions of dairy cows.

Introduction: Although breast and prostate cancers arise in different organs and are more frequent in the opposite sex, multiple studies have reported an association between their family history. Analysis of single nucleotide polymorphism data, based on distant relatives, has revealed a small positive genetic correlation between these cancers explained by common variants. The estimate of genetic correlation based on close relatives reveals the extent to which shared genetic risks are explained by both common and rare variants. This estimate is unknown for breast and prostate cancer. Method: We estimated the relative risks, heritability, and genetic correlation of breast cancer and prostate cancer, based on the Minnesota Breast and Prostate Cancer Study, a family study of 141 families ascertained for breast cancer. Results: Heritability of breast cancer was 0.34 (95% credible interval: 0.23-0.49) and 0.65 (95% credible interval: 0.36-0.97) for prostate cancer, and the genetic correlation was 0.23. In terms of odds ratios, these values correspond to a 1.3 times higher odds of breast cancer among probands, given that the brother has prostate cancer. Conclusion: This study shows the inherent relation between prostate cancer and breast cancer; an incident of one in a family increases the risk of developing the other. The large difference between estimates of genetic correlation from distant and close relatives, if replicated, suggests that rare variants contribute to the shared genetic risk of breast and prostate cancer. However, the difference could steam from genotype-by-family effects shared between the two types of cancers. ©; 2021 The Author(s).

Phenotypic variability of a genotype is relevant both in natural and domestic populations. In the past two decades, variability has been studied as a heritable quantitative genetic trait in its own right, often referred to as inherited variability or environmental canalization. So far, studies on inherited variability have only considered genetic effects of the focal individual, that is, direct genetic effects on inherited variability. Observations from aquaculture populations and some plants, however, suggest that an additional source of genetic variation in inherited variability may be generated through competition. Social interactions, such as competition, are often a source of Indirect Genetic Effects (IGE). An IGE is a heritable effect of an individual on the trait value of another individual. IGEs may substantially affect heritable variation underlying the trait, and the direction and magnitude of response to selection. To understand the contribution of IGEs to evolution of environmental canalization in natural populations, and to exploit such inherited variability in animal and plant breeding, we need statistical models to capture this effect. To our knowledge, it is unknown to what extent the current statistical models commonly used for IGE and inherited variability capture the effect of competition on inherited variability. Here, we investigate the potential of current statistical models for inherited variability and trait values, to capture the direct and indirect genetic effects of competition on variability. Our results show that a direct model of inherited variability almost entirely captures the genetic sensitivity of individuals to competition, whereas an indirect model of inherited variability captures the cooperative genetic effects of individuals on their partners. Models for trait levels, however, capture only a small part of the genetic effects of competition. The estimation of direct and indirect genetic effects of competition, therefore, is pos

The cerebellum has a highly conserved neural structure across species but varies widely in size. The wide variation in cerebellar size (both absolute and in proportion to the rest of the brain) among species and populations suggests that functional specialization is linked to its size. There is increasing recognition that the cerebellum contributes to cognitive processing and emotional control in addition to its role in motor coordination. However, to what extent cerebellum size reflects variation in these behavioral processes within species remains largely unknown. By using a unique intercross chicken population based on parental lines with high divergence in cerebellum size, we compared the behavior of individuals repeatedly exposed to the same fear test (emergence test) early in life and after sexual maturity (eight trials per age group) with proportional cerebellum size and cerebellum neural density. While proportional cerebellum size did not predict the initial fear response of the individuals (trial 1), it did increasingly predict adult individuals response as the trials progressed. Our results suggest that proportional cerebellum size does not necessarily predict an individuals fear response, but rather the habituation process to a fearful stimulus. Cerebellum neuronal density did not predict fear behavior in the individuals which suggests that these effects do not result from changes in neuronal density but due to other variables linked to proportional cerebellum size which might underlie fear habituation., Funding Agencies|Swedish FORMASSwedish Research Council Formas [2019-01508]; European Research CouncilEuropean Research Council (ERC)European Commission [FERALGEN 772874]; Swedish Research CouncilSwedish Research CouncilEuropean Commission; Carl Tryggers Stiftelse; Linkoeping University Neuro-network

Position data from real-time indoor positioning systems are increasingly used for studying individual cow behavior and social behavior in dairy herds. However, missing data challenges achieving reliable continuous activity monitoring and behavior studies. This study investigates the pattern of missing data and alternative interpolation methods in ultra-wideband based real-time indoor positioning systems in a free-stall barn. We collected 3 months of position data from a Swedish farm with around 200 cows. Data sampled for 6 days from 69 cows were used in subsequent analyzes to determine the location and duration of missing data. Data from 20 cows with the most reliable tags were selected to compare the effects of four different interpolation methods (previous, linear interpolation, cubic spline data interpolation and modified Akima interpolation). By comparing the observed data with the interpolations of the simulated missing data, the mean error distance varied from around 55 cm, using the previously last observed position, to around 17 cm for modified Akima. Modified Akima interpolation has the lowest error distance for all investigated activities (rest, walking, standing, feeding). Larger error distances were found in areas where the cows walk and turn, such as the corner between feeding and cubicles. Modified Akima interpolation is expected to be useful in the subsequent analyses of data gathered using real-time indoor positioning systems.

Phenotypic variability of a genotype is relevant both in natural and domestic populations. In the past two decades, variability has been studied as a heritable quantitative genetic trait in its own right, often referred to as inherited variability or environmental canalization. So far, studies on inherited variability have only considered genetic effects of the focal individual, that is, direct genetic effects on inherited variability. Observations from aquaculture populations and some plants, however, suggest that an additional source of genetic variation in inherited variability may be generated through competition. Social interactions, such as competition, are often a source of Indirect Genetic Effects (IGE). An IGE is a heritable effect of an individual on the trait value of another individual. IGEs may substantially affect heritable variation underlying the trait, and the direction and magnitude of response to selection. To understand the contribution of IGEs to evolution of environmental canalization in natural populations, and to exploit such inherited variability in animal and plant breeding, we need statistical models to capture this effect. To our knowledge, it is unknown to what extent the current statistical models commonly used for IGE and inherited variability capture the effect of competition on inherited variability. Here, we investigate the potential of current statistical models for inherited variability and trait values, to capture the direct and indirect genetic effects of competition on variability. Our results show that a direct model of inherited variability almost entirely captures the genetic sensitivity of individuals to competition, whereas an indirect model of inherited variability captures the cooperative genetic effects of individuals on their partners. Models for trait levels, however, capture only a small part of the genetic effects of competition. The estimation of direct and indirect genetic effects of competition, therefore, is pos

To ensure milk quality and detect cows with signs of mastitis, visual inspection of milk by prestripping quarters before milking is recommended in many countries. An objective method to find milk changed in homogeneity (i.e., with clots) is to use commercially available inline filters to inspect the milk. Due to the required manual labor, this method is not applicable in automatic milking systems (AMS). We investigated the possibility of detecting and predicting changes in milk homogeneity using data generated by AMS. In total, 21,335 quarter-level milk inspections were performed on 5,424 milkings of 624 unique cows on 4 farms by applying visual inspection of inline filters that assembled clots from the separate quarters during milking. Images of the filters with clots were scored for density, resulting in 892 observations with signs of clots for analysis (77% traces or mild cases, 15% moderate cases, and 8% heavy cases). The quarter density scores were combined into 1 score indicating the presence of clots during a single cow milking and into 2 scores summarizing the density scores in cow milkings during a 30-h sampling period. Data generated from the AMS, such as milk yield, milk flow, conductivity, and online somatic cell counts, were used as input to 4 multilayer perceptron models to detect or predict single milkings with clots and to detect milking periods with clots. All models resulted in high specificity (98-100%), showing that the models correctly classified cow milkings or cow milking periods with no clots observed. The ability to successfully classify cow milkings or cow periods with observed clots had a low sensitivity. The highest sensitivity (26%) was obtained by the model that detected clots in a single milking. The prevalence of clots in the data was low (2.4%), which was reflected in the results. The positive predictive value depends on the prevalence and was relatively high, with the highest positive predictive value (72%) reached in the model that

Milk quality and clinical mastitis in dairy cows are monitored by detecting visually abnormal milk. A standardized method to evaluate clots in milk and studies of the incidence and dynamics of clots in milk at the quarter level are lacking. We validated a method to score clot density in quarter milk samples and describe the prevalence and dynamics of the density scores between consecutive samplings and periods in 4 farms with automatic milking systems. Using in-line filters, we collected quarter milk samples at each milking during 3 periods of 30 h each in each farm. Clot density was scored based on coverage of the filter area as 0 (negative), 1 (trace), 2 (mild), 3 (moderate), 4 (heavy), and 5 (very heavy). The score for a specific quarter and milking is referred to as the quarter milking score (QMS). Three assessors independently scored 902 images of filter samples with a Fleiss kappa value of 0.72. In total, 21,202 quarter milk samples from 5,398 milkings of 621 cows were collected. Of the quarter filter samples, 2.4% had visible clots, distributed as mild (1.4%), moderate (0.6%), heavy (0.3%), and very heavy (<0.1%, n = 8). Cases with a cow period sum of QMS ≥ 4, corresponding to 9.4% of all periods, harbored 86% and 94% of all QMS of 2 to 5 and 3 to 5, respectively. Of these cases, cows sampled in all 3 periods and clots in only 1 period had a quarter period sum score ≥ 1 in 1.8 different quarters in average. Corresponding numbers for the cows with clots or traces in 2 or 3 periods were 2.2 and 2.5 different quarters, respectively. A QMS of 2 to 5 in the preceding milking increased the chance of a QMS >1 in the following milking, with an average chance of 38%. The probability of a QMS > 1 increased with increasing previous QMS, a higher sum of QMS during the milking period, longer milking interval, and higher lactation number, but decreased with increasing days in milk. Our study showed that the method of clot-density scoring is feasible to perform and

The increasing amount of available biological information on the markers can be used to inform the models applied for genomic selection to improve predictions. The objective of this study was to propose a general model for genomic selection using a link function approach within the hierarchical generalized linear model framework (hglm) that can include external information on the markers. These models can be fitted using the well-established hglm package in R. We also present an R package (CodataGS) to fit these models, which is significantly faster than the hglm package. Simulated data was used to validate the proposed model. We tested categorical, continuous and combination models where the external information on the markers was related to 1) the location of the QTLs on the genome with varying degree of uncertainty, 2) the relationship of the markers with the QTLs calculated as the LD between them, and 3) a combination of both. The proposed models showed improved accuracies from 3.8% up to 23.2% compared to the SNP-BLUP method in a simulated population derived from a base population with 100 individuals. Moreover, the proposed categorical model was tested on a dairy cattle dataset for two traits (Milk Yield and Fat Percentage). These results also showed improved accuracy compared to SNP-BLUP, especially for the Fat% trait. The performance of the proposed models depended on the genetic architecture of the trait, as traits that deviate from the infinitesimal model benefited more from the external information. Also, the gain in accuracy depended on the degree of uncertainty of the external information provided to the model. The usefulness of these type of models is expected to increase with time as more accurate information on the markers becomes available.

eSSENCE

The depredation of semi-domesticated reindeer by large carnivores reflects an important human-wildlife conflict in Fennoscandia. Recent studies have revealed that brown bears (Ursus arctos) may kill substantial numbers of reindeer calves (Rangifer tarandus tarandus) in forest areas in Sweden. Several authors have suggested that predation risk is an important driver of habitat selection in wild Rangifer populations where predation is a limiting factor, but little is known about these mechanisms in semi-domesticated populations. We examined the habitat selection of female reindeer in relation to spatial and temporal variations in brown bear predation risk on the reindeer calving grounds and evaluated the simultaneous responses of brown bears and reindeer to landscape characteristics. We used GPS data from 110 reindeer years (97 individuals) and 29 brown bear years (19 individuals), from two reindeer herding districts in the forest area of northern Sweden. Our results did not indicate that reindeer alter their behavior in response to spatiotemporal variation in brown bear predation risk, on the scale of the calving range. Instead, we suggest that spatiotemporal behavioral adjustments by brown bears were the main driver of prey-predator interactions in our study system. Contrasting responses by brown bears and reindeer to clear-cuts and young forest indicate that forestry can influence species interactions and possibly yield negative consequences for the reindeer herd. Even if clear-cuts may be beneficial in terms of calf survival, logging activity will eventually cause greater abundance of young regenerating forest, reducing available reindeer habitats and increasing habitat preferred by brown bears. Domestication may have made semi-domesticated reindeer in Fennoscandia less adapted to cope with predators. Areal restrictions, limiting the opportunity for dispersion and escape, possibly make the calves more susceptible to predation. Also, a generally higher population d

Over the last decade, we have seen a massive increase in the construction of wind farms in northern Fennoscandia. Wind farms comprising hundreds of wind turbines are being built, with little knowledge of the possible cumulative adverse effects on the habitat use and migration of semi-domesticated free-ranging reindeer. We assessed how reindeer responded to wind farm construction in an already fragmented landscape, with specific reference to the effects on use of movement corridors and reindeer habitat selection. We used GPS-data from reindeer during calving and post-calving in the MalAyen reindeer herding community in Sweden. We analysed data from the pre-development years compared to the construction years of two relatively small wind farms. During construction of the wind farms, use of original migration routes and movement corridors within 2 km of development declined by 76 %. This decline in use corresponded to an increase in activity of the reindeer measured by increased step lengths within 0-5 km. The step length was highest nearest the development and declining with distance, as animals moved towards migration corridors and turned around or were observed in holding patterns while not crossing. During construction, reindeer avoided the wind farms at both regional and landscape scale of selection. The combined construction activities associated with even a few wind turbines combined with power lines and roads in or close to central movement corridors caused a reduction in the use of such corridors and grazing habitat and increased the fragmentation of the reindeer calving ranges.

Generalized linear mixed models are flexible tools for modeling non-normal data and are useful for accommodating overdispersion in Poisson regression models with random effects. Their main difficulty resides in the parameter estimation because there is no analytic solution for the maximization of the marginal likelihood. Many methods have been proposed for this purpose and many of them are implemented in software packages. The purpose of this study is to compare the performance of three different statistical principles - marginal likelihood, extended likelihood, Bayesian analysis-via simulation studies. Real data on contact wrestling are used for illustration.

We present a new version (> 2.0) of the hglm package for fitting hierarchical generalized linear models (HGLMs) with spatially correlated random effects. CAR() and SAR() families for conditional and simultaneous autoregressive random effects were implemented. Eigen decomposition of the matrix describing the spatial structure (e.g., the neighborhood matrix) was used to transform the CAR/SAR random effects into an independent, but eteroscedastic, Gaussian random effect. A linear predictor is fitted for the random effect variance to estimate the parameters in the CAR and SAR models. This gives a computationally efficient algorithm for moderately sized problems.

Generalized linear mixed models are flexible tools for modeling non-normal data and are usefulfor accommodating overdispersion in Poisson regression models with random effects. Theirmain difficulty resides in the parameter estimation because there is no analytic solution for themaximization of the marginal likelihood. Many methods have been proposed for this purpose andmany of them are implemented in software packages. The purpose of this study is to comparethe performance of three different statistical principles –marginal likelihood, extended likelihood,Bayesian analysis – via simulation studies. Real data on contact wrestling are used for illustration., Peer Reviewed