Single-step genomic best linear unbiased predictor genetic parameter estimations and genome-wide associations for milk fatty acid profiles, interval from calving to first insemination, and ketosis in Holstein cattle

Genetic evaluations and genome-wide association studies for specific digital dermatitis diagnoses in dairy cows considering genotype × housing system interactions

The present study aimed to use detailed phenotyping for the claw disorder digital dermatitis (DD) considering specific DD stages in 2 housing systems (conventional cubicle barns [CON] and compost-bedded pack barns [CBPB]) to infer possible genotype × housing system interactions. The DD stages included 2,980 observations for the 3 traits DD-sick, DD-acute, and DD-chronic from 1,311 Holstein-Friesian and 399 Fleckvieh-Simmental cows. Selection of the 5 CBPB and 5 CON herds was based on a specific protocol to achieve a high level of herd similarity with regard to climate, feeding, milking system, and location, but with pronounced housing-system differences. Five other farms had a "mixed system" with 2 subherds, one representing CBPB and the other one CON. The CBPB system was represented by 899 cows (1,530 observations), and 811 cows (1,450 observations) represented the CON system. The average disease prevalence was 20.47% for DD-sick, 13.88% for DD-acute, and 5.34% for DD-chronic, with a higher prevalence in CON than in CBPB. After quality control of 50K genotypes, 38,495 SNPs from 926 cows remained for the ongoing genomic analyses. Genetic parameters for DD-sick, DD-acute, and DD-chronic were estimated by applying single-step approaches for single-trait repeatability animal models considering the whole dataset, and separately for the CON and CBPB subsets. Genetic correlations between same DD traits from different housing systems, and between DD-sick, DD-chronic, and DD-acute, were estimated via bivariate animal models. Heritabilities based on the whole dataset were 0.16 for DD-sick, 0.14 for DD-acute, and 0.11 for DD-chronic. A slight increase of heritabilities and genetic variances was observed in CON compared with the "well-being" CBPB system, indicating a stronger genetic differentiation of diseases in a more challenging environment. Genetic correlations between same DD traits recorded in CON or CBPB were close to 0.80, disproving obvious genotype × housing system interactions. Genetic correlations among DD-sick, DD-acute and DD-chronic ranged from 0.58 to 0.81. SNP main effects and SNP × housing system interaction effects were estimated simultaneously via GWAS, considering only the phenotypes from genotyped cows. Ongoing annotations of potential candidate genes focused on chromosomal segments 100 kb upstream and downstream from the significantly associated candidate SNP. GWAS for main effects indicated heterogeneous Manhattan plots especially for DD-acute and DD-chronic, indicating particularities in disease pathogenesis. Nevertheless, a few shared annotated potential candidate genes, that is, METTL25, AFF3, PRKG1, and TENM4 for DD-sick and DD-acute, were identified. These genes have direct or indirect effects on disease resistance or immunology. For the SNP × housing system interaction, the annotated genes ASXL1 and NOL4L on BTA 13 were relevant for DD-sick and DD-acute. Overall, the very similar genetic parameters for the same traits in different environments and negligible genotype × housing system interactions indicate only minor effects on genetic evaluations for DD due to housing-system particularities.

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.

Detection of Candidate Genes Associated with Fecundity through Genome-Wide Selection Signatures of Katahdin Ewes

One of the strategies to genetically improve reproductive traits, despite their low inheritability, has been the identification of candidate genes. Therefore, the objective of this study was to detect candidate genes associated with fecundity through the fixation index (F_ST) and runs of homozygosity (ROH) of selection signatures in Katahdin ewes. Productive and reproductive records from three years were used and the genotypes (OvineSNP50K) of 48 Katahdin ewes. Two groups of ewes were identified to carry out the genetic comparison: with high fecundity (1.3 ± 0.03) and with low fecundity (1.1 ± 0.06). This study shows for the first time evidence of the influence of the CNOT11, GLUD1, GRID1, MAPK8, and CCL28 genes in the fecundity of Katahdin ewes; in addition, new candidate genes were detected for fecundity that were not reported previously in ewes but that were detected for other species: ANK2 (sow), ARHGAP22 (cow and buffalo cow), GHITM (cow), HERC6 (cow), DPF2 (cow), and TRNAC-GCA (buffalo cow, bull). These new candidate genes in ewes seem to have a high expression in reproduction. Therefore, future studies are needed focused on describing the physiological basis of changes in the reproductive behavior influenced by these genes.