Genome-wide fine-mapping identifies pleiotropic and functional variants that predict many traits across global cattle populations

Functional impact of splicing variants in the elaboration of complex traits in cattle

GWAS conducted directly on imputed whole genome sequence have led to the identification of numerous genetic variants associated with agronomic traits in cattle. However, such variants are often simply markers in linkage disequilibrium with the actual causal variants, which is a limiting factor for the development of accurate genomic predictions. It is possible to identify causal variants by integrating information on how variants impact gene expression into GWAS output. RNA splicing plays a major role in regulating gene expression. Thus, assessing the effect of variants on RNA splicing may explain their function. Here, we use a high-throughput strategy to functionally analyse putative splice-disrupting variants in the bovine genome. Using GWAS, massively parallel reporter assay and deep learning algorithms designed to predict splice-disrupting variants, we identify 38 splice-disrupting variants associated with complex traits in cattle, three of which could be classified as causal. Our results indicate that splice-disrupting variants are widely found in the quantitative trait loci related to these phenotypes. Using our combined approach, we also assess the validity of splicing predictors originally developed to analyse human variants in the context of the bovine genome.

Comparison of genomic prediction accuracies in dairy cattle lactation traits using five classes of functional variants versus generic SNP

BACKGROUND

Genomic selection, typically employing genetic markers from SNP chips, is routine in modern dairy cattle breeding. This study assessed the impact of functional sequence variants on genomic prediction accuracy relative to 50 k SNP chip markers for fat percent, protein percent, milk volume, fat yield, and protein yield in lactating dairy cattle. The functional variants were identified through GWAS, RNA-seq, Histone modification ChIP-seq, ATAC-seq, or were coding variants. The genomic prediction accuracy obtained using each class of functional variants was compared with matched numbers of SNPs randomly selected from the Illumina 50 k SNP chip.

RESULTS

The investigation revealed that variants identified by GWAS or RNA-seq, significantly improved the prediction accuracy across all five traits. Contributions from ChIP-seq, ATAC-seq, and coding variants varied. Some variants identified using ChIP-seq showed marked improvements, while others reduced accuracy in protein yield predictions. Relative to a matched number of 32,595 SNPs from the SNP chip, pooling all the functional variants demonstrated prediction accuracy increases of 1.76% for fat percent, 2.97% for protein percent, 0.51% for milk volume, and 0.26% for fat yield, but with a slight decrease of 0.43% in protein yield.

CONCLUSION

The study demonstrates that functional variants can improve prediction accuracy relative to equivalent numbers of variants from a generic SNP panel, with percent traits showing more significant gains than yield traits. The main advantage of using functional variants for genomic prediction was achievement of comparable accuracy using a smaller, more selective set of loci. This is particularly evident in trait-specific scenarios. Our findings indicate that specific combinations of functional variants comprising 16 k variants can achieve genomic prediction accuracy comparable to employing a standard panel of twice the size (32.6 k), especially for percent traits. This highlights the potential for the development of more efficient, trait-focused SNP panels utilizing functional variants.

Unravelling novel and closely linked association signals for fat-related traits in pigs using prioritised variants from whole-genome sequence data

For most production traits, the largest proportions of genetic variance remain unmapped. Dense whole-genome sequence (WGS) data enable the possibility of discovering novel associations as well as unravelling closely linked association signals with a resolution that marker arrays cannot reach. However, the identification of variants from WGS data that are causal of the variation of complex traits is hindered by the high dimensionality and linkage disequilibrium. Thus, at best, we can narrow the circle around the causal variants to prioritise a set of variants for their posterior validation. In this study, we assessed the utility of WGS data for uncovering associations of weaker effects using, as a model, fat content and composition traits in a Duroc pig population where we previously described major effects of the LEPR and SCD genes. We genotyped 971 pigs for a set of 182 variants from 154 candidate genes that were prioritised from amongst the WGS variants discovered in 205 sequenced individuals. These variants were prioritised conditional to LEPR and SCD. The association of the prioritised variants with the target traits was then tested in the confirmation set of 971 pigs. A total of 17 potentially independent quantitative trait loci (8.4% of the total number of studied genes) were significantly associated (q-value < 0.05) with at least one of the studied traits. We identified novel associations attributable to genes such as ABCC2, MOGAT2, or PLPP1 for backfat thickness, myristic acid content, and monounsaturated fatty acid content, respectively. Our results also revealed a finer granularity of weaker genetic effects in loci such as those around the DGAT2 and FADS2 genes, which may mask the effects of closely located genes like MOGAT2 and DAGLA, respectively. To refine the prioritisation of variants for validation studies, especially when targeting those of weaker effects, we recommend larger and more diverse discovery sets, more precise and complete functional gene annotation, and the integration of other omics data.

Erosion of estimated genomic breeding values with generations is due to long distance associations between markers and QTL

BACKGROUND

Most validation studies of genomic evaluations on candidates (prior to observing phenotypes) present inflation of their predicted breeding values, i.e., regression coefficients of their later observed phenotypes on the early predictions are smaller than one. The aim of this study was to show that this inflation pattern reflects at least partly long-distance associations between markers and quantitative trait loci (QTL) in the reference population and to propose methods to estimate the corresponding "erosion" coefficient.

RESULTS

Across-chromosome linkage disequilibrium (LD) is observed in different dairy cattle breeds, being a result from limited effective population size and from relationships within the reference population. Due to this long distance LD, the estimated SNP effects capture non-zero contributions from distant QTLs, some located on other chromosomes than the SNP itself. Therefore, corresponding SNP effects are partly lost in the next generations and we refer to this loss as "erosion". With the concept of QTL contribution to SNP effects derived from mixed model equations, we show with simulation that this long range LD explains 6-25% of the variance of the estimated genomic breeding values, a proportion that is unchanged when the evaluation model includes a residual polygenic effect. Two methods are proposed to predict this erosion factor assuming known simulated QTL effects. In Method 1, one generation of progeny is simulated from the reference population and the GEBV of these progeny based on SNP effects estimated in this newly simulated generation are regressed on the GEBV of the same progeny based on SNP effects estimated in the reference population. In Method 2 all the QTL contributions to SNP effects are regressed based on SNP-QTL recombination rates and summed to predict the GEBV at the next generation. The regression coefficient of the GEBV based on eroded contributions on the raw GEBV is also an estimate of erosion. An illustration is given with the French Normande female reference bovine population in 2021, showing erosion factors ranging from 0.84 to 0.87.

CONCLUSION

Accounting for erosion is important to avoid inflation and biased predictions. The ways to both reduce inflation and to correct for it in the prediction are discussed.

Evaluation of genomic selection models using whole genome sequence data and functional annotation in Belgian Blue cattle

BACKGROUND

The availability of large cohorts of whole-genome sequenced individuals, combined with functional annotation, is expected to provide opportunities to improve the accuracy of genomic selection (GS). However, such benefits have not often been observed in initial applications. The reference population for GS in Belgian Blue Cattle (BBC) continues to grow. Combined with the availability of reference panels of sequenced individuals, it provides an opportunity to evaluate GS models using whole genome sequence (WGS) data and functional annotation.

RESULTS

Here, we used data from 16,508 cows, with phenotypes for five muscular development traits and imputed at the WGS level, in combination with in silico functional annotation and catalogs of putative regulatory variants obtained from experimental data. We evaluated first GS models using the entire WGS data, with or without functional annotation. At this marker density, we were able to run two approaches, assuming either a highly polygenic architecture (GBLUP) or allowing some variants to have larger effects (BayesRR-RC, a Bayesian mixture model), and observed an increased reliability compared to the official GBLUP model at medium marker density (on average 0.016 and 0.018 for GBLUP and BayesRR-RC, respectively). When functional annotation was used, we observed slightly higher reliabilities with an extension of GBLUP that included multiple polygenic terms (one per functional group), while reliabilities decreased with BayesRR-RC. We then used large subsets of variants selected based on functional information or with a linkage disequilibrium (LD) pruning approach, which allowed us to evaluate two additional approaches, BayesCπ and Bayesian Sparse Linear Mixed Model (BSLMM). Reliabilities were higher for these panels than for the WGS data, with the highest accuracies obtained when markers were selected based on functional information. In our setting, BSLMM systematically achieved higher reliabilities than other methods.

CONCLUSIONS

GS with large panels of functional variants selected from WGS data allowed a significant increase in reliability compared to the official genomic evaluation approach. However, the benefits of using WGS and functional data remained modest, indicating that there is still room for improvement, for example by further refining the functional annotation in the BBC breed.

Mendelian randomisation to uncover causal associations between conformation, metabolism, and production as potential exposure to reproduction in German Holstein dairy cattle

BACKGROUND

Reproduction is vital to welfare, health, and economics in animal husbandry and breeding. Health and reproduction are increasingly being considered because of the observed genetic correlations between reproduction, health, conformation, and performance traits in dairy cattle. Understanding the detailed genetic architecture underlying these traits would represent a major step in comprehending their interplay. Identifying known, putative or novel associations in genomics could improve animal health, welfare, and performance while allowing further adjustments in animal breeding.

RESULTS

We conducted genome-wide association studies for 25 different traits belonging to four different complexes, namely reproduction (n = 13), conformation (n = 6), production (n = 3), and metabolism (n = 3), using a cohort of over 235,000 dairy cows. As a result, we identified genome-wide significant signals for all the studied traits. The obtained summary statistics collected served as the input for a Mendelian randomisation approach (GSMR) to infer causal associations between putative exposure and reproduction traits. The study considered conformation, production, and metabolism as exposure and reproduction as outcome. A range of 139 to 252 genome-wide significant SNPs per combination were identified as instrumental variables (IVs). Out of 156 trait combinations, 135 demonstrated statistically significant effects, thereby enabling the identification of the responsible IVs. Combinations of traits related to metabolism (38 out of 39), conformation (68 out of 78), or production (29 out of 39) were found to have significant effects on reproduction. These relationships were partially non-linear. Moreover, a separate variance component estimation supported these findings, strongly correlating with the GSMR results and offering suggestions for improvement. Downstream analyses of selected representative traits per complex resulted in identifying and investigating potential physiological mechanisms. Notably, we identified both trait-specific SNPs and genes that appeared to influence specific traits per complex, as well as more general SNPs that were common between exposure and outcome traits.

CONCLUSIONS

Our study confirms the known genetic associations between reproduction traits and the three complexes tested. It provides new insights into causality, indicating a non-linear relationship between conformation and reproduction. In addition, the downstream analyses have identified several clustered genes that may mediate this association.

Multiple strategies association revealed functional candidate FASN gene for fatty acid composition in cattle

Fatty acid composition (FA) is an important indicator of meat quality in beef cattle. We investigated potential functional candidate genes for FA in beef cattle by integrating genomic and transcriptomic dataset through multiple strategies. In this study, we observed 65 SNPs overlapping with five candidate genes (CCDC57, FASN, HDAC11, ALG14, and ZMAT4) using two steps association based on the imputed sequencing variants. Using multiple traits GWAS, we further identified three significant SNPs located in the upstream of FASN and one SNP (chr19:50779529) was embedded in FASN. Of those, two SNPs were further identified as the cis-eQTL based on transcriptomic analysis of muscle tissues. Moreover, the knockdown of FASN yielded a significant reduction in intracellular triglyceride content in preadipocytes and impeded lipid droplet accumulation in adipocytes. RNA-seq analysis of preadipocytes with FASN interference revealed that the differentially expressed genes were enriched in cell differentiation and lipid metabolic pathway. Our study underscored the indispensable role of FASN in orchestrating adipocyte differentiation, and fatty acid metabolism. The integrative analysis with multiple strategies may contribute to the understanding of the genetic architecture of FA in farm animals.

Exploration and Enrichment Analysis of the QTLome for Important Traits in Livestock Species

Background: Quantitative trait loci (QTL) are genomic regions that influence essential traits in livestock. Understanding QTL distribution and density across species' genomes is crucial for animal genetics research. Objectives: This study explored the QTLome of cattle, pigs, sheep, and chickens by analyzing QTL distribution and evaluating the correlation between QTL, gene density, and chromosome size with the aim to identify QTL-enriched genomic regions. Methods: Data from 211,715 QTL (1994-2021) were retrieved from the AnimalQTLdb and analyzed using R software v4.2.1. Unique QTL annotations were identified, and redundant or inconsistent data were removed. Statistical analyses included Pearson correlations and binomial, hypergeometric, and bootstrap-based enrichment tests. Results: QTL densities per Mbp were 10 for bovine, 4 for pig, 1 for sheep, and 3 for chicken genomes. Analysis of QTL distribution across chromosomes revealed uneven patterns, with certain regions enriched for QTL. Correlation analysis revealed a strong positive relationship between QTL and gene density/chromosome size across all species (p < 0.05). Enrichment analysis identified pleiotropic regions, where QTL affect multiple traits, often aligning with known candidate and major genes. Significant QTL-enriched windows (p < 0.05) were detected, with 699 (187), 355 (68), 50 (15), and 38 (17) genomic windows for cattle, pigs, sheep, and chickens, respectively, associated with overall traits (and specific phenotypic categories). Conclusions: This study provides critical insights into QTL distribution and its correlation with gene density, offering valuable data for advancing genetic research in livestock species. The identification of QTL-enriched regions also highlights key areas for future exploration in trait improvement programs.

Empirical versus estimated accuracy of imputation: optimising filtering thresholds for sequence imputation

BACKGROUND

Genotype imputation is a cost-effective method for obtaining sequence genotypes for downstream analyses such as genome-wide association studies (GWAS). However, low imputation accuracy can increase the risk of false positives, so it is important to pre-filter data or at least assess the potential limitations due to imputation accuracy. In this study, we benchmarked three different imputation programs (Beagle 5.2, Minimac4 and IMPUTE5) and compared the empirical accuracy of imputation with the software estimated accuracy of imputation (Rsqsoft). We also tested the accuracy of imputation in cattle for autosomal and X chromosomes, SNP and INDEL, when imputing from either low-density or high-density genotypes.

RESULTS

The accuracy of imputing sequence variants from real high-density genotypes was higher than from low-density genotypes. In our software benchmark, all programs performed well with only minor differences in accuracy. While there was a close relationship between empirical imputation accuracy and the imputation Rsqsoft, this differed considerably for Minimac4 compared to Beagle 5.2 and IMPUTE5. We found that the Rsqsoft threshold for removing poorly imputed variants must be customised according to the software and this should be accounted for when merging data from multiple studies, such as in meta-GWAS studies. We also found that imposing an Rsqsoft filter has a positive impact on genomic regions with poor imputation accuracy due to large segmental duplications that are susceptible to error-prone alignment. Overall, our results showed that on average the imputation accuracy for INDEL was approximately 6% lower than SNP for all software programs. Importantly, the imputation accuracy for the non-PAR (non-Pseudo-Autosomal Region) of the X chromosome was comparable to autosomal imputation accuracy, while for the PAR it was substantially lower, particularly when starting from low-density genotypes.

CONCLUSIONS

This study provides an empirically derived approach to apply customised software-specific Rsqsoft thresholds for downstream analyses of imputed variants, such as needed for a meta-GWAS. The very poor empirical imputation accuracy for variants on the PAR when starting from low density genotypes demonstrates that this region should be imputed starting from a higher density of real genotypes.

Genetic and genomic analysis of reproduction traits in holstein cattle using SNP chip data and imputed sequence level genotypes

BACKGROUND

Reproductive performance plays an important role in animal welfare, health and profitability in animal husbandry and breeding. It is well established that there is a negative correlation between performance and reproduction in dairy cattle. This relationship is being increasingly considered in breeding programs. By elucidating the genetic architecture of underlying reproduction traits, it will be possible to make a more detailed contribution to this. Our study followed two approaches to elucidate this area; in a first part, variance components were estimated for 14 different calving and fertility traits, and then genome-wide association studies were performed for 13 reproduction traits on imputed sequence-level genotypes with subsequent enrichment analyses.

RESULTS

Variance components analyses showed a low to moderate heritability (h2) for the traits analysed, ranging from 0.014 for endometritis up to 0.271 for stillbirth, indicating variable degrees of variation within the reproduction traits. For genome-wide association studies, we were able to detect genome-wide significant association signals for nine out of 13 analysed traits after Bonferroni correction on chromosome 6, 18 and the X chromosome. In total, we detected over 2700 associated SNPs encircling more than 90 different genes using the imputed whole-genome sequence data. Functional associations were reviewed so far known and potential candidate regions in the proximity of reproduction events were hypothesised.

CONCLUSION

Our results confirm previous findings of other authors in a comprehensive cohort including 13 different traits at the same time. Additionally, we identified new candidate genes involved in dairy cattle reproduction and made initial suggestions regarding their potential impact, with special regard to the X chromosome as a putative information source for further research. This work can make a contribution to reveal the genetic architecture of reproduction traits in context of trait specific interactions.

Genome-wide fine-mapping improves identification of causal variants

Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 600 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.

Genome-wide fine-mapping improves identification of causal variants

Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 599 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.

Archetypal clustering reveals physiological mechanisms linking milk yield and fertility in dairy cattle

Fertility in dairy cattle has declined as an unintended consequence of single-trait selection for high milk yield. The unfavorable genetic correlation between milk yield and fertility is now well documented; however, the underlying physiological mechanisms are still uncertain. To understand the relationship between these traits, we developed a method that clusters variants with similar patterns of effects and, after the integration of gene expression data, identifies the genes through which they are likely to act. Biological processes that are enriched in the genes of each cluster were then identified. We identified several clusters with unique patterns of effects. One of the clusters included variants associated with increased milk yield and decreased fertility, where the "archetypal" variant (i.e., the one with the largest effect) was associated with the GC gene, whereas others were associated with TRIM32, LRRK2, and U6-associated snRNA. These genes have been linked to transcription and alternative splicing, suggesting that these processes are likely contributors to the unfavorable relationship between the 2 traits. Another cluster, with archetypal variant near DGAT1 and including variants associated with CDH2, BTRC, SFRP2, ZFHX3, and SLITRK5, appeared to affect milk yield but have little effect on fertility. These genes have been linked to insulin, adipose tissue, and energy metabolism. A third cluster with archetypal variant near ZNF613 and including variants associated with ROBO1, EFNA5, PALLD, GPC6, and PTPRT were associated with fertility but not milk yield. These genes have been linked to GnRH neuronal migration, embryonic development, or ovarian function. The use of archetypal clustering to group variants with similar patterns of effects may assist in identifying the biological processes underlying correlated traits. The method is hypothesis generating and requires experimental confirmation. However, we have uncovered several novel mechanisms potentially affecting milk production and fertility such as GnRH neuronal migration. We anticipate our method to be a starting point for experimental research into novel pathways, which have been previously unexplored within the context of dairy production.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes

BACKGROUND

Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected.

RESULTS

A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3.

CONCLUSIONS

The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Genome-wide association study of age at puberty and its (co)variances with fertility and stature in growing and lactating Holstein-Friesian dairy cattle

Reproductive performance is a key determinant of cow longevity in a pasture-based, seasonal dairy system. Unfortunately, direct fertility phenotypes such as intercalving interval or pregnancy rate tend to have low heritabilities and occur relatively late in an animal's life. In contrast, age at puberty (AGEP) is a moderately heritable, early-in-life trait that may be estimated using an animal's age at first measured elevation in blood plasma progesterone (AGEP4) concentrations. Understanding the genetic architecture of AGEP4 in addition to genetic relationships between AGEP4 and fertility traits in lactating cows is important, as is its relationship with body size in the growing animal. Thus, the objectives of this research were 3-fold. First, to estimate the genetic and phenotypic (co)variances between AGEP4 and subsequent fertility during first and second lactations. Second, to quantify the associations between AGEP4 and height, length, and BW measured when animals were approximately 11 mo old (standard deviation = 0.5). Third, to identify genomic regions that are likely to be associated with variation in AGEP4. We measured AGEP4, height, length, and BW in approximately 5,000 Holstein-Friesian or Holstein-Friesian × Jersey crossbred yearling heifers across 54 pasture-based herds managed in seasonal calving farm systems. We also obtained calving rate (CR42, success or failure to calve within the first 42 d of the seasonal calving period), breeding rate (PB21, success or failure to be presented for breeding within the first 21 d of the seasonal breeding period) and pregnancy rate (PR42, success or failure to become pregnant within the first 42 d of the seasonal breeding period) phenotypes from their first and second lactations. The animals were genotyped using the Weatherby's Versa 50K SNP array (Illumina, San Diego, CA). The estimated heritabilities of AGEP4, height, length, and BW were 0.34 (90% credibility interval [CRI]: 0.30, 0.37), 0.28 (90% CRI: 0.25, 0.31), 0.21 (90% CRI: 0.18, 0.23), and 0.33 (90% CRI: 0.30, 0.36), respectively. In contrast, the heritabilities of CR42, PB21 and PR42 were all <0.05 in both first and second lactations. The genetic correlations between AGEP4 and these fertility traits were generally moderate, ranging from 0.11 to 0.60, whereas genetic correlations between AGEP4 and yearling body-conformation traits ranged from 0.02 to 0.28. Our GWAS highlighted a genomic window on chromosome 5 that was strongly associated with variation in AGEP4. We also identified 4 regions, located on chromosomes 14, 6, 1, and 11 (in order of decreasing importance), that exhibited suggestive associations with AGEP4. Our results show that AGEP4 is a reasonable predictor of estimated breeding values for fertility traits in lactating cows. Although the GWAS provided insights into genetic mechanisms underpinning AGEP4, further work is required to test genomic predictions of fertility that use this information.

Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits

BACKGROUND

Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS).

RESULTS

In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits.

CONCLUSIONS

The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits

BACKGROUND

Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected.

RESULTS

All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count.

CONCLUSIONS

Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows

BACKGROUND

Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM.

RESULTS

In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect.

CONCLUSION

Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Unraveling the genomic background of resilience based on variability in milk yield and milk production levels in North American Holstein cattle through genome-wide association study and Mendelian randomization analyses

Breeding more resilient animals will benefit the dairy cattle industry in the long term, especially as global climate changes become more severe. Previous studies have reported genetic parameters for various milk yield-based resilience indicators, but the underlying genomic background of these traits remain unknown. In this study, we conducted GWAS of 62,029 SNPs with 4 milk yield-based resilience indicators, including the weighted occurrence frequency (wfPert) and accumulated milk losses (dPert) of milk yield perturbations, and log-transformed variance (LnVar) and lag-1 autocorrelation (rauto) of daily yield residuals. These variables were previously derived from 5.6 million daily milk yield records from 21,350 lactations (parities 1-3) of 11,787 North American Holstein cows. The average daily milk yield (ADMY) throughout lactation was also included to compare the shared genetic background of resilience indicators with milk yield. The differential genetic background of these indicators was first revealed by the significant genomic regions identified and significantly enriched biological pathways of positional candidate genes, which confirmed the genetic difference among resilience indicators. Interestingly, the functional analyses of candidate genes suggested that the regulation of intestinal homeostasis is most likely affecting resilience derived based on variability in milk yield. Based on Mendelian randomization analyses of multiple instrumental SNPs, we further found an unfavorable causal association of ADMY with LnVar. In conclusion, the resilience indicators evaluated are genetically different traits, and there are causal associations of milk yield with some of the resilience indicators evaluated. In addition to providing biological insights into the molecular regulation mechanisms of resilience derived based on variability in milk yield, this study also indicates the need for developing selection indexes combining multiple indicator traits and taking into account their genetic relationship for breeding more resilient dairy cattle.

The hazards of genotype imputation when mapping disease susceptibility variants

BACKGROUND

The cost-free increase in statistical power of using imputation to infer missing genotypes is undoubtedly appealing, but is it hazard-free? This case study of three type-2 diabetes (T2D) loci demonstrates that it is not; it sheds light on why this is so and raises concerns as to the shortcomings of imputation at disease loci, where haplotypes differ between cases and reference panel.

RESULTS

T2D-associated variants were previously identified using targeted sequencing. We removed these significantly associated SNPs and used neighbouring SNPs to infer them by imputation. We compared imputed with observed genotypes, examined the altered pattern of T2D-SNP association, and investigated the cause of imputation errors by studying haplotype structure. Most T2D variants were incorrectly imputed with a low density of scaffold SNPs, but the majority failed to impute even at high density, despite obtaining high certainty scores. Missing and discordant imputation errors, which were observed disproportionately for the risk alleles, produced monomorphic genotype calls or false-negative associations. We show that haplotypes carrying risk alleles are considerably more common in the T2D cases than the reference panel, for all loci.

CONCLUSIONS

Imputation is not a panacea for fine mapping, nor for meta-analysing multiple GWAS based on different arrays and different populations. A total of 80% of the SNPs we have tested are not included in array platforms, explaining why these and other such associated variants may previously have been missed. Regardless of the choice of software and reference haplotypes, imputation drives genotype inference towards the reference panel, introducing errors at disease loci.

A compendium of genetic regulatory effects across pig tissues

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

Genome-wide association study of anogenital distance and its (co)variances with fertility in growing and lactating Holstein-Friesian dairy cattle

Anogenital distance (AGD) is a moderately heritable trait that can be measured at a young age that may provide an opportunity to indirectly select for improved fertility in dairy cattle. In this study, we characterized AGD and its genetic and phenotypic relationships with a range of body stature and fertility traits. We measured AGD, shoulder height, body length, and body weight in a population of 5,010 Holstein-Friesian and Holstein-Friesian × Jersey crossbred heifers at approximately 11 mo of age (AGD1). These animals were born in 2018 across 54 seasonal calving, pasture-based dairy herds. A second measure of AGD was collected in a subset of herds (n = 17; 1,956 animals) when the animals averaged 29 mo of age (AGD2). Fertility measures included age at puberty (AGEP), then time of calving, breeding, and pregnancy during the first and second lactations. We constructed binary traits reflecting the animal's ability to calve during the first 42 d of their herd's seasonal calving period (CR42), be presented for breeding during the first 21 d of the seasonal breeding period (PB21) and become pregnant during the first 42 d of the seasonal breeding period (PR42). The posterior mean of sampled heritabilities for AGD1 was 0.23, with 90% of samples falling within a credibility interval (90% CRI) of 0.20 to 0.26, whereas the heritability of AGD2 was 0.29 (90% CRI 0.24 to 0.34). The relationship between AGD1 and AGD2 was highly positive, with a genetic correlation of 0.89 (90% CRI 0.82 to 0.94). Using a GWAS analysis of 2,460 genomic windows based on 50k genotype data, we detected a region on chromosome 20 that was highly associated with variation in AGD1, and a second region on chromosome 13 that was moderately associated with variation in AGD1. We did not detect any genomic regions associated with AGD2 which was measured in fewer animals. The genetic correlation between AGD1 and AGEP was 0.10 (90% CRI 0.00 to 0.19), whereas the genetic correlation between AGD2 and AGEP was 0.30 (90% CRI 0.15 to 0.44). The timing of calving, breeding, and pregnancy (CR42, PB21, and PR42) during first or second lactations exhibited moderate genetic relationships with AGD1 (0.19 to 0.52) and AGD2 (0.46 to 0.63). Genetic correlations between AGD and body stature traits were weak (≤0.16). We conclude that AGD is a moderately heritable trait, which may have value as an early-in-life genetic predictor for reproductive success during lactation.

Gene expression and RNA splicing explain large proportions of the heritability for complex traits in cattle

Many quantitative trait loci (QTLs) are in non-coding regions. Therefore, QTLs are assumed to affect gene regulation. Gene expression and RNA splicing are primary steps of transcription, so DNA variants changing gene expression (eVariants) or RNA splicing (sVariants) are expected to significantly affect phenotypes. We quantify the contribution of eVariants and sVariants detected from 16 tissues (n = 4,725) to 37 traits of ∼120,000 cattle (average magnitude of genetic correlation between traits = 0.13). Analyzed in Bayesian mixture models, averaged across 37 traits, cis and trans eVariants and sVariants detected from 16 tissues jointly explain 69.2% (SE = 0.5%) of heritability, 44% more than expected from the same number of random variants. This 69.2% includes an average of 24% from trans e-/sVariants (14% more than expected). Averaged across 56 lipidomic traits, multi-tissue cis and trans e-/sVariants also explain 71.5% (SE = 0.3%) of heritability, demonstrating the essential role of proximal and distal regulatory variants in shaping mammalian phenotypes.

Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences

High-quality whole-genome resequencing in large-scale pig populations with pedigree structure and multiple breeds would enable accurate construction of haplotype and robust selection-signature detection. Here, we sequence 740 pigs, combine with 149 of our previously published resequencing data, retrieve 207 resequencing datasets, and form a panel of worldwide distributed wild boars, aboriginal and highly selected pigs with pedigree structures, amounting to 1096 genomes from 43 breeds. Combining with their haplotype-informative reads and pedigree structure, we accurately construct a panel of 1874 haploid genomes with 41,964,356 genetic variants. We further demonstrate its valuable applications in GWAS by identifying five novel loci for intramuscular fat content, and in genomic selection by increasing the accuracy of estimated breeding value by 36.7%. In evolutionary selection, we detect MUC13 gene under a long-term balancing selection, as well as NPR3 gene under positive selection for pig stature. Our study provides abundant genomic variations for robust selection-signature detection and accurate haplotypes for deciphering complex traits in pigs.

Androgen receptor binding sites enabling genetic prediction of mortality due to prostate cancer in cancer-free subjects

Prostate cancer (PrCa) is the second most common cancer worldwide in males. While strongly warranted, the prediction of mortality risk due to PrCa, especially before its development, is challenging. Here, we address this issue by maximizing the statistical power of genetic data with multi-ancestry meta-analysis and focusing on binding sites of the androgen receptor (AR), which has a critical role in PrCa. Taking advantage of large Japanese samples ever, a multi-ancestry meta-analysis comprising more than 300,000 subjects in total identifies 9 unreported loci including ZFHX3, a tumor suppressor gene, and successfully narrows down the statistically finemapped variants compared to European-only studies, and these variants strongly enrich in AR binding sites. A polygenic risk scores (PRS) analysis restricting to statistically finemapped variants in AR binding sites shows among cancer-free subjects, individuals with a PRS in the top 10% have a strongly higher risk of the future death of PrCa (HR: 5.57, P = 4.2 × 10-10). Our findings demonstrate the potential utility of leveraging large-scale genetic data and advanced analytical methods in predicting the mortality of PrCa.

A genome-wide association study for loin depth and muscle pH in pigs from intensely selected purebred lines

BACKGROUND

Genome-wide association studies (GWAS) aim at identifying genomic regions involved in phenotype expression, but identifying causative variants is difficult. Pig Combined Annotation Dependent Depletion (pCADD) scores provide a measure of the predicted consequences of genetic variants. Incorporating pCADD into the GWAS pipeline may help their identification. Our objective was to identify genomic regions associated with loin depth and muscle pH, and identify regions of interest for fine-mapping and further experimental work. Genotypes for ~ 40,000 single nucleotide morphisms (SNPs) were used to perform GWAS for these two traits, using de-regressed breeding values (dEBV) for 329,964 pigs from four commercial lines. Imputed sequence data was used to identify SNPs in strong ([Formula: see text] 0.80) linkage disequilibrium with lead GWAS SNPs with the highest pCADD scores.

RESULTS

Fifteen distinct regions were associated with loin depth and one with loin pH at genome-wide significance. Regions on chromosomes 1, 2, 5, 7, and 16, explained between 0.06 and 3.55% of the additive genetic variance and were strongly associated with loin depth. Only a small part of the additive genetic variance in muscle pH was attributed to SNPs. The results of our pCADD analysis suggests that high-scoring pCADD variants are enriched for missense mutations. Two close but distinct regions on SSC1 were associated with loin depth, and pCADD identified the previously identified missense variant within the MC4R gene for one of the lines. For loin pH, pCADD identified a synonymous variant in the RNF25 gene (SSC15) as the most likely candidate for the muscle pH association. The missense mutation in the PRKAG3 gene known to affect glycogen content was not prioritised by pCADD for loin pH.

CONCLUSIONS

For loin depth, we identified several strong candidate regions for further statistical fine-mapping that are supported in the literature, and two novel regions. For loin muscle pH, we identified one previously identified associated region. We found mixed evidence for the utility of pCADD as an extension of heuristic fine-mapping. The next step is to perform more sophisticated fine-mapping and expression quantitative trait loci (eQTL) analysis, and then interrogate candidate variants in vitro by perturbation-CRISPR assays.

ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure

The present investigation aimed to identify genome wide SNPs and to carry out diversity and population structure study using ddRAD-seq based genotyping of 58 individuals of six indigenous milch cattle breeds (Bos indicus) such as Sahiwal, Gir, Rathi, Tharparkar, Red Sindhi and Kankrej of India. A high percentage of reads (94.53%) were mapped to the Bos taurus (ARS-UCD1.2) reference genome assembly. Following filtration criteria, a total of 84,027 high quality SNPs were identified across the genome of 6 cattle breeds with the highest number of SNPs observed in Gir (34,743), followed by Red Sindhi (13,092), Kankrej (12,812), Sahiwal (8956), Tharparkar (7356) and Rathi (7068). Most of these SNPs were distributed in the intronic regions (53.87%) followed by intergenic regions (34.94%) while only 1.23% were located in the exonic regions. Together with analysis of nucleotide diversity (π = 0.373), Tajima's D (D value ranging from - 0.295 to 0.214), observed heterozygosity (HO ranging from 0.464 to 0.551), inbreeding coefficient (FIS ranging from - 0.253 to 0.0513) suggested for the presence of sufficient within breed diversity in the 6 major milch breeds of India. The phylogenetic based structuring, principal component and admixture analysis revealed genetic distinctness as well as purity of almost all of the 6 cattle breeds. Overall, our strategy has successfully identified thousands of high-quality genome wide SNPs that will further enrich the Bos indicus representation basic information about genetic diversity and structure of 6 major Indian milch cattle breeds which should have implications for better management and conservation of valuable indicine cattle diversity.

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle

BACKGROUND

A detailed understanding of genetic variants that affect beef merit helps maximize the efficiency of breeding for improved production merit in beef cattle. To prioritize the putative variants and genes, we ran a comprehensive genome-wide association studies (GWAS) analysis for 21 agronomic traits using imputed whole-genome variants in Simmental beef cattle. Then, we applied expression quantitative trait loci (eQTL) mapping between the genotype variants and transcriptome of three tissues (longissimus dorsi muscle, backfat, and liver) in 120 cattle.

RESULTS

We identified 1,580 association signals for 21 beef agronomic traits using GWAS. We then illuminated 854,498 cis-eQTLs for 6,017 genes and 46,970 trans-eQTLs for 1,903 genes in three tissues and built a synergistic network by integrating transcriptomics with agronomic traits. These cis-eQTLs were preferentially close to the transcription start site and enriched in functional regulatory regions. We observed an average of 43.5% improvement in cis-eQTL discovery using multi-tissue eQTL mapping. Fine-mapping analysis revealed that 111, 192, and 194 variants were most likely to be causative to regulate gene expression in backfat, liver, and muscle, respectively. The transcriptome-wide association studies identified 722 genes significantly associated with 11 agronomic traits. Via the colocalization and Mendelian randomization analyses, we found that eQTLs of several genes were associated with the GWAS signals of agronomic traits in three tissues, which included genes, such as NADSYN1, NDUFS3, LTF and KIFC2 in liver, GRAMD1C, TMTC2 and ZNF613 in backfat, as well as TIGAR, NDUFS3 and L3HYPDH in muscle that could serve as the candidate genes for economic traits.

CONCLUSIONS

The extensive atlas of GWAS, eQTL, fine-mapping, and transcriptome-wide association studies aid in the suggestion of potentially functional variants and genes in cattle agronomic traits and will be an invaluable source for genomics and breeding in beef cattle.

Improving the accuracy of genomic prediction for meat quality traits using whole genome sequence data in pigs

BACKGROUND

Pork quality can directly affect customer purchase tendency and meat quality traits have become valuable in modern pork production. However, genetic improvement has been slow due to high phenotyping costs. In this study, whole genome sequence (WGS) data was used to evaluate the prediction accuracy of genomic best linear unbiased prediction (GBLUP) for meat quality in large-scale crossbred commercial pigs.

RESULTS

We produced WGS data (18,695,907 SNPs and 2,106,902 INDELs exceed quality control) from 1,469 sequenced Duroc × (Landrace × Yorkshire) pigs and developed a reference panel for meat quality including meat color score, marbling score, L* (lightness), a* (redness), and b* (yellowness) of genomic prediction. The prediction accuracy was defined as the Pearson correlation coefficient between adjusted phenotypes and genomic estimated breeding values in the validation population. Using different marker density panels derived from WGS data, accuracy differed substantially among meat quality traits, varied from 0.08 to 0.47. Results showed that MultiBLUP outperform GBLUP and yielded accuracy increases ranging from 17.39% to 75%. We optimized the marker density and found medium- and high-density marker panels are beneficial for the estimation of heritability for meat quality. Moreover, we conducted genotype imputation from 50K chip to WGS level in the same population and found average concordance rate to exceed 95% and r² = 0.81.

CONCLUSIONS

Overall, estimation of heritability for meat quality traits can benefit from the use of WGS data. This study showed the superiority of using WGS data to genetically improve pork quality in genomic prediction.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Better safe than sorry-Whole-genome sequencing indicates that missense variants are significant in susceptibility to COVID-19

Undoubtedly, genetic factors play an important role in susceptibility and resistance to COVID-19. In this study, we conducted the GWAS analysis. Out of 15,489,173 SNPs, we identified 18,191 significant SNPs for severe and 11,799 SNPs for resistant phenotype, showing that a great number of loci were significant in different COVID-19 representations. The majority of variants were synonymous (60.56% for severe, 58.46% for resistant phenotype) or located in introns (55.77% for severe, 59.83% for resistant phenotype). We identified the most significant SNPs for a severe outcome (in AJAP1 intron) and for COVID resistance (in FIG4 intron). We found no missense variants with a potential causal function on resistance to COVID-19; however, two missense variants were determined as significant a severe phenotype (in PM20D1 and LRP4 exons). None of the aforementioned SNPs and missense variants found in this study have been previously associated with COVID-19.

Effects of common full-sib families on accuracy of genomic prediction for tagging weight in striped catfish Pangasianodon hypophthalmus

Common full-sib families (c ² ) make up a substantial proportion of total phenotypic variation in traits of commercial importance in aquaculture species and omission or inclusion of the c ² resulted in possible changes in genetic parameter estimates and re-ranking of estimated breeding values. However, the impacts of common full-sib families on accuracy of genomic prediction for commercial traits of economic importance are not well known in many species, including aquatic animals. This research explored the impacts of common full-sib families on accuracy of genomic prediction for tagging weight in a population of striped catfish comprising 11,918 fish traced back to the base population (four generations), in which 560 individuals had genotype records of 14,154 SNPs. Our single step genomic best linear unbiased prediction (ssGLBUP) showed that the accuracy of genomic prediction for tagging weight was reduced by 96.5%-130.3% when the common full-sib families were included in statistical models. The reduction in the prediction accuracy was to a smaller extent in multivariate analysis than in univariate models. Imputation of missing genotypes somewhat reduced the upward biases in the prediction accuracy for tagging weight. It is therefore suggested that genomic evaluation models for traits recorded during the early phase of growth development should account for the common full-sib families to minimise possible biases in the accuracy of genomic prediction and hence, selection response.

Multi-breed genomic predictions and functional variants for fertility of tropical bulls

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Utility of multi-omics data to inform genomic prediction of heifer fertility traits

Biologically informed single nucleotide polymorphisms (SNPs) impact genomic prediction accuracy of the target traits. Our previous genomics, proteomics, and transcriptomics work identified candidate genes related to puberty and fertility in Brahman heifers. We aimed to test this biological information for capturing heritability and predicting heifer fertility traits in another breed i.e., Tropical Composite. The SNP from the identified genes including 10 kilobases (kb) region on either side were selected as biologically informed SNP set. The SNP from the rest of the Bos taurus genes including 10-kb region on either side were selected as biologically uninformed SNP set. Bovine high-density (HD) complete SNP set (628,323 SNP) was used as a control. Two populations-Tropical Composites (N = 1331) and Brahman (N = 2310)-had records for three traits: pregnancy after first mating season (PREG1, binary), first conception score (FCS, score 1 to 3), and rebreeding score (REB, score 1 to 3.5). Using the best linear unbiased prediction method, effectiveness of each SNP set to predict the traits was tested in two scenarios: a 5-fold cross-validation within Tropical Composites using biological information from Brahman studies, and application of prediction equations from one breed to the other. The accuracy of prediction was calculated as the correlation between genomic estimated breeding values and adjusted phenotypes. Results show that biologically informed SNP set estimated heritabilities not significantly better than the control HD complete SNP set in Tropical Composites; however, it captured all the observed genetic variance in PREG1 and FCS when modeled together with the biologically uninformed SNP set. In 5-fold cross-validation within Tropical Composites, the biologically informed SNP set performed marginally better (statistically insignificant) in terms of prediction accuracies (PREG1: 0.20, FCS: 0.13, and REB: 0.12) as compared to HD complete SNP set (PREG1: 0.17, FCS: 0.10, and REB: 0.11), and biologically uninformed SNP set (PREG1: 0.16, FCS: 0.10, and REB: 0.11). Across-breed use of prediction equations still remained a challenge: accuracies by all SNP sets dropped to around zero for all traits. The performance of biologically informed SNP was not significantly better than other sets in Tropical Composites. However, results indicate that biological information obtained from Brahman was successful to predict the fertility traits in Tropical Composite population.

Animal-SNPAtlas: a comprehensive SNP database for multiple animals

Single-nucleotide polymorphisms (SNPs) as the most important type of genetic variation are widely used in describing population characteristics and play vital roles in animal genetics and breeding. Large amounts of population genetic variation resources and tools have been developed in human, which provided solid support for human genetic studies. However, compared with human, the development of animal genetic variation databases was relatively slow, which limits the genetic researches in these animals. To fill this gap, we systematically identified ∼ 499 million high-quality SNPs from 4784 samples of 20 types of animals. On that basis, we annotated the functions of SNPs, constructed high-density reference panels and calculated genome-wide linkage disequilibrium (LD) matrixes. We further developed Animal-SNPAtlas, a user-friendly database (http://gong_lab.hzau.edu.cn/Animal_SNPAtlas/) which includes high-quality SNP datasets and several support tools for multiple animals. In Animal-SNPAtlas, users can search the functional annotation of SNPs, perform online genotype imputation, explore and visualize LD information, browse variant information using the genome browser and download SNP datasets for each species. With the massive SNP datasets and useful tools, Animal-SNPAtlas will be an important fundamental resource for the animal genomics, genetics and breeding community.

Genomic prediction with whole-genome sequence data in intensely selected pig lines

Background

Early simulations indicated that whole-genome sequence data (WGS) could improve the accuracy of genomic predictions within and across breeds. However, empirical results have been ambiguous so far. Large datasets that capture most of the genomic diversity in a population must be assembled so that allele substitution effects are estimated with high accuracy. The objectives of this study were to use a large pig dataset from seven intensely selected lines to assess the benefits of using WGS for genomic prediction compared to using commercial marker arrays and to identify scenarios in which WGS provides the largest advantage.

Methods

We sequenced 6931 individuals from seven commercial pig lines with different numerical sizes. Genotypes of 32.8 million variants were imputed for 396,100 individuals (17,224 to 104,661 per line). We used BayesR to perform genomic prediction for eight complex traits. Genomic predictions were performed using either data from a standard marker array or variants preselected from WGS based on association tests.

Results

The accuracies of genomic predictions based on preselected WGS variants were not robust across traits and lines and the improvements in prediction accuracy that we achieved so far with WGS compared to standard marker arrays were generally small. The most favourable results for WGS were obtained when the largest training sets were available and standard marker arrays were augmented with preselected variants with statistically significant associations to the trait. With this method and training sets of around 80k individuals, the accuracy of within-line genomic predictions was on average improved by 0.025. With multi-line training sets, improvements of 0.04 compared to marker arrays could be expected.

Conclusions

Our results showed that WGS has limited potential to improve the accuracy of genomic predictions compared to marker arrays in intensely selected pig lines. Thus, although we expect that larger improvements in accuracy from the use of WGS are possible with a combination of larger training sets and optimised pipelines for generating and analysing such datasets, the use of WGS in the current implementations of genomic prediction should be carefully evaluated against the cost of large-scale WGS data on a case-by-case basis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00756-0.

A multi-tissue atlas of regulatory variants in cattle

Characterization of genetic regulatory variants acting on livestock gene expression is essential for interpreting the molecular mechanisms underlying traits of economic value and for increasing the rate of genetic gain through artificial selection. Here we build a Cattle Genotype-Tissue Expression atlas (CattleGTEx) as part of the pilot phase of the Farm animal GTEx (FarmGTEx) project for the research community based on 7,180 publicly available RNA-sequencing (RNA-seq) samples. We describe the transcriptomic landscape of more than 100 tissues/cell types and report hundreds of thousands of genetic associations with gene expression and alternative splicing for 23 distinct tissues. We evaluate the tissue-sharing patterns of these genetic regulatory effects, and functionally annotate them using multiomics data. Finally, we link gene expression in different tissues to 43 economically important traits using both transcriptome-wide association and colocalization analyses to decipher the molecular regulatory mechanisms underpinning such agronomic traits in cattle.

Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity

Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.

Development and validation of a horse reference panel for genotype imputation

Background

Genotype imputation is a cost-effective method to generate sequence-level genotypes for a large number of animals. Its application can improve the power of genomic studies, provided that the accuracy of imputation is sufficiently high. The purpose of this study was to develop an optimal strategy for genotype imputation from genotyping array data to sequence level in German warmblood horses, and to investigate the effect of different factors on the accuracy of imputation. Publicly available whole-genome sequence data from 317 horses of 46 breeds was used to conduct the analyses.

Results

Depending on the size and composition of the reference panel, the accuracy of imputation from medium marker density (60K) to sequence level using the software Beagle 5.1 ranged from 0.64 to 0.70 for horse chromosome 3. Generally, imputation accuracy increased as the size of the reference panel increased, but if genetically distant individuals were included in the panel, the accuracy dropped. Imputation was most precise when using a reference panel of multiple but related breeds and the software Beagle 5.1, which outperformed the other two tested computer programs, Impute 5 and Minimac 4. Genome-wide imputation for this scenario resulted in a mean accuracy of 0.66. Stepwise imputation from 60K to 670K markers and subsequently to sequence level did not improve the accuracy of imputation. However, imputation from higher density (670K) was considerably more accurate (about 0.90) than from medium density. Likewise, imputation in genomic regions with a low marker coverage resulted in a reduced accuracy of imputation.

Conclusions

The accuracy of imputation in horses was influenced by the size and composition of the reference panel, the marker density of the genotyping array, and the imputation software. Genotype imputation can be used to extend the limited amount of available sequence-level data from horses in order to boost the power of downstream analyses, such as genome-wide association studies, or the detection of embryonic lethal variants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00740-8.

Identification of Haplotypes Associated With Resistance to Bacterial Cold Water Disease in Rainbow Trout Using Whole-Genome Resequencing

Bacterial cold water disease (BCWD) is an important disease in rainbow trout aquaculture. Previously, we have identified and validated two major QTL (quantitative trait loci) for BCWD resistance, located on chromosomes Omy08 and Omy25, in the odd-year Troutlodge May spawning population. We also demonstrated that marker-assisted selection (MAS) for BCWD resistance using the favorable haplotypes associated with the two major QTL is feasible. However, each favorable haplotype spans a large genomic region of 1.3–1.6 Mb. Recombination events within the haplotype regions will result in new haplotypes associated with BCWD resistance, which will reduce the accuracy of MAS for BCWD resistance over time. The objectives of this study were 1) to identify additional SNPs (single nucleotide polymorphisms) associated with BCWD resistance using whole-genome sequencing (WGS); 2) to validate the SNPs associated with BCWD resistance using family-based association mapping; 3) to refine the haplotypes associated with BCWD resistance; and 4) to evaluate MAS for BCWD resistance using the refined QTL haplotypes. Four consecutive generations of the Troutlodge May spawning population were evaluated for BCWD resistance. Parents and offspring were sequenced as individuals and in pools based on their BCWD phenotypes. Over 12 million SNPs were identified by mapping the sequences from the individuals and pools to the reference genome. SNPs with significantly different allele frequencies between the two BCWD phenotype groups were selected to develop SNP assays for family-based association mapping in three consecutive generations of the Troutlodge May spawning population. Among the 78 SNPs derived from WGS, 77 SNPs were associated with BCWD resistance in at least one of the three consecutive generations. The additional SNPs associated with BCWD resistance allowed us to reduce the physical sizes of haplotypes associated with BCWD resistance to less than 0.5 Mb. We also demonstrated that the refined QTL haplotypes can be used for MAS in the Troutlodge May spawning population. Therefore, the SNPs and haplotypes reported in this study provide additional resources for improvement of BCWD resistance in rainbow trout.

Improving Genomic Selection for Heat Tolerance in Dairy Cattle: Current Opportunities and Future Directions

Heat tolerance is the ability of an animal to maintain production and reproduction levels under hot and humid conditions and is now a trait of economic relevance in dairy systems worldwide because of an escalating warming climate. The Australian dairy population is one of the excellent study models for enhancing our understanding of the biology of heat tolerance because they are predominantly kept outdoors on pastures where they experience direct effects of weather elements (e.g., solar radiation). In this article, we focus on evidence from recent studies in Australia that leveraged large a dataset [∼40,000 animals with phenotypes and 15 million whole-genome sequence variants] to elucidate the genetic basis of thermal stress as a critical part of the strategy to breed cattle adapted to warmer environments. Genotype-by-environment interaction (i.e., G × E) due to temperature and humidity variation is increasing, meaning animals are becoming less adapted (i.e., more sensitive) to changing environments. There are opportunities to reverse this trend and accelerate adaptation to warming climate by 1) selecting robust or heat-resilient animals and 2) including resilience indicators in breeding goals. Candidate causal variants related to the nervous system and metabolic functions are relevant for heat tolerance and, therefore, key for improving this trait. This could include adding these variants in the custom SNP panels used for routine genomic evaluations or as the basis to design specific agonist or antagonist compounds for lowering core body temperature under heat stress conditions. Indeed, it was encouraging to see that adding prioritized functionally relevant variants into the 50k SNP panel (i.e., the industry panel used for genomic evaluation in Australia) increased the prediction accuracy of heat tolerance by up to 10% units. This gain in accuracy is critical because genetic improvement has a linear relationship with prediction accuracy. Overall, while this article used data mainly from Australia, this could benefit other countries that aim to develop breeding values for heat tolerance, considering that the warming climate is becoming a topical issue worldwide.

Multi-omics approaches for comprehensive analysis and understanding of the immune response in the miniature pig breed

The porcine immune system has an important role in pre-clinical studies together with understanding the biological response mechanisms before entering into clinical trials. The size distribution of the Korean minipig is an important feature that make this breed ideal for biomedical research and safe practice in post clinical studies. The extremely tiny (ET) minipig serves as an excellent model for various biomedical research studies, but the comparatively frail and vulnerable immune response to the environment over its Large (L) size minipig breed leads to additional after born care. To overcome this pitfall, comparative analysis of the genomic regions under selection in the L type breed could provide a better understanding at the molecular level and lead to the development of an enhanced variety of ET type minipig. In this study, we utilized whole genome sequencing (WGS) to identify traces of artificial selection and integrated them with transcriptome data generated from blood samples to find strongly selected and differentially expressed genes of interest. We identified a total of 35 common genes among which 7 were differentially expressed and showed selective sweep in the L type over the ET type minipig breed. The stabilization of these genes were further confirmed using nucleotide diversity analysis, and these genes could serve as potential biomarkers for the development of a better variety of ET type pig breed.

Application of Bayesian genomic prediction methods to genome-wide association analyses

Background

Bayesian genomic prediction methods were developed to simultaneously fit all genotyped markers to a set of available phenotypes for prediction of breeding values for quantitative traits, allowing for differences in the genetic architecture (distribution of marker effects) of traits. These methods also provide a flexible and reliable framework for genome-wide association (GWA) studies. The objective here was to review developments in Bayesian hierarchical and variable selection models for GWA analyses.

Results

By fitting all genotyped markers simultaneously, Bayesian GWA methods implicitly account for population structure and the multiple-testing problem of classical single-marker GWA. Implemented using Markov chain Monte Carlo methods, Bayesian GWA methods allow for control of error rates using probabilities obtained from posterior distributions. Power of GWA studies using Bayesian methods can be enhanced by using informative priors based on previous association studies, gene expression analyses, or functional annotation information. Applied to multiple traits, Bayesian GWA analyses can give insight into pleiotropic effects by multi-trait, structural equation, or graphical models. Bayesian methods can also be used to combine genomic, transcriptomic, proteomic, and other -omics data to infer causal genotype to phenotype relationships and to suggest external interventions that can improve performance.

Conclusions

Bayesian hierarchical and variable selection methods provide a unified and powerful framework for genomic prediction, GWA, integration of prior information, and integration of information from other -omics platforms to identify causal mutations for complex quantitative traits.

Weighted Single-Step Genome-Wide Association Study Uncovers Known and Novel Candidate Genomic Regions for Milk Production Traits and Somatic Cell Score in Valle del Belice Dairy Sheep

Simple Summary

Milk production is the most economically crucial dairy sheep trait and constitutes the major genetic enhancement purpose via selective breeding. Also, mastitis is one of the most frequently encountered diseases, having a significant impact on animal welfare, milk yield, and quality. The aim of this study was to identify genomic region(s) associated with the milk production traits and somatic cell score (SCS) in Valle del Belice sheep using single-step genome-wide association (ssGWA) and genotyping data from medium density SNP panels. We identified several genomic regions (OAR1, OAR2, OAR3, OAR4, OAR6, OAR9, and OAR25) and candidate genes implicated in milk production traits and SCS. Our findings offer new insights into the genetic basis of milk production traits and SCS in dairy sheep.

Abstract

The objective of this study was to uncover genomic regions explaining a substantial proportion of the genetic variance in milk production traits and somatic cell score in a Valle del Belice dairy sheep. Weighted single-step genome-wide association studies (WssGWAS) were conducted for milk yield (MY), fat yield (FY), fat percentage (FAT%), protein yield (PY), protein percentage (PROT%), and somatic cell score (SCS). In addition, our aim was also to identify candidate genes within genomic regions that explained the highest proportions of genetic variance. Overall, the full pedigree consists of 5534 animals, of which 1813 ewes had milk data (15,008 records), and 481 ewes were genotyped with a 50 K single nucleotide polymorphism (SNP) array. The effects of markers and the genomic estimated breeding values (GEBV) of the animals were obtained by five iterations of WssGBLUP. We considered the top 10 genomic regions in terms of their explained genomic variants as candidate window regions for each trait. The results showed that top ranked genomic windows (1 Mb windows) explained 3.49, 4.04, 5.37, 4.09, 3.80, and 5.24% of the genetic variances for MY, FY, FAT%, PY, PROT%, and total SCS, respectively. Among the candidate genes found, some known associations were confirmed, while several novel candidate genes were also revealed, including PPARGC1A, LYPLA1, LEP, and MYH9 for MY; CACNA1C, PTPN1, ROBO2, CHRM3, and ERCC6 for FY and FAT%; PCSK5 and ANGPT1 for PY and PROT%; and IL26, IFNG, PEX26, NEGR1, LAP3, and MED28 for SCS. These findings increase our understanding of the genetic architecture of six examined traits and provide guidance for subsequent genetic improvement through genome selection.

Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes

Background

Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic.

Results

There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits.

Conclusions

These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08555-z.

Genome-Wide Genomic and Functional Association Study for Workability and Calving Traits in Holstein Cattle

The goal of our study was to identify the SNPs, metabolic pathways (KEGG), and gene ontology (GO) terms significantly associated with calving and workability traits in dairy cattle. We analysed direct (DCE) and maternal (MCE) calving ease, direct (DSB) and maternal (MSB) stillbirth, milking speed (MSP), and temperament (TEM) based on a Holstein-Friesian dairy cattle population consisting of 35,203 individuals. The number of animals, depending on the trait, ranged from 22,301 bulls for TEM to 30,603 for DCE. We estimated the SNP effects (based on 46,216 polymorphisms from Illumina BovineSNP50 BeadChip Version 2) using a multi-SNP mixed model. The SNP positions were mapped to genes and the GO terms/KEGG pathways of the corresponding genes were assigned. The estimation of the GO term/KEGG pathway effects was based on a mixed model using the SNP effects as dependent variables. The number of significant SNPs comprised 59 for DCE, 25 for DSB and MSP, 17 for MCE and MSB, and 7 for TEM. Significant KEGG pathways were found for MSB (2), TEM (2), and MSP (1) and 11 GO terms were significant for MSP, 10 for DCE, 8 for DSB and TEM, 5 for MCE, and 3 for MSB. From the perspective of a better understanding of the genomic background of the phenotypes, traits with low heritabilities suggest that the focus should be moved from single genes to the metabolic pathways or gene ontologies significant for the phenotype.

Recovery of mitogenomes from whole genome sequences to infer maternal diversity in 1883 modern taurine and indicine cattle

Maternal diversity based on a sub-region of mitochondrial genome or variants were commonly used to understand past demographic events in livestock. Additionally, there is growing evidence of direct association of mitochondrial genetic variants with a range of phenotypes. Therefore, this study used complete bovine mitogenomes from a large sequence database to explore the full spectrum of maternal diversity. Mitogenome diversity was evaluated among 1883 animals representing 156 globally important cattle breeds. Overall, the mitogenomes were diverse: presenting 11 major haplogroups, expanding to 1309 unique haplotypes, with nucleotide diversity 0.011 and haplotype diversity 0.999. A small proportion of African taurine (3.5%) and indicine (1.3%) haplogroups were found among the European taurine breeds and composites. The haplogrouping was largely consistent with the population structure derived from alternate clustering methods (e.g. PCA and hierarchical clustering). Further, we present evidence confirming a new indicine subgroup (I1a, 64 animals) mainly consisting of breeds originating from China and characterised by two private mutations within the I1 haplogroup. The total genetic variation was attributed mainly to within-breed variance (96.9%). The accuracy of the imputation of missing genotypes was high (99.8%) except for the relatively rare heteroplasmic genotypes, suggesting the potential for trait association studies within a breed.