Identification of genetic markers associated with hyperketonemia patterns in early lactation Holstein cows.

Ketosis, evidenced by hyperketonemia with elevated blood ß-hydroxybutyrate (BHB) levels, is a significant metabolic disorder of dairy cattle, typically diagnosed within the first 6 weeks post-calving when high energy levels are essential to milk production. Our study aimed to identify genetic markers linked to hyperketonemia (HYK) patterns in Holstein cows during early lactation and compare these to HYK-negative cows. We screened 964 cows for HYK using a threshold of BHB ≥1.2 mmol/L during the first 2 weeks postpartum (screening period, SP). Cows that tested negative initially were retested the following week. Cows were deemed HYK-negative (CON group) if BHB levels were below 1.2 mmol/L in both tests, while those with BHB levels exceeding this threshold at any test were treated and classified as HYK-positive (HYK+). Post-treatment, HYK+ cows were monitored for two-week follow-up period (FP) and classified based on their recovery: cured (CUR; consistently low BHB), recurrent (REC; fluctuating BHB levels), severe (SEV; high initial BHB that decreased), or chronic (CHR; persistently high BHB). Using 489 cows that were genotyped, a GWAS was conducted using GCTA software, revealing significant associations of several SNPs across different HYK patterns when compared to the CON group. These SNPs were primarily linked to genes affecting milk traits and were enriched in biological pathways relevant to protein glycosylation, inflammatory response, glucose homeostasis, and fatty acid synthesis. Our findings highlight genomic regions, potential candidate genes, and biological pathways related to ketosis, underscoring potential targets for improving health management in dairy cattle. These insights could lead to better strategies for managing ketosis through genetic selection, ultimately enhancing dairy cattle welfare and productivity. Further research with a larger number of cows is recommended to validate these findings and help confirm the implicated SNPs and genes.

Genome-wide association study in thoroughbred horses naturally infected with cyathostomins.

Cyathostomins are considered one of the most important parasites of horses. A group of horses within a herd can be responsible for eliminating the majority of parasite eggs. This phenotype might be explained by genetic factors. This study aimed to identify genomic regions associated with fecal egg count (FEC) and hematological parameters by performing a genomic-wide association study (GWAS) in Thoroughbred horses naturally infected with cyathostomins. Packed cell volume (PCV), differential leukocyte, and FEC were determined from 90 horses. All animals were genotyped using the Illumina Equine 70 K BeadChip panel containing 65,157 SNP markers. The five genomic windows that have explained the highest percentage of the additive genetic variance of a specific trait (top 5) were further explored to identify candidate genes. A total of 33, 21, 30, 21, and 19 genes were identified for FEC, PCV, eosinophils, neutrophils, and lymphocyte count, respectively. The top 5 marker regions explained 2.86, 2.56, 2.73, 2.33, and 2.37% of the additive genetic variation of FEC, PCV, eosinophils, neutrophils, and lymphocytes count, respectively. This is the first study correlating phenotypic horse health traits to GWAS analysis, which may be used for animal breeding activities, reducing losses due to parasite infections.

Accounting for population structure in selective cow genotyping strategies.

The objective of the present study was to investigate the impact of considering population structure in cow genotyping strategies over the accuracy and bias of genomic predictions. A small dairy cattle population was simulated to address these objectives. Based on four main traditional designs (random, top-yield, extreme-yield and top-accuracy cows), different numbers (1,000; 2,000 and 5,000) of cows were sampled and included in the reference population. Traditional designs were replicated considering or not population structure and compared among and with a reference population containing only bulls. The inclusion of cows increased accuracy in all scenarios compared with using only bulls. Scenarios accounting for population structure when choosing cows to the reference population slightly outperformed their traditional versions by yielding higher accuracy and lower bias in genomic predictions. Building a cow-based reference population from groups of related individuals considering the frequency of individuals from those same groups in the validation population yielded promising results with applications on selection for expensive- or difficult-to-measure traits. Methods here presented may be easily implemented in both new or already established breeding programs, as they improved prediction and reduced bias in genomic evaluations while demanding no additional costs.

Genotype imputation in a tropical crossbred dairy cattle population.

The objective of this study was to investigate different strategies for genotype imputation in a population of crossbred Girolando (Gyr × Holstein) dairy cattle. The data set consisted of 478 Girolando, 583 Gyr, and 1,198 Holstein sires genotyped at high density with the Illumina BovineHD (Illumina, San Diego, CA) panel, which includes ∼777K markers. The accuracy of imputation from low (20K) and medium densities (50K and 70K) to the HD panel density and from low to 50K density were investigated. Seven scenarios using different reference populations (RPop) considering Girolando, Gyr, and Holstein breeds separately or combinations of animals of these breeds were tested for imputing genotypes of 166 randomly chosen Girolando animals. The population genotype imputation were performed using FImpute. Imputation accuracy was measured as the correlation between observed and imputed genotypes (CORR) and also as the proportion of genotypes that were imputed correctly (CR). This is the first paper on imputation accuracy in a Girolando population. The sample-specific imputation accuracies ranged from 0.38 to 0.97 (CORR) and from 0.49 to 0.96 (CR) imputing from low and medium densities to HD, and 0.41 to 0.95 (CORR) and from 0.50 to 0.94 (CR) for imputation from 20K to 50K. The CORRanim exceeded 0.96 (for 50K and 70K panels) when only Girolando animals were included in RPop (S1). We found smaller CORRanim when Gyr (S2) was used instead of Holstein (S3) as RPop. The same behavior was observed between S4 (Gyr + Girolando) and S5 (Holstein + Girolando) because the target animals were more related to the Holstein population than to the Gyr population. The highest imputation accuracies were observed for scenarios including Girolando animals in the reference population, whereas using only Gyr animals resulted in low imputation accuracies, suggesting that the haplotypes segregating in the Girolando population had a greater effect on accuracy than the purebred haplotypes. All chromosomes had similar imputation accuracies (CORRsnp) within each scenario. Crossbred animals (Girolando) must be included in the reference population to provide the best imputation accuracies.

Genome-wide association analysis of feed intake and residual feed intake in Nellore cattle.

BACKGROUND: Feed intake plays an important economic role in beef cattle, and is related with feed efficiency, weight gain and carcass traits. However, the phenotypes collected for dry matter intake and feed efficiency are scarce when compared with other measures such as weight gain and carcass traits. The use of genomic information can improve the power of inference of studies on these measures, identifying genomic regions that affect these phenotypes. This work performed the genome-wide association study (GWAS) for dry matter intake (DMI) and residual feed intake (RFI) of 720 Nellore cattle (Bos taurus indicus). RESULTS: In general, no genomic region extremely associated with both phenotypic traits was observed, as expected for the variables that have their regulation controlled by many genes. Three SNPs surpassed the threshold for the Bonferroni multiple test for DMI and two SNPs for RFI. These markers are located on chromosomes 4, 8, 14 and 21 in regions near genes regulating appetite and ion transport and close to important QTL as previously reported to RFI and DMI, thus corroborating the literature that points these two processes as important in the physiological regulation of intake and feed efficiency. CONCLUSIONS: This study showed the first GWAS of DMI to identify genomic regions associated with feed intake and efficiency in Nellore cattle. Some genes and QTLs previously described for DMI and RFI, in other subspecies (Bos taurus taurus), that influences these phenotypes are confirmed in this study.

Genotype by Prenatal Environment Interaction for Postnatal Growth of Nelore Beef Cattle Raised under Tropical Grazing Conditions.

The prenatal environment is recognized as crucial for the postnatal performance in cattle. In tropical regions, pregnant beef cows commonly experience nutritional restriction during the second half of the gestation period. Thus, the present study was designed to analyze the genotype by prenatal environment interaction (G × Epn) and to identify genomic regions associated with the level and response in growth and reproduction-related traits of beef cattle to changes in the prenatal environment. A reaction norm model was applied to data from two Nelore herds using the solutions of contemporary groups for birth weight as a descriptor variable of the gestational environment quality. A better gestational environment favored weights until weaning, scrotal circumference at yearling, and days to first calving of the offspring. The G × Epn was strong enough to result in heterogeneity of variance components and genetic parameters in addition to reranking of estimated breeding values and SNPs effects. Several genomic regions associated with the level of performance and specific responses of the animals to variations in the gestational environment were revealed, which harbor QTLs and can be exploited for selection purposes. Therefore, genetic evaluation models considering G × Epn and special management and nutrition care for pregnant cows are recommended.

Fine mapping of genomic regions associated with female fertility in Nellore beef cattle based on sequence variants from segregating sires.

BACKGROUND: Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers. METHODS: The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries. RESULTS: A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles' QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG "Circadian rhythm" and "Neurotrophin signaling pathway" were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones. CONCLUSIONS: The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.