A Random Model Approach to Interval Mapping of Quantitative Trait Loci

Mapping quantitative trait loci in outbred populations is important because many populations of organisms are noninbred. Unfortunately, information about the genetic architecture of the trait may not be available in outbred populations. Thus, the allelic effects of genes can not be estimated with ease. In addition, under linkage equilibrium, marker genotypes provide no information about the genotype of a QTL (our terminology for a single quantitative trait locus is QTL while multiple loci are referred to as QTLs). To circumvent this problem, an interval mapping procedure based on a random model approach is described. Under a random model, instead of estimating the effects, segregating variances of QTLs are estimated by a maximum likelihood method. Estimation of the variance component of a QTL depends on the proportion of genes identical-by-descent (IBD) shared by relatives at the locus, which is predicted by the IBD of two markers flanking the QTL. The marker IBD shared by two relatives are inferred from the observed marker genotypes. The procedure offers an advantage over the regression interval mapping in terms of high power and small estimation errors and provides flexibility for large sibships, irregular pedigree relationships and incorporation of common environmental and fixed effects.     

Genome-wide association mapping of agronomic traits in sugar beet

Recent results indicate that association mapping in populations from applied plant breeding is a powerful tool to detect QTL which are of direct relevance for breeding. The focus of this study was to unravel the genetic architecture of six agronomic traits in sugar beet. To this end, we employed an association mapping approach, based on a very large population of 924 elite sugar beet lines from applied plant breeding, fingerprinted with 677 single nucleotide polymorphism (SNP) markers covering the entire genome. We show that in this population linkage disequilibrium decays within a short genetic distance and is sufficient for the detection of QTL with a large effect size. To increase the QTL detection power and the mapping resolution a much higher number of SNPs is required. We found that for QTL detection, the mixed model including only the kinship matrix performed best, even in the presence of a considerable population structure. In genome-wide scans, main effect QTL and epistatic QTL were detected for all six traits. Our full two-dimensional epistasis scan revealed that for complex traits there appear to be epistatic master regulators, loci which are involved in a large number of epistatic interactions throughout the genome.     

Congruence in QTL for adventitious rooting in Pinus elliottii × Pinus caribaea hybrids resolves between and within-species effects

Targeting between-species effects for improvement in synthetic hybrid populations derived from outcrossing parental tree species may be one way to increase the efficacy and predictability of hybrid breeding. We present a comparative analysis of the quantitative trait loci (QTL) which resolved between from within-species effects for adventitious rooting in two populations of hybrids between Pinus elliottii and P. caribaea, an outbred F₁ (n=287) and an inbred-like F₂ family (n=357). Most small to moderate effect QTL (each explaining 2–5% of phenotypic variation, PV) were congruent (3 out of 4 QTL in each family) and therefore considered within-species effects as they segregated in both families. A single large effect QTL (40% PV) was detected uniquely in the F₂ family and assumed to be due to a between-species effect, resulting from a genetic locus with contrasting alleles in each parental species. Oligogenic as opposed to polygenic architecture was supported in both families (60% and 20% PV explained by 4 QTL in the F₂ and F₁ respectively). The importance of adventitious rooting for adaptation to survive water-logged environments was thought in part to explain oligogenic architecture of what is believed to be a complex trait controlled by many hundreds of genes.     

Multipoint genetic mapping of quantitative trait loci with dominant markers in outbred populations

We present a multipoint algorithm for mapping quantitative trait loci (QTLs) using dominant markers. The algorithm is designed for outbred populations and is particularly suited for large families. The algorithm works with either codominant or dominant markers, either of which may be interspersed within the same linkage map. Concurrently, the algorithm also partitions dominance variance at the QTL. Computer simulations show that with large families, QTL mapping with dominant markers can be almost as powerful as with bi-allelic, codominant markers. Yet despite this, other situations show a large standard deviation in the estimate of the QTL position, thus making QTL mapping with dominant markers in outbred populations a useful detection tool, albeit limited in its resolution. This revised version was published online in July 2006 with corrections to the Cover Date.     

Detection of QTL for flowering time in multiple families of elite maize

Flowering time is a fundamental quantitative trait in maize that has played a key role in the postdomestication process and the adaptation to a wide range of climatic conditions. Flowering time has been intensively studied and recent QTL mapping results based on diverse founders suggest that the genetic architecture underlying this trait is mainly based on numerous small-effect QTL. Here, we used a population of 684 progenies from five connected families to investigate the genetic architecture of flowering time in elite maize. We used a joint analysis and identified nine main effect QTL explaining approximately 50?% of the genotypic variation of the trait. The QTL effects were small compared with the observed phenotypic variation and showed strong differences between families. We detected no epistasis with the genetic background but four digenic epistatic interactions in a full 2-dimensional genome scan. Our results suggest that flowering time in elite maize is mainly controlled by main effect QTL with rather small effects but that epistasis may also contribute to the genetic architecture of the trait.     

Quantitative trait loci affecting life span in replicated populations of Drosophila melanogaster. II. Response to selection

Three selection experiments were used to identify chromosome regions that contain QTL affecting late-life and early-life fitness in Drosophila melanogaster. The selection experiments were initiated by crossing pairs of inbred lines that had been derived from outbred laboratory populations that had different mean life spans. QTL regions were located by association with microsatellite markers that showed significant selection responses. Regions between recombination map positions 54 and 81 on chromosome 2, between 0 and 30 on chromosome 3, and near locations 49 and 81 on chromosome 3 had the strongest support as locations of life-span QTL. There was good general agreement between the life-span QTL regions that were identified by selection and those that were identified in a companion recombination mapping experiment that used the same fly stocks. Many marker loci responded in opposite directions to selection for late- and early-life fitness, indicating negative genetic correlations or trade-offs between those traits. Indirect evidence suggested that some negative genetic correlations were due to antagonistic pleiotropy.     

A genome-wide association study reveals loci influencing height and other conformation traits in horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ～200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ～40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ～18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.     

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.     

The genetic basis of divergent pigment patterns in juvenile threespine sticklebacks

Animal pigment patterns are important for a range of functions, including camouflage and communication. Repeating pigment patterns, such as stripes, bars and spots have been of particular interest to developmental and theoretical biologists, but the genetic basis of natural variation in such patterns is largely unexplored. In this study, we identify a difference in a periodic pigment pattern among juvenile threespine sticklebacks (Gasterosteus aculeatus) from different environments. Freshwater sticklebacks exhibit prominent vertical bars that visually break up the body shape, but sticklebacks from marine populations do not. We hypothesize that these distinct pigment patterns are tuned to provide crypsis in different habitats. This phenotypic difference is widespread and appears in most of the freshwater populations that we sampled. We used quantitative trait locus (QTL) mapping in freshwater-marine F2 hybrids to elucidate the genetic architecture underlying divergence in this pigmentation pattern. We identified two QTL that were significantly associated with variation in barring. Interestingly, these QTL were associated with two distinct aspects of the pigment pattern: melanophore number and overall pigment level. We compared the QTL locations with positions of known pigment candidate genes in the stickleback genome. We also identified two major QTL for juvenile body size, providing new insights into the genetic basis of juvenile growth rates in natural populations. In summary, although there is a growing literature describing simple genetic bases for adaptive coloration differences, this study emphasizes that pigment patterns can also possess a more complex genetic architecture.     

High resolution mapping and identification of new quantitative trait loci (QTL) affecting susceptibility to Marek's disease

Marek's disease (MD) is a lymphoproliferative disease of chickens that costs the poultry industry approximately $1 billion annually. Genetic resistance to MD is gaining increased attention to augment vaccinal control as disease outbreaks occur more frequently. Previously, analysis of a 272 F2 White Leghorn resource population measured for many MD traits and genotyped for 78 microsatellite markers revealed two and four quantitative trait loci (QTL) with significant and suggestive association, respectively, to one or more MD associated traits. Additional genetic markers have since been scored on the MD resource population to increase QTL resolution and genome coverage. Saturation of four of the QTL regions with 17 markers revealed five new QTL while 32 markers extended the genome coverage by 400 + CM and uncovered three more QTL. QTL analysis by single-point and interval mapping algorithms agreed well when marker saturation was approximately 20 CM or less. Currently 127 genetic markers cover approximately 68% of the genome that contain up to 14 MD QTL associated to one or more MD trait; seven at the significant level and seven at the suggestive level. Individually each QTL accounts for 2-10% of the variation and, in general, resistance was dominant although the resistant allele may come from either parental line. This study suggests that a limited number of genomic regions play a major role in the genetic control of MD resistance. Markers linked to these loci may be useful for selection of MD resistant stock by the poultry industry following verification of the association within their breeding populations.     

Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in Mimulus guttatus

Evolutionary biologists seek to understand the genetic basis for multivariate phenotypic divergence. We constructed an F2 mapping population (N = 539) between two distinct populations of Mimulus guttatus. We measured 20 floral, vegetative, and life-history characters on parents and F1 and F2 hybrids in a common garden experiment. We employed multitrait composite interval mapping to determine the number, effect, and degree of pleiotropy in quantitative trait loci (QTL) affecting divergence in floral, vegetative, and life-history characters. We detected 16 QTL affecting floral traits; 7 affecting vegetative traits; and 5 affecting selected floral, vegetative, and life-history traits. Floral and vegetative traits are clearly polygenic. We detected a few major QTL, with all remaining QTL of small effect. Most detected QTL are pleiotropic, implying that the evolutionary shift between these annual and perennial populations is constrained. We also compared the genetic architecture controlling floral trait divergence both within (our intraspecific study) and between species, on the basis of a previously published analysis of M. guttatus and M. nasutus. Eleven of our 16 floral QTL map to approximately the same location in the interspecific map based on shared, collinear markers, implying that there may be a shared genetic basis for floral divergence within and among species of Mimulus.     

QTL mapping reveals a striking coincidence in the positions of genomic regions associated with adaptive variation in body size in parallel clines of Drosophila melanogaster on different continents

Latitudinal genetic clines in body size are common in many ectotherm species and are attributed to climatic adaptation. Here, we use Quantitative Trait Loci (QTL) mapping to identify genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from extreme ends of a cline in South America. Our results show that there is a significant association between the positions of QTL with strong effects on wing area in South America and those previously reported in a QTL mapping study of Australian cline end populations (P < 0.05). In both continents, the right arm of the third chromosome is associated with QTL with the strongest effect on wing area. We also show that QTL peaks for wing area and thorax length are associated with the same genomic regions, indicating that the clinal variation in the body size traits may have a similar genetic basis. The consistency of the results found for the South American and Australian cline end populations indicate that the genetic basis of the two clines may be similar and future efforts to identify the genes producing the response to selection should be focused on the genomic regions highlighted by the present work.     

Tests of candidate genes in breed cross populations for QTL mapping in livestock

In recent years, several F₂ crosses between outbred lines of livestock have been developed to identify quantitative trait loci (QTL). These populations are valuable for further genetic analysis, including positional candidate gene loci (CGL). Analysis of CGL in F₂ populations is, however, hindered by extensive between-breed linkage disequilibrium (LD). The objectives here were to develop and evaluate three tests for CGL in simulated F₂ breed-cross populations. 1) A standard association test, based on the fixed effect of CGL genotype. This test was significant for CGL at considerable distances from the QTL. 2) A marker-assisted association test, based on a test at the CGL of the fixed effect of CGL genotype in a breed-cross QTL interval mapping model. This removed the impact of between-breed LD, but was not powerful in detecting CGL closely linked to the QTL, unless the CGL was the QTL. 3) An F-drop test, comparing F ratios for a QTL at the CGL with and without the CGL included as fixed effect. It had low power to distinguish close from distant CGL. Power to distinguish two CGL within 10 cM from the QTL was limited and little improved by including QTL effects associated with markers to remove between-breed LD, although the power was greater when one of the CGL was the causative mutation. Therefore, while we conclude that candidate gene tests in QTL mapping populations must be interpreted with caution, we now have a clearer picture of the value of candidate gene tests in these populations.     

Do large effect QTL fractionate? A case study at the maize domestication QTL teosinte branched1

Quantitative trait loci (QTL) mapping is a valuable tool for studying the genetic architecture of trait variation. Despite the large number of QTL studies reported in the literature, the identified QTL are rarely mapped to the underlying genes and it is usually unclear whether a QTL corresponds to one or multiple linked genes. Similarly, when QTL for several traits colocalize, it is usually unclear whether this is due to the pleiotropic action of a single gene or multiple linked genes, each affecting one trait. The domestication gene teosinte branched1 (tb1) was previously identified as a major domestication QTL with large effects on the differences in plant and ear architecture between maize and teosinte. Here we present the results of two experiments that were performed to determine whether the single gene tb1 explains all trait variation for its genomic region or whether the domestication QTL at tb1 fractionates into multiple linked QTL. For traits measuring plant architecture, we detected only one QTL per trait and these QTL all mapped to tb1. These results indicate that tb1 is the sole gene for plant architecture traits that segregates in our QTL mapping populations. For most traits related to ear morphology, we detected multiple QTL per trait in the tb1 genomic region, including a large effect QTL at tb1 itself plus one or two additional linked QTL. tb1 is epistatic to two of these additional QTL for ear traits. Overall, these results provide examples for both a major QTL that maps to a single gene, as well as a case in which a QTL fractionates into multiple linked QTL.     

Use of a novel outbred by inbred F1 cross to detect genetic markers for growth

A unique outbred by inbred F1 resource population was established. The population structure facilitated the unique opportunity of examining gene by genetic background interaction through crossing two modern broiler sires with dams from two unrelated inbred lines, with no selection for growth rate, to produce about 600 F1 chicks. Pools of DNA were generated from the phenotypic extremes (20% high and low) for 8-week body weight for each of the four combinations of sire and dam line. For one sire family, pools were also separately generated for each sex. The pools were genoyped with 25 informative (segregating) microsatellites. This unique F1 cross between outbred and inbred populations allowed use of the inbred alleles as an 'internal control' for polymerase chain reaction amplification quality in DNA pools. Ten microsatellites showed marked differences (P < 0.05) in allele frequencies between high and low pools, suggesting an association between marker and quantitative trait loci (QTL). These differences were verified using selective genotyping. For many markers, differences in allele frequencies between the high and the low pools, or marker effect, varied between the two dam lines and the two sexes, suggesting an interaction between some genes and the genetic background as represented by different dam lines or sexes. The suggestive marker-QTL associations identified in this F1 population demonstrate the efficiency of this population design while different QTL effects in different genetic line crosses and sexes highlight the importance of gene by genetic background interaction in QTL detection.     

A Further Look at Porcine Chromosome 7 Reveals VRTN Variants Associated with Vertebral Number in Chinese and Western Pigs

The number of vertebrae is an economically important trait that affects carcass length and meat production in pigs. A major quantitative trait locus (QTL) for thoracic vertebral number has been repeatedly identified on pig chromosome (SSC) 7. To dissect the genetic basis of the major locus, we herein genotyped a large sample of animals from 3 experimental populations of Chinese and Western origins using 60K DNA chips. Genome-wide association studies consistently identified the locus across the 3 populations and mapped the locus to a 947-Kb region on SSC7. An identical-by-descent sharing assay refined the locus to a 100-Kb segment that harbors only two genes including VRTN and SYNDIG1L. Of them, VRNT has been proposed as a strong candidate of the major locus in Western modern breeds. Further, we resequenced the VRTN gene using DNA samples of 35 parental animals with known QTL genotypes by progeny testing. Concordance tests revealed 4 candidate causal variants as their genotypes showed the perfect segregation with QTL genotypes of the tested animals. An integrative analysis of evolutional constraints and functional elements supported two VRTN variants in a complete linkage disequilibrium phase as the most likely causal mutations. The promising variants significantly affect the number of thoracic vertebrae (one vertebra) in large scale outbred animals, and are segregating at rather high frequencies in Western pigs and at relatively low frequencies in a number of Chinese breeds. Altogether, we show that VRTN variants are significantly associated with the number of thoracic vertebrae in both Chinese and Western pigs. The finding advances our understanding of the genetic architecture of the vertebral number in pigs. Furthermore, our finding is of economical importance as it provides a robust breeding tool for the improvement of vertebral number and meat production in both Chinese indigenous pigs and Western present-day commercial pigs.     

On the efficiency of composite interval mapping in an outbred population

Composite interval mapping (CIM) has been successfully applied to the detection of QTL in experimental animals and plants. However, practical analyses based on CIM have been reported mainly for populations derived from cross between inbred lines. There are few studies on QTL analyses with CIM in outbred populations. To evaluate the applicability of CIM to outbred populations is prerequisite for the fine mapping of QTL in industrial animals such as pig and chicken. Some markers are usually not fully informative in outbred populations. In application of CIM to outbred populations, the influence of inclusion of such uninformative markers used as covariates on the efficiency of CIM should be investigated. In this paper a least-squares method for CIM was formalized in an F(2) population derived by crossing two outbred lines. The efficiencies of CIM were evaluated for outbred populations in comparison with simple interval mapping (SIM) for several cases of marker informativeness using simulations. By incorporating markers linked to a tested position as well as those unlinked, CIM showed a higher efficiency to separate two linked QTL over SIM. The efficiency of dissection was enhanced as the marker informativeness was increased. The power of CIM to detect an isolated QTL was improved by excluding markers linked to a tested position from covariates and higher than SIM regardless of marker informativeness. In conclusion, CIM is a useful procedure for the analysis of QTL in outbred populations even under low marker informativeness.     

Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse

The extent to which epistasis contributes to adaptation, population differentiation, and speciation is a long-standing and important problem in evolutionary genetics. Using recombinant inbred (RI) lines of Arabidopsis thaliana grown under natural field conditions, we have examined the genetic architecture of fitness-correlated traits with respect to epistasis; we identified both single-locus additive and two-locus epistatic QTL for natural variation in fruit number, germination, and seed length and width. For fruit number, we found seven significant epistatic interactions, but only two additive QTL. For seed germination, length, and width, there were from two to four additive QTL and from five to eight epistatic interactions. The epistatic interactions were both positive and negative. In each case, the magnitude of the epistatic effects was roughly double that of the effects of the additive QTL, varying from -41% to +29% for fruit number and from -5% to +4% for seed germination, length, and width. A number of the QTL that we describe participate in more than one epistatic interaction, and some loci identified as additive also may participate in an epistatic interaction; the genetic architecture for fitness traits may be a network of additive and epistatic effects. We compared the map positions of the additive and epistatic QTL for germination, seed width, and seed length from plants grown in both the field and the greenhouse. While the total number of significant additive and epistatic QTL was similar under the two growth conditions, the map locations were largely different. We found a small number of significant epistatic QTL x environment effects when we tested directly for them. Our results support the idea that epistatic interactions are an important part of natural genetic variation and reinforce the need for caution in comparing results from greenhouse-grown and field-grown plants.