Clinal variation in phyB2, a candidate gene for day-length-induced growth cessation and bud set, across a latitudinal gradient in European aspen (Populus tremula)

The initiation of growth cessation and dormancy represents a critical ecological and evolutionary trade-off between survival and growth in most forest trees. The most important environmental cue regulating the initiation of dormancy is a shortening of the photoperiod and phytochrome genes have been implicated in short-day-induced bud set and growth cessation in Populus. We characterized patterns of DNA sequence variation at the putative candidate gene phyB2 in 4 populations of European aspen (Populus tremula) and scored single-nucleotide polymorphisms in an additional 12 populations collected along a latitudinal gradient in Sweden. We also measured bud set from a subset of these trees in a growth chamber experiment. Buds set showed significant clinal variation with latitude, explaining approximately 90% of the population variation in bud set. A sliding-window scan of phyB2 identified six putative regions with enhanced population differentiation and four SNPs showed significant clinal variation. The clinal variation at individual SNPs is suggestive of an adaptive response in phyB2 to local photoperiodic conditions. Three of four SNPs showing clinal variation were located in regions with excessive genetic differentiation, demonstrating that searching for regions of high genetic differentiation can be useful for identifying sites putatively involved in local adaptation.     

Nucleotide polymorphism and linkage disequilibrium within and among natural populations of European aspen (Populus tremula L., Salicaceae)

Populus is an important model organism in forest biology, but levels of nucleotide polymorphisms and linkage disequilibrium have never been investigated in natural populations. Here I present a study on levels of nucleotide polymorphism, haplotype structure, and population subdivision in five nuclear genes in the European aspen Populus tremula. Results show substantial levels of genetic variation. Levels of silent site polymorphisms, pi(s), averaged 0.016 across the five genes. Linkage disequilibrium was generally low, extending only a few hundred base pairs, suggesting that rates of recombination are high in this obligate outcrossing species. Significant genetic differentiation was found at all five genes, with an average estimate of F(ST) = 0.116. Levels of polymorphism in P. tremula are 2- to 10-fold higher than those in other woody, long-lived perennial plants, such as Pinus and Cryptomeria. The high levels of nucleotide polymorphism and low linkage disequilibrium suggest that it may be possible to map functional variation to very fine scales in P. tremula using association-mapping approaches.     

An excess of nonsynonymous polymorphism and extensive haplotype structure at the PtABI1B locus in European aspen (Populus tremula): a case of balancing selection in an obligately outcrossing plant?

Here, we describe an unusually pronounced haplotype structure at the PtABI1B locus in the obligately outcrossing tree Populus tremula. Both nucleotide diversity and divergence at PtABI1B was low compared to other P. tremula genes suggesting that the gene is located in a region with a low mutation rate. Despite this, PtABI1B shows a very marked excess of nonsynonymous polymorphisms across the entire coding region and linkage disequilibrium (LD) extending across the entire PtABI1B region of approximately 2.6 kb. Such extensive LD is normally not seen in P. tremula. The extensive LD at PtABI1B is caused by the presence of two distinct haplotypes. The haplotype structure is not caused by a lack of recombination in the region, because evidence of recombination can be detected. In addition, several statistical tests strongly reject neutrality for the PtABI1B region, suggesting that the unusual haplotype structure could be actively maintained by balancing selection.     

Genetics of sputum gene expression in chronic obstructive pulmonary disease

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.     

Use of Ecotilling as an efficient SNP discovery tool to survey genetic variation in wild populations of Populus trichocarpa

Abstract Ecotilling was used as a simple nucleotide polymorphism (SNP) discovery tool to examine DNA variation in natural populations of the western black cottonwood, Populus trichocarpa, and was found to be more efficient than sequencing for large-scale studies of genetic variation in this tree. A publicly available, live reference collection of P. trichocarpa from the University of British Columbia Botanical Garden was used in this study to survey variation in nine different genes among individuals from 41 different populations. A large amount of genetic variation was detected, but the level of variation appears to be less than in the related species, Populus tremula, based on reported statistics for that tree. Genes examined varied considerably in their level of variation, from PoptrTB1 which had a single SNP, to PoptrLFY which had more than 23 in the 1000-bp region examined. Overall nucleotide diversity, measured as (Total), was relatively low at 0.00184. Linkage disequilibrium, on the other hand, was higher than reported for some woody plant species, with mean r2 equal to 0.34. This study reveals the potential of Ecotilling as a rapid genotype discovery method to explore and utilize the large pool of genetic variation in tree species.     

Multilocus patterns of nucleotide polymorphism and the demographic history of Populus tremula

I have studied nucleotide polymorphism and linkage disequilibrium using multilocus data from 77 fragments, with an average length of fragments of 550 bp, in the deciduous tree Populus tremula (Salicaceae). The frequency spectrum across loci showed a modest excess of mutations segregating at low frequency and a marked excess of high-frequency derived mutations at silent sites, relative to neutral expectations. These excesses were also seen at replacement sites, but were not so pronounced for high-frequency derived mutations. There was a marked excess of low-frequency mutations at replacement sites, likely indicating deleterious amino acid-changing mutations that segregate at low frequencies in P. tremula. I used approximate Bayesian computation (ABC) to evaluate a number of different demographic scenarios and to estimate parameters for the best-fitting model. The data were found to be consistent with a historical reduction in the effective population size of P. tremula through a bottleneck. The timing inferred for this bottleneck is largely consistent with geological data and with data from several other long-lived plant species. The results show that P. tremula harbors substantial levels of nucleotide polymorphism with the posterior mode of the scaled mutation rate, = 0.0177 across loci. The ABC analyses also provided an estimate of the scaled recombination rate that indicates that recombination rates in P. tremula are likely to be 2-10 times higher than the mutation rate. This study reinforces the notion that linkage disequilibrium is low and decays to negligible levels within a few hundred base pairs in P. tremula.     

The adaptive potential of Populus balsamifera L. to phenology requirements in a warmer global climate

The manner in which organisms adapt to climate change informs a broader understanding of the evolution of biodiversity as well as conservation and mitigation plans. We apply common garden and association mapping approaches to quantify genetic variance and identify loci affecting bud flush and bud set, traits that define a tree's season for height growth, in the boreal forest tree Populus balsamifera L. (balsam poplar). Using data from 478 genotypes grown in each of two common gardens, one near the southern edge and another near the northern edge of P. balsamifera's range, we found that broad‐sense heritability for bud flush and bud set was generally high (H² > 0.5 in most cases), suggesting that abundant genetic variation exists for phenological response to changes in the length of the growing season. To identify the molecular genetic basis of this variation, we genotyped trees for 346 candidate single nucleotide polymorphisms (SNPs) from 27 candidate genes for the CO/FT pathway in poplar. Mixed‐model analyses of variance identified SNPs in 10 genes to be associated with variation in either bud flush or bud set. Multiple SNPs within FRIGIDA were associated with bud flush, whereas multiple SNPs in LEAFY and GIGANTEA 5 were associated with bud set. Although there was strong population structure in stem phenology, the geographic distribution of multilocus association SNP genotypes was widespread except at the most northern populations, indicating that geographic regions may harbour sufficient diversity in functional genes to facilitate adaption to future climatic conditions in many sites.     

Comparison of the power between microsatellite and single-nucleotide polymorphism markers for linkage and linkage disequilibrium mapping of an electrophysiological phenotype

We performed linkage and linkage disequilibrium (LD) mapping analyses to compare the power between microsatellite and single nucleotide polymorphism (SNP) markers. Chromosome-wide analyses were performed for a quantitative electrophysiological phenotype, ttth1, on chromosome 7. Multipoint analysis of microsatellite markers using the variance component (VC) method showed the highest LOD score of 4.20 at 162 cM, near D7S509 (163.7 cM). Two-point analysis of SNPs using the VC method yielded the highest LOD score of 3.98 in the Illumina SNP data and 3.45 in the Affymetrix SNP data around 152-153 cM. In family-based single SNP and SNP haplotype LD analysis, we identified seven SNPs associated with ttth1. We searched for any potential candidate genes in the location of the seven SNPs. The SNPs rs1476640 and rs768055 are located in the FLJ40852 gene (a hypothetical protein), and SNP rs1859646 is located in the TAS2R5 gene (a taste receptor). The other four SNPs are not located in any known or annotated genes. We found the high density SNP scan to be superior to microsatellites because it is effective in downstream fine mapping due to a better defined linkage region. Our study proves the utility of high density SNP in genome-wide mapping studies.     

Admixture in European Populus hybrid zones makes feasible the mapping of loci that contribute to reproductive isolation and trait differences

The use of admixed human populations to scan the genome for chromosomal segments affecting complex phenotypic traits has proved a powerful analytical tool. However, its potential in other organisms has not yet been evaluated. Here, we use DNA microsatellites to assess the feasibility of this approach in hybrid zones between two members of the 'model tree' genus Populus: Populus alba (white poplar) and Populus tremula (European aspen). We analyzed samples of both species and a Central European hybrid zone (N=544 chromosomes) for a genome-wide set of 19 polymorphic DNA microsatellites. Our results indicate that allele frequency differentials between the two species are substantial (mean delta=0.619+/-0.067). Background linkage disequilibrium (LD) in samples of the parental gene pools is moderate and should respond to sampling schemes that minimize drift and account for rare alleles. LD in hybrids decays with increasing number of backcross generations as expected from theory and approaches background levels of the parental gene pools in advanced generation backcrosses. Introgression from P. tremula into P. alba varies strongly across marker loci. For several markers, alleles from P. tremula are slightly over-represented relative to neutral expectations, whereas a single locus exhibits evidence of selection against P. tremula genotypes. We interpret our results in terms of the potential for admixture mapping in these two ecologically divergent Populus species, and we validate a modified approach of studying genotypic clines in 'mosaic' hybrid zones.     

Comparative LD mapping using single SNPs and haplotypes identifies QTL for plant height and biomass as secondary traits of drought tolerance in maize

Drought often delays developmental events so that plant height and above-ground biomass are reduced, resulting in yield loss due to inadequate photosynthate. In this study, plant height and biomass measured by the Normalized Difference Vegetation Index (NDVI) were used as criteria for drought tolerance. A total of 305 lines representing temperate, tropical and subtropical maize germplasm were genotyped using two single nucleotide polymorphism (SNP) chips each containing 1536 markers, from which 2052 informative SNPs and 386 haplotypes each constructed with two or more SNPs were used for linkage disequilibrium (LD) or association mapping. Single SNP- and haplotype-based LD mapping identified two significant SNPs and three haplotype loci [a total of four quantitative trait loci (QTL)] for plant height under well-watered and water-stressed conditions. For biomass, 32 SNPs and 12 haplotype loci (30 QTL) were identified using NDVIs measured at seven stages under the two water regimes. Some significant SNP and haplotype loci for NDVI were shared by different stages. Comparing significant loci identified by single SNP- and haplotype-based LD mapping, we found that six out of the 14 chromosomal regions defined by haplotype loci each included at least one significant SNP for the same trait. Significant SNP haplotype loci explained much higher phenotypic variation than individual SNPs. Moreover, we found that two significant SNPs (two QTL) and one haplotype locus were shared by plant height and NDVI. The results indicate the power of comparative LD mapping using single SNPs and SNP haplotypes with QTL shared by plant height and biomass as secondary traits for drought tolerance in maize.     

No evidence for an association between common nonsynonymous polymorphisms in delta and bristle number variation in natural and laboratory populations of Drosophila melanogaster

We test the hypothesis that naturally occurring nonsynonymous variants in the Delta ligand of the Notch signaling pathway contribute to standing variation in sternopleural and/or abdominal bristle number in Drosophila melanogaster, for both a large cohort of wild-caught flies and previously described laboratory lines. We sequenced the transcribed region of Delta for 16 naturally occurring chromosomes and 65 SNPs, including 7 nonsynonymous SNPs (nsSNPs), were observed. Identified nsSNPs and 6 additional common SNPs, all located in exon 6 and the 3' UTR, were genotyped in 2060 wild-caught flies using an OLA-based methodology and genotyped in 38 additional natural chromosomes via DNA sequencing. None of the genotyped nsSNPs were significantly associated with natural variation in bristle number as assessed by a permutation test. A 95% upper bound on the additive genetic variance attributable to each genotyped SNP in the large natural cohort is <2% of the total phenotypic variation. Results suggest that two previously detected genotype/phenotype associations between bristle number and variants in the introns of Delta cannot be explained by linkage disequilibrium between these variants and nearby nonsynonymous variants. Unidentified regulatory variants more parsimoniously explain previous observations.     

Comparative Nucleotide Diversity Across North American and European Populus Species

Nucleotide polymorphisms in two North American balsam poplars (Populus trichocarpa Torr. & Gray and P. balsamifera L.; section Tacamahaca), and one Eurasian aspen (P. tremula L.; section Populus) were compared using nine loci involved in defense, stress response, photoperiodism, freezing tolerance, and housekeeping. Nucleotide diversity varied among species and was highest for P. tremula (θ(w) = 0.005, π(T) = 0.007) as compared to P. balsamifera (θ(w) = 0.004, π(T) = 0.005) or P. trichocarpa (θ(w) = 0.002, π(T) = 0.003). Across species, the defense and the stress response loci accounted for the majority of the observed level of nucleotide diversity. In general, the studied loci did not deviate from neutral expectation either at the individual locus (non-significant normalized Fay and Wu's H) or at the multi-locus level (non-significant HKA test). Using molecular clock analysis, section Tacamahaca probably shared a common ancestor with section Populus approximately 4.5 million year ago. Divergence between the two closely related balsam poplars was about 0.8 million years ago, a pattern consistent with an isolation-with-migration (IM) model. As expected, P. tremula showed a five-fold higher substitution rate (2 × 10(-8) substitution/site/year) compared to the North American species (0.4 × 10(-8) substitution/site/year), probably reflecting its complex demographic history. Linkage disequilibrium (LD) varied among species with a more rapid decay in the North American species (<400 bp) in comparison to P. tremula (?400 bp). The similarities in nucleotide diversity pattern and LD decay of the two balsam poplar species likely reflects the recent time of their divergence.     

Nucleotide Diversity and Haplotype Structure of the Human Angiotensinogen Gene in Two Populations

Variation in the angiotensinogen gene, AGT, has been associated with variation in plasma angiotensinogen levels. In addition, the T235M polymorphism in the AGT product is associated with an increased risk of essential hypertension in multiple populations, making AGT a good example of a quantitative-trait locus underlying susceptibility to a common disease. To better understand genetic variation in AGT, we sequenced a 14.4-kb genomic region spanning the entire AGT and identified 44 single-nucleotide polymorphisms (SNPs). Forty-two SNPs were observed both in 88 white and in 77 Japanese unselected subjects. Six major haplotypes accounted for most of the variation in this region, indicating less allelic complexity than in many other genomic regions. Although the two populations were found to share all of the major AGT haplotypes, there were substantial differences in haplotype frequencies. Pairwise linkage disequilibrium (LD), measured by the D', r(2), and d(2) statistics, demonstrated a general pattern of decline with increasing distance, but, as expected in a small genomic region, individual LD values were highly variable. LD between T235M and each of the other 39 SNPs was assessed in order to model the usefulness of LD to detect a disease-associated mutation. Among the Japanese subjects, 13 (33%) of the SNPs had r(2) values >0.1, whereas this statistic was substantially higher for the white subjects (occurring in 35/39 [90%]). LD between a hypertension-associated promoter mutation, A-6G, and 39 SNPs was also measured. Similar results were obtained, with 33% of the SNPs showing r(2)>0.1 in the Japanese subjects and 92% of the SNPs showing r(2)>0.1 in the white subjects. This difference, which occurs despite an overall similarity in LD patterns in the two populations, reflects a much higher frequency of the M235-associated haplotype in the white sample. These results have important implications for the usefulness of LD approaches in the mapping of genes underlying susceptibility to complex diseases.     

High diversity of the chicken growth hormone gene and effects on growth and carcass traits

The chicken growth hormone (cGH) gene plays a crucial role in controlling growth and metabolism, leading to potential correlations between cGH polymorphisms and economic traits. In this study, DNA from four divergent chicken breeds were screened for single nucleotide polymorphisms (SNPs) in the cGH gene using denaturing high-performance liquid chromatography and sequencing. A total of 46 SNPs were identified, of which 4 were in the 5' untranslated region, 1 in the 3' untranslated region, 5 in exons (two of which are nonsynonymous), with the remaining 36 in introns. The nucleotide diversity in the cGH gene ( theta = 2.7 x 10(-3)) was higher than that reported for other chicken genes, even within the same breeds. The associations of five of these SNPs and their haplotypes with chicken growth and carcass traits were determined using polymerase chain reaction-restriction fragment length polymorphism analysis in a F2 resource population cross of two of the four chicken breeds (White Recessive Rock and Xinghua). This analysis shows that, among other correlations, G+1705A was significantly associated with body weight at all ages measured, shank length at three of four ages measured, and average daily gain within weeks 0 to 4. Thus, this cGH polymorphism, or another polymorphism that is in linkage disequilibrium with G+1705A, appears to correspond to a significant growth-related quantitative trait locus difference between the two breeds used to construct the resource population.     

Evidence that the phytochrome gene family in black cottonwood has one PHYA locus and two PHYB loci but lacks members of the PHYC/F and PHYE subfamilies

The phytochrome photoreceptors play important roles in the photoperiodic control of vegetative bud set, growth cessation, dormancy induction, and cold-hardiness in trees. Interestingly, ecotypic differences in photoperiodic responses are observed in many temperate- zone tree species. Northern and southern ecotypes of black cottonwood (Populus trichocarpa Torr. & Gray), for example, exhibit marked differences in the timing of short-day-induced bud set and growth cessation, and these responses are controlled by phytochrome. Therefore, as a first step toward determining the molecular genetic basis of photoperiodic ecotypes in trees, we characterized the phytochrome gene (PHY) family in black cottonwood. We recovered fragments of one PHYA and two PHYB using PCR-based cloning and by screening a genomic library. Results from Southern analyses confirmed that black cottonwood has one PHYA locus and two PHYB loci, which we arbitrarily designated PHYB1 and PHYB2. Phylogenetic analyses which included PHY from black cottonwood, Arabidopsis thaliana and tomato (Solanum lycopersicum) suggest that the PHYB/D duplications in these species occurred independently. When Southern blots were probed with PHYC, PHYE, and PHYE heterologous probes, the strongest bands that we detected were those of black cottonwood PHYA and/or PHYB. These results suggest that black cottonwood lacks members of the PHYC/F and PHYE subfamilies. Although black cottonwood could contain additional PHY that are distantly related to known angiosperm PHY, our results imply that the PHY family of black cottonwood is less complex than that of other well-characterized dicot species such as Arabidopsis and tomato. Based on Southern analyses of five black cottonwood genotypes representing three photoperiodic ecotypes, substantial polymorphism was detected for at least one of the PHYB loci but not for the PHYA locus. The novel character of the PHY family in black cottonwood, as well as the differences in polymorphism we observed between the PHYA and PHYB subfamilies, indicates that a number of fundamental macro- and microevolutionary questions remain to be answered about the PHY family in dicots.      

EphB2 SNPs and sporadic prostate cancer risk in African American men

The EphB2 gene has been implicated as a tumor suppressor gene somatically altered in both prostate cancer (PC) and colorectal cancer. We have previously shown an association between an EphB2 germline nonsense variant and risk of familial prostate cancer among African American Men (AAM). Here we set out to test the hypothesis that common variation within the EphB2 locus is associated with increased risk of sporadic PC in AAM. We genotyped a set of 341 single nucleotide polymorphisms (SNPs) encompassing the EphB2 locus, including known and novel coding and noncoding variants, in 490 AA sporadic PC cases and 567 matched controls. Single marker-based logistical regression analyses revealed seven EphB2 SNPs showing statistically significant association with prostate cancer risk in our population. The most significant association was achieved for a novel synonymous coding SNP, TGen-624, (Odds Ratio (OR) =?0.22; 95% Confidence Interval (CI) 0.08-0.66, p?=?1×10(-5)). Two other SNPs also show significant associations toward a protective effect rs10465543 and rs12090415 (p?=?1×10(-4)), OR?=?0.49 and 0.7, respectively. Two additional SNPs revealed trends towards an increase in risk of prostate cancer, rs4612601 and rs4263970 (p?=?0.001), OR?=?1.35 and 1.31, respectively. Furthermore, haplotype analysis revealed low levels of linkage disequilibrium within the region, with two blocks being associated with prostate cancer risk among our population. These data suggest that genetic variation at the EphB2 locus may increase risk of sporadic PC among AAM.     

Nucleotide diversity and linkage disequilibrium in cold-hardiness- and wood quality-related candidate genes in Douglas fir

Nuclear sequence variation and linkage disequilibrium (LD) were studied in 15 cold-hardiness- and 3 wood quality-related candidate genes in Douglas fir [Pseudotsuga menziesii (Mirb.) Franco]. This set of genes was selected on the basis of its function in other plants and collocation with cold-hardiness-related quantitative trait loci (QTL). The single-nucleotide polymorphism (SNP) discovery panel represented 24 different trees from six regions in Washington and Oregon plus parents of a segregating population used in the QTL study. The frequency of SNPs was one SNP per 46 bp across coding and noncoding regions on average. Haplotype and nucleotide diversities were also moderately high with H(d) = 0.827 +/- 0.043 and pi = 0.00655 +/- 0.00082 on average, respectively. The nonsynonymous (replacement) nucleotide substitutions were almost five times less frequent than synonymous ones and substitutions in noncoding regions. LD decayed relatively slowly but steadily within genes. Haploblock analysis was used to define haplotype tag SNPs (htSNPs). These data will help to select SNPs for association mapping, which is already in progress.     

Canine polydactyl mutations with heterogeneous origin in the conserved intronic sequence of LMBR1

Canine preaxial polydactyly (PPD) in the hind limb is a developmental trait that restores the first digit lost during canine evolution. Using a linkage analysis, we previously demonstrated that the affected gene in a Korean breed is located on canine chromosome 16. The candidate locus was further limited to a linkage disequilibrium (LD) block of <213 kb composing the single gene, LMBR1, by LD mapping with single nucleotide polymorphisms (SNPs) for affected individuals from both Korean and Western breeds. The ZPA regulatory sequence (ZRS) in intron 5 of LMBR1 was implicated in mammalian polydactyly. An analysis of the LD haplotypes around the ZRS for various dog breeds revealed that only a subset is assigned to Western breeds. Furthermore, two distinct affected haplotypes for Asian and Western breeds were found, each containing different single-base changes in the upstream sequence (pZRS) of the ZRS. Unlike the previously characterized cases of PPD identified in the mouse and human ZRS regions, the canine mutations in pZRS lacked the ectopic expression of sonic hedgehog in the anterior limb bud, distinguishing its role in limb development from that of the ZRS.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.