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.

Spatial dynamics of the Girolando breed in Brazil: analysis of genetic integration and environmental factors.

Brazil is one of the world's largest milk producers. Several scientific studies have been developed related to landscape analyses that combine genetic with landscape structure data. In the present study, we aimed to analyze the relationship between genetic, environmental, and socioeconomic aspects of production in Girolando cattle in Brazil, as well as verify the spatial patterns of its genetic diversity. Genetic values and accuracy of 46,289 animals were used as well as information from DNA of 310 Girolando animals. Canonic, discriminant, and cluster analyses were conducted in SAS® and K-means method in ArcGIS 10.3 software. The relationship between genetic and geographic distance was analyzed using different methods in software Alleles in Space®. Clusters with animals with higher genetic values for milk production are located in municipalities with lower gross domestic product, fewer family-based establishments, and lower human development index. These clusters are associated with regions with higher area planted with crops, lower percentage of pastures that were less degraded, higher humidity, lower temperature range, and lower normalized difference vegetation index (NDVI) values. The greater the geographical distance between groups of animals, the greater the genetic distance between them with a significant distinction over 504 km. There is high genetic heterogeneity among animals. From these results, it will be possible to develop methodologies for better evaluation of the animals within the production systems.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data.

BACKGROUND: The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics. RESULTS: A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds. CONCLUSION: The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

Assessment of runs of homozygosity islands and estimates of genomic inbreeding in Gyr (Bos indicus) dairy cattle.

BACKGROUND: Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (FROH), genomic relationship matrix approach (FGRM) and based on the observed versus expected number of homozygous genotypes (FHOM), and from pedigree-based coefficient (FPED). RESULTS: ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH1-2 Mb) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of FROH, FGRM, FHOM, and FPED approaches. FPED estimates ranged from 0.00 to 0.327 and FROH from 0.001 to 0.201. Low to moderate correlations were observed between FPED-FROH and FGRM-FROH, with values ranging from -0.11 to 0.51. Low to high correlations were observed between FROH-FHOM and moderate between FPED-FHOM and FGRM-FHOM. Correlations between FROH from different lengths and FPED gradually increased with ROH length. CONCLUSIONS: Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low FPED-FROH correlations for small segments indicate that FPED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that FROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.