Search

Whole-genome resequencing, alignment, and annotation analyses were undertaken for ten sires representing Gyr, Girolando, and Holstein cattle breeds to detect and make publicly available genome-wide single nucleotide variations (SNVs) and insertions/deletions (InDels). A total of approximately 3.5 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. Sequencing reads were mapped to the UMD 3.1 bovine genome assembly using the Burrows-Wheeler Aligner tool (v. 0.7.10-r789) with default parameters. Picard tools (v. 1.54) were used to eliminate PCR duplicates, and the variant calling was conducted with FreeBayes. The resulting variant lists obtained for each animal were filtered by vcffilter in order to remove variants with quality scores lower than 30 or coverage lower than 7. A total of 23,743,970 SNVs and 3,171,711 InDels were detected in the samples, of which approximately 2,224,013 SNVs and 798,8910 [Editor: Please check number.] InDels were novel. The high number of genetic variants identified for each animal within the same breed shows the importance of resequencing to identify novel variants for monitoring genetic diversity in the cattle breeds and for developing strategies to prevent some eventual loss of genetic variability. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous nonsynonymous SNVs and frameshift InDels, which could affect phenotypic variation. Functional enrichment analysis was performed by the David (v. 6.8) tool and revealed that variants in the olfactory transduction pathway were overrepresented in all three cattle breeds, while the ECM-receptor interaction pathway was overrepresented only in the Girolando breed, the ABC transporters pathway was overrepresented only in the Holstein breed, and the metabolic pathways were overrepresented only in the Gyr breed. Therefore, the genetic variants discovered in this study provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, and Holstein breeding programs. [ABSTRACT FROM AUTHOR]

The aim of this study was to detect, by whole-genome resequencing, SNVs (single nucleotide variants) and indels (insertion/deletions) and annotate them into functional categories in order to identify functionally relevant variants in three important cattle breeds in Brazil: Gyr and Girolando (dairy production) and Guzerat (dual-purpose). A total of approximately 2.7 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.9- to 16.4-fold genome coverage. Sequencing reads were mapped to the UMD 3.1 bovine genome assembly using the Burrows-Wheeler Aligner tool (v.0.7.10-r789). Picard tools (v.1.54) were used to eliminate PCR duplicates, and the variant calling was conducted with FreeBayes. The resulting variant lists were filtered by vcffilter in order to remove SNVs and indels with quality scores lower than 30 or coverage lower than 7. A total of 25,020,024 SNVs (16,743,392 from Guzerat, 15,941,804 from Gyr, and 13,286,669 from Girolando) and 3,249,148 indels (1,975,563 from Guzerat, 1,833,387 from Gyr, and 1,413,047 from Girolando) were detected in the DNA samples. A total of 13,253, 13,817, and 12,480 genes showed genetic variation in Guzerat, Gyr, and Girolando, respectively, of which 64.84% of genes with SNVs and 31.29% of genes with indels were shared among all breeds. The functional enrichment analysis by the DAVID (v.6.8) tool revealed 24, 27, and 28 enriched KEGG pathways (FDR < 10%) in Guzerat, Gyr, and Girolando, respectively, of which 14 pathways were common to all breeds, three pathways were significantly overrepresented only in Guzerat (arachidonic acid metabolism, Fc gamma R-mediated phagocytosis, and aldosterone-regulated sodium reabsorption), six pathways only in Gyr (alanine, aspartate and glutamate metabolism, inflammatory mediator regulation of TRP channels, thyroid hormone synthesis, pancreatic secretion, central carbon metabolism in cancer, and choline metabolism in cancer), and eight pathways only in Girolando (amino sugar and nucleotide sugar metabolism, inositol phosphate metabolism, vascular smooth muscle contraction, tight junction, regulation of actin cytoskeleton, amoebiasis, small cell lung cancer, and dilated cardiomyopathy). Although the genetic variants identified were distributed throughout the genomes, a large number of novel variants were clustered in specific genes. A total of 61 genes were identified with novel variants common in all breeds, while 349, 404, and 206 genes were identified as enriched with novel variants exclusively in the Guzerat, Gyr, and Girolando breeds, respectively. These genes are related to many biological processes, providing valuable information about genomic variants that may be responsible for variation in economically important traits among these breeds. [ABSTRACT FROM AUTHOR]

Genomic selection may accelerate genetic progress in breeding programs of indicine breeds when compared with traditional selection methods. We present results of genomic predictions in Gyr (Bos indicus) dairy cattle of Brazil for milk yield (MY), fat yield (FY), protein yield (PY), and age at first calving using information from bulls and cows. Four different single nucleotide polymorphism (SNP) chips were studied. Additionally, the effect of the use of imputed data on genomic prediction accuracy was studied. A total of 474 bulls and 1,688 cows were genotyped with the Illumina BovineHD (HD; San Diego, CA) and BovineSNP50 (50K) chip, respectively. Genotypes of cows were imputed to HD using FImpute v2.2. After quality check of data, 496,606 markers remained. The HD markers present on the GeneSeek SGGP-20Ki (15,727; Lincoln, NE), 50K (22,152), and GeneSeek GGP-75Ki (65,018) were subset and used to assess the effect of lower SNP density on accuracy of prediction. Deregressed breeding values were used as pseudophenotypes for model training. Data were split into reference and validation to mimic a forward prediction scheme. The reference population consisted of animals whose birth year was ≤2004 and consisted of either only bulls (TR1) or a combination of bulls and dams (TR2), whereas the validation set consisted of younger bulls (born after 2004). Genomic BLUP was used to estimate genomic breeding values (GEBV) and reliability of GEBV (R²PEV) was based on the prediction error variance approach. Reliability of GEBV ranged from ~0.46 (FY and PY) to 0.56 (MY) with TR1 and from 0.51 (PY) to 0.65 (My) with TR2. When averaged across all traits, R²PEV were substantially higher (R²PEV of TR1 = 0.50 and TR2 = 0.57) compared with reliabilities of parent averages (0.35) computed from pedigree data and based on diagonals of the coefficient matrix (prediction error variance approach). Reliability was similar for all the 4 marker panels using either TR1 or TR2, except that imputed HD cow data set led to an inflation of reliability. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information. A reduced panel of ~15K markers resulted in reliabilities similar to using HD markers. Reliability of GEBV could be increased by enlarging the limited bull reference population with cow information. [ABSTRACT FROM AUTHOR]

Genotype imputation is widely used as a cost-effective strategy in genomic evaluation of cattle. Key determinants of imputation accuracies, such as linkage disequilibrium patterns, marker densities, and ascertainment bias, differ between Bos indicus and Bos taurus breeds. Consequently, there is a need to investigate effectiveness of genotype imputation in indicine breeds. Thus, the objective of the study was to investigate strategies and factors affecting the accuracy of genotype imputation in Gyr (Bos indicus) dairy cattle. Four imputation scenarios were studied using 471 sires and 1,644 dams genotyped on Illumina BovineHD (HD-777K; San Diego, CA) and BovineSNP50 (50K) chips, respectively. Scenarios were based on which reference high-density single nucleotide polymorphism (SNP) panel (HDP) should be adopted [HD-777K, 50K, and GeneSeek GGP-75Ki (Lincoln, NE)]. Depending on the scenario, validation animals had their genotypes masked for one of the lower-density panels: Illumina (3K, 7K, and 50K) and GeneSeek (SGGP-20Ki and GGP-75Ki). We randomly selected 171 sires as reference and 300 as validation for all the scenarios. Additionally, all sires were used as reference and the 1,644 dams were imputed for validation. Genotypes of 98 individuals with 4 and more offspring were completely masked and imputed. Imputation algorithms FImpute and Beagle v3.3 and v4 were used. Imputation accuracies were measured using the correlation and allelic correct rate. FImpute resulted in highest accuracies, whereas Beagle 3.3 gave the least-accurate imputations. Accuracies evaluated as correlation (allelic correct rate) ranged from 0.910 (0.942) to 0.961 (0.974) using 50K as HDP and with 3K (7K) as low-density panels. With GGP-75Ki as HDP, accuracies were moderate for 3K, 7K, and 50K, but high for SGGP-20Ki. The use of HD-777K as HDP resulted in accuracies of 0.888 (3K), 0.941 (7K), 0.980 (SGGP-20Ki), 0.982 (50K), and 0.993 (GGP-75Ki). Ungenotyped individuals were imputed with an average accuracy of 0.970. The average top 5 kinship coefficients between reference and imputed individuals was a strong predictor of imputation accuracy. FImpute was faster and used less memory than Beagle v4. Beagle v4 outperformed Beagle v3.3 in accuracy and speed of computation. A genotyping strategy that uses the HD-777K SNP chip as a reference panel and SGGP-20Ki as the lower-density SNP panel should be adopted as accuracy was high and similar to that of the 50K. However, the effect of using imputed HD-777K genotypes from the SGGP-20Ki on genomic evaluation is yet to be studied. [ABSTRACT FROM AUTHOR]

The aim of this study was to elucidate the differential gene expression in the RNA sequencing transcriptome of isolated perfused udders collected from 4 slaughtered Holstein × Zebu crossbred dairy cows experimentally inoculated with Streptococcus agalactiae. We studied 3 different statistical tools (edgeR, baySeq, and Cuffdiff 2). In summary, 2 quarters of each udder were experimentally inoculated with Strep. agalactiae and the other 2 were used as a control. Mammary tissue biopsies were collected at times 0 and 3 h after infection. The total RNA was extracted and sequenced on an Illumina HiSeq 2000 (Illumina Inc., San Diego, CA). Transcripts were assembled from the reads aligned to the bovine UMD 3.1 reference genome, and the statistical analyses were performed using the previously mentioned tools (edgeR, baySeq, and Cuffdiff 2). Finally, the identified genes were submitted to pathway enrichment analysis. A total of 1,756, 1,161, and 3,389 genes with differential gene expression were identified when using edgeR, baySeq, and Cuffdiff 2, respectively. A total of 122 genes were identified by the overlapping of the 3 methods; however, only the platelet activation presented a significantly enriched pathway. From the results, we suggest the FCER1G, GNAI2, ORAI1, and VASP genes shared among the 3 methods in this pathway for posterior biological validation., (Copyright © 2019 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Fourteen Brazilian dairy Gyr sire families with 657 daughters were analyzed for quantitative trait loci (QTL) on chromosome 6 by using a daughter design for 5 economic traits: milk, fat, and protein production, fat and protein percentage. The cows and sires were genotyped for 27 microsatellites with average spacing between markers of 4.9 cM. In the analyses across 14 families, for the largest significant families, and within family, a QTL was located for milk yield and fat yield close to marker BMS2508 at the 5% chromosome-wide significance level across families and 1% chromosome-wide within families. For fat percentage, a QTL near DIK4482 was identified at the 5% chromosome-wide significance level when all families were analyzed together and at the 1% chromosome-wide significance level within the largest significant families. The different analyses yielded results that were generally consistent for milk yield, fat yield, and fat percentage. The order of the markers in the derived map was consistent with that in the consensus map. Some QTL and candidate genes in dairy cattle for milk production traits are probably preserved in Bos taurus and Bos indicus., (Copyright © 2011 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)