Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Accuracy of marker-assisted selection with single markers and marker haplotypes in cattle

A key question for the implementation of marker-assisted selection (MAS) using markers in linkage disequilibrium with quantitative trait loci (QTLs) is how many markers surrounding each QTL should be used to ensure the marker or marker haplotypes are in sufficient linkage disequilibrium (LD) with the QTL. In this paper we compare the accuracy of MAS using either single markers or marker haplotypes in an Angus cattle data set consisting of 9323 genome-wide single nucleotide polymorphisms (SNPs) genotyped in 379 Angus cattle. The extent of LD in the data set was such that the average marker-marker r2 was 0.2 at 200 kb. The accuracy of MAS increased as the number of markers in the haplotype surrounding the QTL increased, although only when the number of markers in the haplotype was 4 or greater did the accuracy exceed that achieved when the SNP in the highest LD with the QTL was used. A large number of phenotypic records (>1000) were required to accurately estimate the effects of the haplotypes.     

Genome-wide scan for bovine twinning rate QTL using linkage disequilibrium

Twinning is a complex trait with negative impacts on health and reproduction, which cause economic loss in dairy production. Several twinning rate quantitative trait loci (QTL) have been detected in previous studies, but confidence intervals for QTL location are broad and many QTL are unreplicated. To identify genomic regions or genes associated with twinning rate, QTL analysis based on linkage combined with linkage disequilibrium (LLD) and individual marker associations was conducted across the genome using high-throughput single nucleotide polymorphism (SNP) genotypes. A total of 9919 SNP markers were genotyped with 200 sires and sons in 19 half-sib North American Holstein dairy cattle families. After SNPs were genotyped, informative markers were selected for genome-wide association tests and QTL searches. Evidence for twinning rate QTL was found throughout the genome. Thirteen markers significantly associated with twinning rate were detected on chromosomes 2, 5 and 14 ( P  < 2.3 × 10⁻⁵). Twenty-six regions on fourteen chromosomes were identified by LLD analysis at P  < 0.0007. Seven previously reported ovulation or twinning rate QTL were supported by results of single marker association or LLD analyses. Single marker association analysis and LLD mapping were complementary tools for the identification of putative QTL in this genome scan.     

Validated quantitative trait loci for eggshell quality in experimental and commercial laying hens

Compromised eggshell quality causes considerable economic losses for the egg industry. Breeding for improved eggshell quality has been very challenging. Eggshell quality is a trait that would greatly benefit from marker‐assisted selection, which would allow the selection of sires for their direct contribution to the trait and would also allow implementation of measurements integrating a number of shell parameters that are difficult to measure. In this study, we selected the most promising autosomal quantitative trait loci (QTL) affecting eggshell quality on chromosomes 2, 3, 6 and 14 from earlier experiments and we extended the F₂ population to include 1599 F₂ females. The study was repeated on two commercial populations: Lohmann Tierzucht Rhode Island Red line (n = 692 females) and a Hy‐Line White Plymouth Rock line (n = 290 progeny tested males). We analyzed the selected autosomal QTL regions on the three populations with SNP markers at 4–13 SNPs/Mb density. QTL for eggshell quality were replicated on all studied regions in the F₂ population. New QTL were detected for eggshell color on chromosomes 3 and 6. Marker associations with eggshell quality traits were validated in the tested commercial lines on chromosomes 2, 3 and 6, thus paving the way for marker‐assisted selection for improved eggshell quality.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

A genome‐wide association study to detect genetic variation for postpartum dysgalactia syndrome in five commercial pig breeding lines

Postpartum dysgalactia syndrome (PDS) in sows is an important disease after parturition with a relevant economic impact, affecting the health and welfare of both sows and piglets. The genetic background of this disease has been discussed and its heritability estimated, but further genetic analyses are lacking in detail. The aim of the current study was to detect loci affecting the susceptibility to PDS through a genome‐wide association approach. The study was designed as a family‐based association study with matched sampling of affected sows and healthy half‐ or full‐sib control sows on six farms. For the study, 597 sows (322 affected vs. 275 healthy control sows) were genotyped on 62 163 single nucleotide polymorphisms (SNPs) using the Illumina PorcineSNP60 BeadChip. After quality control, 585 sows (314 affected vs. 271 healthy control sows) and 49 740 SNPs remained for further analysis. Statistics were performed mainly with the r package genabel and included a principal component analysis. A statistically significant genome‐wide associated SNP was identified on porcine chromosome (SSC) 17. Further promising results with moderate significance were detected on SSC 13 and on an unplaced scaffold with an older annotation on SSC 15. The PRICKLE2 and NRP2 genes were identified as candidate genes near associated SNPs. Several quantitative trait loci (QTL) have been previously described in these genomic regions, including QTL for mammary gland condition, as teat number and non‐functional nipples QTL, as well as QTL for body temperature and gestation length.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

A Multi-Trait,Meta-analysis for Detecting Pleiotropic Polymorphisms for Stature,Fatness and Reproduction in Beef Cattle

Polymorphisms that affect complex traits or quantitative trait loci (QTL) often affect multiple traits. We describe two novel methods (1) for finding single nucleotide polymorphisms (SNPs) significantly associated with one or more traits using a multi-trait, meta-analysis, and (2) for distinguishing between a single pleiotropic QTL and multiple linked QTL. The meta-analysis uses the effect of each SNP on each of n traits, estimated in single trait genome wide association studies (GWAS). These effects are expressed as a vector of signed t-values (t) and the error covariance matrix of these t values is approximated by the correlation matrix of t-values among the traits calculated across the SNP (V). Consequently, t''V⁻¹t is approximately distributed as a chi-squared with n degrees of freedom. An attractive feature of the meta-analysis is that it uses estimated effects of SNPs from single trait GWAS, so it can be applied to published data where individual records are not available. We demonstrate that the multi-trait method can be used to increase the power (numbers of SNPs validated in an independent population) of GWAS in a beef cattle data set including 10,191 animals genotyped for 729,068 SNPs with 32 traits recorded, including growth and reproduction traits. We can distinguish between a single pleiotropic QTL and multiple linked QTL because multiple SNPs tagging the same QTL show the same pattern of effects across traits. We confirm this finding by demonstrating that when one SNP is included in the statistical model the other SNPs have a non-significant effect. In the beef cattle data set, cluster analysis yielded four groups of QTL with similar patterns of effects across traits within a group. A linear index was used to validate SNPs having effects on multiple traits and to identify additional SNPs belonging to these four groups.     

Fatty acid synthase effects on bovine adipose fat and milk fat

A quantitative trait locus (QTL) was identified by linkage analysis on bovine Chromosome 19 that affects the fatty acid, myristic acid (C14:0), in subcutaneous adipose tissue of pasture-fed beef cattle (99% level: experiment-wise significance). The QTL was also shown to have significant effects on ten fatty acids in the milk fat of pasture-fed dairy cattle. A positional candidate gene for this QTL was identified as fatty acid synthase (FASN), which is a multifunctional enzyme with a central role in the metabolism of lipids. Five single nucleotide polymorphisms (SNPs) were identified in the bovine FASN gene, and animals were genotyped for FASN SNPs in three different cattle resource populations. Linkage and association mapping results using these SNPs were consistent with FASN being the gene underlying the QTL. SNP substitution effects for C14:0 percentage were found to have an effect in the opposite direction in adipose fat to that in milk fat. It is concluded that SNPs in the bovine FASN gene are associated with variation in the fatty acid composition of adipose fat and milk fat.     

Effects of genetic sex and genomic background on epistasis in rainbow trout (Oncorhynchus mykiss)

Epistasis among quantitative trait loci (QTL) for survival (upper thermal tolerance, UTT) and morphological (fork length, FL and condition factor, K) traits was detected in purestrain and interstrain rainbow trout (Oncorhynchus mykiss) families. One sex-linked (OmyFGT19TUF) and three autosomal (Omy325UoG, Ssa14DU and Ssa20.19NUIG; linkage groups B, D and S, respectively) microsatellite loci linked to UTT QTL in this species were used. Within half sib families, significant effects of full sib family on epistasis involving Omy325UoG and OmyFGT19TUF were detected at a rate significantly higher than expected for UTT (p < 0.001*) and FL (p < 0.01*), using results significant at comparisonwise significance thresholds derived from permutational analysis. Measured across half sib families, the phenotype of female genotypic classes was more divergent from the family trait mean than that of males where epistasis involved the sex-linked locus OmyFGT19TUF (p = 0.0176*), and also for means over all families (p = 0.0355*). Female genotypic classes were also more divergent (p = 0.0011**) from the full sib trait mean where three-way interaction between OmyFGT19TUF, one of the autosomal loci and full sib family was significant, and marginally more divergent for trait means of genotypic classes across all full sib families (p = 0.0856). There was no evidence that these effects were more pronounced in hybrid F₁ families than purestrains.     

Association analysis of adiponectin and somatostatin polymorphisms on BTA1 with growth and carcass traits in Angus cattle

This study tested positional candidate genes adiponectin (ADIPOQ) and somatostatin (SST) for effects on carcass traits in a commercially relevant cattle population. Both genes are located within a region of BTA1 previously reported to harbour quantitative trait loci (QTL) that affect marbling, quality grade, yield grade, ribeye area and weaning weight in Bos taurus x Bos indicus crosses. Except for the first intron of ADIPOQ, both genes, including over 2 kb upstream of the promoters, were sequenced in five registered Angus sires to identify polymorphisms. A variable copy duplication and three single nucleotide polymorphisms (SNPs) in ADIPOQ and one SNP in SST were genotyped and tested for association with 19 traits in a 14-generation pedigree of 1697 registered Angus artificial insemination sires representing all the major USA lineages of the breed. Linear models that parameterized predicted genetic merits in terms of allele substitution effects were fit by weighted least squares, and goodness-of-fit tests were employed to differentiate causal mutations or polymorphisms in strong linkage disequilibrium (LD) with causal mutations from markers in weak LD with QTL. We confirmed the presence of QTL affecting marbling, ribeye muscle area and fat thickness in the vicinity of SST and ADIPOQ on BTA1 in Angus; excluded SST as underlying the ribeye muscle area QTL; and excluded ADIPOQ as underlying the marbling score QTL. However, association analysis provides very limited information about QTL location and has little intrinsic value when performed in the absence of linkage or LD analysis using flanking marker data to localize the QTL effect relative to positional candidate genes.     

Congenital Sensorineural Deafness in Dalmatian Dogs Associated with Quantitative Trait Loci

A genome-wide association study (GWAS) was performed for 235 Dalmatian dogs using the canine Illumina high density bead chip to identify quantitative trait loci (QTL) associated with canine congenital sensorineural deafness (CCSD). Data analysis was performed for all Dalmatian dogs and in addition, separately for brown-eyed and blue-eyed dogs because of the significant influence of eye colour on CCSD in Dalmatian dogs. Mixed linear model analysis (MLM) revealed seven QTL with experiment-wide significant associations (-log₁₀P>5.0) for CCSD in all Dalmatian dogs. Six QTL with experiment-wide significant associations for CCSD were found in brown-eyed Dalmatian dogs and in blue-eyed Dalmatian dogs, four experiment-wide significant QTL were detected. The experiment-wide CCSD-associated SNPs explained 82% of the phenotypic variance of CCSD. Five CCSD-loci on dog chromosomes (CFA) 6, 14, 27, 29 and 31 were in close vicinity of genes shown as causative for hearing loss in human and/or mouse.     

Detection of quantitative trait loci for carcass traits in the pig by using AFLP

For evaluation of the suitability of Amplified Fragment Length Polymorphism (AFLP) for detection of quantitative trait loci in farm animals, a combination of AFLP and selective genotyping has been applied as a rapid screening method for marker–QTL associations. Focusing on loci affecting eye muscle area, six extreme discordant sib pairs were selected from a Duroc × Berlin Miniature Pig F₂ experimental cross and examined by using 48 AFLP primer combinations. Two prominent AFLP markers were converted into simple codominant PCR markers (STS-Bo1 and STS-Bo3) and assigned to Sscr4 by physical and linkage mapping. Single marker analysis indicated association of the STS markers with a putative QTL influencing eye muscle area. Interval mapping confirmed the presence of a significant QTL for eye muscle area (P_genomewide < 0.01) on the Sscr4, with STS-Bo1 being the closer marker. At the same location, significant effects (P_genomewide < 0.01) on carcass length and backfat thickness were also detected. Our results demonstrate the capability of the combination of AFLP analysis and selective genotyping as a method for detection of genome regions containing QTL in livestock.     

Genome‐wide association analysis for quantitative trait loci influencing Warner–Bratzler shear force in five taurine cattle breeds

We performed a genome‐wide association study for Warner–Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single‐nucleotide polymorphisms (SNPs) within μ‐calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within‐ and across‐breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across‐breed analysis were moderately correlated (0.31–0.66) with those from the individual within‐breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across‐breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within‐breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across‐breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome‐wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine‐map the CAPN1 causal mutation to a 4581‐bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively.     

Mixed model association scans of multi-environmental trial data reveal major loci controlling yield and yield related traits in Hordeum vulgare in Mediterranean environments

An association panel consisting of 185 accessions representative of the barley germplasm cultivated in the Mediterranean basin was used to localise quantitative trait loci (QTL) controlling grain yield and yield related traits. The germplasm set was genotyped with 1,536 SNP markers and tested for associations with phenotypic data gathered over 2?years for a total of 24?year?×?location combinations under a broad range of environmental conditions. Analysis of multi-environmental trial (MET) data by fitting a mixed model with kinship estimates detected from two to seven QTL for the major components of yield including 1000 kernel weight, grains per spike and spikes per m², as well as heading date, harvest index and plant height. Several of the associations involved SNPs tightly linked to known major genes determining spike morphology in barley (vrs1 and int-c). Similarly, the largest QTL for heading date co-locates with SNPs linked with eam6, a major locus for heading date in barley for autumn sown conditions. Co-localization of several QTL related to yield components traits suggest that major developmental loci may be linked to most of the associations. This study highlights the potential of association genetics to identify genetic variants controlling complex traits.     

Identification of candidate genes and mutations in QTL regions for immune responses in chicken

There are two categories of immune responses – innate and adaptive immunity – both having polygenic backgrounds and a significant environmental component. In our study, adaptive immunity was represented by the specific antibody response toward keyhole limpet hemocyanin (KLH); innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Defining genetic bases of immune responses leads from defining quantitative trait loci (QTL) toward a single mutation responsible for variation in the phenotypic trait. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chicken by performing an association study of SNPs located in candidate genes defined in QTL regions. Candidate genes and SNPs in QTL regions were selected in silico. SNP association was based on a custom SNP panel, GoldenGate genotyping assay (Illumina) and two statistical models: random mixed model and CAR score. The most significant SNP for immune response toward KLH was located in the JMJD6 gene located on GGA18. Four SNPs in candidate genes FOXJ1 (GGA18), EPHB1 (GGA9), PTGER4 (GGAZ) and PRKCB (GGA14) showed association with natural antibodies for LPS. A single SNP in ITGB4 (GGA18) was associated with natural antibodies for LTA. All associated SNPs mentioned above showed additive effects.     

Quantitative trait locus analysis of a recombinant inbred line population derived from a Lycopersicon esculentum x Lycopersicon cheesmanii cross

Quantitative trait loci influencing fruit traits were identified by restriction fragment length polymorphism (RFLP) analysis in a population of recombinant inbred lines (RIL) derived from a cross of the cultivated tomato, Lycopersicon esculentum with a related wild species Lycopersicon cheesmanii. One hundred thirty-two polymorphic RFLP loci spaced throughout the tomato genome were scored for 97 F₈ RIL families. Fruit weight and soluble solids were measured in replicated trials during 1991 and 1992. Seed weight was measured in 1992. Significant (P<0.01 level) quantitative trait locus (QTL) associations of marker loci were identified for each trait. A total of 73 significant marker locus-trait associations were detected for the three traits measured. Fifty-three of these associations were for fruit weight and soluble solids, many of which involved marker loci signficantly associated with both traits. QTL with large effects on all three traits were detected on chromosome 6. Greater homozygosity at many loci in the RIL population as compared to F₂ populations and greater genomic coverage resulted in increased precision in the estimation of QTL effects, and large proportions of the total phenotypic variance were explained by marker class variation at significant marker loci for many traits. The RIL population was effective in detecting and discriminating among QTL for these traits previously identified in other investigations despite skewed segregation ratios at many marker loci. Large additive effects were measured at significant marker loci. Lower fruit weight, higher soluble solids, and lower seed weight were generally associated with RFLP alleles from theL. cheesmanii parent.     

Identification of polymorphisms influencing feed intake and efficiency in beef cattle

Feed efficiency is an economically important trait in beef cattle. Net feed efficiency, measured as residual feed intake (RFI), is the difference between actual feed intake and the predicted feed intake required for maintenance and gain of the animal. SNPs that show associations with RFI may be useful quantitative trait nucleotides for marker-assisted selection. This study identified associations between SNPs underlying five RFI QTL on five bovine chromosomes (BTA2, 5, 10, 20 and 29) with measures of dry matter intake (DMI), RFI and feed conversion ratio (FCR) in beef cattle. Six SNPs were found to have effects on RFI (P < 0.05). The largest single SNP allele substitution effect for RFI was -0.25 kg/day located on BTA2. The combined effects of the SNPs found significant in this experiment explained 6.9% of the phenotypic variation of RFI. Not all the RFI SNPs showed associations with DMI and FCR even though these traits are highly correlated with RFI (r = 0.77 and r = 0.62 respectively). This shows that these SNPs may be affecting the underlying biological mechanisms of feed efficiency beyond feed intake control and weight gain efficiency. These SNPs can be used in marker-assisted selection but first it will be important to verify these effects in independent populations of cattle.     

Mapping of a Quantitative Trait Locus for Beef Marbling on Bovine Chromosome 9 in Purebred Japanese Black Cattle

The goal of this study is to detect quantitative trait loci (QTL) for carcass traits applicable for a DNA-based breeding system in a Japanese Black cattle population. A purebred paternal half-sib family from a commercial line composed of 65 steers was initially analyzed using 188 informative microsatellites giving a 16-cM average interval covering 29 autosomes. A significant QTL for marbling was detected in the centromeric portion of bovine chromosome (BTA) 9. After additional marker genotyping across a larger sample size composed of 169 individuals, this locus was refined to a 20-cM confidence interval between microsatellites BM1227 (24 cM) and DIK2741 (50 cM) at a 1% chromosome-wise threshold. The allele substitution effect between Q and q for a beef marbling standard score (1 to 12 range) on BTA9 was 1.0 (5.7% of total phenotypic variance in QTL contribution in this family). This result provides a primary platform for a marker-assisted selection system of the beef marbling trait within the Japanese Black (Wagyu) cattle population.     

Integrating Milk Metabolite Profile Information for the Prediction of Traditional Milk Traits Based on SNP Information for Holstein Cows

In this study the benefit of metabolome level analysis for the prediction of genetic value of three traditional milk traits was investigated. Our proposed approach consists of three steps: First, milk metabolite profiles are used to predict three traditional milk traits of 1,305 Holstein cows. Two regression methods, both enabling variable selection, are applied to identify important milk metabolites in this step. Second, the prediction of these important milk metabolite from single nucleotide polymorphisms (SNPs) enables the detection of SNPs with significant genetic effects. Finally, these SNPs are used to predict milk traits. The observed precision of predicted genetic values was compared to the results observed for the classical genotype-phenotype prediction using all SNPs or a reduced SNP subset (reduced classical approach). To enable a comparison between SNP subsets, a special invariable evaluation design was implemented. SNPs close to or within known quantitative trait loci (QTL) were determined. This enabled us to determine if detected important SNP subsets were enriched in these regions. The results show that our approach can lead to genetic value prediction, but requires less than 1% of the total amount of (40,317) SNPs., significantly more important SNPs in known QTL regions were detected using our approach compared to the reduced classical approach. Concluding, our approach allows a deeper insight into the associations between the different levels of the genotype-phenotype map (genotype-metabolome, metabolome-phenotype, genotype-phenotype).