Single-step genome-wide association studies and post-GWAS analyses for the number of oocytes and embryos in Gir cattle.

Genome-Wide Association Studies (GWAS) are used for identification of quantitate trait loci (QTL) and genes associated with several traits. We aimed to identify genomic regions, genes, and biological processes associated with number of total and viable oocytes, and number of embryos in Gir dairy cattle. A dataset with 17,526 follicular aspirations, including the following traits: number of viable oocytes (VO), number of total oocytes (TO), and number of embryos (EMBR) from 1641 Gir donors was provided by five different stock farms. A genotype file with 2093 animals and 395,524 SNP markers was used to perform a single-step GWAS analysis for each trait. The top 10 windows with the highest percentage of additive genetic variance explained by 100 adjacent SNPs were selected. The genomic regions identified in our work were overlapped with QTLs from QTL database on chromosomes 1, 2, 5, 6, 7, 8, 9, 13, 17, 18, 20, 21, 22, 24, and 29. These QTLs were classified as External, Health, Meat and carcass, Production or Reproduction traits, and about 38% were related to Reproduction. In total, 117 genes were identified, of which 111 were protein-coding genes. Exclusively associations were observed for 42 genes with EMBR, and 1 with TO. Also, 42 genes were in common between VO and TO, 28 between VO and EMBR, and four genes were in common among all traits. In conclusion, great part of the identified genes plays a functional role in initial embryo development or general cell functions. The protein-coding genes ARNT, EGR1, HIF1A, AHR, and PAX2 are good markers for the production of oocytes and embryos in Gir cattle.

Imputation accuracy for genomic selection using embryo biopsy samples in Gir.

The aim of this study was to establish a platform for genomic selection of in vitro-fertilized (IVF) Gir embryos. Multiple displacement amplification (MDA)-based embryo biopsy samples were genotyped, and genomic estimated breeding values (GEBV) for milk yield (305MY) were calculated. The concordance of GEBV and accuracy between embryo biopsies and the respective liveborn were assessed. Imputation was performed using two panels (Z-Chip and Bovine HD, Illumina) based on a database of 73,110 lactating cow's database and pedigree files from 147,131 animals. Biopsied embryos had similar pregnancy rates (39% vs 40%), pregnancy loss rates (18% vs 20%), and pregnancy length compared to Control embryos. After genotyping, low call rate means were detected for biopsy samples compared to the respective calf samples (0.80 vs 0.98). Imputation presented 0.83 (Z-Chip) and 0.96 (HD) accuracy (CORRanim). Embryo GEBV accuracy levels were higher in BovineHD imputation (0.82) than Z-Chip imputation (0.55) or no imputation (0.62), and the correlation between embryo/calf pairs' accuracy was 0.85 for BovineHD imputation, 0.11 for Z-Chip imputation, and 0.02 for no imputation. GEVB estimates correlation between embryo/calf pairs was 0.87 for BovineHD imputation, 0.80 for Z-Chip imputation, and 0.41 before imputation. The call rate of embryo samples did not affect the correlation between embryo/calf pairs for accuracy and GEBV before and after BovineHD imputation. Embryos obtained on the same farm presented GEBV 305MY differences of up to 800 kg, emphasizing the expected impact of embryo genomic selection for the Gir breed.

Runs of homozygosity and signatures of selection for number of oocytes and embryos in the Gir Indicine cattle.

Runs of homozygosity (ROH) and signatures of selection are the results of selection processes in livestock species that have been shown to affect several traits in cattle. The aim of the current work was to verify the profile of ROH and inbreeding depression in the number of total (TO) and viable oocytes (VO) and the number of embryos (EMBR) in Gir Indicine cattle. In addition, we aim to identify signatures of selection, genes, and enriched regions between Gir subpopulations sorted by breeding value for these traits. The genotype file contained 2093 animals and 420,718 SNP markers. Breeding values used to sort Gir animals were previously obtained. ROH and signature of selection analyses were performed using PLINK software, followed by ROH-based (FROH) and pedigree-based inbreeding (Fped) and a search for genes and their functions. An average of 50 ± 8.59 ROHs were found per animal. ROHs were separated into classes according to size, ranging from 1 to 2 Mb (ROH1-2Mb: 58.17%), representing ancient inbreeding, ROH2-4Mb (22.74%), ROH4-8Mb (11.34%), ROH8-16Mb (5.51%), and ROH>16Mb (2.24%). Combining our results, we conclude that the increase in general FROH and Fped significantly decreases TO and VO; however, in different chromosomes traits can increase or decrease with FROH. In the analysis for signatures of selection, we identified 15 genes from 47 significant genomic regions, indicating differences in populations with high and low breeding value for the three traits.

Pedigree reconstruction and population structure using SNP markers in Gir cattle.

Our objective was to establish a SNPs panel for pedigree reconstruction using microarrays of different densities and evaluate the genomic relationship coefficient of the inferred pedigree, in addition to analyzing the population structure based on genomic analyses in Gir cattle. For parentage analysis and genomic relationship, 16,205 genotyped Gir animals (14,458 females and 1747 males) and 1810 common markers to the four SNP microarrays were used. For population structure analyses, including linkage disequilibrium, effective population size, and runs of homozygosity (ROH), genotypes from 21,656 animals were imputed. Likelihood ratio (LR) approach was used to reconstruct the pedigree, deepening the pedigree and showing it is well established in terms of recent information. Coefficients for each relationship category of the inferred pedigree were adequate. Linkage disequilibrium showed rapid decay. We detected a decrease in the effective population size over the last 50 generations, with the average generation interval around 9.08 years. Higher ROH-based inbreeding coefficient in a class of short ROH segments, with moderate to high values, was also detected, suggesting bottlenecks in the Gir genome. Breeding strategies to minimize inbreeding and avoid massive use of few proven sires with high genetic value are suggested to maintain genetic variability in future generations. In addition, we recommend reducing the generation interval to maximize genetic progress and increase effective population size.

Fine Mapping of a Major Backfat QTL Reveals a Causal Regulatory Variant Affecting the CCND2 Gene.

Backfat is an important trait in pork production, and it has been included in the breeding objectives of genetic companies for decades. Although adipose tissue is a good energy storage, excessive fat results in reduced efficiency and economical losses. A large QTL for backfat thickness on chromosome 5 is still segregating in different commercial pig breeds. We fine mapped this QTL region using a genome-wide association analysis (GWAS) with 133,358 genotyped animals from five commercial populations (Landrace, Pietrain, Large White, Synthetic, and Duroc) imputed to the porcine 660K SNP chip. The lead SNP was located at 5:66103958 (G/A) within the third intron of the CCND2 gene, with the G allele associated with more backfat, while the A allele is associated with less backfat. We further phased the QTL region to discover a core haplotype of five SNPs associated with low backfat across three breeds. Linkage disequilibrium analysis using whole-genome sequence data revealed three candidate causal variants within intronic regions and downstream of the CCND2 gene, including the lead SNP. We evaluated the association of the lead SNP with the expression of the genes in the QTL region (including CCND2) in a large cohort of 100 crossbred samples, sequenced in four different tissues (lung, spleen, liver, muscle). Results show that the A allele increases the expression of CCND2 in an additive way in three out of four tissues. Our findings indicate that the causal variant for this QTL region is a regulatory variant within the third intron of the CCND2 gene affecting the expression of CCND2.

Genome-Wide Analyses Reveal the Genetic Architecture and Candidate Genes of Indicine, Taurine, Synthetic Crossbreds, and Locally Adapted Cattle in Brazil.

Cattle population history, breeding systems, and geographic subdivision may be reflected in runs of homozygosity (ROH), effective population size (N e), and linkage disequilibrium (LD) patterns. Thus, the assessment of this information has become essential to the implementation of genomic selection on purebred and crossbred cattle breeding programs. In this way, we assessed the genotype of 19 cattle breeds raised in Brazil belonging to taurine, indicine, synthetic crossbreds, and Iberian-derived locally adapted ancestries to evaluate the overall LD decay patterns, N e, ROH, and breed composition. We were able to obtain a general overview of the genomic architecture of cattle breeds currently raised in Brazil and other tropical countries. We found that, among the evaluated breeds, different marker densities should be used to improve the genomic prediction accuracy and power of genome-wide association studies. Breeds showing low N e values indicate a recent inbreeding, also reflected by the occurrence of longer ROH, which demand special attention in the matting schemes to avoid extensive inbreeding. Candidate genes (e.g., ABCA7, PENK, SPP1, IFNAR1, IFNAR2, SPEF2, PRLR, LRRTM1, and LRRTM4) located in the identified ROH islands were evaluated, highlighting biological processes involved with milk production, behavior, rusticity, and fertility. Furthermore, we were successful in obtaining the breed composition regarding the taurine and indicine composition using single-nucleotide polymorphism (SNP) data. Our results were able to observe in detail the genomic backgrounds that are present in each breed and allowed to better understand the various contributions of ancestor breeds to the modern breed composition to the Brazilian cattle.