Single-step genome-wide association analyses of claw horn lesions in Holstein cattle using linear and threshold models

BACKGROUND

Lameness in dairy cattle is primarily caused by foot lesions including the claw horn lesions (CHL) of sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL). This study investigated the genetic architecture of the three CHL based on detailed animal phenotypes of CHL susceptibility and severity. Estimation of genetic parameters and breeding values, single-step genome-wide association analyses, and functional enrichment analyses were performed.

RESULTS

The studied traits were under genetic control with a low to moderate heritability. Heritability estimates of SH and SU susceptibility on the liability scale were 0.29 and 0.35, respectively. Heritability of SH and SU severity were 0.12 and 0.07, respectively. Heritability of WL was relatively lower, indicating stronger environmental influence on the presence and development of WL than the other two CHL. Genetic correlations between SH and SU were high (0.98 for lesion susceptibility and 0.59 for lesion severity), whereas genetic correlations of SH and SU with WL also tended to be positive. Candidate quantitative trait loci (QTL) were identified for all CHL, including some on Bos taurus chromosome (BTA) 3 and 18 with potential pleiotropic effects associated with multiple foot lesion traits. A genomic window of 0.65 Mb on BTA3 explained 0.41, 0.50, 0.38, and 0.49% of the genetic variance for SH susceptibility, SH severity, WL susceptibility, and WL severity, respectively. Another window on BTA18 explained 0.66, 0.41, and 0.70% of the genetic variance for SH susceptibility, SU susceptibility, and SU severity, respectively. The candidate genomic regions associated with CHL harbour annotated genes that are linked to immune system function and inflammation responses, lipid metabolism, calcium ion activities, and neuronal excitability.

CONCLUSIONS

The studied CHL are complex traits with a polygenic mode of inheritance. Most traits exhibited genetic variation suggesting that animal resistance to CHL can be improved with breeding. The CHL traits were positively correlated, which will facilitate genetic improvement for resistance to CHL as a whole. Candidate genomic regions associated with lesion susceptibility and severity of SH, SU, and WL provide insights into a global profile of the genetic background underlying CHL and inform genetic improvement programmes aiming at enhancing foot health in dairy cattle.

The bovine foot skin microbiota is associated with host genotype and the development of infectious digital dermatitis lesions

BACKGROUND

Bovine Digital Dermatitis (BDD) is a prevalent infectious disease, causing painful foot skin lesions and lameness in cattle. We describe herein the bovine foot skin microbiota and its associations with BDD using 16S rRNA gene amplicon and shotgun metagenomic sequencing on samples from 259 dairy cows from three UK dairy farms.

RESULTS

We show evidence of dysbiosis, and differences in taxonomy and functional profiles in the bovine foot skin microbiome of clinically healthy animals that subsequently develop BDD lesions, compared to those that do not. Our results suggest that taxonomical and functional differences together with alterations in ecological interactions between bacteria in the normal foot skin microbiome may predispose an animal to develop BDD lesions. Using genome-wide association and regional heritability mapping approaches, we provide first evidence for interactions between host genotype and certain members of the foot skin microbiota. We show the existence of significant genetic variation in the relative abundance of Treponema spp. and Peptoclostridium spp. and identify regions in the bovine genome that explain a significant proportion of this variation.

CONCLUSIONS

Collectively this work shows early changes in taxonomic and functional profiles of the bovine foot-skin microbiota in clinically healthy animals which are associated with subsequent development of BDD and could be relevant to prevention of disease. The description of host genetic control of members of the foot skin microbiota, combined with the association of the latter with BDD development offer new insights into a complex relationship that can be exploited in selective breeding programmes. Video Abstract.

Genomic analysis of the rare British Lop pig and identification of distinctive genomic markers

Concentration of production on a few commercial pig breeds has led to the marginalization of many native, numerically small breeds, increasing their risk of endangerment. In the UK, one such rare breed is the British Lop, a lop-eared breed, of similar origin to the Welsh breed. The objective of the present study was to address the genomic status of the British Lop and its relationship with other breeds and identify a small set of genomic markers that uniquely characterize and distinguish British Lop animals. Results have shown that the British Lop is a relatively distinct population with reduced genomic diversity and effective size consistent with its status as a rare breed. Furthermore, we demonstrated the genetic closeness of the British Lop to phenotypically similar breeds such as Landrace and Welsh as well Large White, Middle White and Pietrain. Finally, a set of 75 Single Nucleotide Polymorphisms distributed across multiple chromosomes were identified and validated as markers that can consistently distinguish British Lops from other closely related breeds. Results may inform breeding and management strategies aiming to enhance diversity as well as the development of a breed purity test.

A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep

Background

The management of farm animal genetic resources and the adaptation of animals to climate change will probably have major effects on the long-term sustainability of the livestock sector. Genomic data harbour useful relevant information that needs to be harnessed for effectively managing genetic resources. In this paper, we report the genome characterization of the highly productive Mediterranean Chios dairy sheep and focus on genetic diversity measures related with local adaptation and selection and the genetic architecture of animal resilience to weather fluctuations as a novel adaptative trait linked to climate change.

Results

We detected runs of homozygosity (ROH) and heterozygosity (ROHet) that revealed multiple highly homozygous and heterozygous hotspots across the Chios sheep genome. A particularly highly homozygous region was identified on chromosome 13 as a candidate of directional genetic selection associated with milk traits, which includes annotated genes that were previously shown to be linked to local adaptation to harsh environmental conditions. Favourable heterozygosity related with a potentially protective role against livestock diseases and enhanced overall fitness was revealed in heterozygous-rich regions on sheep chromosomes 3, 10, 13 and 19. Furthermore, genomic analyses were conducted on sheep resilience phenotypes that display changes in milk production in response to weather variation. Sheep resilience to heat stress was a significantly heritable trait (h² = 0.26) and genetically antagonistic to milk production. Genome-wide association and regional heritability mapping analyses revealed novel genomic markers and regions on chromosome 5 that were significantly associated with sheep resilience to climate change. Subsequently, an annotation analysis detected a set of genes on chromosome 5 that were associated with olfactory receptor complexes that could participate in heat stress mitigation through changes in respiration rate and respiratory evaporation. Other genes were grouped in previously reported biological processes relevant to livestock heat dissipation, including stress and immune response.

Conclusions

Our results may contribute to the optimal management of sheep genetic resources and inform modern selective breeding programmes that aim at mitigating future environmental challenges towards sustainable farming, while better balancing animal adaptation and productivity. Our results are directly relevant to the studied breed and the respective environmental conditions; however, the methodology may be extended to other livestock species of interest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00682-7.

Breeding Strategies for Weather Resilience in Small Ruminants in Atlantic and Mediterranean Climates

Many efforts are being made to cope with negative consequences of climate change (CC) on livestock. Among them, selective breeding of resilient animals to CC is presented as an opportunity to maintain high levels of performance regardless of variation in weather. In the present work, we proposed a set of breeding strategies to improve weather resilience in dairy goats raised in north-western European Atlantic conditions and dairy sheep raised in Mediterranean conditions while improving production efficiency at the same time. Breeding strategies differed in the selection emphasis placed on resilience traits, ranging from 0 to 40% in the index. Simulations were carried out mimicking real breeding programs including: milk yield, length of productive life, age at first kidding and mastitis incidence in dairy goats and milk, fat and protein yields, and fertility for dairy sheep. Considering the particular climatic conditions in the two regions, the predicted future climate scenarios, and genetic correlations among breeding objectives, resilience was defined as stability to weather changes for dairy goats and as the ability to improve performance under heat stress for dairy sheep. A strategy giving a selection weight of 10 and 20% for goat and sheep resilience, respectively, resulted in the best overall genetic response in terms of both, production and resilience ability. Not considering resilience in breeding programs could lead to a major production loss in future climate scenarios, whereas putting too much emphasis on resilience would result in a limited progress in milk production.

Understanding the seasonality of performance resilience to climate volatility in Mediterranean dairy sheep

As future climate challenges become increasingly evident, enhancing performance resilience of farm animals may contribute to mitigation against adverse weather and seasonal variation, and underpin livestock farming sustainability. In the present study, we develop novel seasonal resilience phenotypes reflecting milk production changes to fluctuating weather. We evaluate the impact of calendar season (autumn, winter and spring) on animal performance resilience by analysing 420,534 milk records of 36,908 milking ewes of the Chios breed together with relevant meteorological data from eastern Mediterranean. We reveal substantial seasonal effects on resilience and significant heritable trait variation (h² = 0.03–0.17). Resilience to cold weather (10 °C) of animals that start producing milk in spring was under different genetic control compared to autumn and winter as exemplified by negative genetic correlations (− 0.09 to − 0.27). Animal resilience to hot weather (25 °C) was partially under the same genetic control with genetic correlations between seasons ranging from 0.43 to 0.86. We report both favourable and antagonistic associations between animal resilience and lifetime milk production, depending on calendar season and the desirable direction of genetic selection. Concluding, we emphasise on seasonal adaptation of animals to climate and the need to incorporate the novel seasonal traits in future selective breeding programmes.

Breeding strategies for animal resilience to weather variation in meat sheep

BACKGROUND

The alteration in weather patterns expected due to climate change will affect farm animal performance, probably resulting in lower quantity and quality of available products. A potential mitigation strategy would be to breed selected animals for enhanced resilience to climate change. In this context, resilience would reflect stable animal performance in spite of weather variation. The objectives of this study were to (i) derive and characterise novel animal resilience phenotypes, (ii) investigate their genetic profiles and (iii) assess the impact of integrating them in breeding strategies for genetic improvement in meat sheep.

RESULTS

Random regression models were used to jointly analyse live body weight measured in different time points throughout the growth phases of 4469 Scottish Blackface sheep and weather variables during the same period to derive novel resilience phenotypes. The genetic analysis of these phenotypes revealed significant genetic variance and heritability, and an antagonistic genetic correlation with some animal performance traits. Simulated breeding strategies demonstrated that a relative emphasis of 10% on resilience compared to other traits would enhance performance stability against weather volatility without compromising animal growth.

CONCLUSIONS

Novel resilience traits exhibited sufficient genetic variation to be amenable to genetic improvement with selective breeding and are recommended to be included in future breeding goals.

Genetic dissection of complex behaviour traits in German Shepherd dogs

A favourable genetic structure and diversity of behavioural features highlights the potential of dogs for studying the genetic architecture of behaviour traits. However, behaviours are complex traits, which have been shown to be influenced by numerous genetic and non-genetic factors, complicating their analysis. In this study, the genetic contribution to behaviour variation in German Shepherd dogs (GSDs) was analysed using genomic approaches. GSDs were phenotyped for behaviour traits using the established Canine Behavioural Assessment and Research Questionnaire (C-BARQ). Genome-wide association study (GWAS) and regional heritability mapping (RHM) approaches were employed to identify associations between behaviour traits and genetic variants, while accounting for relevant non-genetic factors. By combining these complementary methods we endeavoured to increase the power to detect loci with small effects. Several behavioural traits exhibited moderate heritabilities, with the highest identified for Human-directed playfulness, a trait characterised by positive interactions with humans. We identified several genomic regions associated with one or more of the analysed behaviour traits. Some candidate genes located in these regions were previously linked to behavioural disorders in humans, suggesting a new context for their influence on behaviour characteristics. Overall, the results support dogs as a valuable resource to dissect the genetic architecture of behaviour traits and also highlight the value of focusing on a single breed in order to control for background genetic effects and thus avoid limitations of between-breed analyses.

Genetic analysis of novel phenotypes for farm animal resilience to weather variability

BACKGROUND

Climate change is expected to have a negative impact on food availability. While most efforts have been directed to reducing greenhouse gas emissions, complementary strategies are necessary to control the detrimental effects of climate change on farm animal performance. The objective of this study was to develop novel animal resilience phenotypes using reaction norm slopes, and examine their genetic and genomic parameters. A closely monitored dairy goat population was used for this purpose.

RESULTS

Individual animals differed in their response to changing atmospheric temperature and a temperature-humidity index. Significant genetic variance and heritability estimates were derived for these animal resilience phenotypes. Furthermore, some resilience traits had a significant unfavourable genetic correlation with animal performance. Genome-wide association analyses identified several candidate genes related to animal resilience to environment change.

CONCLUSIONS

Heritable variation exists among dairy goats in their production response to fluctuating weather variables. Results may inform future breeding programmes aimed to ensure efficient animal performance under changing climatic conditions.

Quantitative Trait Loci Mapping for Lameness Associated Phenotypes in Holstein-Friesian Dairy Cattle

Lameness represents a significant challenge for the dairy cattle industry, resulting in economic losses and reduced animal health and welfare. The existence of underlying genomic variation for lameness associated traits has the potential to improve selection strategies by using genomic markers. Therefore, the aim of this study was to identify genomic regions and potential candidate genes associated with lameness traits. Lameness related lesions and digital cushion thickness were studied using records collected by our research team, farm records, and a combination of both. Genome-wide analyses were performed to identify significant genomic effects, and a combination of single SNP association analysis and regional heritability mapping was used to identify associated genomic regions. Significant genomic effects were identified for several lameness related traits: Two genomic regions were identified on chromosome 3 associated with digital dermatitis and interdigital hyperplasia, one genomic region on chromosome 23 associated with interdigital hyperplasia, and one genomic region on chromosome 2 associated with sole haemorrhage. Candidate genes in those regions are mainly related to immune response and fibroblast proliferation. Quantitative trait loci (QTL) identified in this study could enlighten the understanding of lameness pathogenesis, providing an opportunity to improve health and welfare in dairy cattle with the addition of these regions into selection programs.

Impact of Genetic Selection for Increased Cattle Resistance to Bovine Tuberculosis on Disease Transmission Dynamics

Bovine tuberculosis (bTB) poses a challenge to animal health and welfare worldwide. Presence of genetic variation in host resistance to Mycobacterium bovis infection makes the trait amenable to improvement with genetic selection. Genetic evaluations for resistance to infection in dairy cattle are currently available in the United Kingdom (UK), enabling genetic selection of more resistant animals. However, the extent to which genetic selection could contribute to bTB eradication is unknown. The objective of this study was to quantify the impact of genetic selection for bTB resistance on cattle-to-cattle disease transmission dynamics and prevalence by developing a stochastic genetic epidemiological model. The model was used to implement genetic selection in a simulated cattle population. The model considered various levels of selection intensity over 20 generations assuming genetic heterogeneity in host resistance to infection. Our model attempted to represent the dairy cattle population structure and current bTB control strategies in the UK, and was informed by genetic and epidemiological parameters inferred from data collected from UK bTB infected dairy herds. The risk of a bTB breakdown was modeled as the percentage of herds where initially infected cows (index cases) generated secondary cases by infecting herd-mates. The model predicted that this risk would be reduced by half after 4, 6, 9, and 15 generations for selection intensities corresponding to genetic selection of the 10, 25, 50, and 70% most resistant sires, respectively. In herds undergoing bTB breakdowns, genetic selection reduced the severity of breakdowns over generations by reducing both the percentage of secondary cases and the duration over which new secondary cases were detected. Selection of the 10, 25, 50, and 70% most resistant sires reduced the percentage of secondary cases to <1% in 4, 5, 7, and 11 generations, respectively. Similarly, the proportion of long breakdowns (breakdowns in which secondary cases were detected for more than 365 days) was reduced by half in 2, 2, 3, and 4 generations, respectively. Collectively, results suggest that genetic selection could be a viable tool that can complement existing management and surveillance methods to control and ultimately eradicate bTB.

Joint Genomic Prediction of Canine Hip Dysplasia in UK and US Labrador Retrievers

Canine hip dysplasia, a debilitating orthopedic disorder that leads to osteoarthritis and cartilage degeneration, is common in several large-sized dog breeds and shows moderate heritability suggesting that selection can reduce prevalence. Estimating genomic breeding values require large reference populations, which are expensive to genotype for development of genomic prediction tools. Combining datasets from different countries could be an option to help build larger reference datasets without incurring extra genotyping costs. Our objective was to evaluate genomic prediction based on a combination of UK and US datasets of genotyped dogs with records of Norberg angle scores, related to canine hip dysplasia. Prediction accuracies using a single population were 0.179 and 0.290 for 1,179 and 242 UK and US Labrador Retrievers, respectively. Prediction accuracies changed to 0.189 and 0.260, with an increased bias of genomic breeding values when using a joint training set (biased upwards for the US population and downwards for the UK population). Our results show that in this study of canine hip dysplasia, little or no benefit was gained from using a joint training set as compared to using a single population as training set. We attribute this to differences in the genetic background of the two populations as well as the small sample size of the US dataset.

The Genomic Architecture of Fowl Typhoid Resistance in Commercial Layers

Salmonella enterica serovar Gallinarum causes devastating outbreaks of fowl typhoid across the globe, especially in developing countries. With the use of antimicrobial agents being reduced due to legislation and the absence of licensed vaccines in some parts of the world, an attractive complementary control strategy is to breed chickens for increased resistance to Salmonella. The potential for genetic control of salmonellosis has been demonstrated by experimental challenge of inbred populations. Quantitative trait loci (QTL) associated with resistance have been identified in many genomic regions. A major QTL associated with systemic salmonellosis has been identified in a region termed SAL1. In the present study, two outbreaks of fowl typhoid in 2007 and 2012 in the United Kingdom were used to investigate the genetic architecture of Salmonella resistance in commercial laying hens. In the first outbreak 100 resistant and 150 susceptible layers were genotyped using 11 single nucleotide polymorphism (SNP) and 3 microsatellite markers located in the previously identified SAL1 region on chromosome 5. From the second outbreak 100 resistant and 200 susceptible layers, belonging to a different line, were genotyped with a high-density (600 K) genome-wide SNP array. Substantial heritability estimates were obtained in both populations (h ² = 0.22 and 0.26, for the layers in the first and second outbreak, respectively). Significant associations with three markers on chromosome 5 located close to AKT1 and SIVA1 genes, coding for RAC-alpha serine/threonine protein kinase, and the CD27-binding protein SIVA1, respectively, were identified in the first outbreak. From analysis of the second outbreak, eight genome-wide significant associations with Salmonella resistance were identified on chromosomes 1, 6, 7, 11, 23, 24, 26, 28 and several others with suggestive genome-wide significance were found. Pathway and network analysis revealed the presence of many innate immune pathways related to Salmonella resistance. Although, significant associations with SNPs located in the SAL1 locus were not identified by the genome-wide scan for layers from the second outbreak, pathway analysis revealed P13K/AKT signaling as the most significant pathway. In summary, resistance to fowl typhoid is a heritable polygenic trait that could possibly be enhanced through selective breeding.

Genomic data illuminates demography, genetic structure and selection of a popular dog breed

BACKGROUND

Genomic methods have proved to be important tools in the analysis of genetic diversity across the range of species and can be used to reveal processes underlying both short- and long-term evolutionary change. This study applied genomic methods to investigate population structure and inbreeding in a common UK dog breed, the Labrador Retriever.

RESULTS

We found substantial within-breed genetic differentiation, which was associated with the role of the dog (i.e. working, pet, show) and also with coat colour (i.e. black, yellow, brown). There was little evidence of geographical differentiation. Highly differentiated genomic regions contained genes and markers associated with skull shape, suggesting that at least some of the differentiation is related to human-imposed selection on this trait. We also found that the total length of homozygous segments (runs of homozygosity, ROHs) was highly correlated with inbreeding coefficient.

CONCLUSIONS

This study demonstrates that high-density genomic data can be used to quantify genetic diversity and to decipher demographic and selection processes. Analysis of genetically differentiated regions in the UK Labrador Retriever population suggests the possibility of human-imposed selection on craniofacial characteristics. The high correlation between estimates of inbreeding from genomic and pedigree data for this breed demonstrates that genomic approaches can be used to quantify inbreeding levels in dogs, which will be particularly useful where pedigree information is missing.

Genomic regions underlying susceptibility to bovine tuberculosis in Holstein-Friesian cattle

BACKGROUND

The significant social and economic loss as a result of bovine tuberculosis (bTB) presents a continuous challenge to cattle industries in the UK and worldwide. However, host genetic variation in cattle susceptibility to bTB provides an opportunity to select for resistant animals and further understand the genetic mechanisms underlying disease dynamics.

METHODS

The present study identified genomic regions associated with susceptibility to bTB using genome-wide association (GWA), regional heritability mapping (RHM) and chromosome association approaches. Phenotypes comprised de-regressed estimated breeding values of 804 Holstein-Friesian sires and pertained to three bTB indicator traits: i) positive reactors to the skin test with positive post-mortem examination results (phenotype 1); ii) positive reactors to the skin test regardless of post-mortem examination results (phenotype 2) and iii) as in (ii) plus non-reactors and inconclusive reactors to the skin tests with positive post-mortem examination results (phenotype 3). Genotypes based on the 50 K SNP DNA array were available and a total of 34,874 SNPs remained per animal after quality control.

RESULTS

The estimated polygenic heritability for susceptibility to bTB was 0.26, 0.37 and 0.34 for phenotypes 1, 2 and 3, respectively. GWA analysis identified a putative SNP on Bos taurus autosomes (BTA) 2 associated with phenotype 1, and another on BTA 23 associated with phenotype 2. Genomic regions encompassing these SNPs were found to harbour potentially relevant annotated genes. RHM confirmed the effect of these genomic regions and identified new regions on BTA 18 for phenotype 1 and BTA 3 for phenotypes 2 and 3. Heritabilities of the genomic regions ranged between 0.05 and 0.08 across the three phenotypes. Chromosome association analysis indicated a major role of BTA 23 on susceptibility to bTB.

CONCLUSION

Genomic regions and candidate genes identified in the present study provide an opportunity to further understand pathways critical to cattle susceptibility to bTB and enhance genetic improvement programmes aiming at controlling and eradicating the disease.

A practical approach to detect ancestral haplotypes in livestock populations

Background

The effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle.

Results

The identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content.

Conclusion

In agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0405-2) contains supplementary material, which is available to authorized users.

Genomic-Based Optimum Contribution in Conservation and Genetic Improvement Programs with Antagonistic Fitness and Productivity Traits

Animal selection for genetic improvement of productivity may lead to an increase in inbreeding through the use of techniques that enhance the reproductive capability of selected animals. Therefore, breeding strategies aim to balance maintaining genetic variability and acceptable fitness levels with increasing productivity. The present study demonstrates the effectiveness of genomic-based optimum contribution strategies at addressing this objective when fitness and productivity are genetically antagonistic traits. Strategies are evaluated in directional selection (increasing productivity) or conservation (maintaining fitness) scenarios. In the former case, substantial rates of genetic gain can be achieved while greatly constraining the rate of increase in inbreeding. Under a conservation approach, inbreeding depression can be effectively halted while also achieving a modest rate of genetic gain for productivity. Furthermore, the use of optimum contribution strategies when combined with a simple non-random mating scheme (minimum kinship method) showed an additional delay in the increase of inbreeding in the short term. In conclusion, genomic-based optimum contribution methods can be effectively used to control inbreeding and inbreeding depression, and still allow genetic gain for productivity traits even when fitness and productivity are antagonistically correlated.

Genomic prediction of traits related to canine hip dysplasia

Increased concern for the welfare of pedigree dogs has led to development of selection programs against inherited diseases. An example is canine hip dysplasia (CHD), which has a moderate heritability and a high prevalence in some large-sized breeds. To date, selection using phenotypes has led to only modest improvement, and alternative strategies such as genomic selection (GS) may prove more effective. The primary aims of this study were to compare the performance of pedigree- and genomic-based breeding against CHD in the UK Labrador retriever population and to evaluate the performance of different GS methods. A sample of 1179 Labrador Retrievers evaluated for CHD according to the UK scoring method (hip score, HS) was genotyped with the Illumina CanineHD BeadChip. Twelve functions of HS and its component traits were analyzed using different statistical methods (GBLUP, Bayes C and Single-Step methods), and results were compared with a pedigree-based approach (BLUP) using cross-validation. Genomic methods resulted in similar or higher accuracies than pedigree-based methods with training sets of 944 individuals for all but the untransformed HS, suggesting that GS is an effective strategy. GBLUP and Bayes C gave similar prediction accuracies for HS and related traits, indicating a polygenic architecture. This conclusion was also supported by the low accuracies obtained in additional GBLUP analyses performed using only the SNPs with highest test statistics, also indicating that marker-assisted selection (MAS) would not be as effective as GS. A Single-Step method that combines genomic and pedigree information also showed higher accuracy than GBLUP and Bayes C for the log-transformed HS, which is currently used for pedigree based evaluations in UK. In conclusion, GS is a promising alternative to pedigree-based selection against CHD, requiring more phenotypes with genomic data to improve further the accuracy of prediction.

Quantitative trait loci mapping for canine hip dysplasia and its related traits in UK Labrador Retrievers

Background

Canine hip dysplasia (CHD) is characterised by a malformation of the hip joint, leading to osteoarthritis and lameness. Current breeding schemes against CHD have resulted in measurable but moderate responses. The application of marker-assisted selection, incorporating specific markers associated with the disease, or genomic selection, incorporating genome-wide markers, has the potential to dramatically improve results of breeding schemes. Our aims were to identify regions associated with hip dysplasia or its related traits using genome and chromosome-wide analysis, study the linkage disequilibrium (LD) in these regions and provide plausible gene candidates. This study is focused on the UK Labrador Retriever population, which has a high prevalence of the disease and participates in a recording program led by the British Veterinary Association (BVA) and The Kennel Club (KC).

Results

Two genome-wide and several chromosome-wide QTLs affecting CHD and its related traits were identified, indicating regions related to hip dysplasia.

Conclusion

Consistent with previous studies, the genetic architecture of CHD appears to be based on many genes with small or moderate effect, suggesting that genomic selection rather than marker-assisted selection may be an appropriate strategy for reducing this disease.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-833) contains supplementary material, which is available to authorized users.

Assessing the impact of genomic selection against hip dysplasia in the Labrador Retriever dog

Many purebred dogs exhibit a higher prevalence of inherited diseases compared with non-purebred dogs. One of the most popular breeds in the UK is the Labrador Retriever, which has a high prevalence of hip dysplasia resulting in high costs for surgical operations and impaired animal welfare. Considering the many complications of highly managed populations, mainly due to breeder's conventions and the resulting population structure, is of great importance for the proper development of a strategy against the disease. In this study, we have compared the utilities and performances of both genomic and phenotypic selection against hip dysplasia in a simulated population with the characteristics of the British Veterinary Association and Kennel Club (BV /KC) hip dysplasia scheme. The results confirm the potential benefits of genomic selection by showing a moderate increase of 1.15-fold (assuming a realistic accuracy of r(2) = 0.5) in response to selection due to the higher accuracy (between 0.96- and 1.32-fold, considering 0.35 ≤ r(2) ≤ 0.7) and more than a threefold increase when all the offspring in each litter are tested (between 3.25- and 4.55-fold, again considering 0.35 ≤ r(2) ≤ 0.7).

Efficiency of conservation management methods for subdivided populations under local adaptation

Computer simulations were used to investigate the efficiency of management methods for the conservation of a structured population when local adaptation exists. A subdivided population, with subpopulations adapted to different optima for a quantitative trait under stabilizing selection, was managed in order to maintain the highest genetic diversity in a 10-generation period. Two procedures were compared. For the first, minimum coancestry contributions were carried out independently for each subpopulation, and random migration of individuals was accomplished thereafter. For the second, minimum coancestry contributions from individuals were globally implemented, including an optimal migration design. This optimal method can be adjusted to control local inbreeding to different extents. Adaptation to local optima implies a reduction in the efficiency of the management methods because of the effective failure in the established migrations. For strong selection, the optimal design can be very inefficient, even more than the random migration scheme because the intended migrants have usually low fitness in the recipient subpopulations. However, for more realistic moderate or weak selection, the optimal method is more efficient than random migration, especially if inbreeding depression on fitness is also taken into account. It is concluded that the optimal management method can be recommended in conservation programs with local adaptation of subpopulations, but this issue should be accounted for when designing the management strategies.

Detection of growth-related QTL in turbot (Scophthalmus maximus)

Background

The turbot (Scophthalmus maximus) is a highly appreciated European aquaculture species. Growth related traits constitute the main goal of the ongoing genetic breeding programs of this species. The recent construction of a consensus linkage map in this species has allowed the selection of a panel of 100 homogeneously distributed markers covering the 26 linkage groups (LG) suitable for QTL search. In this study we addressed the detection of QTL with effect on body weight, length and Fulton's condition factor.

Results

Eight families from two genetic breeding programs comprising 814 individuals were used to search for growth related QTL using the panel of microsatellites available for QTL screening. Two different approaches, maximum likelihood and regression interval mapping, were used in order to search for QTL. Up to eleven significant QTL were detected with both methods in at least one family: four for weight on LGs 5, 14, 15 and 16; five for length on LGs 5, 6, 12, 14 and 15; and two for Fulton's condition factor on LGs 3 and 16. In these LGs an association analysis was performed to ascertain the microsatellite marker with the highest apparent effect on the trait, in order to test the possibility of using them for marker assisted selection.

Conclusions

The use of regression interval mapping and maximum likelihood methods for QTL detection provided consistent results in many cases, although the high variation observed for traits mean among families made it difficult to evaluate QTL effects. Finer mapping of detected QTL, looking for tightly linked markers to the causative mutation, and comparative genomics are suggested to deepen in the analysis of QTL in turbot so they can be applied in marker assisted selection programs.

The build up of mutation-selection- drift balance in laboratory Drosophila populations

The build up of an equilibrium between mutation, selection, and drift in populations of moderate size is an important evolutionary issue, and can be critical in the conservation of endangered populations. We studied this process in two Drosophila melanogaster populations initially lacking genetic variability (C1 and C2) that were subsequently maintained during 431 or 165 generations with effective population size N(e) approximately 500 (estimated by lethal complementation analysis). Each population originated synchronously to a companion set of full-sib mutation accumulation (MA) lines, C1 and MA1 were derived from an isogenic origin and C2 and MA2 from a single MA1 line at generation 265. The results suggest that both C1 and C2 populations were close to the mutation-selection-drift balance for viability and bristle traits, and are consistent with a 2.5-fold increase of the mutation rate in C2 and MA2. Despite this increase, the average panmictic viability in C2 was only slightly below that of C1, indicating that the expressed loads due to segregating deleterious mutation were small, in agreement with the low deleterious mutation rate (0.015-0.045) previously reported for the MA1 lines. In C1, the nonlethal inbreeding depression rate for viability was 30% of that usually estimated in segregating populations. The genetic variance for bristles regenerated in C1 and C2 was moderately smaller than the average value reported for natural populations, implying that they have accumulated a substantial adaptive potential. In light of neutral and selective predictions, these results suggest that bristle additive variance was predominantly due to segregation of mutations with deleterious effects of the order of 10(-3), and is consistent with relatively weak causal stabilizing selection (V(s) approximately 30).