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

Inheritance of genomic regions and genes associated with number of oocytes and embryos in Gir cattle through daughter design

Over the past decades, daughter designs, including genotyped sires and their genotyped daughters, have been used as an approach to identify QTL related to economic traits. The aim of this study was to identify genomic regions inherited by Gir sire families and genes associated with number of viable oocytes (VO), total number of oocytes (TO), and number of embryos (EMBR) based on a daughter design approach. In total, 15 Gir sire families were selected. The number of daughters per family ranged from 26 to 395, which were genotyped with different SNP panels and imputed to the Illumina BovineHD BeadChip (777K) and had phenotypes for oocyte and embryo production. Daughters had phenotypic data for VO, TO, and EMBR. The search for QTL was performed through GWAS based on GBLUP. The QTL were found for each trait among and within families based on the top 10 genomic windows with the greatest genetic variance. For EMBR, genomic windows identified among families were located on BTA4, BTA5, BTA6, BTA7, BTA8, BTA13, BTA16, and BTA17, and they were most frequent on BTA7 within families. For VO, genomic windows were located on BTA2, BTA4, BTA5, BTA7, BTA17, BTA21, BTA22, BTA23, and BTA27 among families, being most frequent on BTA8 within families. For TO, the top 10 genomic windows were identified on BTA2, BTA4, BTA5, BTA7, BTA17, BTA21, BTA22, BTA26, and BTA27, being most frequent on BTA7 and BTA8 within families. Considering all results, the greatest number of genomic windows was found on BTA7, where the VCAN, XRCC4, TRNAC-ACA, HAPLN1, and EDIL3 genes were identified in the common regions. In conclusion, 15 Gir sire families with 26 to 395 daughters per family with phenotypes for oocyte and embryo production helped to identify the inheritance of several genomic regions, especially on BTA7, where the EDIL3, HAPLN1, and VCAN candidate genes were associated with number of oocytes and embryos in Gir cattle families.

Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle

Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.

Sustainable Intensification of Beef Production in the Tropics: The Role of Genetically Improving Sexual Precocity of Heifers

Simple Summary

Tropical pasture-based beef production systems play a vital role in global food security. The importance of promoting sustainable intensification of such systems has been debated worldwide. Demand for beef is growing together with concerns over the impact of its production on the environment. Implementing sustainable livestock intensification programs relies on animal genetic improvement. In tropical areas, the lack of sexual precocity is a bottleneck for cattle efficiency, directly impacting the sustainability of production systems. In the present review we present and discuss the state of the art of genetic evaluation for sexual precocity in Bos indicus beef cattle, covering the definition of measurable traits, genetic parameter estimates, genomic analyses, and a case study of selection for sexual precocity in Nellore breeding programs.

Abstract

Increasing productivity through continued animal genetic improvement is a crucial part of implementing sustainable livestock intensification programs. In Zebu cattle, the lack of sexual precocity is one of the main obstacles to improving beef production efficiency. Puberty-related traits are complex, but large-scale data sets from different “omics” have provided information on specific genes and biological processes with major effects on the expression of such traits, which can greatly increase animal genetic evaluation. In addition, genetic parameter estimates and genomic predictions involving sexual precocity indicator traits and productive, reproductive, and feed-efficiency related traits highlighted the feasibility and importance of direct selection for anticipating heifer reproductive life. Indeed, the case study of selection for sexual precocity in Nellore breeding programs presented here show that, in 12 years of selection for female early precocity and improved management practices, the phenotypic means of age at first calving showed a strong decreasing trend, changing from nearly 34 to less than 28 months, with a genetic trend of almost −2 days/year. In this period, the percentage of early pregnancy in the herds changed from around 10% to more than 60%, showing that the genetic improvement of heifer’s sexual precocity allows optimizing the productive cycle by reducing the number of unproductive animals in the herd. It has a direct impact on sustainability by better use of resources. Genomic selection breeding programs accounting for genotype by environment interaction represent promising tools for accelerating genetic progress for sexual precocity in tropical beef cattle.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Investigating anogenital distance and antral follicle count as novel markers of fertility within a herd of cows with positive or negative genetic merit for fertility traits

Using early-in-life markers of reproductive characteristics may enhance the speed and success of genetic improvement in fertility. We investigated 2 phenotypes that can be measured early in life and are moderately heritable to determine their association with traditional measures of reproductive success and genetic variation under a seasonal-calving, pasture-based system. Cows were bred to be divergent in the New Zealand Fertility Breeding Value, which estimates genetic merit for fertility. Cows consisted of 2 groups with an average positive (+5%) or negative (-5%) genetic merit for fertility traits and were expected to have large diversity in reproductive outcomes. Calves were genotyped at 41 ± 3.1 d of age (mean ± SD; n = 538), and antral follicle counts (AFC) were measured when they were postpubertal heifers before their first breeding (416 ± 15 d old; 92-d range; n = 520). The anogenital distance (AGD) was measured in 478 primiparous cows of this same population 50 to 60 d after the breeding start date when they were 881 ± 25 d old (145-d range). The AGD was shorter in animals with a positive genetic merit for fertility traits (based on parent averages). An indicator of herd reproductive success in a seasonal-calving system (recalving by 6 wk in lactation 2) was chosen for logistic regression with cross-validation, and if significant, a cut-off was calculated that categorized animals into groups. Both linear and quadratic regression was undertaken, and the model with the greatest sensitivity for detection of nonpregnant cows used. The AGD linear model was significant with a sensitivity of 64% and a specificity of 48%. This model resulted in a cut-off of 102 mm, which we used to classify cows as short (≤102 mm) or long (>102 mm) AGD animals. Primiparous cows with a short AGD were more likely to be pregnant within the first 3 and 6 wk of mating, and become pregnant as a primiparous cow, than those with a long AGD. The time from calving to conception was 20 d earlier in short AGD compared with long-AGD cows. None of the models tested for AFC were significant; therefore cows were categorized into 3 groups based on previous work in seasonal systems. However, associations between fertility phenotypes and AFC group were limited. Genomic regions of interest for AGD and AFC did not overlap, indicating phenotypes were genetically independent. Overall, AGD appears as a promising early marker of fertility in seasonal grazing systems.