Genetic analysis of Verticillium wilt resistance in a backcross inbred line population and a meta-analysis of quantitative trait loci for disease resistance in cotton

Background

Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding.

Results

This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated.

Conclusions

Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Electronic supplementary material

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

Limits on the reproducibility of marker associations with southern leaf blight resistance in the maize nested association mapping population

Background

A previous study reported a comprehensive quantitative trait locus (QTL) and genome wide association study (GWAS) of southern leaf blight (SLB) resistance in the maize Nested Association Mapping (NAM) panel. Since that time, the genomic resources available for such analyses have improved substantially. An updated NAM genetic linkage map has a nearly six-fold greater marker density than the previous map and the combined SNPs and read-depth variants (RDVs) from maize HapMaps 1 and 2 provided 28.5 M genomic variants for association analysis, 17 fold more than HapMap 1. In addition, phenotypic values of the NAM RILs were re-estimated to account for environment-specific flowering time covariates and a small proportion of lines were dropped due to genotypic data quality problems. Comparisons of original and updated QTL and GWAS results confound the effects of linkage map density, GWAS marker density, population sample size, and phenotype estimates. Therefore, we evaluated the effects of changing each of these parameters individually and in combination to determine their relative impact on marker-trait associations in original and updated analyses.

Results

Of the four parameters varied, map density caused the largest changes in QTL and GWAS results. The updated QTL model had better cross-validation prediction accuracy than the previous model. Whereas joint linkage QTL positions were relatively stable to input changes, the residual values derived from those QTL models (used as inputs to GWAS) were more sensitive, resulting in substantial differences between GWAS results. The updated NAM GWAS identified several candidate genes consistent with previous QTL fine-mapping results.

Conclusions

The highly polygenic nature of resistance to SLB complicates the identification of causal genes. Joint linkage QTL are relatively stable to perturbations of data inputs, but their resolution is generally on the order of tens or more Mbp. GWAS associations have higher resolution, but lower power due to stringent thresholds designed to minimize false positive associations, resulting in variability of detection across studies. The updated higher density linkage map improves QTL estimation and, along with a much denser SNP HapMap, greatly increases the likelihood of detecting SNPs in linkage with causal variants. We recommend use of the updated genetic resources and results but emphasize the limited repeatability of small-effect associations.

Electronic supplementary material

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

A Genome-Wide Association Study Reveals That Genes with Functions for Bone Development Are Associated with Body Conformation in Catfish

Body conformation is of great scientific and commercial interest for aquaculture fish species because it affects biological adaptation of the organism to environments, and is of economic importance to the aquaculture industry considering its direct effect on fillet yield. Catfish is the primary aquaculture species in the USA. Two major species used in the aquaculture industry, channel catfish and blue catfish, differ in body shape and therefore the backcross progenies serve as a good model for quantitative trait locus (QTL) analysis. Here, a genome-wide association study (GWAS) with hybrid catfish was conducted to identify the QTL for body conformation, including deheaded body length (DBL), body length (BL), body depth (BD), and body breadth (BB), which were all standardized by cubic root of body weight. Overall, the results indicate that the traits are polygenic. For DBL, linkage group (LG) 2 and LG 24 contain significant QTL, and LG 13 and LG 26 contain suggestively associated QTL (?log₁₀(P value) > 4.5). Compared with DBL, additional SNPs were identified to be associated with body length on LG 2, LG 7, and LG 18. Although no significant QTL for body depth was found, three suggestively associated QTLs were identified on LG 5, LG 13, and LG 14. No SNP for body breadth reached the threshold for suggestive association. Genes close to the associated SNPs were determined, many of which are known to be involved in bone development. This work therefore provides the basis for future identification of causal genes for the control of body conformation.     

Quantitative trait analysis of the development of pulmonary tolerance to inhaled zinc oxide in mice

Background

Individuals may develop tolerance to the induction of adverse pulmonary effects following repeated exposures to inhaled toxicants. Previously, we demonstrated that genetic background plays an important role in the development of pulmonary tolerance to inhaled zinc oxide (ZnO) in inbred mouse strains, as assessed by polymorphonuclear leukocytes (PMNs), macrophages, and total protein in bronchoalveolar lavage (BAL) phenotypes. The BALB/cByJ (CBy) and DBA/2J (D2) strains were identified as tolerant and non-tolerant, respectively. The present study was designed to identify candidate genes that control the development of pulmonary tolerance to inhaled ZnO.

Methods

Genome-wide linkage analyses were performed on a CByD2F2 mouse cohort phenotyped for BAL protein, PMNs, and macrophages following 5 consecutive days of exposure to 1.0 mg/m³ inhaled ZnO for 3 hours/day. A haplotype analysis was carried out to determine the contribution of each quantitative trait locus (QTL) and QTL combination to the overall BAL protein phenotype. Candidate genes were identified within each QTL interval using the positional candidate gene approach.

Results

A significant quantitative trait locus (QTL) on chromosome 1, as well as suggestive QTLs on chromosomes 4 and 5, for the BAL protein phenotype, was established. Suggestive QTLs for the BAL PMN and macrophage phenotypes were also identified on chromosomes 1 and 5, respectively. Analysis of specific haplotypes supports the combined effect of three QTLs in the overall protein phenotype. Toll-like receptor 5 (Tlr5) was identified as an interesting candidate gene within the significant QTL for BAL protein on chromosome 1. Wild-derived Tlr5-mutant MOLF/Ei mice were tolerant to BAL protein following repeated ZnO exposure.

Conclusion

Genetic background is an important influence in the acquisition of pulmonary tolerance to BAL protein, PMNs, and macrophages following ZnO exposure. Promising candidate genes exist within the identified QTL intervals that would be good targets for additional studies, including Tlr5. The implications of tolerance to health risks in humans are numerous, and this study furthers the understanding of gene-environment interactions that are likely to be important factors from person-to-person in regulating the development of pulmonary tolerance to inhaled toxicants.     

Large-effect pleiotropic or closely linked QTL segregate within and across ten US cattle breeds

Background

The availability of high-density SNP assays including the BovineSNP50 (50 K) enables the identification of novel quantitative trait loci (QTL) and improvement of the resolution of the locations of previously mapped QTL. We performed a series of genome-wide association studies (GWAS) using 50 K genotypes scored in 18,274 animals from 10 US beef cattle breeds with observations for twelve body weights, calving ease and carcass traits.

Results

A total of 159 large-effects QTL (defined as 1-Mb genome windows explaining more than 1% of additive genetic variance) were identified. In general, more QTL were identified in analyses with bigger sample sizes. Four large-effect pleiotropic or closely linked QTLs located on BTA6 at 37–42 Mb (primarily at 38 Mb), on BTA7 at 93 Mb, on BTA14 at 23–26 Mb (primarily at 25 Mb) and on BTA20 at 4 Mb were identified in more than one breed. Several breed-specific large-effect pleiotropic or closely linked QTL were also identified. Some identified QTL regions harbor genes known to have large effects on a variety of traits in cattle such as PLAG1 and MSTN and others harbor promising candidate genes including NCAPG, ARRDC3, ERGIC1, SH3PXD2B, HMGA2, MSRB3, LEMD3, TIGAR, SEPT7, and KIRREL3. Gene ontology analysis revealed that genes involved in ossification and in adipose tissue development were over-represented in the identified pleiotropic QTL. Also, the MAPK signaling pathway was identified as a common pathway affected by the genes located near the pleiotropic QTL.

Conclusions

This largest GWAS ever performed in beef cattle, led us to discover several novel across-breed and breed-specific large-effect pleiotropic QTL that cumulatively account for a significant percentage of additive genetic variance (e.g. more than a third of additive genetic variance of birth and mature weights; and calving ease direct in Hereford). These results will improve our understanding of the biology of growth and body composition in cattle.

Electronic supplementary material

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

The genetics of domestication of rice bean, Vigna umbellata

Background and Aims

The Asian genus Vigna, to which four cultivated species (rice bean, azuki bean, mung bean and black gram) belong, is suitable for comparative genomics. The aims were to construct a genetic linkage map of rice bean, to identify the genomic regions associated with domestication in rice bean, and to compare these regions with those in azuki bean.

Methods

A genetic linkage map was constructed by using simple sequence repeat and amplified fragment length polymorphism markers in the BC₁F₁ population derived from a cross between cultivated and wild rice bean. Using this map, 31 domestication-related traits were dissected into quantitative trait loci (QTLs). The genetic linkage map and QTLs of rice bean were compared with those of azuki bean.

Key Results

A total of 326 markers converged into 11 linkage groups (LGs), corresponding to the haploid number of rice bean chromosomes. The domestication-related traits in rice bean associated with a few major QTLs distributed as clusters on LGs 2, 4 and 7. A high level of co-linearity in marker order between the rice bean and azuki bean linkage maps was observed. Major QTLs in rice bean were found on LG4, whereas major QTLs in azuki bean were found on LG9.

Conclusions

This is the first report of a genetic linkage map and QTLs for domestication-related traits in rice bean. The inheritance of domestication-related traits was so simple that a few major QTLs explained the phenotypic variation between cultivated and wild rice bean. The high level of genomic synteny between rice bean and azuki bean facilitates QTL comparison between species. These results provide a genetic foundation for improvement of rice bean; interchange of major QTLs between rice bean and azuki bean might be useful for broadening the genetic variation of both species.     

Jasmonate-dependent plant defense restricts thrips performance and preference

Background

Worldwide, diseases are important reducers of peanut (Arachis hypogaea) yield. Sources of resistance against many diseases are available in cultivated peanut genotypes, although often not in farmer preferred varieties. Wild species generally harbor greater levels of resistance and even apparent immunity, although the linkage of agronomically un-adapted wild alleles with wild disease resistance genes is inevitable. Marker-assisted selection has the potential to facilitate the combination of both cultivated and wild resistance loci with agronomically adapted alleles. However, in peanut there is an almost complete lack of knowledge of the regions of the Arachis genome that control disease resistance.

Results

In this work we identified candidate genome regions that control disease resistance. For this we placed candidate disease resistance genes and QTLs against late leaf spot disease on the genetic map of the A-genome of Arachis, which is based on microsatellite markers and legume anchor markers. These marker types are transferable within the genus Arachis and to other legumes respectively, enabling this map to be aligned to other Arachis maps and to maps of other legume crops including those with sequenced genomes. In total, 34 sequence-confirmed candidate disease resistance genes and five QTLs were mapped.

Conclusion

Candidate genes and QTLs were distributed on all linkage groups except for the smallest, but the distribution was not even. Groupings of candidate genes and QTLs for late leaf spot resistance were apparent on the upper region of linkage group 4 and the lower region of linkage group 2, indicating that these regions are likely to control disease resistance.     

High-resolution association mapping of number of teats in pigs reveals regions controlling vertebral development

Background

Selection pressure on the number of teats has been applied to be able to provide enough teats for the increase in litter size in pigs. Although many QTL were reported, they cover large chromosomal regions and the functional mutations and their underlying biological mechanisms have not yet been identified. To gain a better insight in the genetic architecture of the trait number of teats, we performed a genome-wide association study by genotyping 936 Large White pigs using the Illumina PorcineSNP60 Beadchip. The analysis is based on deregressed breeding values to account for the dense family structure and a Bayesian approach for estimation of the SNP effects.

Results

The genome-wide association study resulted in 212 significant SNPs. In total, 39 QTL regions were defined including 170 SNPs on 13 Sus scrofa chromosomes (SSC) of which 5 regions on SSC7, 9, 10, 12 and 14 were highly significant. All significantly associated regions together explain 9.5% of the genetic variance where a QTL on SSC7 explains the most genetic variance (2.5%). For the five highly significant QTL regions, a search for candidate genes was performed. The most convincing candidate genes were VRTN and Prox2 on SSC7, MPP7, ARMC4, and MKX on SSC10, and vertebrae δ-EF1 on SSC12. All three QTL contain candidate genes which are known to be associated with vertebral development. In the new QTL regions on SSC9 and SSC14, no obvious candidate genes were identified.

Conclusions

Five major QTL were found at high resolution on SSC7, 9, 10, 12, and 14 of which the QTL on SSC9 and SSC14 are the first ones to be reported on these chromosomes. The significant SNPs found in this study could be used in selection to increase number of teats in pigs, so that the increasing number of live-born piglets can be nursed by the sow. This study points to common genetic mechanisms regulating number of vertebrae and number of teats.

Electronic supplementary material

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

QTL analysis of heterostyly in Primula sieboldii and its application for morph identification in wild populations

Background and Aims

Primula sieboldii is a perennial clonal herb that is distributed around the Sea of Japan and is endangered in Japan. Its breeding system is characterized by heteromorphic self-incompatibility, and the morph ratio within a population is very important for reproductive success. The aims of this study were to construct a linkage map, map the S locus as a qualitative trait and quantitative trait loci (QTLs) for floral morphological traits related to heterostyly, and predict the morph type in wild populations by using molecular markers for devising a conservation strategy.

Methods

A linkage map was constructed with 126 markers. The QTLs for four floral traits and the S locus were mapped. Using the genotypes of loci that were located near both the S locus and the QTLs with large effects, morphs of 59 wild genets were predicted.

Key Results

The linkage map consisted of 14 linkage groups (LGs). The S locus was mapped to LG 7. Major QTLs for stigma and anther heights were detected in the same region as the S locus. These QTLs exhibited high logarithm of the odds scores and explained a high percentage of the phenotypic variance (>85 %). By analysing these two traits within each morph, additional QTLs for each trait were detected. Using the four loci linked to the S locus, the morphs of 43 genets in three wild populations could be predicted.

Conclusions

This is the first report of a linkage map and QTL analysis for floral morphology related to heterostyly in P. sieboldii. Floral morphologies related to heterostyly are controlled by the S locus in LG 7 and by several QTLs in other LGs. Additionally, this study showed that molecular markers are effective tools for investigating morph ratios in a population containing the non-flowering individuals or during the non-flowering seasons.     

QTL analysis of root traits as related to phosphorus efficiency in soybean

Background and Aims

Low phosphorus (P) availability is a major constraint to soybean growth and production, especially in tropical and subtropical areas. Root traits have been shown to play critical roles in P efficiency in crops. Identification of the quantitative trait loci (QTLs) conferring superior root systems could significantly enhance genetic improvement in soybean P efficiency.

Methods

A population of 106 F₉ recombinant inbred lines (RILs) derived from a cross between BD2 and BX10, which contrast in both P efficiency and root architecture, was used for mapping and QTL analysis. Twelve traits were examined in acid soils. A linkage map was constructed using 296 simple sequence repeat (SSR) markers with the Kosambi function, and the QTLs associated with these traits were detected by composite interval mapping and multiple-QTL mapping.

Key Results

The first soybean genetic map based on field data from parental genotypes contrasting both in P efficiency and root architecture was constructed. Thirty-one putative QTLs were detected on five linkage groups, with corresponding contribution ratios of 9·1–31·1 %. Thirteen putative QTLs were found for root traits, five for P content, five for biomass and five for yield traits. Three clusters of QTLs associated with the traits for root and P efficiency at low P were located on the B1 linkage group close to SSR markers Satt519 and Satt519-Sat_128, and on the D2 group close to Satt458; and one cluster was on the B1 linkage group close to Satt519 at high P.

Conclusions

Most root traits in soybean were conditioned by more than two minor QTLs. The region closer to Satt519 on the B1 linkage group might have great potential for future genetic improvement for soybean P efficiency through root selection.     

Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc?×?(Landrace × Yorkshire) commercial population

Background

Understanding the genetic mechanisms that underlie meat quality traits is essential to improve pork quality. To date, most quantitative trait loci (QTL) analyses have been performed on F₂ crosses between outbred pig strains and have led to the identification of numerous QTL. However, because linkage disequilibrium is high in such crosses, QTL mapping precision is unsatisfactory and only a few QTL have been found to segregate within outbred strains, which limits their use to improve animal performance. To detect QTL in outbred pig populations of Chinese and Western origins, we performed genome-wide association studies (GWAS) for meat quality traits in Chinese purebred Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) (DLY) commercial population.

Methods

Three hundred and thirty six Chinese Erhualian and 610 DLY pigs were genotyped using the Illumina PorcineSNP60K Beadchip and evaluated for 20 meat quality traits. After quality control, 35 985 and 56 216 single nucleotide polymorphisms (SNPs) were available for the Chinese Erhualian and DLY datasets, respectively, and were used to perform two separate GWAS. We also performed a meta-analysis that combined P-values and effects of 29 516 SNPs that were common to Erhualian, DLY, F₂ and Sutai pig populations.

Results

We detected 28 and nine suggestive SNPs that surpassed the significance level for meat quality in Erhualian and DLY pigs, respectively. Among these SNPs, ss131261254 on pig chromosome 4 (SSC4) was the most significant (P = 7.97E-09) and was associated with drip loss in Erhualian pigs. Our results suggested that at least two QTL on SSC12 and on SSC15 may have pleiotropic effects on several related traits. All the QTL that were detected by GWAS were population-specific, including 12 novel regions. However, the meta-analysis revealed seven novel QTL for meat characteristics, which suggests the existence of common underlying variants that may differ in frequency across populations. These QTL regions contain several relevant candidate genes.

Conclusions

These findings provide valuable insights into the molecular basis of convergent evolution of meat quality traits in Chinese and Western breeds that show divergent phenotypes. They may contribute to genetic improvement of purebreds for crossbred performance.

Electronic supplementary material

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

Quantitative Trait Locus (QTL) meta-analysis and comparative genomics for candidate gene prediction in perennial ryegrass (<Emphasis Type="Italic">Lolium perenne</Emphasis> L.)

Background

In crop species, QTL analysis is commonly used for identification of factors contributing to variation of agronomically important traits. As an important pasture species, a large number of QTLs have been reported for perennial ryegrass based on analysis of biparental mapping populations. Further characterisation of those QTLs is, however, essential for utilisation in varietal improvement programs.

Results

A bibliographic survey of perennial ryegrass trait-dissection studies identified a total of 560 QTLs from previously published papers, of which 189, 270 and 101 were classified as morphology-, physiology- and resistance/tolerance-related loci, respectively. The collected dataset permitted a subsequent meta-QTL study and implementation of a cross-species candidate gene identification approach. A meta-QTL analysis based on use of the BioMercator software was performed to identify two consensus regions for pathogen resistance traits. Genes that are candidates for causal polymorphism underpinning perennial ryegrass QTLs were identified through in silico comparative mapping using rice databases, and 7 genes were assigned to the p150/112 reference map. Markers linked to the Lp DGL1, Lp Ph1 and Lp PIPK1 genes were located close to plant size, leaf extension time and heading date-related QTLs, respectively, suggesting that these genes may be functionally associated with important agronomic traits in perennial ryegrass.

Conclusions

Functional markers are valuable for QTL meta-analysis and comparative genomics. Enrichment of such genetic markers may permit further detailed characterisation of QTLs. The outcomes of QTL meta-analysis and comparative genomics studies may be useful for accelerated development of novel perennial ryegrass cultivars with desirable traits.

     

Unlocking new alleles for leaf rust resistance in the Vavilov wheat collection

Key message

Thirteen potentially new leaf rust resistance loci were identified in a Vavilov wheat diversity panel. We demonstrated the potential of allele stacking to strengthen resistance against this important pathogen.

Abstract

Leaf rust (LR) caused by Puccinia triticina is an important disease of wheat (Triticum aestivum L.), and the deployment of genetically resistant cultivars is the most viable strategy to minimise yield losses. In this study, we evaluated a diversity panel of 295 bread wheat accessions from the N. I. Vavilov Institute of Plant Genetic Resources (St Petersburg, Russia) for LR resistance and performed genome-wide association studies (GWAS) using 10,748 polymorphic DArT-seq markers. The diversity panel was evaluated at seedling and adult plant growth stages using three P. triticina pathotypes prevalent in Australia. GWAS was applied to 11 phenotypic data sets which identified a total of 52 significant marker–trait associations representing 31 quantitative trait loci (QTL). Among them, 29 QTL were associated with adult plant resistance (APR). Of the 31 QTL, 13 were considered potentially new loci, whereas 4 co-located with previously catalogued Lr genes and 14 aligned to regions reported in other GWAS and genomic prediction studies. One seedling LR resistance QTL located on chromosome 3A showed pronounced levels of linkage disequilibrium among markers (r ² = 0.7), suggested a high allelic fixation. Subsequent haplotype analysis for this region found seven haplotype variants, of which two were strongly associated with LR resistance at seedling stage. Similarly, analysis of an APR QTL on chromosome 7B revealed 22 variants, of which 4 were associated with resistance at the adult plant stage. Furthermore, most of the tested lines in the diversity panel carried 10 or more combined resistance-associated marker alleles, highlighting the potential of allele stacking for long-lasting resistance.

     

Identification of QTLs for behavioral reactivity to social separation and humans in sheep using the OvineSNP50 BeadChip

Background

Current trends in sheep farming practices rely on animals with a greater level of behavioral autonomy than before, a phenotype that actively contributes to the sustainability of animal production. Social reactivity and reactivity to humans are relevant behavioral traits in sheep, known for their strong gregariousness and weak tolerance to handling, which have previously been reported with moderate to high heritabilities. To identify loci underlying such behaviors, we performed a genome study in Romane lambs.

Results

The experiment was carried out on 934 male and female lambs allocated into 9 half-sib families (average of 103 lambs per family) and reared outside. After weaning, all the lambs were individually exposed to 4 standardized behavioral tests combining social isolation, exposure to humans or handling, confinement and novelty (i.e. arena test, corridor test, isolation box test, shearing test). A broad range of behaviors including vocalizations, locomotion, vigilance and flight distance, as well as the cortisol response to handling, were collected. All lambs were genotyped using the Illumina OvineSNP50 BeadChip. QTL detection was performed by linkage, association and joint linkage and association analyses using the QTLmap software. Five main QTL regions were identified on sheep chromosomes (Ovis Aries Region, OAR) 12, 16, 19, 21 and 23 among many other QTLs with small to moderate effects. The QTLs on OAR12, 16 and 21 showed significant associations with social reactivity. The QTLs on OAR19 and 23 were found to be associated with reactivity to humans. No overlapping QTLs were identified for the different traits measured in the behavioral tests, supporting the hypothesis that different genetic factors influence social reactivity and tolerance to humans.

Conclusion

The results of this study using ovine SNP data suggest that in domestic sheep the behavioral responses to social separation and exposure to humans are under polygenic influence. The most relevant QTLs reported in the present study contain interesting candidate genes previously described to be associated with various emotional and social behaviors in mammals.

Electronic supplementary material

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

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

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