Genome‐wide association study uncovers four QTL predisposing to supernumerary teats in cattle

Supernumerary teats (hyperthelia, SNTs) are a common abnormality of the bovine udder with a medium to high heritability and a postulated oligogenic or polygenic inheritance pattern. SNTs not only negatively affect machine milking ability but also act as a reservoir for bacteria. A genome‐wide association study was carried out to identify genes involved in the development of SNTs in the dual‐purpose Fleckvieh breed. A total of 2467 progeny‐tested bulls were genotyped at 43 698 single nucleotide polymorphisms, and daughter yield deviations (DYDs) for ‘udder clearness' (UC) were used as high‐heritability phenotypes. Massive structuring of the study population was accounted for by principal components analysis‐based and mixed model‐based approaches. Four loci on BTA5, BTA6, BTA11 and BTA17 were significantly associated with the UC DYD. Three associated regions contain genes of the highly conserved Wnt signalling pathway. The four QTL together account for 10.7% of the variance of the UC DYD, whereas the major fraction of the DYD variance is attributable to chromosomes with no identified QTL. Our results support both an oligogenic and a polygenic inheritance pattern of SNTs in cattle. The identified candidate genes permit insights into the genetic architecture of teat malformations in cattle and provide clues to unravel the molecular mechanisms of mammary gland alterations in cattle and other species.     

A genome-screen experiment to detect quantitative trait loci affecting resistance to facial eczema disease in sheep

Facial eczema (FE) is a secondary photosensitization disease arising from liver cirrhosis caused by the mycotoxin sporidesmin. The disease affects sheep, cattle, deer and goats, and costs the New Zealand sheep industry alone an estimated NZ$63M annually. A long-term sustainable solution to this century-old FE problem is to breed for disease-resistant animals by marker-assisted selection. As a step towards finding a diagnostic DNA test for FE sensitivity, we have conducted a genome-scan experiment to screen for quantitative trait loci (QTL) affecting this trait in Romney sheep. Four F₁ sires, obtained from reciprocal matings of FE resistant and susceptible selection-line animals, were used to generate four outcross families. The resulting half-sib progeny were artificially challenged with sporidesmin to phenotype their FE traits measured in terms of their serum levels of liver-specific enzymes, namely gamma-glutamyl transferase and glutamate dehydrogenase. In a primary screen using selective genotyping on extreme progeny of each family, a total of 244 DNA markers uniformly distributed over all 26 ovine autosomes (with an autosomal genome coverage of 79–91%) were tested for linkage to the FE traits. Data were analysed using Haley–Knott regression. The primary screen detected one significant and one suggestive QTL on chromosomes 3 and 8 respectively. Both the significant and suggestive QTL were followed up in a secondary screen where all progeny were genotyped and analysed; the QTL on chromosome 3 was significant in this analysis.     

Mapping quantitative trait Loci underlying fitness-related traits in a free-living sheep population

We searched for quantitative trait loci (QTL) underlying fitness-related traits in a free-living pedigree of 588 Soay sheep in which a genetic map using 251 markers with an average spacing of 15 cM had been established previously. Traits examined included birth date and weight, considered both as maternal and offspring traits, foreleg length, hindleg length, and body weight measured on animals in August and jaw length and metacarpal length measured on cleaned skeletal material. In some cases the data were split to consider different age classes separately, yielding a total of 15 traits studied. Genetic and environmental components of phenotypic variance were estimated for each trait and, for those traits showing nonzero heritability (N= 12), a QTL search was conducted by comparing a polygenic model with a model including a putative QTL. Support for a QTL at genome-wide significance was found on chromosome 11 for jaw length; suggestive QTL were found on chromosomes 2 and 5 (for birth date as a trait of the lamb), 8 (birth weight as a trait of the lamb), and 15 (adult hindleg length). We discuss the prospects for refining estimates of QTL position and effect size in the study population, and for QTL searches in free-living pedigrees in general.     

Identification of an ovulation rate QTL in cattle on BTA14 using selective DNA pooling and interval mapping

Increased twinning incidence in beef cattle has the potential to improve production efficiency. However, phenotypic selection for twinning rate is difficult because of the trait's low heritability and the long time interval necessary to collect phenotypic records. Therefore, this trait and the correlated trait of ovulation rate are ideal candidates for marker-assisted selection. The objective of this study was to complete a genome-wide search for ovulation rate quantitative trait loci (QTL) in two related sire families. The families (paternal halfsib sires 839802 and 839803) were from a population of cattle selected for ovulation rate at the USDA Meat Animal Research Center, Clay Center, Nebraska. Putative ovulation rate QTL have previously been identified in the 839802 family on chromosomes 7 and 19; however, marker coverage in the original scan was not complete. This study fills the gaps in marker coverage of the earlier study by adding approximately 60 informative microsatellites to each sire family. Each family was genotyped using selective DNA pooling. Sons and daughters were included in either the high or low pool based on their estimated breeding value deviations from the mid-parent average (EBVMD) for ovulation rate. Approximately 40% (839802) and 26% (839803) of available progeny comprised the high and low pools combined. Pooled typing revealed possible associations (nominal P < 0.05) between ovulation rate and marker genotype for 11 and 15 microsatellites in the 839802 and 839803 families, respectively. Subsequent interval mapping strengthened support for the presence of an ovulation rate QTL on BTA14 (chromosome-wise P < 0.02).     

Mendelian Factors Underlying Quantitative Traits in Tomato: Comparison across Species, Generations, and Environments

As part of ongoing studies regarding the genetic basis of quantitative variation in phenotype, we have determined the chromosomal locations of quantitative trait loci (QTLs) affecting fruit size, soluble solids concentration, and pH, in a cross between the domestic tomato (Lycopersicon esculentum Mill.) and a closely-related wild species, L. cheesmanii. Using a RFLP map of the tomato genome, we compared the inheritance patterns of polymorphisms in 350 F2 individuals with phenotypes scored in three different ways: (1) from the F2 progeny themselves, grown near Davis, California; (2) from F3 families obtained by selfing each F2 individual, grown near Gilroy, California (F3-CA); and (3) from equivalent F3 families grown near Rehovot, Israel (F3-IS). Maximum likelihood methods were used to estimate the approximate chromosomal locations, phenotypic effects (both additive effects and dominance deviations), and gene action of QTLs underlying phenotypic variation in each of these three environments. A total of 29 putative QTLs were detected in the three environments. These QTLs were distributed over 11 of the 12 chromosomes, accounted for 4.7-42.0% of the phenotypic variance in a trait, and showed different types of gene action. Among these 29 QTLs, 4 were detected in all three environments, 10 in two environments, and 15 in only a single environment. The two California environments were most similar, sharing 11/25 (44%) QTLs, while the Israel environment was quite different, sharing 7/20 (35%) and 5/26 (19%) QTLs with the respective California environments. One major goal of QTL mapping is to predict, with maximum accuracy, which individuals will produce progeny showing particular phenotypes. Traditionally, the phenotype of an individual alone has been used to predict the phenotype of its progeny. Our results suggested that, for a trait with low heritability (soluble solids), the phenotype of F3 progeny could be predicted more accurately from the genotype of the F2 parent at QTLs than from the phenotype of the F2 individual. For a trait with intermediate heritability (fruit pH), QTL genotype and observed phenotype were about equally effective at predicting progeny phenotype. For a trait with high heritability (mass per fruit), knowing the QTL genotype of an individual added little if any predictive value, to simply knowing the phenotype. The QTLs mapped in the L. esculentum X L. cheesmanii F2 appear to be at similar locations to many of those mapped in a previous cross with a different wild tomato (L. chmielewskii).(ABSTRACT TRUNCATED AT 400 WORDS)     

Genetic aspects of twinning in cattle

Twinning in cattle ranges from about 1% for beef breeds to about 4% for dairy breeds. The incidence of double births may have both positive and negative effects, which mainly depends on the purpose for which the cattle are raised. Because of freemartinism, as well as management problems connected e.g. with a greater risk of dystocia and retained placenta, it is an undesirable trait in dairy herds. In beef cattle, however, twinning can considerably increase the efficiency of production. Low heritability, a long generation interval for progeny testing, sex-limited expression and an unfavourable correlation with milk yield make twinning difficult to control by selection. Hence, it is the type of trait for which the identification of the genetic marker - quantitative trait loci (QTL) linkage and the implementation of marker-assisted selection in breeding strategies are expected to be especially beneficial. Searching for QTL influencing the reproductive rate in cattle was performed mainly in the US Meat Animal Research Center twinning herd and in the commercial Norwegian cattle population. Among several genome regions that appear to control twinning and ovulation rates, the most interesting seem to be chromosomes 5, 7, 19 and 23.     

A QTL study of cattle behavioral traits in embryo transfer families.

Two behavioral traits, temperament and habituation, were measured in 130 calves from 17 full-sib families which comprise the Canadian Beef Cattle Reference Herd. Using variance components, heritability was calculated as 0.36 for temperament and 0.46 for habituation. Genotyping of 162 microsatellites at approximately 20 cM intervals allowed the detection of six quantitative trait loci (QTL) for behavior traits on cattle chromosomes 1, 5, 9, 11, 14, 15.     

Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

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.     

Epistasis affecting litter size in mice

Litter size is an important reproductive trait as it makes a major contribution to fitness. Generally, traits closely related to fitness show low heritability perhaps because of the corrosive effects of directional natural selection on the additive genetic variance. Nonetheless, low heritability does not imply, necessarily, a complete absence of genetic variation because genetic interactions (epistasis and dominance) contribute to variation in traits displaying strong heterosis in crosses, such as litter size. In our study, we investigated the genetic architecture of litter size in 166 females from an F2 intercross of the SM/J and LG/J inbred mouse strains. Litter size had a low heritability (h2 = 12%) and a low repeatability (r = 33%). Using interval-mapping methods, we located two quantitative trait loci (QTL) affecting litter size at locations D7Mit21 + 0 cM and D12Mit6 + 8 cM, on chromosomes 7 and 12 respectively. These QTL accounted for 12.6% of the variance in litter size. In a two-way genome-wide epistasis scan we found eight QTL interacting epistatically involving chromosomes 2, 4, 5, 11, 14, 15 and 18. Taken together, the QTL and their interactions explain nearly 49% (39.5% adjusted multiple r2) of the phenotypic variation for litter size in this cross, an increase of 36% over the direct effects of the QTL. This indicates the importance of epistasis as a component of the genetic architecture of litter size and fitness in our intercross population.     

Identification and Validation of Quantitative Trait Loci (QTL) for Canine Hip Dysplasia (CHD) in German Shepherd Dogs

Canine hip dysplasia (CHD) is the most common hereditary skeletal disorder in dogs. To identify common alleles associated with CHD, we genotyped 96 German Shepherd Dogs affected by mild, moderate and severe CHD and 96 breed, sex, age and birth year matched controls using the Affymetrix canine high density SNP chip. A mixed linear model analysis identified five SNPs associated with CHD scores on dog chromosomes (CFA) 19, 24, 26 and 34. These five SNPs were validated in a by sex, age, birth year and coancestry stratified sample of 843 German Shepherd Dogs including 277 unaffected dogs and 566 CHD-affected dogs. Mean coancestry coefficients among and within cases and controls were <0.1%. Genotype effects of these SNPs explained 20–32% of the phenotypic variance of CHD in German Shepherd Dogs employed for validation. Genome-wide significance in the validation data set could be shown for each one CHD-associated SNP on CFA24, 26 and 34. These SNPs are located within or in close proximity of genes involved in bone formation and related through a joint network. The present study validated positional candidate genes within two previously known quantitative trait loci (QTL) and a novel QTL for CHD in German Shepherd Dogs.     

Multiple loci associated with canine hip dysplasia (CHD) in German shepherd dogs

Canine hip dysplasia (CHD) is the most common hereditary skeletal disorder in dogs. To identify common alleles associated with CHD, we developed 37 informative single nucleotide polymorphisms (SNPs) within 13 quantitative trait loci (QTL) previously identified for German shepherd dogs. These SNPs were genotyped in 95 German shepherd dogs affected by CHD and 95 breed, sex, and birth year-matched controls. A total of ten SNPs significant at a nominal P value of 0.05 were validated in 843 German shepherd dogs including 277 unaffected dogs and 566 CHD-affected dogs. Cases and controls were sampled from the whole German shepherd dog population in Germany in such a way that mean coancestry coefficients were below 0.1 % within cases and controls as well as among cases and controls. We identified nine SNPs significantly associated with CHD within five QTL on dog chromosomes (CFA) 3, 9, 26, 33, and 34. Genotype effects of these nine SNPs explained between 22 and 34 % of the phenotypic variance of hip dysplasia in German shepherd dogs. The strongest associated SNPs were located on CFA33 and 34 within the candidate genes PNCP, TRIO, and SLC6A3. Thus, the present study validated positional candidate genes within five QTL for CHD.     

Mapping of quantitative trait loci controlling tick [Riphicephalus (Boophilus) microplus] resistance on bovine chromosomes 5, 7 and 14

Differences in domestication and selection processes have contributed to considerable phenotypic and genotypic differences between Bos taurus and Bos indicus cattle breeds. Of particular interest in tropical and subtropical production environments are those genetic differences between subspecies that underlie the phenotypic extremes in tolerance and susceptibility to parasite infection. In general, B. taurus cattle are more susceptible to ectoparasites than B. indicus cattle in tropical environments, and much of this difference is under genetic control. To identify genomic regions involved in tick resistance, we developed a B. taurus x B. indicus F(2) experimental population to map quantitative trait loci (QTL) for resistance to the Riphicephalus (Boophilus) microplus tick. About 300 individuals were measured for parasite load in two seasons (rainy and dry) and genotyped for 23 microsatellite markers covering chromosomes 5, 7 and 14. We mapped a suggestive chromosome-wide QTL for tick load in the rainy season (P < 0.05) on chromosome 5. For the dry season, suggestive (P < 0.10) chromosome-wide QTL were mapped on chromosomes 7 and 14. The additive effect of the QTL on chromosome 14 corresponds to 3.18% of the total observed phenotypic variance. Our QTL-mapping study has identified different genomic regions controlling tick resistance; these QTL were dependent upon the season in which the ticks were counted, suggesting that the QTL in question may depend on environmental factors.     

Mapping Quantitative Trait Loci Controlling Milk Production in Dairy Cattle by Exploiting Progeny Testing

We have exploited ``progeny testing' to map quantitative trait loci (QTL) underlying the genetic variation of milk production in a selected dairy cattle population. A total of 1,518 sires, with progeny tests based on the milking performances of > 150,000 daughters jointly, was genotyped for 159 autosomal microsatellites bracketing 1645 centimorgan or approximately two thirds of the bovine genome. Using a maximum likelihood multilocus linkage analysis accounting for variance heterogeneity of the phenotypes, we identified five chromosomes giving very strong evidence (LOD score >/= 3) for the presence of a QTL controlling milk production: chromosomes 1, 6, 9, 10 and 20. These findings demonstrate that loci with considerable effects on milk production are still segregating in highly selected populations and pave the way toward marker-assisted selection in dairy cattle breeding.     

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.     

A genome scan for quantitative trait loci influencing carcass,post‐natal growth and reproductive traits in commercial Angus cattle

To gain insight into the number of loci of large effect that underlie variation in cattle, a quantitative trait locus (QTL) scan for 14 economically important traits was performed in two commercial Angus populations using 390 microsatellites, 11 single nucleotide polymorphisms (SNPs) and one duplication loci. The first population comprised 1769 registered Angus bulls born between 1955 and 2003, with Expected Progeny Differences computed by the American Angus Association. The second comprised 38 half‐sib families containing 1622 steers with six post‐natal growth and carcass phenotypes. Linkage analysis was performed by half‐sib least squares regression with gridqtl or Bayesian Markov chain Monte Carlo analysis of complex pedigrees with loki . Of the 673 detected QTL, only 118 have previously been reported, reflecting both the conservative approach to QTL reporting in the literature, and the more liberal approach taken in this study. From 33 to 71% of the genetic variance and 35 to 56% of the phenotypic variance in each trait was explained by the detected QTL. To analyse the effects of 11 SNPs and one duplication locus within candidate genes on each trait, a single marker analysis was performed by fitting an additive allele substitution model in both mapping populations. There were 53 associations detected between the SNP/duplication loci and traits with ?log₁₀P_nominal≥ 4.0, where each association explained 0.92% to 4.4% of the genetic variance and 0.01% to 1.86% of the phenotypic variance. Of these associations, only six SNP/duplication loci were located within 8 cM of a QTL peak for the trait, with two being located at the QTL peak: SST_DG156121:c.362A>G for ribeye muscle area and TG_X05380:c.422C>T for calving ease. Strong associations between several SNP/duplication loci and trait variation were obtained in the absence of any detected linked QTL. However, we reject the causality of several commercialized DNA tests, including an association between TG_X05380:c.422C>T and marbling in Angus cattle.     

Mapping of fertility traits in Finnish Ayrshire by genome-wide association analysis

A whole-genome scan using single marker association was used to detect chromosome regions associated with seven female fertility traits in Finnish Ayrshire dairy cattle. The phenotypic data consisted of de-regressed estimated breeding values for 340 bulls which were estimated using a single trait model. Genotypes were obtained with the Illumina BovineSNP50 panel and a total of 35 630 informative, high-quality single nucleotide polymorphism (SNP) markers were used. The association analysis was performed using a mixed-model approach which fitted a fixed effect for each SNP and a random polygenic effect. We detected eleven genome-wide significant associations on eight different chromosomes. With at least chromosome-wise significance after Bonferroni correction, sixteen SNPs on nine chromosomes showed significant associations with one or more fertility traits. The results confirmed quantitative trait loci on three chromosomes (1, 2 and 20) for fertility traits previously reported for the same breed and one on chromosome four previously detected in Holstein cattle.     

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.     

Genome-wide association study identifies two major loci affecting calving ease and growth-related traits in cattle

Identifying quantitative trait loci (QTL) underlying complex, low-heritability traits is notoriously difficult. Prototypical for such traits, calving ease is an important breeding objective of cattle (Bos taurus)-improving programs. To identify QTL underlying calving ease, we performed a genome-wide association study using estimated breeding values (EBVs) as highly heritable phenotypes for paternal calving ease (pCE) and related traits. The massively structured study population consisted of 1800 bulls of the German Fleckvieh (FV) breed. Two pCE-associated regions on bovine chromosomes (BTA) 14 and 21 (P = 5.72 × 10(-15) and P = 2.27 × 10(-8), respectively) were identified using principal components analysis to correct for population stratification. The two most significantly associated SNPs explain 10% of the EBV variation. Since marker alleles with negative effect on pCE have positive effects on growth-related traits, the QTL may exert their effects on the birthing process through fetal growth traits. The QTL region on BTA14 corresponds to a human chromosome (HSA) region that is associated with growth characteristics. The HSA region corresponding to the BTA21 pCE QTL is maternally imprinted and involved in the Prader-Willi and Angelman syndromes. Resequencing of positional candidate genes on BTA14 revealed a highly significantly (P = 1.96 × 10(-14)) associated polymorphism ablating a polyadenylation signal of the gene encoding ribosomal protein S20 (RPS20). Our study demonstrates the leverage potential of EBVs in unraveling the genetic architecture of lowly heritable traits.