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, π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 FST = 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.
QTL mapping is currently the key tool for identifying the genetic basis of quantitative traits. However, an approach that has gained in popularity recently is the use of natural populations to map traits by means of association analysis. Association analysis, or linkage disequilibrium (LD) mapping, has been used extensively to dissect traits in species where traditional QTL-mapping approaches are not feasible, such as in humans (Rafalski and Morgante 2004). Association methods have also been extended to plants and have helped increase the resolution considerably compared to standard mapping populations (Flint-Garcia et al. 2003). In principle, association studies can identify variation down to the single-nucleotide substitutions that are responsible for variation in phenotypes (quantitative trait nucleotides, QTNs). However, a move from a traditional QTL-mapping approach to association-based population surveys requires detailed knowledge about basic population genetic parameters, such as levels of genetic variation and the extent of linkage disequilibrium and population structure, and also how these parameters vary across the genome of the species in question.
The long generation times of most forest trees have led to slow progress in elucidating the genetic architecture of complex traits through traditional QTL-mapping approaches, and suggestions have been made that linkage-disequilibrium mapping approaches might be more fruitful (Brunner et al. 2004; Neale and Savolainen 2004). However, in plants, the majority of polymorphism data at the nucleotide level come from a few, well-studied model species such as Arabidopsis, maize, rice, and barley and from only a few studies that have estimated levels of nucleotide polymorphism and linkage disequilibrium in long-lived outcrossing plants, such as woody perennials (Dvornyk et al. 2002; Garcia-Gil et al. 2003; Kado et al. 2003; Semerikov and Lascoux 2003).
It is hard to know to what extent results from short-lived, selfing, or domesticated species generalize to plants with different life histories, mating systems, and/or domestication histories. For instance, even though levels of nucleotide polymorphism at the alcohol dehydrogenase (Adh) locus were comparable in selfing Arabidopsis thaliana and outcrossing A. lyrata, the distribution of segregating sites across the gene was very different in the two species (Savolainen et al. 2000). Furthermore, in predominantly selfing species like selfing A. thaliana and rice (Oryza sativa) linkage disequilibrium extends over large physical distances [>150 kb in Arabidopsis (Nordborg et al. 2000) and ∼100 kb in rice (Garris et al. 2003)] whereas in outcrossing maize, linkage disequilibrium declines to negligible levels in <1 kb (Remington et al. 2001).
The genus Populus has a long tradition as a model system for tree physiology and Populus has recently also been established as the de facto model species for tree genomics (Wullschleger et al. 2002; Brunner et al. 2004). The hope is that Populus will complement the “classical” plant model systems, such as Arabidopsis and rice, in elucidating areas such as wood formation, the transition from juvenility to maturity, and adaptation to seasonal and other environmental changes (Wullschleger et al. 2002; Brunner et al. 2004). Populus is dioecious, and hence outcrossing, and wind pollinated, and the expectation is that levels of linkage disequilibrium and population structure are generally low across the genome in Populus (Brunner et al. 2004). However, despite the long history of Populus as an experimental system for tree physiology and the recent public release of the complete genome sequence for a species of Populus (Populus trichocarpa), remarkably little is known about levels of naturally occurring genetic variation and this is especially true for polymorphisms at the nucleotide level. In this article I therefore provide the first characterization of levels of nucleotide polymorphism in a species of Populus, the European aspen (P. tremula), using data from five nuclear genes. I also characterize linkage disequilibrium and population structure by sampling from four natural populations across Europe. The results show substantial levels of nucleotide polymorphism in all five genes, comparable to levels in worldwide samples of Arabidopsis and maize. Linkage disequilibrium is generally low, extending only a few hundred base pairs. Estimates of population subdivision are low, although significant genetic differentiation among populations occurs at all five loci included in the study. Because of the high polymorphism levels and low linkage disequilibrium genome-wide linkage disequilibrium mapping approaches may not be feasible in Populus. However, very fine-scale mapping is possible if candidate gene approaches are used. The presence of significant, although modest, population structure and variation between loci in the decline of linkage disequilibrium within local populations highlights the care that has to be taken when sampling for association studies.
MATERIALS AND METHODS
Leaf material was sampled from 24 different trees of P. tremula from four different sites throughout Europe (five to seven plants per site) in June/July 2002 and July 2003. Samples were taken within a few kilometers of Besancon in eastern France (FRA), Klagenfurt in southern Austria (AUT), Färjestaden in southeastern Sweden (SWE S), and Umeå in northern Sweden (SWE N). Three to four young, undamaged leaves were collected from each tree, dried in silica gel, and stored at room temperature until DNA extraction.
Five loci were included in the study. The Alcohol dehydrogenase 1 locus (Adh1, EC 184.108.40.206) is a key enzyme for coping with low oxygen stress in higher plants. The Adh gene family has been extensively studied both at the functional and at the molecular level in a number of plant species and there are thus ample data to compare. Gibberellin 20-oxidase (GA20ox1, EC 1.14.11.-) is a multifunctional 2-oxoglutarate-dependent dioxygenase that is a key enzyme in controlling the synthesis of the plant hormone gibberellin (GAs). GAs are a group of powerful growth promoters in plants and have been implicated in such processes as seed germination, elongation growth, fruit development, and flowering (Eriksson and Moritz 2002). Glyceraldehyde-3-phosphate dehydrogenase (Gapdh, 220.127.116.11) is a tetrameric nicotinamide adenine dinucleotide-binding enzyme that plays an important role in glycolysis and glyconeogenesis. CI-1 is a cysteine protease inhibitor and TI-3 is a member of a small family of wound-inducible trypsin inhibitors that have recently been isolated in Populus (Haruta et al. 2001).
DNA extraction, PCR amplification, and sequencing:
Total genomic DNA was extracted from dried leaf tissue using the DNeasy plant mini-prep kit (QIAGEN, Valencia, CA). Primers to amplify P. tremula Adh1 and CI-1 genes were designed on the basis of BLAST searches in the publicly available Populus expressed sequence tag collection (http://poppel.fysbot.umu.se) using homologs from A. thaliana (Adh1, AF110456) and Malus x domestica (CI-1, AY173139). Primers to amplify GA20ox1 (Eriksson and Moritz 2002) and a TI-3 (Haruta et al. 2001) were designed from published sequences in GenBank (accession nos. AJ001326 and AF349443, respectively). Finally, part of the Gapdh gene was amplified in Populus using primers GPDX7F and GPDX9R from Strand et al. (1997). Primer sequences are available as supplementary material (supplementary Table 1).
PCR products were cloned into the pCR2.1 vector using a TA-cloning kit from Invitrogen (Carlsbad, CA). Clones were sequenced using BigDye chemistry (Applied Biosystems, Foster City, CA) on an ABI377 automated sequencer (Applied Biosystems) at the Umeå Plant Science Center sequencing facility. Four to 10 different clones of each fragment were sequenced to identify the presence of multiple haplotypes within individuals and to control for Taq polymerase errors. All sequences used in this article have been deposited in the EMBL/GenBank databases (accession nos. AJ842873, AJ842874, AJ842875, AJ842876, AJ842877, AJ842878, AJ842879, AJ842880, AJ842881, AJ842882, AJ842883, AJ842884, AJ842885, AJ842886, AJ842887, AJ842888, AJ842889, AJ842890, AJ842891, AJ842892, AJ842893, AJ842894, AJ842895, AJ842896, AJ842897, AJ842898, AJ842899, AJ842900, AJ842901, AJ842902, AJ842903, AJ842904, AJ842905, AJ842906, AJ842907, AJ842908, AJ842909, AJ842910, AJ842911, AJ842912, AJ842913, AJ842914, AJ842915, AJ842916, AJ842917, AJ842918, AJ842919, AJ842920, AJ842921, AJ842922, AJ842923, AJ842924, AJ842925, AJ842926, AJ842927, AJ842928, AJ842929, AJ842930, AJ842931, AJ842932, AJ842933, AJ842934, AJ842935, AJ842936, AJ842937, AJ842938, AJ842939, AJ842940, AJ842941, AJ842942, AJ842943, AJ842944, AJ842945, AJ842946, AJ842947, AJ842948, AJ842949, AJ842950, AJ842951, AJ842952 and AJ843576, AJ843713). Homologous sequences of Adh, CI-1, Gapdh, GA20-ox1, and TI-3 region were also obtained from P. trichocarpa by manually assembling contigs identified through BLAST searches in the Populus Genome Project database at Joint Genome Institute (JGI, http://genome.jgi-psf.org/).
Sequences were verified manually and contigs were assembled using the computer program Sequencher v4.0 (Gene Codes, Ann Arbor, MI). Multiple sequence alignments were made using ClustalW (Thompson et al. 1994) and adjusted manually using BioEdit (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). Estimates of nucleotide polymorphism (segregating sites, S, nucleotide diversity, π), the scaled mutation rate (θ = 4Nμ), and between-species divergences were obtained using the computer program DnaSP v4.00.5 (http://www.ub.es/dnasp/). DnaSP was also used to calculate a statistical test of neutrality, Tajima's D (Tajima 1989). The significance of D was evaluated using coalescent simulations of an island model with an average genetic differentiation equal to that estimated for each locus (see below). All simulations were run using a fixed number of segregating sites, equal to the observed number of segregating sites at each locus, and without recombination, making all tests conservative.
The decay of linkage disequilibrium with physical distance was estimated using nonlinear regression of linkage disequilibrium between polymorphic sites vs. the distance, in base pairs, between sites (Remington et al. 2001). This analysis was done both within populations and for the complete data set at each locus. Linkage disequilibrium was scored between pairs of polymorphic sites using the squared allele frequency correlations, r2 (Weir 1990). Under a simple drift-recombination model the expected value of r2 is E(r2) = 1/(1 + r), where r = 4Nc is the scaled recombination rate for the gene region. In the presence of mutations, the expectation changes to 1
(Hill and Weir 1988), where n is number of haplotypes sampled. Equation 1 was fitted using the statistical package R (http://cran.us.r-project.org/). The nonlinear regression yields a least-squares estimate of r per base pair that may be biased because of the nonindependence between different sites (Remington et al. 2001). However, it may still be useful as a representation of the rate of decay of linkage disequilibrium with physical distance. Singletons were excluded in all linkage disequilibrium analyses.
Nucleotide polymorphism and divergence:
Sequences were obtained from 17–24 individuals for the five genes studied. The length of the regions analyzed varied between ∼650 and ∼2200 bp [including insertions/deletions (indels), Table 1]. All regions contained both coding and noncoding sites (introns and/or flanking regions), although the amount of noncoding sites varied from as little as 30 bp in TI-3 to almost 1300 bp in Adh1 (Table 1).
Both indel polymorphisms and single-nucleotide polymorphisms (SNPs) were present in the sequenced regions. There were a total of nine indel polymorphisms, all in introns or flanking regions. In particular the 5′-UTR region of GA20ox1 was indel rich, with several short microsatellites that were highly polymorphic. There were no indels in the TI-3 region, which included only 30 noncoding sites. Indels were excluded in all further analyses.
SNP polymorphisms were substantially higher than indel polymorphisms; a total of 384 unique SNPs were found in the ∼6.2 kb sequenced from the five genes (Table 2) or ∼1 SNP every 60 bp. Total (πT), silent (synonymous sites and sites in introns, πsil), synonymous (πsyn), and nonsynonymous (πnonsyn) nucleotide diversities for the five genes are summarized in Table 2. Levels of polymorphism were generally high; a weighted average of silent site diversity across the five genes was 0.0160 and πsil varied between 9.4 × 10−3 at GA20ox1 and 22.9 × 10−3 at TI-3 (Table 2). Nonsynonymous diversities were generally lower than silent variation and the ratio of nonsynonymous to synonymous diversities ranged from 0.1 to 0.23 for Adh1, Gapdh, and GA20ox1, indicating strong purifying selection at most codons in these regions. CI-1 and TI-3 had about three- to fourfold higher levels of nonsynonymous diversity than the other three loci (Table 2), suggesting that the level of selective constraint varies between the genes.
There was also abundant nucleotide variation within populations at all five loci (Table 2), the only real exception being Adh1 samples from the FRA population, which had less than half the diversity of other populations, πsil = 0.0045. However, this is also the population from which the fewest samples were taken (N = 6), so this could be just a sampling artifact.
Silent site divergence, using P. trichocarpa as outgroup, was fairly constant across the five genes and averaged 0.049 (Table 2). However, the nonsynonymous substitution rate (Ka) varied >30-fold, from 0.0017 in G3PDH to 0.0463 in TI-3, suggesting that selective constraints and/or the history of adaptive evolution vary between loci. The low level of nonsynonymous divergence at Gapdh is consistent with strong selection for conserved amino acid sequence in this gene that plays a crucial role in glycolysis. There was a tendency across loci for an excess of singleton mutations, shown by negative values of Tajima's D, and this excess was significant in four out of five genes (Table 2).
The nonlinear regression model for analyzing the decay of linkage disequilibrium with distance showed that LD decays quite rapidly with distance when total samples were used; the expected value of r2 declined to <0.05 in around a few hundred bases (Figure 1). However, LD extends on average two to five times further within local populations, although there is substantial variation among populations. It should be noted that estimates of within-population rate of decay of LD are subject to much larger standard errors, due to the smaller number of sites that were polymorphic within populations. Despite the rapid decline of LD, several sites in Adh1 and GA20ox1 show extensive linkage disequilibrium over distances that approach the length of the sequenced region (Figure 1).
All five genes showed significant genetic differentiation among populations (Table 3) and using the combined data yielded an estimate of FST = 0.117, which is surprisingly high given that Populus is outbreeding, wind pollinated, and has seeds that are adapted to long-distance wind dispersal (Hamrick and Godt 1996). There was substantial variation across the five genes in levels of genetic differentiation, with the lowest estimate of FST = 0.04 for TI-3 and the highest estimate FST = 0.161 for Adh1. Most of the pairwise estimates of FST between populations were also significant, although in a few cases pairwise FST estimates were negative, probably reflecting low levels of genetic differentiation and/or stochastic variation associated with small sample sizes (supplementary Table 2 at http://www.genetics.org/supplemental/).
To test whether there was evidence for any heterogeneity in the estimates of FST across loci, I ran a coalescent simulation (with 106 replicates). The simulations assumed an island model structure with a total of 20 populations and with a migration parameter, M = 4Nm, equal to 7.57. This produces an expected FST that matches the average FST observed from the five genes (FST = 0.117). A total of 34 chromosomes were simulated, distributed among 4 of the 20 populations using the same sample configuration as in the real samples (Table 2). The simulation was standardized to a locus 1 kb in length, with θ = 4Nμ matching the weighted average across loci. The approximate 95% confidence region for FST obtained from the simulation is 0.012–0.276 and all five loci have FST values that fall with this region, suggesting that the among-locus variation in FST across loci is not unusually high in this data set.
Nucleotide polymorphism and divergence:
There are substantial levels of nucleotide polymorphisms in all five genes included in this study, both species-wide and within populations (Table 2). Levels of silent site diversity within populations averaged 0.0145, although diversity varied almost threefold between loci, from 0.0089 in GA20ox1 to 0.022 in TI-3. These values are comparable to within-population diversity in outcrossing A. lyrata (pW = 0.014, Wright et al. 2003) and to the average within-population diversity seen in Drosophila melanogaster (πW = 0.014, Moriyama and Powell 1996). It is also about two to fivefold higher than nucleotide polymorphism in the BpMADS2 gene in Betula pendula, another long-lived woody perennial, which was sampled from two locations in Finland (πW = 0.0047, Järvinen et al. 2002).
Species-wide levels of silent polymorphism in P. tremula (πsil = 0.016) are also comparable to levels of polymorphism in several plant species. In selfing A. thaliana, πsil = 0.0114, whereas in the two closely related, outcrossing species A. lyrata ssp. petrea and A. halleri species-wide diversity at silent sites averaged 0.023 and 0.015, respectively (Aguade 2001; Wright et al. 2003; Ramos-Onsins et al. 2004). In two wild species of maize, Zea diploperennis and Z. perennis, silent site diversity in four genes averaged 0.012 and 0.013 (Tiffin and Gaut 2001). It is not surprising to see such similar levels of species-wide diversities, despite differences in breeding system, life histories, and/or demography since these factors are known to have far less effect on species-wide levels of polymorphism than on levels of polymorphism within populations (Ingvarsson 2002; Charlesworth 2003).
It thus appears that P. tremula harbors substantial genetic diversity, both within populations and across the species range. This is in agreement with earlier studies of allozyme variation that have shown that Populus has a substantially higher proportion of polymorphic loci and higher gene diversity than the average among long-lived, woody perennials (Jelinski and Cheliak 1992; Liu and Furnier 1993; Hamrick and Godt 1996).
This is in stark contrast to recent studies of several species of conifers, which appear to have substantially lower genome-wide levels of polymorphism. For instance, levels of silent polymorphism range from 0.0004 to 0.0049 in P. sylvestris (Dvornyk et al. 2002; Garcia-Gil et al. 2003) and from 0.00017 to 0.0081 in Cryptomeria japonica (Kado et al. 2003). Similarily, in Pinus taeda, silent site diversity averaged 0.0064 across 63 genes involved in wood formation and pathogen and drought resistance (Neale and Savolainen 2004).
Polymorphism in P. tremula is thus ∼2- to 10-fold higher than that in Pinus and Cryptomeria. P. tremula has, together with P. sylvestris, among the widest-known geographic distributions of any tree species and the samples analyzed in this article cover only a small fraction of the entire range of P. tremula, so species-wide diversities are, if anything, underestimated in P. tremula. This suggests that the low nucleotide diversities seen in Pinus and Cryptomeria are not a general phenomenon in long-lived tree species but rather appear to be unique for conifers.
Species-wide samples show an excess of low-frequency polymorphisms that is remarkably consistent across loci (Table 2). Selection seems like an unlikely cause of negative values of these statistics because the excess of low-frequency polymorphisms is consistent across loci. Also, sequences were obtained through cloning of PCR products and although care was taken to verify all singleton mutations in multiple clones, it is possible that some of the remaining singletons are the result of Taq polymerase errors. However, given the error rate of Taq of 1 × 10−4–2 × 10−5 (Cline et al. 1996), PCR errors would account for only between one and five singletons at each locus and thus do not appear to explain the significant excess of singleton mutations observed. Therefore, it is more likely that the departures of the frequency spectra from the standard neutral model are a product of demographic processes, such as postglacial expansion (Petit et al. 2003) or effects of population subdivision. Although population structure is expected to shift the frequency spectrum toward an excess of sites at intermediate frequencies (i.e., positive values of DTajima) under some scenarios, it is possible that population subdivision can also lead to negative Tajima's D, depending on strength of population subdivision and the actual sample configuration. Thus the reason for the overall excess of low-frequency mutation in Populus must be investigated in more detail.
Silent site divergence from P. trichocarpa averaged 6.1% and, using a silent substitution rate in the range 5.0–8.0 × 10−9 (Wolfe et al. 1987; Gaut 1998), yields an estimated time of divergence between P. tremula and P. trichocarpa between 3.8 and 6.2 million years. P. tremula and P. trichocarpa belong to different sections within the genus Populus (Leuce, aspens, and Tacamahaca poplars, respectively) and are morphologically quite distinct. It thus appears that the radiation that produced the ∼300 species contained in the genus Populus has been rapid and relatively recent. This view is also supported by other molecular data (Leskinen and Alström-Rapaport 1999) and by the prevalence of both intra- and intersectional hybridization in Populus (Brunner et al. 2004).
LD declines rapidly with distance in all five genes. In general, average linkage disequilibrium declined to negligible levels (r2 < 0.05) in <500 bp, although in some cases LD extends much further within local populations (up to 1 kb or more). Levels of LD observed in Populus are on par with those observed in maize, where LD also declines rapidly. Remington et al. (2001) showed that LD declined to negligible levels in ∼1 kb in 5 of 6 genes studied in a set of 102 inbred lines of maize. The sixth gene studied, su1, is part of the starch production pathway in maize and has been a target of selection during the domestication of maize (Whitt et al. 2002) and this artificial selection has presumably contributed to the extensive LD seen in su1. A study of a more limited set of elite maize lines showed that LD can extend over great distances also in maize, up to 100 kb, and that there was virtually no decline in r2 over distances of ∼500 bp in the collection of 18 genes they surveyed, consistent with a strong effect of breeding-induced bottlenecks and selection in producing the elite germplasm pool (Ching et al. 2002).
These results from Populus and maize, both outcrossers, are in stark contrast to the patterns of LD seen in primarily self-fertilizing species, such as A. thaliana and rice (Oryza spp.). Selfing dramatically reduced the effective recombination (Nordborg 2000), and in Arabidopsis LD generally extends over distances up to 250 kb (Nordborg et al. 2000) while in rice LD extends ∼100 kb (Garris et al. 2003). However, in wild barley (Hordeum vulgare ssp. spontaneum), another highly selfing species, some loci show high levels of LD, as expected in a selfing species, where as other loci show a surprisingly rapid decline of LD with distance (Lin et al. 2002). Although the reasons behind the different patterns of LD among loci in barley and also differences in LD patterns between Arabidopsis and barley are not fully understood, it is clear that mating system alone does not explain differences in LD and that other factors, such as the demographic history of a species, are important as well (Lin et al. 2002).
All five genes showed significant genetic differentiation among populations, although estimates of FST were moderate (Table 3). The estimates of FST are still surprisingly high given that Populus is outbreeding, wind pollinated, and has seeds that are adapted to long-distance wind dispersal (Hamrick and Godt 1996). There was a greater than fourfold variation among loci in estimates of FST, but a coalescent simulation suggested that this dispersion is no greater than expected. A recent study of genetic differentiation in maternally inherited cpDNA markers also showed moderate levels of genetic differentiation across Europe in P. tremula (FST = 0.11, Petit et al. 2003). Interestingly, the cpDNA estimate of FST is close to the average of the five nuclear genes studied here. This is somewhat surprising, since maternally inherited cpDNA markers are haploid and their lower effective population size in a dioecious species, like P. tremula, should result in an increased genetic differentiation, compared to nuclear genes [Ne(cpDNA)/Ne(nuclear) ≈ 0.25 with equal sex ratios, Laporte and Charlesworth 2002]. Sex ratios of Scandinavian P. tremula populations are also predominantly male biased (Blumenthal 1942), reinforcing the differences in effective population sizes between nuclear and organelle loci. More loci are needed to establish whether this discrepancy in genetic differentiation between organelle and nuclear markers is a real phenomenon or just a sampling artifact.
This is the first study to quantify levels of nucleotide polymorphism and linkage disequilibrium from a member of the genus Populus. The results presented here show that P. tremula harbors significant amounts of genetic variation both within populations and across the species range. Since estimates of population subdivision are low (albeit significant), a substantial fraction of the species-wide diversity can be sampled from single populations. In the species-wide sample LD declines rapidly with physical distance, although it may extend much further within local populations.
The rate of decay of LD has implications for the potential utility of association approaches in mapping QTL as it effectively determines whether genome-wide scans are feasible or whether a candidate gene approach has to be taken (Gaut and Long 2003). The rapid decay of LD in P. tremula suggests that the number of markers needed to ensure adequate genome coverage may be simply too high for genome-wide scans to be feasible in this species and a more realistic approach might therefore be to search for functional variation in putative candidate loci. On the other hand, the rapid decay of LD suggests that it may be possible to map functional variation to very fine scales in P. tremula.
Population structure is low, but significant for all loci investigated. This structure creates additional linkage disequilibrium between loci that have different allelic frequencies in different subpopulations. Population structure may therefore interfere with association studies by producing spurious associations. However, methods to properly account for population structure have been developed (Pritchard et al. 2000) and applied successfully in other plant species (Thornsberry et al. 2001). Nevertheless, the presence of significant population structure and variation in rate of decay of LD in different populations emphasizes the care that has to be taken when sampling for association studies.
I thank Carin Olofsson for assistance in obtaining the DNA sequences presented in this article. Outi Savolainen and two anonymous reviewers provided comments that notably improved the article. Funding from the Swedish Research Council (Vetenskapsrådet) is also acknowledged.
Communicating editor: O. Savolainen
- Received August 12, 2004.
- Accepted October 25, 2004.
- Genetics Society of America