Genome‐wide association study of insect bite hypersensitivity in Dutch Shetland pony mares

Insect bite hypersensitivity (IBH) is the most common allergic disease present in horses worldwide. It has been shown that IBH is under genetic control, but the knowledge of associated genes is limited. We conducted a genome‐wide association study to identify and quantify genomic regions contributing to IBH in the Dutch Shetland pony population. A total of 97 cases and 91 controls were selected and matched on withers height, coat colour and pedigree to minimise the population stratification. A blood sample was collected from participating Shetland pony mares, their IBH phenotype was scored and the owner filled in a questionnaire. A total of 40 021 single‐nucleotide polymorphisms (SNPs) were fitted in a univariable logistic model fitting an additive effect. Analysis revealed no effects of population stratification. Significant associations with IBH were detected for 24 SNPs on 12 chromosomes [?log₁₀(P‐value) > 2.5]. Odds ratios of allele substitution effects of the unfavourable allele were between 1.94 and 5.95. The most significant SNP was found on chromosome 27, with an odds ratio of 2.31 and with an allele frequency of the unfavourable allele of 0.72 in cases and 0.53 in controls. Genome‐wide association studies on additional horse populations are desired to validate the identified associations, to identify the genes involved in IBH and to develop genomic tools to decrease IBH prevalence.     

Genome-wide association study for calving performance using high-density genotypes in dairy and beef cattle

Background

Calving difficulty and perinatal mortality are prevalent in modern-day cattle production systems. It is well-established that there is a genetic component to both traits, yet little is known about their underlying genomic architecture, particularly in beef breeds. Therefore, we performed a genome-wide association study using high-density genotypes to elucidate the genomic architecture of these traits and to identify regions of the bovine genome associated with them.

Results

Genomic regions associated with calving difficulty (direct and maternal) and perinatal mortality were detected using two statistical approaches: (1) single-SNP (single nucleotide polymorphism) regression and (2) a Bayesian approach. Data included high-density genotypes on 770 Holstein-Friesian, 927 Charolais and 963 Limousin bulls. Several novel or previously identified genomic regions were detected but associations differed by breed. For example, two genomic associations, one each on chromosomes 18 and 2 explained 2.49 % and 3.13 % of the genetic variance in direct calving difficulty in the Holstein-Friesian and Charolais populations, respectively. Imputed Holstein-Friesian sequence data was used to refine the genomic regions responsible for significant associations. Several candidate genes on chromosome 18 were identified and four highly significant missense variants were detected within three of these genes (SIGLEC12, CTU1, and ZNF615). Nevertheless, only CTU1 contained a missense variant with a putative impact on direct calving difficulty based on SIFT (0.06) and Polyphen (0.95) scores. Using imputed sequence data, we refined a genomic region on chromosome 4 associated with maternal calving difficulty in the Holstein-Friesian population and found the strongest association with an intronic variant in the PCLO gene. A meta-analysis was performed across the three breeds for each calving performance trait to identify common variants associated with these traits in the three breeds. Our results suggest that a portion of the genetic variation in calving performance is common to all three breeds.

Conclusion

The genomic architecture of calving performance is complex and mainly influenced by many polymorphisms of small effect. We identified several associations of moderate effect size but the majority were breed-specific, indicating that breed-specific alleles exist for calving performance or that the linkage phase between genotyped allele and causal mutation varies between breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0126-4) contains supplementary material, which is available to authorized users.     

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

Background

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F_ST. The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

Results

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F_ST, potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F_ST signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F_ST values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

Conclusions

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Electronic supplementary material

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

Short Copy Number Variations Potentially Associated with Tonic Immobility Responses in Newly Hatched Chicks

Introduction

Tonic immobility (TI) is fear-induced freezing that animals may undergo when confronted by a threat. It is principally observed in prey species as defence mechanisms. In our preliminary research, we detected large inter-individual variations in the frequency and duration of freezing behavior among newly hatched domestic chicks (Gallus gallus). In this study we aim to identify the copy number variations (CNVs) in the genome of chicks as genetic candidates that underlie the behavioral plasticity to fearful stimuli.

Methods

A total of 110 domestic chicks were used for an association study between TI responses and copy number polymorphisms. Array comparative genomic hybridization (aCGH) was conducted between chicks with high and low TI scores using an Agilent 4×180 custom microarray. We specifically focused on 3 genomic regions (>60 Mb) of chromosome 1 where previous quantitative trait loci (QTL) analysis showed significant F-values for fearful responses.

Results

ACGH successfully detected short CNVs within the regions overlapping 3 QTL peaks. Eleven of these identified loci were validated by real-time quantitative polymerase chain reaction (qPCR) as copy number polymorphisms. Although there wkas no significant p value in the correlation analysis between TI scores and the relative copy number within each breed, several CNV loci showed significant differences in the relative copy number between 2 breeds of chicken (White Leghorn and Nagoya) which had different quantitative characteristics of fear-induced responses.

Conclusion

Our data shows the potential CNVs that may be responsible for innate fear response in domestic chicks.     

Identification and functional characterization of copy number variations in diverse chicken breeds

Background

The detection and functional characterization of genomic structural variations are important for understanding the landscape of genetic variation in the chicken. A recently recognized aspect of genomic structural variation, called copy number variation (CNV), is gaining interest in chicken genomic studies. The aim of the present study was to investigate the pattern and functional characterization of CNVs in five characteristic chicken breeds, which will be important for future studies associating phenotype with chicken genome architecture.

Results

Using a commercial 385 K array-based comparative genomic hybridization (aCGH) genome array, we performed CNV discovery using 10 chicken samples from four local Chinese breeds and the French breed Houdan chicken. The female Anka broiler was used as a reference. A total of 281 copy number variation regions (CNVR) were identified, covering 12.8 Mb of polymorphic sequences or 1.07% of the entire chicken genome. The functional annotation of CNVRs indicated that these regions completely or partially overlapped with 231 genes and 1032 quantitative traits loci, suggesting these CNVs have important functions and might be promising resources for exploring differences among various breeds. In addition, we employed quantitative PCR (qPCR) to further validate several copy number variable genes, such as prolactin receptor, endothelin 3 (EDN3), suppressor of cytokine signaling 2, CD8a molecule, with important functions, and the results suggested that EDN3 might be a molecular marker for the selection of dark skin color in poultry production. Moreover, we also identified a new CNVR (chr24: 3484617–3512275), encoding the sortilin-related receptor gene, with copy number changes in only black-bone chicken.

Conclusions

Here, we report a genome-wide analysis of the CNVs in five chicken breeds using aCGH. The association between EDN3 and melanoblast proliferation was further confirmed using qPCR. These results provide additional information for understanding genomic variation and related phenotypic characteristics.

Electronic supplementary material

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

Relative extended haplotype homozygosity signals across breeds reveal dairy and beef specific signatures of selection

Background

A number of methods are available to scan a genome for selection signatures by evaluating patterns of diversity within and between breeds. Among these, “extended haplotype homozygosity” (EHH) is a reliable approach to detect genome regions under recent selective pressure. The objective of this study was to use this approach to identify regions that are under recent positive selection and shared by the most representative Italian dairy and beef cattle breeds.

Results

A total of 3220 animals from Italian Holstein (2179), Italian Brown (775), Simmental (493), Marchigiana (485) and Piedmontese (379) breeds were genotyped with the Illumina BovineSNP50 BeadChip v.1. After standard quality control procedures, genotypes were phased and core haplotypes were identified. The decay of linkage disequilibrium (LD) for each core haplotype was assessed by measuring the EHH. Since accurate estimates of local recombination rates were not available, relative EHH (rEHH) was calculated for each core haplotype. Genomic regions that carry frequent core haplotypes and with significant rEHH values were considered as candidates for recent positive selection. Candidate regions were aligned across to identify signals shared by dairy or beef cattle breeds. Overall, 82 and 87 common regions were detected among dairy and beef cattle breeds, respectively. Bioinformatic analysis identified 244 and 232 genes in these common genomic regions. Gene annotation and pathway analysis showed that these genes are involved in molecular functions that are biologically related to milk or meat production.

Conclusions

Our results suggest that a multi-breed approach can lead to the identification of genomic signatures in breeds of cattle that are selected for the same production goal and thus to the localisation of genomic regions of interest in dairy and beef production.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0113-9) contains supplementary material, which is available to authorized users.     

Next generation sequencing gives an insight into the characteristics of highly selected breeds versus non-breed horses in the course of domestication

Background

Domestication has shaped the horse and lead to a group of many different types. Some have been under strong human selection while others developed in close relationship with nature. The aim of our study was to perform next generation sequencing of breed and non-breed horses to provide an insight into genetic influences on selective forces.

Results

Whole genome sequencing of five horses of four different populations revealed 10,193,421 single nucleotide polymorphisms (SNPs) and 1,361,948 insertion/deletion polymorphisms (indels). In comparison to horse variant databases and previous reports, we were able to identify 3,394,883 novel SNPs and 868,525 novel indels. We analyzed the distribution of individual variants and found significant enrichment of private mutations in coding regions of genes involved in primary metabolic processes, anatomical structures, morphogenesis and cellular components in non-breed horses and in contrast to that private mutations in genes affecting cell communication, lipid metabolic process, neurological system process, muscle contraction, ion transport, developmental processes of the nervous system and ectoderm in breed horses.

Conclusions

Our next generation sequencing data constitute an important first step for the characterization of non-breed in comparison to breed horses and provide a large number of novel variants for future analyses. Functional annotations suggest specific variants that could play a role for the characterization of breed or non-breed horses.

Electronic supplementary material

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

Detection of selection signatures in dairy and beef cattle using high-density genomic information

Background

Artificial selection for economically important traits in cattle is expected to have left distinctive selection signatures on the genome. Access to high-density genotypes facilitates the accurate identification of genomic regions that have undergone positive selection. These findings help to better elucidate the mechanisms of selection and to identify candidate genes of interest to breeding programs.

Results

Information on 705 243 autosomal single nucleotide polymorphisms (SNPs) in 3122 dairy and beef male animals from seven cattle breeds (Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental) were used to detect selection signatures by applying two complementary methods, integrated haplotype score (iHS) and global fixation index (F_ST). To control for false positive results, we used false discovery rate (FDR) adjustment to calculate adjusted iHS within each breed and the genome-wide significance level was about 0.003. Using the iHS method, 83, 92, 91, 101, 85, 101 and 86 significant genomic regions were detected for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental cattle, respectively. None of these regions was common to all seven breeds. Using the F_ST approach, 704 individual SNPs were detected across breeds. Annotation of the regions of the genome that showed selection signatures revealed several interesting candidate genes i.e. DGAT1, ABCG2, MSTN, CAPN3, FABP3, CHCHD7, PLAG1, JAZF1, PRKG2, ACTC1, TBC1D1, GHR, BMP2, TSG1, LYN, KIT and MC1R that play a role in milk production, reproduction, body size, muscle formation or coat color. Fifty-seven common candidate genes were found by both the iHS and global F_ST methods across the seven breeds. Moreover, many novel genomic regions and genes were detected within the regions that showed selection signatures; for some candidate genes, signatures of positive selection exist in the human genome. Multilevel bioinformatic analyses of the detected candidate genes suggested that the PPAR pathway may have been subjected to positive selection.

Conclusions

This study provides a high-resolution bovine genomic map of positive selection signatures that are either specific to one breed or common to a subset of the seven breeds analyzed. Our results will contribute to the detection of functional candidate genes that have undergone positive selection in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0127-3) contains supplementary material, which is available to authorized users.     

Heart rate and heart rate variability in pregnant warmblood and Shetland mares as well as their fetuses

Heart rate (HR) is an important parameter of fetal well-being. In horses, HR and heart rate variability (HRV) can be determined by fetomaternal electrocardiography (ECG) from mid-pregnancy to foaling. Normal values for physiological parameters in larger breeds are often used as reference values in ponies. However, HR increases with decreasing size of the animal and in ponies is higher than in warmblood horses. It is not known if fetal HR is affected by breed and if values obtained in larger breeds can be used to assess Shetland fetuses. We have determined fetomaternal beat-to-beat (RR) interval (inversely correlated to HR) and HRV in warmblood (n=6) and Shetland pregnancies (n=7) at days 280 and 300 of gestation by ECG. Maternal RR interval was lower in pony than in warmblood mares (day 280: Shetland: 958±110, warmblood: 1489±126ms, p<0.01) The SDRR (standard deviation of RR interval) and the RMSSD (root mean square of successive RR differences) did not differ between breeds at any time. Also RR interval as well as HRV did not differ between warmblood and pony fetuses (RR interval day 280: Shetland: 606±39, warmblood: 589±38ms). In conclusion, although maternal RR interval is clearly higher in Shetland than in warmblood mares, fetal RR interval in the two breeds is on the same level.     

Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes

Background

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Results

Divergent lines of Romney and Perendale sheep, selectively bred for high and low faecal nematode egg count, were genotyped using the Illumina® Ovine SNP50 BeadChip. The resulting genome-wide SNP data were analysed for selective sweeps on loci associated with resistance or susceptibility to gastrointestinal nematode infection. Population differentiation using F_ST and Peddrift revealed sixteen regions, which included candidate genes involved in chitinase activity and the cytokine response. Two of the sixteen regions identified were contained within previously identified QTLs associated with nematode resistance.

Conclusions

In this study we identified fourteen novel regions associated with resistance or susceptibility to gastrointestinal nematodes. Results from this study support the hypothesis that host resistance to internal nematode parasites is likely to be controlled by a number of loci of moderate to small effects.

Electronic supplementary material

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

Differential Gene Repertoire in Mycobacterium ulcerans Identifies Candidate Genes for Patho-Adaptation

Background

Based on large genomic sequence polymorphisms, several haplotypes belonging to two major lineages of the human pathogen Mycobacterium ulcerans could be distinguished among patient isolates from various geographic origins. However, the biological relevance of insertional/deletional diversity is not understood.

Methodology

Using comparative genomics, we have investigated the genes located in regions of difference recently identified by DNA microarray based hybridisation analysis. The analysed regions of difference comprise ∼7% of the entire M. ulcerans genome.

Principal Findings

Several different mechanisms leading to loss of functional genes were identified, ranging from pseudogenization, caused by frame shift mutations or mobile genetic element interspersing, to large sequence polymorphisms. Four hot spot regions for genetic instability were unveiled. Altogether, 229 coding sequences were found to be differentially inactivated, constituting a repertoire of coding sequence variation in the rather monomorphic M. ulcerans.

Conclusions/Significance

The differential gene inactivation patterns associated with the M. ulcerans haplotypes identified candidate genes that may confer enhanced adaptation upon ablation of expression. A number of gene conversions confined to the classical lineage may contribute to particular virulence of this group comprising isolates from Africa and Australia. Identification of this spectrum of anti-virulence gene candidates expands our understanding of the pathogenicity and ecology of the emerging infectious disease Buruli ulcer.     

Recombination locations and rates in beef cattle assessed from parent-offspring pairs

Background

Recombination events tend to occur in hotspots and vary in number among individuals. The presence of recombination influences the accuracy of haplotype phasing and the imputation of missing genotypes. Genes that influence genome-wide recombination rate have been discovered in mammals, yeast, and plants. Our aim was to investigate the influence of recombination on haplotype phasing, locate recombination hotspots, scan the genome for Quantitative Trait Loci (QTL) and identify candidate genes that influence recombination, and quantify the impact of recombination on the accuracy of genotype imputation in beef cattle.

Methods

2775 Angus and 1485 Limousin parent-verified sire/offspring pairs were genotyped with the Illumina BovineSNP50 chip. Haplotype phasing was performed with DAGPHASE and BEAGLE using UMD3.1 assembly SNP (single nucleotide polymorphism) coordinates. Recombination events were detected by comparing the two reconstructed chromosomal haplotypes inherited by each offspring with those of their sires. Expected crossover probabilities were estimated assuming no interference and a binomial distribution for the frequency of crossovers. The BayesB approach for genome-wide association analysis implemented in the GenSel software was used to identify genomic regions harboring QTL with large effects on recombination. BEAGLE was used to impute Angus genotypes from a 7K subset to the 50K chip.

Results

DAGPHASE was superior to BEAGLE in haplotype phasing, which indicates that linkage information from relatives can improve its accuracy. The estimated genetic length of the 29 bovine autosomes was 3097 cM, with a genome-wide recombination distance averaging 1.23 cM/Mb. 427 and 348 windows containing recombination hotspots were detected in Angus and Limousin, respectively, of which 166 were in common. Several significant SNPs and candidate genes, which influence genome-wide recombination were localized in QTL regions detected in the two breeds. High-recombination rates hinder the accuracy of haplotype phasing and genotype imputation.

Conclusions

Small population sizes, inadequate half-sib family sizes, recombination, gene conversion, genotyping errors, and map errors reduce the accuracy of haplotype phasing and genotype imputation. Candidate regions associated with recombination were identified in both breeds. Recombination analysis may improve the accuracy of haplotype phasing and genotype imputation from low- to high-density SNP panels.     

Whole-genome sequencing identifies a homozygous deletion encompassing exons 17 to 23 of the integrin beta 4 gene in a Charolais calf with junctional epidermolysis bullosa

Background

Since 2010, four Charolais calves with a congenital mechanobullous skin disorder that were born in the same herd from consanguineous matings were reported to us. Clinical and histopathological examination revealed lesions that are compatible with junctional epidermolysis bullosa (JEB).

Results

Fifty-four extended regions of homozygosity (>1 Mb) were identified after analysing the whole-genome sequencing (WGS) data from the only case available for DNA sampling at the beginning of the study. Filtering of variants located in these regions for (i) homozygous polymorphisms observed in the WGS data from eight healthy Charolais animals and (ii) homozygous or heterozygous polymorphisms found in the genomes of 234 animals from different breeds did not reveal any deleterious candidate SNPs (single nucleotide polymorphisms) or small indels. Subsequent screening for structural variants in candidate genes located in the same regions identified a homozygous deletion that includes exons 17 to 23 of the integrin beta 4 (ITGB4), a gene that was previously associated with the same defect in humans. Genotyping of a second case and of six parents of affected calves (two sires and four dams) revealed a perfect association between this mutation and the assumed genotypes of the individuals. Mining of Illumina BovineSNP50 Beadchip genotyping data from 6870 Charolais cattle detected only 44 heterozygous animals for a 5.6-Mb haplotype around ITGB4 that was shared with the carriers of the mutation. Interestingly, none of the 16 animals genotyped for the deletion carried the mutation, which suggests a rather recent origin for the mutation.

Conclusions

In conclusion, we successfully identified the causative mutation for a very rare autosomal recessive mutation with only one case by exploiting the most recent DNA sequencing technologies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0110-z) contains supplementary material, which is available to authorized users.     

Genetic variants in mammary development,prolactin signalling and involution pathways explain considerable variation in bovine milk production and milk composition

Background

The maintenance of lactation in mammals is the result of a balance between competing signals from mammary development, prolactin signalling and involution pathways. Dairy cattle are an interesting case study to investigate the effect of polymorphisms that affect the function of genes in these pathways. In dairy cattle, lactation yields and milk composition (for example protein percentage and fat percentage) are routinely recorded, and these vary greatly between individuals. In this study, we test 8058 single nucleotide polymorphisms in or close to genes in these pathways for association with milk production traits and determine the proportion of variance explained by each pathway, using data on 16 812 dairy cattle, including Holstein-Friesian and Jersey bulls and cows.

Results

Single nucleotide polymorphisms close to genes in the mammary development, prolactin signalling and involution pathways were significantly associated with milk production traits. The involution pathway explained the largest proportion of genetic variation for production traits. The mammary development pathway also explained additional genetic variation for milk volume, fat percentage and protein percentage.

Conclusions

Genetic variants in the involution pathway explained considerably more genetic variation in milk production traits than expected by chance. Many of the associations for single nucleotide polymorphisms in genes in this pathway have not been detected in conventional genome-wide association studies. The pathway approach used here allowed us to identify some novel candidates for further studies that will be aimed at refining the location of associated genomic regions and identifying polymorphisms contributing to variation in lactation volume and milk composition.     

Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA

Background

Artificial selection has caused rapid evolution in domesticated species. The identification of selection footprints across domesticated genomes can contribute to uncover the genetic basis of phenotypic diversity.

Methodology/Main Findings

Genome wide footprints of pig domestication and selection were identified using massive parallel sequencing of pooled reduced representation libraries (RRL) representing ∼2% of the genome from wild boar and four domestic pig breeds (Large White, Landrace, Duroc and Pietrain) which have been under strong selection for muscle development, growth, behavior and coat color. Using specifically developed statistical methods that account for DNA pooling, low mean sequencing depth, and sequencing errors, we provide genome-wide estimates of nucleotide diversity and genetic differentiation in pig. Widespread signals suggestive of positive and balancing selection were found and the strongest signals were observed in Pietrain, one of the breeds most intensively selected for muscle development. Most signals were population-specific but affected genomic regions which harbored genes for common biological categories including coat color, brain development, muscle development, growth, metabolism, olfaction and immunity. Genetic differentiation in regions harboring genes related to muscle development and growth was higher between breeds than between a given breed and the wild boar.

Conclusions/Significance

These results, suggest that although domesticated breeds have experienced similar selective pressures, selection has acted upon different genes. This might reflect the multiple domestication events of European breeds or could be the result of subsequent introgression of Asian alleles. Overall, it was estimated that approximately 7% of the porcine genome has been affected by selection events. This study illustrates that the massive parallel sequencing of genomic pools is a cost-effective approach to identify footprints of selection.     

Improved precision of QTL mapping using a nonlinear Bayesian method in a multi-breed population leads to greater accuracy of across-breed genomic predictions

Background

Genomic selection is increasingly widely practised, particularly in dairy cattle. However, the accuracy of current predictions using GBLUP (genomic best linear unbiased prediction) decays rapidly across generations, and also as selection candidates become less related to the reference population. This is likely caused by the effects of causative mutations being dispersed across many SNPs (single nucleotide polymorphisms) that span large genomic intervals. In this paper, we hypothesise that the use of a nonlinear method (BayesR), combined with a multi-breed (Holstein/Jersey) reference population will map causative mutations with more precision than GBLUP and this, in turn, will increase the accuracy of genomic predictions for selection candidates that are less related to the reference animals.

Results

BayesR improved the across-breed prediction accuracy for Australian Red dairy cattle for five milk yield and composition traits by an average of 7% over the GBLUP approach (Australian Red animals were not included in the reference population). Using the multi-breed reference population with BayesR improved accuracy of prediction in Australian Red cattle by 2 – 5% compared to using BayesR with a single breed reference population. Inclusion of 8478 Holstein and 3917 Jersey cows in the reference population improved accuracy of predictions for these breeds by 4 and 5%. However, predictions for Holstein and Jersey cattle were similar using within-breed and multi-breed reference populations. We propose that the improvement in across-breed prediction achieved by BayesR with the multi-breed reference population is due to more precise mapping of quantitative trait loci (QTL), which was demonstrated for several regions. New candidate genes with functional links to milk synthesis were identified using differential gene expression in the mammary gland.

Conclusions

QTL detection and genomic prediction are usually considered independently but persistence of genomic prediction accuracies across breeds requires accurate estimation of QTL effects. We show that accuracy of across-breed genomic predictions was higher with BayesR than with GBLUP and that BayesR mapped QTL more precisely. Further improvements of across-breed accuracy of genomic predictions and QTL mapping could be achieved by increasing the size of the reference population, including more breeds, and possibly by exploiting pleiotropic effects to improve mapping efficiency for QTL with small effects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0074-4) contains supplementary material, which is available to authorized users.     

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

Background

Signatures of selection are regions in the genome that have been preferentially increased in frequency and fixed in a population because of their functional importance in specific processes. These regions can be detected because of their lower genetic variability and specific regional linkage disequilibrium (LD) patterns.

Methods

By comparing the differences in regional LD variation between dairy and beef cattle types, and between indicine and taurine subspecies, we aim at finding signatures of selection for production and adaptation in cattle breeds. The VarLD method was applied to compare the LD variation in the autosomal genome between breeds, including Angus and Brown Swiss, representing taurine breeds, and Nelore and Gir, representing indicine breeds. Genomic regions containing the top 0.01 and 0.1 percentile of signals were characterized using the UMD3.1 Bos taurus genome assembly to identify genes in those regions and compared with previously reported selection signatures and regions with copy number variation.

Results

For all comparisons, the top 0.01 and 0.1 percentile included 26 and 165 signals and 17 and 125 genes, respectively, including TECRL, BT.23182 or FPPS, CAST, MYOM1, UVRAG and DNAJA1.

Conclusions

The VarLD method is a powerful tool to identify differences in linkage disequilibrium between cattle populations and putative signatures of selection with potential adaptive and productive importance.     

Genomic consequences of selection and genome-wide association mapping in soybean

Background

Crop improvement always involves selection of specific alleles at genes controlling traits of agronomic importance, likely resulting in detectable signatures of selection within the genome of modern soybean (Glycine max L. Merr.). The identification of these signatures of selection is meaningful from the perspective of evolutionary biology and for uncovering the genetic architecture of agronomic traits.

Results

To this end, two populations of soybean, consisting of 342 landraces and 1062 improved lines, were genotyped with the SoySNP50K Illumina BeadChip containing 52,041 single nucleotide polymorphisms (SNPs), and systematically phenotyped for 9 agronomic traits. A cross-population composite likelihood ratio (XP-CLR) method was used to screen the signals of selective sweeps. A total of 125 candidate selection regions were identified, many of which harbored genes potentially involved in crop improvement. To further investigate whether these candidate regions were in fact enriched for genes affected by selection, genome-wide association studies (GWAS) were conducted on 7 selection traits targeted in soybean breeding (grain yield, plant height, lodging, maturity date, seed coat color, seed protein and oil content) and 2 non-selection traits (pubescence and flower color). Major genomic regions associated with selection traits overlapped with candidate selection regions, whereas no overlap of this kind occurred for the non-selection traits, suggesting that the selection sweeps identified are associated with traits of agronomic importance. Multiple novel loci and refined map locations of known loci related to these traits were also identified.

Conclusions

These findings illustrate that comparative genomic analyses, especially when combined with GWAS, are a promising approach to dissect the genetic architecture of complex traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1872-y) contains supplementary material, which is available to authorized users.     

Imputation of sequence level genotypes in the Franches-Montagnes horse breed

Background

A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions.

Methods

We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute.

Results

The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%.

Conclusions

Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0063-7) contains supplementary material, which is available to authorized users.