Effect of genotyping density on the detection of runs of homozygosity and heterozygosity in cattle.

Runs of homozygosity (ROHom) are contiguous stretches of homozygous regions of the genome. In contrast, runs of heterozygosity (ROHet) are heterozygosity-rich regions. The detection of these two types of genomic regions (ROHom and ROHet) is influenced by the parameters involved in their identification and the number of available single nucleotide polymorphisms (SNP). The present study aimed to test the effect of chip density in detecting ROHom and ROHet in the Italian Simmental cattle breed. A sample of 897 animals were genotyped at low-density (50k SNP; 397 individuals), medium-density (140k SNP; 348 individuals), or high-density (800k SNP; 152 individuals). The number of ROHom and ROHet per animal (nROHom and nROHet, respectively) and their average length were calculated. ROHom or ROHet shared by more than one animal and the number of times a particular SNP was inside a run were also computed (SNPROHom and SNPROHet). As the chip density increased, the nROHom increased, whereas their average length decreased. In contrast, the nROHet decreased and the average length increased as the chip density increased. The most repeated ROHom harbored no genes, whereas in the most repeated ROHet four genes (SNRPN, SNURF, UBE3A and ATP10A) previously associated with reproductive traits were found. Across the three datasets, 31 SNP, located on Bos taurus autosome (BTA) 6, and 37 SNP (located on BTA21) exceeded the 99th percentile in the distribution of the SNPROHom and SNPROHet, respectively. The genomic region on BTA6 mapped the SLIT2, PACRGL, and KCNIP4 genes, whereas 19 and 18 genes were mapped on BTA16 and BTA21, respectively. Interestingly, most of genes found through the ROHet analysis were previously reported to be related to health, reproduction, and fitness traits. The results of the present study confirm that the detection of ROHom is more reliable when the chip density increases, whereas the ROHet trend seems to be the opposite. Genes and Quantitative Trait Loci (QTL) mapped in the highlighted regions confirm that ROHet can be due to balancing selection, thus related to fitness traits, health, and reproduction, whereas ROHom are mainly involved in production traits. The results of the present study strengthened the usefulness of these parameters in analyzing the genomes of livestock and their biological meaning.

Characterization of heterozygosity-rich regions in Italian and worldwide goat breeds.

Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north-south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

Investigation of phenotypic, genetic and genomic background of Milk spectra in Sarda dairy sheep.

Aim of this study was to analyse the genetic background of milk Fourier transform infrared (FTIR) spectra in dairy sheep. Individual milk FTIR spectra, with 1060 wavenumbers each, were available for 793 adult Sarda breed ewes genotyped at 45,813 SNP. The absorbance values of each wavenumber was analysed using a linear mixed model that included dim class, parity and lambing month as fixed effects and flock-test date and animal as random effects. The model was applied to estimate variance components and heritability and to perform a genome-wide association study for each wavenumber. Average h2 of wavenumbers absorbance was 0.13 ± 0.08, with the largest values observed in the regions associated with the characteristic bonds of carbonylic and methylenic groups of milk fat (h2 = 0.57 at 1724-1728 cm-1; and h2 = 0.34 at 2811-2834 cm-1, respectively). The absorbance values of wavenumbers were moderately correlated with the estimated heritabilities. After the Bonferroni correction, a total of nine markers were found to be significantly associated with 32 different wavenumbers. Of particular interest was the SNP s63269.1, mapped on chromosome 2, that was found to be associated with 27 wavenumbers. Genes previously found to be related to traits of interest (e.g. disease resistance, milk yield and quality, cheese firmness) are located close to the significant SNP. As expected, the heritability estimated for the absorbance of each wavenumbers seems to be associated with the related milk components.

Toward Sustainability: An Overview of the Use of Green Hydrogen in the Agriculture and Livestock Sector.

The agro-livestock sector produces about one third of global greenhouse gas (GHG) emissions. Since more energy is needed to meet the growing demand for food and the industrial revolution in agriculture, renewable energy sources could improve access to energy resources and energy security, reduce dependence on fossil fuels, and reduce GHG emissions. Hydrogen production is a promising energy technology, but its deployment in the global energy system is lagging. Here, we analyzed the theoretical and practical application of green hydrogen generated by electrolysis of water, powered by renewable energy sources, in the agro-livestock sector. Green hydrogen is at an early stage of development in most applications, and barriers to its large-scale deployment remain. Appropriate policies and financial incentives could make it a profitable technology for the future.

A transgenerational study on the effect of great-granddam birth month on granddaughter EBV for production traits in Italian Simmental cattle.

Heat tolerance is a key feature of resilient animals. Offspring of animals that suffer environmental stress during pregnancy could show physiological, morphological, and metabolic modifications. This is due to a dynamic reprogramming of the epigenetics of the mammalian genome that occurs in the early life cycle. Thus, the aim of this study was to investigate the extent of the transgenerational effect of heat stress during the pregnancy of Italian Simmental cows. The effects of dam and granddam birth months (as indicator of pregnancy period) on their daughter and granddaughter estimated breeding values (EBV) for some dairy traits as well as of the temperature-humidity index (THI) during the pregnancy were tested. A total of 128,437 EBV (milk, fat, and protein yields, and somatic cell score) were provided by the Italian Association of Simmental Breeders. The best birth months (of both dam and granddam) for milk yield and protein yield were May and June, whereas the worst were January and March. Great-granddam pregnancies developed during the winter and spring seasons positively affected the EBV for milk and protein yields of their great-granddaughters; in contrast, pregnancies during summer and autumn had negative effects. These findings were confirmed by the effects of maximum and minimum THI in different parts of the great-granddam pregnancy on the performances of their great-granddaughters. Thus, a negative effect of high temperatures during the pregnancy of female ancestors was observed. Results of the present study suggest a transgenerational epigenetic inheritance in Italian Simmental cattle due to environmental stressors.

Dry Matter Intake Prediction from Milk Spectra in Sarda Dairy Sheep.

Individual dry matter intake (DMI) is a relevant factor for evaluating feed efficiency in livestock. However, the measurement of this trait on a large scale is difficult and expensive. DMI, as well as other phenotypes, can be predicted from milk spectra. The aim of this work was to predict DMI from the milk spectra of 24 lactating Sarda dairy sheep ewes. Three models (Principal Component Regression, Partial Least Squares Regression, and Stepwise Regression) were iteratively applied to three validation schemes: records, ewes, and days. DMI was moderately correlated with the wavenumbers of the milk spectra: the largest correlations (around ±0.30) were observed at ~1100-1330 cm-1 and ~2800-3000 cm-1. The average correlations between real and predicted DMI were 0.33 (validation on records), 0.32 (validation on ewes), and 0.23 (validation on days). The results of this preliminary study, even if based on a small number of animals, demonstrate that DMI can be routinely estimated from the milk spectra.

Analysis of runs of homozygosity of cattle living in different climate zones.

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Domestication and evolution of cattle originated different modern breeds in different places worldwide. The interaction between natural and artificial selection and the adaptation to environment shaped the genome, and the three different types of cattle here considered (taurus, indicus, and composite) may harbor different selection signatures. To study the difference among types and region of origin (tropics, temperate, and continental zones), Regions of Homozygosity (ROH) were used. ROH are continuous homozygous chromosomal segments identical by descendent, which characteristics can give information about inbreeding occurrence and natural and artificial selection. Moreover, it had been investigated which genes were mapped in these regions, and if interesting differences pertaining to environmental adaptation or fitness in general, could be observed.

A Preliminary Study on the Interplay between the Serum Levels of Neurotransmitters and Thyroid Hormones for the Evaluation of the Behavioral Phenotype of Dogs.

A total of 112 dogs (49 males and 63 females) belonging to different breeds (i.e., Boxer, Cirneco dell'Etna, Fonni's Dog, Labrador, Crossbreed, German Shepherd, Pit Bull, Shar-Pei, Yorkshire) were analyzed to compare the serum concentration of serotonin, dopamine, norepinephrine, prolactin, beta-endorphins, thyroxine (T4), triiodothyronine (T3), thyroid-stimulating hormone (TSH), and assess whether these parameters can be correlated with the behavioral phenotype of the investigated breeds. T4 was above or below the threshold in 61% and 14% of dogs, respectively; T3, in contrast, 41% of dogs showed values below the limit, while 26% above it. TSH was within the reference range in 58% of dogs; 94% of the dogs had prolactin in the reference range and only five animals showed values above the limit. For beta-endorphins, 49% of dogs had values above the limit, while 46% had values within the reference range. Serotonin and dopamine values below physiological limits were found in 62% and 70% of dogs, respectively. Finally, 61% of the dogs showed norepinephrine values within the reference range. The study confirmed that the assessment of the serum values of hormones and neurotransmitters in dogs could be useful to better understand the behavioral phenotype of the animal and could be useful for breeders and trainers for the selection of the most suitable subjects for specific tasks.

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds.

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (FPED ), or runs of homozygosity (ROH; FROH ) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of FROH . The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.

Ecotoxicity of Copper(I) Chloride in Grooved Carpet Shell (Ruditapes decussatus).

Copper (Cu) is a ubiquitous trace element in the aquatic environment, and is usually found at low levels. Copper environmental concentrations can be altered as a result anthropogenic activities. Shellfish are useful bioindicators to ensure adequate environmental monitoring. Thus, the aim of the present study was as follows: (a) determine the LC50 of copper(I) chloride in grooved carpet shell (Ruditapes decussatus) collected in the Santa Gilla lagoon (Sardinia, Italy), and (b) analyze the antioxidant biomarkers in digestive gland and gills of same specimens exposed to different concentrations of the above-mentioned metal (0.045, 0.45, and 0.90 mg/L) for 96 h. A withdrawal period of 96 h was considered for the treated clam, carrying out the same biochemical analyses, superoxide dismutase (SOD), catalase (CAT), selenium-dependent glutathione peroxidase (Se-GPx), glutathione S-transferases (GSTs), and total glutathione (GSH+2GSSG) in the two tissues. Different time and dose responses of the antioxidant biomarkers were recorded in the digestive glands and gills. Oxidative stress biomarkers highlighted the ability of Cu to induce oxidative stress in R. decussatus. Clam, following the withdrawal period of 96 h, has not been able to achieve the control levels of all biochemical markers in the digestive gland and gills. R. decussatus can be a suitable model to assess the ecotoxicity of copper in aquatic ecosystems. These findings may advance knowledge on the role and the effects of copper on oxidative stress biomarkers in grooved carpet shell. The metal ecotoxicity response can be useful to perform accurate biomarker-based monitoring programs using this bivalve species.

Using Pedigree and Genomic Data toward Better Management of Inbreeding in Italian Dairy Sheep and Goat Breeds.

The inbreeding coefficient is an important parameter for livestock management. Small ruminant breeders and associations mainly rely on pedigree information, but genomic tools are gaining relevance, overcoming possible pedigree inconsistencies. This study investigates the relationship between pedigree-based and genomic inbreeding in two goat and four sheep dairy breeds. Pedigree and genomic data (medium-density SNPchip) were obtained for 3107 goats and 2511 sheep. We estimated pedigree depth (number of fully traced generations, FullGen) and inbreeding (FPED), as well as two genomic inbreeding indexes, using runs of the homozygosity (FROH) and genomic relationship matrix (FGRM). The correlation between the inbreeding coefficients was assessed. A linear regression model (LRM) was fitted for estimating FPED from FROH. After quality control on genomic data, we retained 5085 animals. Mean inbreeding values were low, with higher FROH than FPED and FGRM. Breed differences can partially depend on different managements. The correlation between FPED and FROH was the highest and directly related to pedigree depth. The best LRM was chosen for FullGen ≥4 and ≥6 for goats and sheep, respectively; after excluding animals with extreme residuals, a new refined regression equation was calculated. Since massive genotyping is not affordable to small ruminant breeders, it is important to understand the distinction and relationship between differently calculated inbreeding coefficients, also in view of the introduction of genomic enhanced breeding values. Our study highlights the importance of accurate pedigree information and, especially if not obtainable, of calculating genomic-based inbreeding coefficients. A better estimation of animals' relatedness contributes to improve animal breeding and conservation.

Assessment of Genetic Diversity, Runs of Homozygosity, and Signatures of Selection in Tropical Milking Criollo Cattle Using Pedigree and Genomic Data.

The objective of this study was to evaluate the genetic diversity of the Tropical Milking Criollo cattle (TMC) breed in Mexico through parameters derived from pedigree and genomic information assessment. The pedigree file consisted of 3780 animals. Seventy-nine bovines were genotyped with the medium-density single nucleotide polymorphism chip and considered a reference population for pedigree analysis. The effective population size and the probability of gene origin used to assess the evolution of genetic diversity were calculated with pedigree information. Inbreeding coefficients were evaluated based on pedigree (FPed), the genomic relationship matrix (FGRM), and runs of homozygosity (FROH) of different length classes. The average inbreeding was 2.82 ± 2.66%, −0.7 ± 3.8%, and 10.9 ± 3.0% for FPED, FGRM, and FROH, respectively. Correlation between FPED and FROH was significant only for runs of homozygosity > 4 Mb, indicating the FPED of a population with an average equivalent complete generation of five only recovers the most recent inbreeding. The parameters of the probability of gene origin indicated the existence of genetic bottlenecks and the loss of genetic diversity in the history of the TMC cattle population; however, pedigree and genomic information revealed the existence of current sufficient genetic diversity to design a sustainable breeding program.

Feedlot pens with greenhouse roofs improve beef cattle performance in temperate weather.

Muddy pens can negatively affect welfare and performance of feedlot beef cattle. In some regions with temperate weather, plastic greenhouse covers, above the entire pens are used to fatten cattle in a clean and dry environment. The objective of this research was to investigate effects of greenhouse roofed pens on beef cattle feedlot performance in temperate weather. Data were collected from a feedlot located in Central Mexico between 2016 and 2019. The study included 1,062 closeouts of pens with 68,305 crossbred bulls fed in pens with or without a greenhouse roof. Feeding ranged from 82 to 210 d, depending on the initial weight of cattle, which ranged from 255 to 511 kg. For each pen, average daily dry matter intake (DMI; kg of DMI·animal-1·d-1), average daily gain (ADG, kg·animal-1·d-1), and feed efficiency (G:F, ADG/DMI) were measured. Factorial analyses were performed to test the interaction and main effects of initial weight grouping (light, medium, and heavy), roof, and season as fixed effects, and year as a random effect. None of the three-way interactions were significant (P > 0.51). There was no initial weight grouping × roof interactions for DMI and ADG (P > 0.31). There was (P = 0.03) an initial weight grouping × roof interaction for G:F, as pens of all initial weight groups had greater (P < 0.01) G:F in pens with greenhouse roof than its counterpart in pens without greenhouse roof, but the advantage was greater for pens with light cattle (0.178 vs. 0.166; P < 0.01). There was no initial weight grouping × season interactions for all variables (P > 0.39). There was no roof × season interaction for DMI (P = 0.47), but there were interactions for ADG and G:F (P < 0.01). The ADG was not different (P > 0.13) during summer and autumn based on the roofing system, but pens with greenhouse roofs had greater ADG during spring (1.70 vs. 1.61) and winter (1.68 vs. 1.64; P ≤ 0.01). The G:F was greater (P < 0.01) in all seasons for pens with a greenhouse roof, with the most prominent advantage during spring (0.173 vs. 0.160). There were main effects for cattle initial weight grouping and roof for all variables (P < 0.01). Season affected DMI and G:F (P < 0.01). Pens with greenhouse roofs had decreased DMI (9.70 vs. 9.86), greater ADG (1.67 vs. 1.63), and increased G:F (0.173 vs. 0.166) compared to pens without greenhouse roofs (P < 0.01). Pens with greenhouse roofs in feedlots located in temperate regions positively affect beef cattle performance.

A Dynamic Model for Estimating the Interaction of ROS-PUFA-Antioxidants in Rabbit.

Defining optimal nutrition in animals and humans remains a main scientific challenge. The objective of the work was to develop a dynamic model of reactive oxygen species (ROS)-polyunsaturated fatty acid (PUFA)-antioxidant homeostasis using the rabbit as a model. The problem entity was to evaluate the main metabolites generated from interactions between traits included in the conceptual model and identified by three main sub-models: (i) ROS generation, (ii) PUFA oxidation and (iii) antioxidant defence. A mathematical model (VENSIM software) that consisted of molecular stocks (INPUTs, OUTPUTs), exchange flows (intermediate OUTPUTs) and process rates was developed. The calibration was performed by using standard experimental data (Experiment 1), whereas the validation was carried out in Experiments 2 and 3 by using supra-nutritional dietary inputs (VIT E+ and PUFA+). The accuracy of the models was measured using 95% confidence intervals. Analytical OUTPUTs (ROS, PUFA, Vit E, Ascorbic acid, Iso-/NeuroProstanes, Aldehydes) were well described by the standard model. There was also good accuracy for the VIT E+ scenario, whereas some compensatory rates (Kc1-Kc4) were added to assess body compensation when high levels of dietary PUFA were administered (Experiment 3). In conclusion, the model can be very useful for predicting the effects of dietary treatments on the redox homeostasis of rabbits.

Identification of differentially expressed genes in early-postmortem Semimembranosus muscle of Italian Large White heavy pigs divergent for glycolytic potential.

Glycolytic potential (GP) is one of the postmortem traits used to predict the quality of the final meat products. Despite that, the knowledge of the molecular and metabolic pathways controlling this trait is still not complete. To add some information on this field we used two pools of Italian Large White heavy pigs divergent for GP to investigate through a microarray approach the differences of gene expressions between the two pools. On the whole, 32 genes were differentially expressed, and among them 31 were overexpressed in low GP pool. These genes were involved in mitochondrial functions and adenosine triphosphate (ATP) biosynthetic processes, in calcium homeostasis, and in lipid metabolism, with Peroxisome proliferator-activated receptor (PPAR) signaling being a possible master regulator of the molecular differences observed between the two pools. The different GP levels between the two pools could have determined in low GP muscles a more rapid occurrence of the molecular cascade related to the events triggering cell death.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity.

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

Occurrence of microplastics in the gastrointestinal tract of benthic by-catches from an eastern Mediterranean deep-sea environment.

Concern about microplastic pollution little is known about levels in deep-sea species; to fill this knowledge gap, levels of microplastics in the gastrointestinal (GI) tracts of 34 fish from eight different deep-sea by-catches: blackmouth catshark, lesser spotted dogfish, and velvet belly, armless snake eel, hollowsnout grenadier, phaeton dragonet, royal flagfin, and slender snipe eel were measured. All were collected at the same site (east Sardinia, Mediterranean Sea; 40°10'12.49″N, 9°44'12.31″E) using a bottom gillnet at depths between -820/250 and -1148 ft./350 m. Microplastics (MPs) were retrieved in 16 out of 34 fish. At least one microplastic item was found in 48% (33%, E. spinax - 75%, G. melastomus) of the samples. The most frequent was polyethylene (PE), with nine items (filaments, films, fragments) found in five specimens. This preliminary study of by-catches adds new data on MPs ingestion by species inhabiting a deep-sea environment of the Mediterranean.

Investigation of ß-hydroxybutyrate in early lactation of Simmental cows: Genetic parameters and genomic predictions.

Genomic information allows for a more accurate calculation of relationships among animals than the pedigree information, leading to an increase in accuracy of breeding values. Here, we used pedigree-based and single-step genomic approaches to estimate variance components and breeding values for ß-hydroxybutyrate milk content (BHB). Additionally, we performed a genome-wide association study (GWAS) to depict its genetic architecture. BHB concentrations within the first 90 days of lactation, estimated from milk medium infrared spectra, were available for 30,461 cows (70,984 records). Genotypes at 42,152 loci were available for 9,123 animals. Low heritabilities were found for BHB using pedigree-based (0.09 ± 0.01) and genomic (0.10 ± 0.01) approaches. Genetic correlation between BHB and milk traits ranged from -0.27 ± 0.06 (BHB and protein percentage) to 0.13 ± 0.07 (BHB and fat-to-protein ratio) using pedigree and from -0.26 ± 0.05 (BHB and protein percentage) to 0.13 ± 0.06 (BHB and fat-to-protein ratio) using genomics. Breeding values were validated for 344 genotyped cows using linear regression method. The genomic EBV (GEBV) had greater accuracy (0.51 vs. 0.45) and regression coefficient (0.98 vs. 0.95) compared to EBV. The correlation between two subsequent evaluations, without and with phenotypes for validation cows, was 0.85 for GEBV and 0.82 for EBV. Predictive ability (correlation between (G)EBV and adjusted phenotypes) was greater when genomic information was used (0.38) than in the pedigree-based approach (0.31). Validation statistics in the pairwise two-trait models (milk yield, fat and protein percentage, urea, fat/protein ratio, lactose and logarithmic transformation of somatic cells count) were very similar to the ones highlighted for the single-trait model. The GWAS allowed discovering four significant markers located on BTA20 (57.5-58.2 Mb), where the ANKH gene is mapped. This gene has been associated with lactose, alpha-lactalbumin and BHB. Results of this study confirmed the usefulness of genomic information to provide more accurate variance components and breeding values, and important insights about the genomic determination of BHB milk content.

Farm Animals Are Long Away from Natural Behavior: Open Questions and Operative Consequences on Animal Welfare.

The concept of welfare applied to farm animals has undergone a remarkable evolution. The growing awareness of citizens pushes farmers to guarantee the highest possible level of welfare to their animals. New perspectives could be opened for animal welfare reasoning around the concept of domestic, especially farm, animals as partial human artifacts. Therefore, it is important to understand how much a particular behavior of a farm animal is far from the natural one of its ancestors. This paper is a contribution to better understand the role of genetics of the farm animals on their behavior. This means that the naïve approach to animal welfare regarding returning animals to their natural state should be challenged and that welfare assessment should be considered.

Genetic Background and Inbreeding Depression in Romosinuano Cattle Breed in Mexico.

The ultimate goal of genetic selection is to improve genetic progress by increasing favorable alleles in the population. However, with selection, homozygosity, and potentially harmful recessive alleles can accumulate, deteriorating genetic variability and hampering continued genetic progress. Such potential adverse side effects of selection are of particular interest in populations with a small effective population size like the Romosinuano beef cattle in Mexico. The objective of this study was to evaluate the genetic background and inbreeding depression in Mexican Romosinuano cattle using pedigree and genomic information. Inbreeding was estimated using pedigree (FPED) and genomic information based on the genomic relationship matrix (FGRM) and runs of homozygosity (FROH) of different length classes. Linkage disequilibrium (LD) was evaluated using the correlation between pairs of loci, and the effective population size (Ne) was calculated based on LD and pedigree information. The pedigree file consisted of 4875 animals born between 1950 and 2019, of which 71 had genotypes. LD decreased with the increase in distance between markers, and Ne estimated using genomic information decreased from 610 to 72 animals (from 109 to 1 generation ago), the Ne estimated using pedigree information was 86.44. The reduction in effective population size implies the existence of genetic bottlenecks and the decline of genetic diversity due to the intensive use of few individuals as parents of the next generations. The number of runs of homozygosity per animal ranged between 18 and 102 segments with an average of 55. The shortest and longest segments were 1.0 and 36.0 Mb long, respectively, reflecting ancient and recent inbreeding. The average inbreeding was 2.98 ± 2.81, 2.98 ± 4.01, and 7.28 ± 3.68% for FPED, FGRM, and FROH, respectively. The correlation between FPED and FGRM was -0.25, and the correlations among FPED and FROH of different length classes were low (from 0.16 to 0.31). The correlations between FGRM and FROH of different length classes were moderate (from 0.44 to 0.58), indicating better agreement. A 1% increase in population inbreeding decreased birth weight by 0.103 kg and weaning weight by 0.685 kg. A strategy such as optimum genetic contributions to maximize selection response and manage the long-term genetic variability and inbreeding could lead to more sustainable breeding programs for the Mexican Romosinuano beef cattle breed.