Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data.

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects.

BACKGROUND: Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs. RESULTS: The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings. CONCLUSIONS: Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Runs of homozygosity and cross-generational inbreeding of Iranian fat-tailed sheep.

The Lori-Bakhtiari fat-tailed sheep is one of the most important heavyweight native breeds of Iran. The breed is robust and well-adapted to semi-arid regions and an important resource for smallholder farms. An established nucleus-based breeding scheme is used to improve their production traits but there is an indication of inbreeding depression and loss of genetic diversity due to selection. Here, we estimated the inbreeding levels and the distribution of runs of homozygosity (ROH) islands in 122 multi-generational female Lori-Bakhtiari from different half-sib families selected from a breeding station that were genotyped on the 50k array. A total of 2404 ROH islands were identified. On average, there were 19.70 ± 1.4 ROH per individual ranging between 6 and 41. The mean length of the ROH was 4.1 ± 0.14 Mb. There were 1999 short ROH of length 1-6 Mb and another 300 in the 6-12 Mb range. Additionally long ROH indicative of inbreeding were found in the ranges of 12-24 Mb (95) and 24-48 Mb (10). The average inbreeding coefficient (FROH) was 0.031 ± 0.003 with estimates varying from 0.006 to 0.083. Across generations, FROH increased from 0.019 ± 0.012 to 0.036 ± 0.007. Signatures of selection were identified on chromosomes 2, 6, and 10, encompassing 55 genes and 23 QTL associated with production traits. Inbreeding coefficients are currently within acceptable levels but across generations, inbreeding is increasing due to selection. The breeding program needs to actively monitor future inbreeding rates and ensure that the breed maintains or improves on its current levels of environmental adaptation.

Genome-wide association studies for epistatic genetic effects on fertility and reproduction traits in Holstein cattle.

Non-additive genetic effects are well known to play an important role in the phenotypic expression of complex traits, such as fertility and reproduction. In this study, a genome scan was performed using 41,640 single nucleotide polymorphism (SNP) markers to identify genomic regions associated with epistatic (additive-by-additive) effects in fertility and reproduction traits in Holstein cattle. Nine fertility and reproduction traits were analysed on 5825 and 6090 Holstein heifers and cows with phenotypes and genotypes, respectively. The Marginal Epistasis Test (MAPIT) was used to identify SNPs with significant marginal epistatic effects at a chromosome-wise 5% and 10% false discovery rate (FDR) level. The -log10 (p) values were adjusted by the genomic inflation factor (λ) to correct for the potential bias on the p-values and minimize the possible effects of population stratification. After adjustments, MAPIT enabled the identification of genomic regions with significant marginal epistatic effects for heifers on BTA5 for age at first insemination, BTA3 and BTA24 for non-return rate (NRR); BTA16 and BTA28 for gestation length (GL); BTA1, BTA4 and BTA17 for stillbirth (SB). For the cow traits, MAPIT enabled the identification of regions on BTA11 for GL, BTA11 and BTA16 for SB and BTA19 for calf size (CZ). An additional approach for mapping epistasis in a genome-wide association study was also proposed, in which the genome scan was performed using estimates of epistatic values as the input pseudo-phenotypes, computed using single-trait animal models. Significant SNPs were identified at the chromosome-wise 5% and 10% FDR levels for all traits. For the heifer traits, significant regions were found on BTA7 for AFS; BTA12 for NRR; BTA14 and BTA19 for GL; BTA19 for calving ease (CE); BTA5, BTA24, BTA25 and in the X chromosome for SB; BTA23 and in the X chromosome for CZ and in the X chromosome for the number of services (NS). For the cow traits, significant regions were found on BTA29 and in the X chromosome for NRR, BTA11, BTA16 and in the X chromosome for SB, BTA2 for GL, BTA28 for CZ, BTA19 for calving to first insemination, and in the X chromosome for NS and first insemination to conception. The results suggest that the epistatic genetic effects are likely due to many loci with a small effect rather than few loci with a large effect and/or a single SNP marker alone do not capture the epistatic effects well. The genomic architecture of fertility and reproduction traits is complex, and these results should be validated in independent dairy cattle populations and using alternative statistical models.

Genome-wide detection of copy number variation in American mink using whole-genome sequencing.

BACKGROUND: Copy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta. RESULTS: A total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3). CONCLUSIONS: This study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.

Imputation accuracy to whole-genome sequence in Nellore cattle.

BACKGROUND: A cost-effective strategy to explore the complete DNA sequence in animals for genetic evaluation purposes is to sequence key ancestors of a population, followed by imputation mechanisms to infer marker genotypes that were not originally reported in a target population of animals genotyped with single nucleotide polymorphism (SNP) panels. The feasibility of this process relies on the accuracy of the genotype imputation in that population, particularly for potential causal mutations which may be at low frequency and either within genes or regulatory regions. The objective of the present study was to investigate the imputation accuracy to the sequence level in a Nellore beef cattle population, including that for variants in annotation classes which are more likely to be functional. METHODS: Information of 151 key sequenced Nellore sires were used to assess the imputation accuracy from bovine HD BeadChip SNP (~ 777 k) to whole-genome sequence. The choice of the sires aimed at optimizing the imputation accuracy of a genotypic database, comprised of about 10,000 genotyped Nellore animals. Genotype imputation was performed using two computational approaches: FImpute3 and Minimac4 (after using Eagle for phasing). The accuracy of the imputation was evaluated using a fivefold cross-validation scheme and measured by the squared correlation between observed and imputed genotypes, calculated by individual and by SNP. SNPs were classified into a range of annotations, and the accuracy of imputation within each annotation classification was also evaluated. RESULTS: High average imputation accuracies per animal were achieved using both FImpute3 (0.94) and Minimac4 (0.95). On average, common variants (minor allele frequency (MAF) > 0.03) were more accurately imputed by Minimac4 and low-frequency variants (MAF ≤ 0.03) were more accurately imputed by FImpute3. The inherent Minimac4 Rsq imputation quality statistic appears to be a good indicator of the empirical Minimac4 imputation accuracy. Both software provided high average SNP-wise imputation accuracy for all classes of biological annotations. CONCLUSIONS: Our results indicate that imputation to whole-genome sequence is feasible in Nellore beef cattle since high imputation accuracies per individual are expected. SNP-wise imputation accuracy is software-dependent, especially for rare variants. The accuracy of imputation appears to be relatively independent of annotation classification.

Signatures of selection analysis using whole-genome sequence data reveals novel candidate genes for pony and light horse types.

Natural selection and domestication have shaped modern horse populations, resulting in a vast range of phenotypically diverse breeds. Horse breeds are classified into three types (pony, light, and draft) generally based on their body type. Understanding the genetic basis of horse type variation and selective pressures related to the evolutionary trend can be particularly important for current selection strategies. Whole-genome sequences were generated for 14 pony and 32 light horses to investigate the genetic signatures of selection of the horse type in pony and light horses. In the overlapping extremes of the fixation index and nucleotide diversity results, we found novel genomic signatures of selective sweeps near key genes previously implicated in body measurements including C4ORF33, CRB1, CPN1, FAM13A, and FGF12 that may influence variation in pony and light horse types. This study contributes to a better understanding of the genetic background of differences between pony and light horse types.

Estimation of genetic parameters for mid-infrared-predicted lactoferrin and milk fat globule size in Holstein cattle.

Lactoferrin (LF) and milk fat globule (MFG) are 2 biologically active components of milk with great economical and nutritional value in the dairy industry. The objectives of this study were to estimate (1) the heritability of mid-infrared (MIR)-predicted LF and MFG size (MFGS) and (2) the genetic correlations between predicted LF and MFGS with milk, fat, and protein yields, fat and protein percentages, and somatic cell score in first-parity Canadian Holstein cattle. A total of 109,029 test-day records from 22,432 cows and 1,572 farms for MIR-predicted LF and 109,212 test-day records from 22,424 cows and 1,559 farms for MIR-predicted MFGS were used in the analyses. Four separate 5-trait random regression test-day models were used. The models included days in milk, herd test date, and a polynomial regression on DIM nested in age-season of calving classes as fixed effects, random polynomial regressions on DIM nested in herd-year of calving, animal additive genetic and permanent environment classes, and a residual effect. Regression curves were modeled using orthogonal Legendre polynomials of order 4 for the fixed age-season of calving effect and of order 5 for the random effects. Moderate overall heritability estimates of 0.34 and 0.46 were estimated for the MIR-predicted LF and MIR-predicted MFGS, respectively. These heritability estimates were similar to the ones estimated for the direct measure of MFGS in a previous study. The genetic correlations between predicted MFGS and fat percentage (0.53) and between predicted LF and protein percentage (0.41) were both moderate and positive. Predicted LF and somatic cell score showed a weaker correlation (0.06) compared with other studies. The moderate genetic correlation between MIR-predicted MFGS and fat percentage and between MIR-predicted LF and protein percentage suggests that MIR predictions of MFGS and LF are not simply a function of the amount of fat and protein percentage, respectively, in the milk (i.e., the prediction equations are not simply predicting fat or protein percentages). Thus, these MIR-predicted values may provide additional information for selecting for fine milk components in Holstein cattle.

Estimation of additive and non-additive genetic effects for fertility and reproduction traits in North American Holstein cattle using genomic information.

Non-additive genetic effects are usually ignored in animal breeding programs due to data structure (e.g., incomplete pedigree), computational limitations and over-parameterization of the models. However, non-additive genetic effects may play an important role in the expression of complex traits in livestock species, such as fertility and reproduction traits. In this study, components of genetic variance for additive and non-additive genetic effects were estimated for a variety of fertility and reproduction traits in Holstein cattle using pedigree and genomic relationship matrices. Four linear models were used: (a) an additive genetic model; (b) a model including both additive and epistatic (additive by additive) genetic effects; (c) a model including both additive and dominance effects; and (d) a full model including additive, epistatic and dominance genetic effects. Nine fertility and reproduction traits were analysed, and models were run separately for heifers (N = 5,825) and cows (N = 6,090). For some traits, a larger proportion of phenotypic variance was explained by non-additive genetic effects compared with additive effects, indicating that epistasis, dominance or a combination thereof is of great importance. Epistatic genetic effects contributed more to the total phenotypic variance than dominance genetic effects. Although these models varied considerably in the partitioning of the components of genetic variance, the models including a non-additive genetic effect did not show a clear advantage over the additive model based on the Akaike information criterion. The partitioning of variance components resulted in a re-ranking of cows based solely on the cows' additive genetic effects between models, indicating that adjusting for non-additive genetic effects could affect selection decisions made in dairy cattle breeding programs. These results suggest that non-additive genetic effects play an important role in some fertility and reproduction traits in Holstein cattle.

The effect of host genetics on in vitro performance of bovine monocyte-derived macrophages.

The dynamic interaction between the host and pathogens, along with environmental factors, influences the regulation of mammalian immune responses. Therefore, comprehensive in vivo immune-phenotyping during an active response to a pathogen can be complex and prone to confounding effects. Evaluating critical fundamental aspects of the immune system at a cellular level is an alternative approach to reduce this complexity. Therefore, the objective of the current study was to examine an in vitro model for functional phenotyping of bovine monocyte-derived macrophages (MDM), cells which play a crucial role at all phases of inflammation, as well influence downstream immune responses. As indicators of MDM function, phagocytosis and nitric oxide (NO-) production were tested in MDM of 16 cows in response to 2 common bacterial pathogens of dairy cows, Escherichia coli and Staphylococcus aureus. Notable functional variations were observed among the individuals (coefficient of variation: 33% for phagocytosis and 70% in the production of NO-). The rank correlation analysis revealed a significant, positive, and strong correlation (rho = 0.92) between NO- production in response to E. coli and S. aureus, and a positive but moderate correlation (rho = 0.58) between phagocytosis of E. coli and S. aureus. To gain further insight into this trait, another 58 cows were evaluated solely for NO- response against E. coli. The pedigree of the tested animals was added to the statistical model and the heritability was estimated to be 0.776. Overall, the finding of this study showed a strong effect of host genetics on the in vitro activities of MDM and the possibility of ranking Holstein cows based on the in vitro functional variation of MDM.

Impact of including information from bulls and their daughters in the training population of multiple-step genomic evaluations in dairy cattle: A simulation study.

The objective of this study was to investigate the impact of accounting for parent average (PA) and genotyped daughters' average (GDA) on the estimation of deregressed estimated breeding values (dEBVs) used as pseudo-phenotypes in multiple-step genomic evaluations. Genomic estimated breeding values (GEBVs) were predicted, in eight different simulated scenarios, using dEBVs calculated based on four methods. These methods included PA and GDA in the dEBV (VR) or only GDA (VRpa) and excluded both PA and GDA from the dEBV with either all information or only information from PA and GDA (JA and NEW, respectively). In general, VR and NEW showed the lowest and highest GEBV reliabilities across scenarios, respectively. Among all deregression methods, VRpa and NEW provided the most consistent bias estimates across the majority of scenarios, and they significantly yielded the least biased GEBVs. Our results indicate that removing PA and GDA information from dEBVs used in multiple-step genomic evaluations can increase the reliability of GEBVs, when both bulls and their daughters are included in the training population.

Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle.

BACKGROUND: While autozygosity as a consequence of selection is well understood, there is limited information on the ability of different methods to measure true inbreeding. In the present study, a gene dropping simulation was performed and inbreeding estimates based on runs of homozygosity (ROH), pedigree, and the genomic relationship matrix were compared to true inbreeding. Inbreeding based on ROH was estimated using SNP1101, PLINK, and BCFtools software with different threshold parameters. The effects of different selection methods on ROH patterns were also compared. Furthermore, inbreeding coefficients were estimated in a sample of genotyped North American Holstein animals born from 1990 to 2016 using 50 k chip data and ROH patterns were assessed before and after genomic selection. RESULTS: Using ROH with a minimum window size of 20 to 50 using SNP1101 provided the closest estimates to true inbreeding in simulation study. Pedigree inbreeding tended to underestimate true inbreeding, and results for genomic inbreeding varied depending on assumptions about base allele frequencies. Using an ROH approach also made it possible to assess the effect of population structure and selection on distribution of runs of autozygosity across the genome. In the simulation, the longest individual ROH and the largest average length of ROH were observed when selection was based on best linear unbiased prediction (BLUP), whereas genomic selection showed the largest number of small ROH compared to BLUP estimated breeding values (BLUP-EBV). In North American Holsteins, the average number of ROH segments of 1 Mb or more per individual increased from 57 in 1990 to 82 in 2016. The rate of increase in the last 5 years was almost double that of previous 5 year periods. Genomic selection results in less autozygosity per generation, but more per year given the reduced generation interval. CONCLUSIONS: This study shows that existing software based on the measurement of ROH can accurately identify autozygosity across the genome, provided appropriate threshold parameters are used. Our results show how different selection strategies affect the distribution of ROH, and how the distribution of ROH has changed in the North American dairy cattle population over the last 25 years.

A genome-wide association study for natural antibodies measured in blood of Canadian Holstein cows.

BACKGROUND: Natural antibodies (NAb) are an important component of the innate immune system, and fight infections as a part of the first line defence. NAb are poly-reactive and can respond non-specifically to antigens. Therefore, NAb may be a key trait when evaluating an animal's potential natural disease resistance. Variation in NAb is caused by both genetic and environmental factors. In this study genetic parameters of NAb were estimated and a genome-wide association study (GWAS) was performed to gain further understanding on the genes that are responsible for the observed genetic variation of NAb in Canadian Holsteins. RESULTS: In total, blood samples of 1327 cows from 64 farms were studied. NAb binding to keyhole limpet hemocyanin (KLH) were determined via indirect ELISA. Immunoglobulin (Ig) isotypes, IgG and IgM, were evaluated. From the sample population, 925 cows were genotyped for 45,187 markers and each individual marker was tested to detect genetic variation in NAb levels. The relationships among animals was accounted for with genomic relationship. Results show heritabilities of 0.27 ± 0.064 (IgG) and 0.31 ± 0.065 (IgM). In total, 23 SNPs were found to be associated with IgG, but no SNPs were associated with IgM (FDR p-value < 0.05). The significant SNPs were located on autosomal chromosomes 1, 20 and 21 of the cow genome. Functional annotation analysis of the positional candidate genes revealed two sets of genes with biologically relevant functions related to NAb. In one set, seven genes with crucial roles in the production of antibody in B cells were associated with the trafficking of vesicles inside the cells between organelles. In the second set, two genes among positional candidate genes were associated with isotype class-switching and somatic hypermutation of B cells. CONCLUSIONS: This study demonstrated the possibility of increasing NAb through selective breeding. In addition, the effects of two candidate pathways are proposed for further investigation of NAb production in Holsteins.

Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers.

BACKGROUND: The detection of signatures of selection has the potential to elucidate the identities of genes and mutations associated with phenotypic traits important for livestock species. It is also very relevant to investigate the levels of genetic diversity of a population, as genetic diversity represents the raw material essential for breeding and has practical implications for implementation of genomic selection. A total of 1151 animals from nine goat populations selected for different breeding goals and genotyped with the Illumina Goat 50K single nucleotide polymorphisms (SNP) Beadchip were included in this investigation. RESULTS: The proportion of polymorphic SNPs ranged from 0.902 (Nubian) to 0.995 (Rangeland). The overall mean HO and HE was 0.374 ± 0.021 and 0.369 ± 0.023, respectively. The average pairwise genetic distance (D) ranged from 0.263 (Toggenburg) to 0.323 (Rangeland). The overall average for the inbreeding measures FEH, FVR, FLEUT, FROH and FPED was 0.129, -0.012, -0.010, 0.038 and 0.030, respectively. Several regions located on 19 chromosomes were potentially under selection in at least one of the goat breeds. The genomic population tree constructed using all SNPs differentiated breeds based on selection purpose, while genomic population tree built using only SNPs in the most significant region showed a great differentiation between LaMancha and the other breeds. We hypothesized that this region is related to ear morphogenesis. Furthermore, we identified genes potentially related to reproduction traits, adult body mass, efficiency of food conversion, abdominal fat deposition, conformation traits, liver fat metabolism, milk fatty acids, somatic cells score, milk protein, thermo-tolerance and ear morphogenesis. CONCLUSIONS: In general, moderate to high levels of genetic variability were observed for all the breeds and a characterization of runs of homozygosity gave insights into the breeds' development history. The information reported here will be useful for the implementation of genomic selection and other genomic studies in goats. We also identified various genome regions under positive selection using smoothed FST and hapFLK statistics and suggested genes, which are potentially under selection. These results can now provide a foundation to formulate biological hypotheses related to selection processes in goats.

Novel methods for genotype imputation to whole-genome sequence and a simple linear model to predict imputation accuracy.

BACKGROUND: Accurate imputation plays a major role in genomic studies of livestock industries, where the number of genotyped or sequenced animals is limited by costs. This study explored methods to create an ideal reference population for imputation to Next Generation Sequencing data in cattle. METHODS: Methods for clustering of animals for imputation were explored, using 1000 Bull Genomes Project sequence data on 1146 animals from a variety of beef and dairy breeds. Imputation from 50 K to 777 K was first carried out to choose an ideal clustering method, using ADMIXTURE or PLINK clustering algorithms with either genotypes or reconstructed haplotypes. RESULTS: Due to efficiency, accuracy and ease of use, clustering with PLINK using haplotypes as quasi-genotypes was chosen as the most advantageous grouping method. It was found that using a clustered population slightly decreased computing time, while maintaining accuracy across the population. Although overall accuracy remained the same, a slight increase in accuracy was observed for groups of animals in some breeds (primarily purebred beef cattle from breeds with fewer sequenced animals) and for other groups, primarily crossbreed animals, a slight decrease in accuracy was observed. However, it was noted that some animals in each breed were poorly imputed across all methods. When imputed sequences were included in the reference population to aid imputation of poorly imputed animals, a small increase in overall accuracy was observed for nearly every individual in the population. Two models were created to predict imputation accuracy, a complete model using all information available including Euclidean distances from genotypes and haplotypes, pedigree information, and clustering groups and a simple model using only breed and an Euclidean distance matrix as predictors. Both models were successful in predicting imputation accuracy, with correlations between predicted and true imputation accuracy as measured by concordance rate of 0.87 and 0.83, respectively. CONCLUSIONS: A clustering methodology can be very useful to subgroup cattle for efficient genotype imputation. In addition, accuracy of genotype imputation from medium to high-density Single Nucleotide Polymorphisms (SNP) chip panels to whole-genome sequence can be predicted well using a simple linear model defined in this study.

Genomic predictions for economically important traits in Brazilian Braford and Hereford beef cattle using true and imputed genotypes.

BACKGROUND: Genomic selection (GS) has played an important role in cattle breeding programs. However, genotyping prices are still a challenge for implementation of GS in beef cattle and there is still a lack of information about the use of low-density Single Nucleotide Polymorphisms (SNP) chip panels for genomic predictions in breeds such as Brazilian Braford and Hereford. Therefore, this study investigated the effect of using imputed genotypes in the accuracy of genomic predictions for twenty economically important traits in Brazilian Braford and Hereford beef cattle. Various scenarios composed by different percentages of animals with imputed genotypes and different sizes of the training population were compared. De-regressed EBVs (estimated breeding values) were used as pseudo-phenotypes in a Genomic Best Linear Unbiased Prediction (GBLUP) model using two different mimicked panels derived from the 50 K (8 K and 15 K SNP panels), which were subsequently imputed to the 50 K panel. In addition, genomic prediction accuracies generated from a 777 K SNP (imputed from the 50 K SNP) were presented as another alternate scenario. RESULTS: The accuracy of genomic breeding values averaged over the twenty traits ranged from 0.38 to 0.40 across the different scenarios. The average losses in expected genomic estimated breeding values (GEBV) accuracy (accuracy obtained from the inverse of the mixed model equations) relative to the true 50 K genotypes ranged from -0.0007 to -0.0012 and from -0.0002 to -0.0005 when using the 50 K imputed from the 8 K or 15 K, respectively. When using the imputed 777 K panel the average losses in expected GEBV accuracy was -0.0021. The average gain in expected EBVs accuracy by including genomic information when compared to simple BLUP was between 0.02 and 0.03 across scenarios and traits. CONCLUSIONS: The percentage of animals with imputed genotypes in the training population did not significantly influence the validation accuracy. However, the size of the training population played a major role in the accuracies of genomic predictions in this population. The losses in the expected accuracies of GEBV due to imputation of genotypes were lower when using the 50 K SNP chip panel imputed from the 15 K compared to the one imputed from the 8 K SNP chip panel.

Genetic diversity, extent of linkage disequilibrium and persistence of gametic phase in Canadian pigs.

BACKGROUND: Knowledge on the levels of linkage disequilibrium (LD) across the genome, persistence of gametic phase between breed pairs, genetic diversity and population structure are important parameters for the successful implementation of genomic selection. Therefore, the objectives of this study were to investigate these parameters in order to assess the feasibility of a multi-herd and multi-breed training population for genomic selection in important purebred and crossbred pig populations in Canada. A total of 3,057 animals, representative of the national populations, were genotyped with the Illumina Porcine SNP60 BeadChip (62,163 markers). RESULTS: The overall LD (r 2) between adjacent SNPs was 0.49, 0.38, 0.40 and 0.31 for Duroc, Landrace, Yorkshire and Crossbred (Landrace x Yorkshire) populations, respectively. The highest correlation of phase (r) across breeds was observed between Crossbred animals and either Landrace or Yorkshire breeds, in which r was approximately 0.80 at 1 Mbp of distance. Landrace and Yorkshire breeds presented r ≥ 0.80 in distances up to 0.1 Mbp, while Duroc breed showed r ≥ 0.80 for distances up to 0.03 Mbp with all other populations. The persistence of phase across herds were strong for all breeds, with r ≥ 0.80 up to 1.81 Mbp for Yorkshire, 1.20 Mbp for Duroc, and 0.70 Mbp for Landrace. The first two principal components clearly discriminate all the breeds. Similar levels of genetic diversity were observed among all breed groups. The current effective population size was equal to 75 for Duroc and 92 for both Landrace and Yorkshire. CONCLUSIONS: An overview of population structure, LD decay, demographic history and inbreeding of important pig breeds in Canada was presented. The rate of LD decay for the three Canadian pig breeds indicates that genomic selection can be successfully implemented within breeds with the current 60 K SNP panel. The use of a multi-breed training population involving Landrace and Yorkshire to estimate the genomic breeding values of crossbred animals (Landrace × Yorkshire) should be further evaluated. The lower correlation of phase at short distances between Duroc and the other breeds indicates that a denser panel may be required for the use of a multi-breed training population including Duroc.

Prediction of genomic breeding values for growth, carcass and meat quality traits in a multi-breed sheep population using a HD SNP chip.

BACKGROUND: New Zealand has some unique Terminal Sire composite sheep breeds, which were developed in the last three decades to meet commercial needs. These composite breeds were developed based on crossing various Terminal Sire and Maternal breeds and, therefore, present high genetic diversity compared to other sheep breeds. Their breeding programs are focused on improving carcass and meat quality traits. There is an interest from the industry to implement genomic selection in this population to increase the rates of genetic gain. Therefore, the main objectives of this study were to determine the accuracy of predicted genomic breeding values for various growth, carcass and meat quality traits using a HD SNP chip and to evaluate alternative genomic relationship matrices, validation designs and genomic prediction scenarios. A large multi-breed population (n = 14,845) was genotyped with the HD SNP chip (600 K) and phenotypes were collected for a variety of traits. RESULTS: The average observed accuracies (± SD) for traits measured in the live animal, carcass, and, meat quality traits ranged from 0.18 ± 0.07 to 0.33 ± 0.10, 0.28 ± 0.09 to 0.55 ± 0.05 and 0.21 ± 0.07 to 0.36 ± 0.08, respectively, depending on the scenario/method used in the genomic predictions. When accounting for population stratification by adjusting for 2, 4 or 6 principal components (PCs) the observed accuracies of molecular breeding values (mBVs) decreased or kept constant for all traits. The mBVs observed accuracies when fitting both G and A matrices were similar to fitting only G matrix. The lowest accuracies were observed for k-means cross-validation and forward validation performed within each k-means cluster. CONCLUSIONS: The accuracies observed in this study support the feasibility of genomic selection for growth, carcass and meat quality traits in New Zealand Terminal Sire breeds using the Ovine HD SNP chip. There was a clear advantage on using a mixed training population instead of performing analyzes per genomic clusters. In order to perform genomic predictions per breed group, genotyping more animals is recommended to increase the size of the training population within each group and the genetic relationship between training and validation populations. The different scenarios evaluated in this study will help geneticists and breeders to make wiser decisions in their breeding programs.

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

BACKGROUND: Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value predictions and followed by gene identification and in silico functional analyses. The association results were further validated using five scenarios with different numbers of SNPs. RESULTS: Seven hundred and eighty-two SNPs were significantly associated with calving performance at a genome-wise false discovery rate (FDR) of 5%. Most of these significant SNPs were on chromosomes 18 (71.9%), 17 (7.4%), 5 (6.8%) and 7 (2.4%) and mapped to 675 genes, among which 142 included at least one significant SNP and 532 were nearby one (100 kbp). For body conformation traits, 607 SNPs were significant at a genome-wise FDR of 5% and most of them were located on chromosomes 5 (30%), 18 (27%), 20 (13%), 6 (6%), 7 (5%), 14 (5%) and 13 (3%). SNP enrichment functional analyses for calving traits at a FDR of 1% suggested potential biological processes including musculoskeletal movement, meiotic cell cycle, oocyte maturation and skeletal muscle contraction. Furthermore, pathway analyses suggested potential pathways associated with calving performance traits including tight junction, oxytocin signaling, and MAPK signaling (P < 0.10). The prediction ability of the 1206 significant SNPs was between 78 and 83% of the prediction ability of the BovineSNP50 SNPs for calving performance traits and between 35 and 79% for body conformation traits. CONCLUSIONS: Various SNPs that are significantly associated with calving performance are located within or nearby genes with potential roles in tight junction, oxytocin signaling, and MAPK signaling. Combining the significant SNPs or SNPs within or nearby gene(s) from the HD panel with the BovineSNP50 panel yielded a marginal increase in the accuracy of prediction of genomic estimated breeding values for all traits compared to the use of the BovineSNP50 panel alone.

A comparison of different algorithms for phasing haplotypes using Holstein cattle genotypes and pedigree data.

Phasing genotypes to haplotypes is becoming increasingly important due to its applications in the study of diseases, population and evolutionary genetics, imputation, and so on. Several studies have focused on the development of computational methods that infer haplotype phase from population genotype data. The aim of this study was to compare phasing algorithms implemented in Beagle, Findhap, FImpute, Impute2, and ShapeIt2 software using 50k and 777k (HD) genotyping data. Six scenarios were considered: no-parents, sire-progeny pairs, sire-dam-progeny trios, each with and without pedigree information in Holstein cattle. Algorithms were compared with respect to their phasing accuracy and computational efficiency. In the studied population, Beagle and FImpute were more accurate than other phasing algorithms. Across scenarios, phasing accuracies for Beagle and FImpute were 99.49-99.90% and 99.44-99.99% for 50k, respectively, and 99.90-99.99% and 99.87-99.99% for HD, respectively. Generally, FImpute resulted in higher accuracy when genotypic information of at least one parent was available. In the absence of parental genotypes and pedigree information, Beagle and Impute2 (with double the default number of states) were slightly more accurate than FImpute. Findhap gave high phasing accuracy when parents' genotypes and pedigree information were available. In terms of computing time, Findhap was the fastest algorithm followed by FImpute. FImpute was 30 to 131, 87 to 786, and 353 to 1,400 times faster across scenarios than Beagle, ShapeIt2, and Impute2, respectively. In summary, FImpute and Beagle were the most accurate phasing algorithms. Moreover, the low computational requirement of FImpute makes it an attractive algorithm for phasing genotypes of large livestock populations.