| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Genetics, Vol. 171, 1247-1256, November 2005, Copyright © 2005
doi:10.1534/genetics.105.041566
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
* Institute for Genomic Diversity, Cornell University, Ithaca, New York 14853 and Plant Genome Mapping Laboratory, University of Georgia, Athens, Georgia 30602
1 Corresponding author: Institute for Genomic Diversity, 157 Biotechnology Bldg., Cornell University, Ithaca, NY 14853.
E-mail: sk20{at}cornell.edu
 |
ABSTRACT |
|---|
 TOP ABSTRACT MATERIALS AND METHODSRESULTSDISCUSSIONACKNOWLEDGEMENTSLITERATURE CITED |
|---|
, however, based on empirical estimates of rates of recombination, mutation, and self-pollination, show that LD is more extensive than expected under a neutral equilibrium model. The disparity between and the population mutation parameter, 
, is similar to that observed in other species whose population history appears to be complex. From a practical standpoint, these results suggest that S. bicolor is well suited for association studies using reasonable numbers of markers, since LD typically extends at least several kilobases but has largely decayed by 15 kb.
, where Ne (effective population size) is a function of the long-term historical size of the population, population structure, and mating system. Similarly, the mutational process that generates associations is summarized in the population mutation parameter 4Neµ, or , where µ is the mutation rate. At equilibrium in a randomly mating population without selection, the extent of allelic associations is simply a function of the relative rates of mutation and recombination, since 4Ne cancels out in the ratio /. Nevertheless, even in this simplest of scenarios, LD decay varies widely in unlinked regions due to the substantial stochasticity of the evolutionary process (NORDBORG and TAVARE 2002).
When a population departs from the panmictic equilibrium model, we can no longer accurately estimate and from empirical data, and observed levels of LD may be inconsistent with empirically determined rates of recombination under simple demographic models. In the case of a partially self-pollinating species, the resulting reduction in heterozygous genotypes means that the effective rate of recombination is lower for a given rate of crossing over, so LD is, on average, expected to be more extensive. This is an equilibrium effect and can be accounted for by appropriate scaling (NORDBORG and DONNELLY 1997; NORDBORG 2000). In other cases, however, mating is random but nonequilibrium population history leads to inconsistencies that are less easily explained. For example, LD in Drosophila is more extensive than expected on the basis of genetic distances and estimates of Ne from (ANDOLFATTO and PRZEWORSKI 2000). In humans, LD can be either more or less extensive than expected, depending on the range of interlocus comparisons (PRITCHARD and PRZEWORSKI 2001). Selection can also generate LD, although the locus-specific effects of selection can be hard to distinguish from the noise generated by neutral processes (HUTTLEY et al. 1999; KIM and NIELSEN 2004).
Our knowledge of the extent of LD in plants is limited (FLINT-GARCIA et al. 2003). Only in Arabidopsis, with a fully sequenced genome and high-density genetic map, has it been possible to conduct analyses comparable to those in Drosophila and humans, namely to evaluate LD over large, defined physical distances in the context of local rates of recombination. Very extensive LD on the order of 250 kb has been observed in this highly self-pollinating species (HAGENBLAD and NORDBORG 2002; NORDBORG et al. 2002), although a recent genome-wide study has shown that LD at most loci decays within 25–50 kb (NORDBORG et al. 2005). In contrast, LD in maize, an outcrosser, decays within a few hundred base pairs in diverse samples (TENAILLON et al. 2001), although the extent of LD increases when narrower samples of germplasm (REMINGTON et al. 2001; CHING et al. 2002; JUNG et al. 2004) or targets of selection (CLARK et al. 2004; PALAISA et al. 2004) are analyzed. Studies of gene-sized regions in both Populus (INGVARSSON 2004) and loblolly pine (BROWN et al. 2004) indicate low levels of LD in these highly outcrossing trees. While these contrasts are frequently explained in terms of mating system, LD in wild barley, with a selfing rate similar to Arabidopsis, has patterns of LD more similar to maize (LIN et al. 2002). Very extensive LD is observed in rice (GARRIS et al. 2003; SEMON et al. 2004) where the effects of mating system are likely confounded with population structure.
The extent of LD is a key issue in the design and feasibility of association mapping methods (FLINT-GARCIA et al. 2003). Association studies (also called LD mapping) have been successful in maize (e.g.,THORNSBERRY et al. 2001), where they may be limited to candidate genes because of the small extent of LD. Sorghum bicolor, a largely (
70%) self-pollinating domesticated grass (ROONEY and SMITH 2000) that is important for human nutrition in semiarid regions of sub-Saharan Africa (FAO 1996), is closely related to maize but has a smaller and less complex genome (DRAYE et al. 2001). In a previous study (HAMBLIN et al. 2004), we reported that LD over very short distances in sorghum was more extensive than in maize, suggesting that sorghum may be suitable for LD mapping of genes underlying complex, agronomically important traits common to both species. However, our limited data did not reveal anything about the scale over which LD dissipates, and we had no information about local relationships between genetic and physical distance for the regions analyzed. A physical map of the sorghum genome is being assembled and integrated with the genetic map (DRAYE et al. 2001; MULLET et al. 2002; BOWERS et al. 2003), which ultimately will allow for estimation of rates of recombination, as well as LD, over fairly large distances, on the scale of centimorgans and megabases.
Meanwhile, genomic regions represented by BAC clones containing genetic markers can be used to sample patterns of LD and the relationship between physical and genetic distance on a scale of tens of kilobases. Six fully sequenced BAC clones, representing five different chromosomes, were used for this purpose. The goals of this study were to examine the pattern of pairwise associations among a large number of single-nucleotide polymorphisms (SNPs) in six large (40–100 kb) unlinked regions and to estimate (the population recombination parameter, 4Ner), (the population mutation parameter, 4Neµ), and r (the rate of recombination per base pair per generation) for those same regions. Analysis of this data set allowed us to assess the contribution of mating system and recombination rate to patterns of LD in sorghum and provides a general picture of those patterns that may prove useful in LD-based methods of genetic analysis.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONACKNOWLEDGEMENTSLITERATURE CITED |
|---|
(see below). DNA was prepared from young leaves of individual plants according to the method of DOYLE and DOYLE (1987).
Choice of regions for resequencing:
Annotation of predicted genes was used to identify putative intronic regions between 500 and 1700 bp in length that could be amplified from PCR primers on the basis of flanking exon sequence. Note that this is a technical consideration only and that the functional status of the sequence has no bearing on the analysis.
PCR products ranged from 700 to 1700 bp in size (Table 1) and were prepared for sequencing by treatment with shrimp alkaline phosphatase and exonuclease I digestion. PCR primers (and internal primers as necessary) were used for cycle sequencing with ABI Big Dye V. 3.1, and sequencing reactions were analyzed on a 3730 or 3700 sequencer at the Bioresources Center at Cornell University.
|
Files of variable sites generated by MEGA (KUMAR et al. 2001) were formatted for LD estimation after removal of sites for which the minor allele had a frequency of <10% (REMINGTON et al. 2001). In a few regions we observed an individual that appeared to result from introgression from a divergent wild relative, as it differed from all other alleles in the sample at many sites. For example, in region 2, of a total of 91 segregating sites, accession LWA22 contributed 51. Such individuals were removed from the sample for those regions only. Note that, because these individuals contribute only singletons, they have no effect on estimates of LD.
In a small number of cases, some accessions produced sequence of insufficient quality for calling all bases, but most of the polymorphic SNPs could nonetheless be reliably scored and were added to the LD analysis of the full sequence data. These sequences were not included in estimates of sequence diversity. These exceptions to the strategy of full resequencing are noted in Table 1. The programs dipdat (kindly provided by Dick Hudson) and maxdip (http://genapps.uchicago.edu/maxdip/index.html) were used to estimate r2 (HILL 1974) and (i.e., 4Ner), respectively, from genotypic (i.e., unphased diploid) data (HUDSON 2001). LD triangle plots were constructed using TASSEL (http://www.maizegenetics.net/bioinformatics/tasselindex.htm) after removal of genotypes containing more than one heterozygous site, since TASSEL requires that phase be known.
We also used the program PHASE (LI and STEPHENS 2003) to obtain estimates of and of the variation in recombination rates across the surveyed regions. Default parameters were used, except that the –X10 option was used to increase the number of iterations in the final run, as suggested in the documentation. Point estimates of and 
(the factor by which the recombination rate in an interval between two loci exceeds the background rate) were obtained by taking the median value from 1000 iterations, as suggested in the documentation. An interval was considered a potential hotspot if the fifth percentile of was >1.0.
Estimation of rates of recombination:
We used a recombinant inbred line (RIL) population, DNAs of which were kindly provided by Tom Hash of the International Crop Research Institute for the Semiarid Tropics. There were 244 lines in the population, which had been self-pollinated and advanced to the F6 generation by single-seed descent. Within each BAC sequence, we identified at least two simple sequence repeats that showed length variation between the parents of the RIL population (BTx623 and IS18551). We designed primers for fragment analysis and scored the markers in the 244 DNAs. Recombination per generation (r) was estimated using the formula: observed recombination fraction = 2r/(1 + 2r) (BURR et al. 1988) and divided by the number of base pairs between markers.
|
RESULTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONACKNOWLEDGEMENTSLITERATURE CITED |
|---|
700–1700 bp spaced irregularly across each region (Figure 1, Table 1), spanning a total distance of between 38 and 103 kb/BAC. Our panel of 32 S. bicolor accessions (see MATERIALS AND METHODS), which includes both cultivated and wild representatives, was chosen to maximize diversity and to capture as much of the evolutionary history of S. bicolor as possible. All 32 lines were fully sequenced for most regions (for exceptions, see MATERIALS AND METHODS). The numbers of SNPs observed in each subregion (i.e., amplicon) vary considerably (Table 1). Only SNPs for which the minor allele was present three or more times in the sample were included in the analyses of LD; all subsequent references to SNPs refer to this subset of 249 of 427 total SNPs observed.
|
|
There is considerable interest in whether LD is "block-like," as strong haplotype structure, if real, may simplify LD-mapping studies (ZHANG et al. 2002). The plots in Figure 2 obscure patterns of LD among blocks of sites, so we have also made triangle plots of the data (Figure S3 at http://www.genetics.org/supplemental/). Again, the patterns are very different among the regions. In region 4, for example, there are only a few isolated associations of any significance (< 0.01) between sites >1 kb apart. Region 5 has associations of P < 0.0001 that extend >15 kb, but those associations are patchy, not block like. In contrast, region 1 has an almost solid block of associations of P < 0.001 extending >13 kb. These patterns are not easy to compare because the spacing of subregions and the numbers of SNPs observed in each are not uniform; nonetheless it is apparent that some regions of the sorghum genome have extended haplotypes while others do not.
Multilocus estimates of LD:
As is evident in Figure 2, pairwise estimates of LD are noisy and often difficult to interpret (NORDBORG and TAVARE 2002). A statistic that summarizes LD over an entire region is , the population recombination parameter (see Introduction). A number of methods have been proposed for estimating , but many have been found to perform poorly, particularly when sample sizes are small (WALL 2000). Due to the nature of our data set (namely, fairly large, discontinuous regions), we chose to use the composite likelihood (CL) method of HUDSON (2001) as well as the product of approximate conditionals (PAC) method of LI and STEPHENS (2003). Consistent with simulation studies that showed that these estimators perform well, the two estimates for each region do not differ by more than about twofold (Table 2). PAC was lower than CL in five of six cases, but the relative order of the values, by region, was similar for the two methods.
|
; however, the impact on this data set is modest. Assuming l = 300 or 500, a model without gene conversion (i.e., f = 0) fit the data best for all regions except regions 2 and 5 (Table 2). For region 2, the likelihood curve was extremely shallow, so that a model without gene conversion fit the data almost as well as several other models that included high rates of gene conversion. Support for a model including gene conversion was strong only for region 5 but only a modest amount of gene conversion (f = 1) was inferred; CL for this region was 7.55 under a model with f = 0.
For both methods, the estimates of (Table 2) vary 12-fold across the six genomic regions, and the 90% C.I.'s for regions 3 and 4 do not overlap. Variation in recombination rate among regions of the genome is a possible factor underlying these differences. To test this hypothesis, we experimentally measured rates of recombination for five of our regions. (Recombination was not measured for region 1 because the entire BAC clone is only 40 kb in size and recombination events were not likely to be observed.) These studies, shown in Table 3, indicate that rates of recombination vary 4-fold from region to region and may explain a portion of the difference in patterns of LD that we observe, in particular the difference between regions 3 and 4. None of the estimates of r, however, are significantly different. Variation in r explains only 25% of the variation in CL, but explains almost 60% of the variation in PAC.
|
, allows rates of recombination to vary along the region and reports a value of for each interval between adjacent SNPs (see MATERIALS AND METHODS). The vast majority (98%) of -values were between 0.5 and 2.0, indicating a fairly uniform rate of recombination. Two exceptions, noted in Table 2, were found. In both cases, point estimates of were 6, and 95% of the iterations produced a -value >1.0. While this does not constitute a test of significance, and does not account for multiple tests (252 intervals were tested), these two intervals are clearly outliers in the data set.
Cultivated sorghum has experienced a domestication bottleneck (ALDRICH and DOEBLEY 1992), which is expected to affect patterns of LD (PRITCHARD and PRZEWORSKI 2001). Furthermore, admixture of accessions from the wild and cultivated populations could result in elevated LD (PRITCHARD and PRZEWORSKI 2001). However, estimates of for the cultivated accessions only are almost all (five of six) lower than that for the total sample (Table S5 at http://www.genetics.org/supplemental/), suggesting that the effect of the bottleneck is more important than that of admixture. This is not surprising, given that the wild accessions are only moderately differentiated from the cultivated ones and that there is little structure among the cultivated accessions (CASA et al 2005). The exception is region 1, for which is considerably higher in the cultivated sample. This is because the cultivated sample had 47 SNPs as compared to 78 in the total sample, and a large number of alleles in strong LD were present only in the wild accessions (most of the SNPs in subregions 1A–1D). When wild accessions were eliminated, the number of SNPs in region 5 dropped from 70 to 62, and for the four other regions there was no difference. In general, most of the polymorphisms that were observed only in the wild accessions were in low frequency and had not been included in the analysis. The greater LD in the cultivated sample is therefore presumably due to the loss of haplotypes that provide evidence of recombination among the same pairs of alleles.
Because the ability to detect evidence of recombination depends on the presence of informative polymorphisms, it is useful to look at the ratio of to in comparing across regions. Table 4 shows the ratio of and for the total sample and the cultivated accessions only. (In this analysis, we use the higher estimate of , usually CL, which is conservative for testing the hypothesis that LD is higher than expected.) This ratio, which ranges from 0.040 to 0.375 in the total sample and from 0.014 to 0.249 in the cultivated sample, indicates that recombination is relatively infrequent relative to mutation.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONACKNOWLEDGEMENTSLITERATURE CITED |
|---|
, based on six large, unlinked, fully resequenced regions for which we have also estimated the local rate of recombination per base pair. We find that the extent of allelic associations in sorghum, as assessed by pairwise measures of LD, is higher than in maize but lower than in rice and Arabidopsis, in qualitative agreement with expectations based on differences in mating system. Multilocus estimates of the population recombination parameter , however, are among the lowest observed in any species, including Arabidopsis (KUITTINEN and AGUADE 2000; HAGENBLAD and NORDBORG 2002; NORDBORG et al. 2005). In attempting to account for these observations, several factors should be considered.
Estimation of :
A number of different methods have been proposed for estimating , many of which do not perform well when tested against simulated data with known values of (WALL 2000; HUDSON 2001): point estimates are often far from the known value, and/or confidence intervals are very large. Likewise, different estimators of may produce very different results for the same empirical data (e.g., TENAILLON et al. 2002). Therefore it is reasonable to ask how much confidence we have in the estimates that we report. True confidence intervals for estimates of are not trivial to obtain, particularly when the data collection scheme is not simple (as in our case). For PAC, however, it is easy to obtain an approximation of the uncertainty from the distribution of sampled -values; each of these 90% credible intervals (Table 2) contains the corresponding CL. In fact, the two estimates are in all cases within about twofold of each other, suggesting that they are not very far from the true value.
Testing an equilibrium model:
Both mating system and rates of recombination affect LD in predictable ways that can be accounted for in an equilibrium model. If accounting for these factors does not explain the data, then we must invoke nonequilibrium phenomena such as population structure and history and selection.
The effects of mating system, recombination rate, and mutation rate on in a partially selfing organism can be described by the equation / = (r/µ)(1 – F), where F is the inbreeding coefficient (HAGENBLAD and NORDBORG 2002). The rate of self-pollination in sorghum is 70% (ROONEY and SMITH 2000), so F = 0.7/(2 – 0.7) = 0.54. On the basis of synonymous substitutions between maize and sorghum at multiple loci (SWIGONOVA et al. 2004), we use an estimate of µ = 1 x 10–8/bp/generation. We therefore expect that r/µ should be on the order of (r /1 x 10–8) x (0.46), a value that is >1 for all the regions in our study. Actual values of / are 5–33 times lower (Table 4), using the higher of our two estimates of . While mutation rates may vary among regions, as reflected by the differences in , it appears that the effects of mating system and recombination rate cannot account for the low values of observed in this sample. Much higher rates of self-pollination (>90%) would be required to fit these data with an equilibrium model. While some cultivated sorghum accessions do have rates in this range, sorghum spent the vast majority of its evolutionary history as a wild species with an outcrossing rate believed to be 
30% (DOGGETT 1988). Rates of recombination much lower (10–9/bp/generation), or mutation much higher (2 x 10–7/bp/generation), could also resolve the discrepancy, but such values are not plausible.
Population structure and history:
Population structure can contribute to elevated LD (PRITCHARD and PRZEWORSKI 2001), so our wide sampling could bias estimates of LD upward if population structure were strong. However, we know from studies of SSR diversity that there is little structure in sorghum populations (CASA et al. 2005). Furthermore, results of WAKELEY and LESSARD (2003) suggest that our sampling strategy may minimize the impact of any population structure on estimates of LD. In samples drawn from just two demes, they show that correlations in histories of alleles are increased. The properties of "scattered" samples like ours, however, where each individual comes from a different deme, approach those of samples from panmictic populations. In any case, estimates of in the cultivated sample are lower than those in the total (i.e., admixed) sample, including wild accessions. These results are consistent with the findings for maize, where samples that capture greater genetic diversity (i.e., more ancestral recombination events) show increasingly less LD (FLINT-GARCIA et al. 2003). Thus our sampling strategy could be considered conservative for testing the hypothesis that LD in sorghum is more extensive than expected.
Nonequilibrium population history can result in discrepancies between and , as has been observed in humans and Drosophila (FRISSE et al. 2001; WALL et al. 2002), and this factor is likely to be important in sorghum, which has experienced a domestication bottleneck. In our data, the cultivated subsample has 58% of the segregating sites of the total sample (Table S6 at http://www.genetics.org/supplemental/). Consistent with the effects of a bottleneck, the frequency spectrum of variation in the cultivated data set shows a strong departure from the neutral expectation: coalescent simulations show that both the mean (0.45) and the variance of Tajima's D (1.42) are significantly too large. We have explored a small range of recent simple bottleneck models to attempt to find one that is consistent with the domestication history of sorghum and with our empirical data. Using = 0.0056 based on variation in wild S. bicolor (unpublished data from our lab), and µ = 1 x 10–8/bp/generation, we estimate ancestral Ne to be 1.4 x 105. A domestication event 5000–6000 years ago (KIMBER 2000) would thus correspond to 0.01(4Ne) generations. Using coalescent simulations, we were unable to find bottleneck parameters for which both the average D and the variance of D were close to what we observe. We are currently performing more exhaustive analyses with a larger data set that will be published elsewhere. Nonetheless, while the details remain to be determined, it seems reasonable to conclude that a nonequilibrium history has perturbed the frequency spectrum and also elevated linkage disequilibrium in this sample.
The weak relationship between association and distance for alleles that are <4 kb apart (Figure 2) is an interesting observation that may also be due to population history. If a bottleneck generated excess LD, that excess LD will have decayed over time more quickly for alleles that are farther apart (i.e., for which r is larger), while closely linked alleles may still retain the signature of that demographic event. A similar pattern is observed at some loci in Arabidopsis; e.g., see Figure 4 in SHEPARD and PURUGGANAN (2003).
Selection:
While the data may be consistent with a strictly neutral, nonequilibrium model, this of course does not preclude that selection may also have influenced the observed patterns of LD. In particular, we might expect that directional selection associated with domestication has played a role in the cultivated subsample. A multilocus HKA test of polymorphism and divergence (HUDSON et al. 1987) in that subsample showed that the data were highly unlikely under a neutral model (P = 0.00001), but there was no convincing evidence of selection at any particular loci. It is possible that the departure is due to demography rather than to selection. Interestingly, subregions A–D of region 1 in the cultivated sample have only 10–15% of the diversity of the total sample, possibly due to selection at the shrunken2 locus near subregion 1A, but this has not resulted in higher estimates of LD for this region in the cultivated sample (see RESULTS). Region 3, with the highest LD, contains the phytochrome A locus, for which a previous study showed no evidence of selection in sorghum (WHITE et al. 2004). Thus, while we cannot rule out that selection may have played a role in shaping observed patterns of LD, there is no clear relationship in the data between the extent of LD and any evidence of a history of selection.
Features of recombination in sorghum:
In maize, where LD is much less extensive, it appears that most recombination occurs within genes rather than in intergenic regions, perhaps because the substantial variation in transposable element complements from chromosome to chromosome disrupts the homology necessary for recombination to occur (FU et al. 2002). Our sequencing strategy was not designed to address this question, and most evidence for recombination occurs between the sequenced subregions in genomic areas that consist of both coding and noncoding sequence. There is evidence of only 12 recombination events [using the four-gamete test of HUDSON and KAPLAN (1985)] within the 27 fully resequenced subregions, almost all of which correspond to introns. However, the small number of bases surveyed represents only 7% of the total length of the genomic regions analyzed (24 of 330 kb). Interestingly, the most gene-rich region, region 3, has the highest LD, while regions 4 and 5, relatively sparse in genes, have the least. Furthermore, rates of recombination, as assessed by the PAC method of LI and STEPHENS (2003), appeared to be quite uniform across each region. These observations, although anecdotal, do not suggest that recombination in sorghum is concentrated primarily within genes.
HAUBOLD et al. (2002) and NORDBORG et al. (2005) concluded that, in Arabidopsis, most "recombination" is in fact caused by gene conversion. High rates of gene conversion have also been implicated in the discrepancy between short- and long-range LD in humans (FRISSE et al. 2001; PRZEWORSKI and WALL 2001; ANDOLFATTO and WALL 2003) and a deficit of LD in short regions of low recombination in Drosophila (ANDOLFATTO and WALL 2003). In sorghum, if there is any discrepancy between short- and long-range LD, it goes in the other direction: short-range LD is more extensive than expected, relative to long-range LD. Consistent with this observation, we saw little evidence for gene conversion, and there is little evidence of recombination within subregions, the scale on which short gene conversion tracts would be observed.
To understand and make use of information about LD, we made empirical estimates of r per base pair. Rates of recombination within and between species vary tremendously, such that a genetic interval defined by two markers 1 cM apart may correspond to 50 kb or 5 Mb of DNA, depending on the organism and the local rate of recombination. In wheat, for example, rates as low as 0.04 cM/Mb and as high as 8.5 cM/Mb have been measured (GILL et al. 1996a,b). In maize, measurement of recombination nodules per cytological distance vary >12-fold across chromosome 1 (TENAILLON et al. 2002). In sorghum, our estimates of r varied only 4-fold in five unlinked regions, in the range of 2–8 cM/Mb, and the differences were not significant. All measured rates were higher than the estimate of 1.5 cM/Mb on the basis of total genome size and genetic map distance. This suggests that other genomic regions, not sampled in this study, may be recombinationally relatively inert.
Conclusions:
The extent of LD in sorghum is greater than that in maize, where it is generally low (TENAILLON et al. 2001), and less than that in Arabidopsis (NORDBORG et al. 2005). This qualitative observation is consistent with the mixed mating system of sorghum producing intermediate levels of effective recombination. Quantitative analyses, however, based on estimated rates of recombination, mutation, and selfing, show that both and the ratio of to are lower than expected under an equilibrium model. These analyses, as well as the frequency spectrum of polymorphism, suggest that a genome-wide departure from equilibrium underlies this phenomenon.
The greater extent of LD in sorghum makes it amenable for association studies using a limited number of markers. Genotyping of a few SNPs per gene in many cases can capture most of the haplotypic variation.
|
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONACKNOWLEDGEMENTSLITERATURE CITED |
|---|
Note added in proof: J.-S. KIM, M. N. ISLAM-FARIDI, P. E. KLEIN, D. M. STELLY, H. J. PRICE, R. R. KLEIN and J. E. MULLET (2005, Comprehensive molecular cytogenetic analysis of sorghum genome architecture: distribution of euchromatin, heterochromatin, genes and recombination in comparison to rice. Genetics 171 (in press)) have recently estimated the average genome-wide rate of recombination for euchromatic regions in S. bicolor. Their estimate is very similar to ours: 0.254 Mbp/cM or 4 x 10–8/bp.
|
FOOTNOTES |
|---|
|
LITERATURE CITED |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONACKNOWLEDGEMENTSLITERATURE CITED |
|---|
ALDRICH, P. R., and J. DOEBLEY, 1992 Restriction fragment variation in the nuclear and chloroplast genomes of cultivated and wild Sorghum bicolor. Theor. Appl. Genet. 85: 293–302.
ANDOLFATTO, P., and M. PRZEWORSKI, 2000 A genome-wide departure from the standard neutral model in natural populations of Drosophila. Genetics 156: 257–268.
ANDOLFATTO, P., and J. D. WALL, 2003 Linkage disequilibrium patterns across a recombination gradient in African Drosophila melanogaster. Genetics 165: 1289–1305.
BOWERS, J. E., C. ABBEY, S. ANDERSON, C. CHANG, X. DRAYE et al., 2003 A high-density genetic recombination map of sequence-tagged sites for Sorghum, as a framework for comparative structural and evolutionary genomics of tropical grains and grasses. Genetics 165: 367–386.
BROWN, G. R., G. P. GILL, R. J. KUNTZ, C. H. LANGLEY and D. B. NEALE, 2004 Nucleotide diversity and linkage disequilibrium in loblolly pine. Proc. Natl. Acad. Sci. USA 101: 15255–15260.
BURR, B., F. A. BURR, K. H. THOMPSON, M. C. ALBERTSON and C. W. STUBER, 1988 Gene mapping with recombinant inbreds in maize. Genetics 118: 519–526.
CASA, A. M., S. E. MITCHELL, M. T. HAMBLIN, H. SUN, J. E. BOWERS et al., 2005 Diversity and selection in sorghum: simultaneous analyses using simple sequence repeats. Theor. Appl. Genet. 111: 23–30.[CrossRef][Medline]
CHING, A. S., K. S. CALDWELL, M. T. JUNG, M. DOLAN, O. H. SMITH et al., 2002 SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines. BMC Genet. 3: 19.[CrossRef][Medline]
CLARK, R. M., E. LINTON, J. MESSING and J. F. DOEBLEY, 2004 Pattern of diversity in the genomic region near the maize domestication gene tb1. Proc. Natl. Acad. Sci. USA 101: 700–707.
DOGGETT, H., 1988 Sorghum. Longman Scientific & Technical, Essex, UK.
DOYLE, J. J., and J. L. DOYLE, 1987 A rapid DNA isolation procedure for small amounts of leaf tissue. Phytochem. Bull. 19: 11–15.
DRAYE, X., Y. R. LIN, X. Y. QIAN, J. E. BOWERS, G. B. BUROW et al., 2001 Toward integration of comparative genetic, physical, diversity, and cytomolecular maps for grasses and grains, using the sorghum genome as a foundation. Plant Physiol. 125: 1325–1341.
FAO, 1996 The World Sorghum and Millet Economies: Facts, Trends and Outlook. Food and Agricultural Organization of the United Nations, Rome.
FEARNHEAD, P., and P. DONNELLY, 2001 Estimating recombination rates from population genetic data. Genetics 159: 1299–1318.
FLINT-GARCIA, S. A., J. M. THORNSBERRY and E. S. BUCKLER, IV, 2003 Structure of linkage disequilibrium in plants. Annu. Rev. Plant Biol. 54: 357–374.[CrossRef][Medline]
FRISSE, L., R. R. HUDSON, A. BARTOSZEWICZ, J. D. WALL, J. DONFACK et al., 2001 Gene conversion and different population histories may explain the contrast between polymorphism and linkage disequilibrium levels. Am. J. Hum. Genet. 69: 831–843.[CrossRef][Medline]
FU, H., Z. ZHENG and H. K. DOONER, 2002 Recombination rates between adjacent genic and retrotransposon regions in maize vary by 2 orders of magnitude. Proc. Natl. Acad. Sci. USA 99: 1082–1087.
GARRIS, A. J., S. R. MCCOUCH and S. KRESOVICH, 2003 Population structure and its effect on haplotype diversity and linkage disequilibrium surrounding the xa5 locus of rice (Oryza sativa L.). Genetics 165: 759–769.
GILL, K. S., B. S. GILL, T. R. ENDO and E. V. BOYKO, 1996a Identification and high-density mapping of gene-rich regions in chromosome group 5 of wheat. Genetics 143: 1001–1012.[Abstract]
GILL, K. S., B. S. GILL, T. R. ENDO and T. TAYLOR, 1996b Identification and high-density mapping of gene-rich regions in chromosome group 1 of wheat. Genetics 144: 1883–1891.[Abstract]
HAGENBLAD, J., and M. NORDBORG, 2002 Sequence variation and haplotype structure surrounding the flowering time locus FRI in Arabidopsis thaliana. Genetics 161: 289–298.
HAMBLIN, M. T., S. E. MITCHELL, G. M. WHITE, J. GALLEGO, R. KUKATLA et al., 2004 Comparative population genetics of the panicoid grasses: sequence polymorphism, linkage disequilibrium and selection in a diverse sample of sorghum bicolor. Genetics 167: 471–483.
HAUBOLD, B., J. KROYMANN, A. RATZKA, T. MITCHELL-OLDS and T. WIEHE, 2002 Recombination and gene conversion in a 170-kb genomic region of Arabidopsis thaliana. Genetics 161: 1269–1278.
HILL, W. G., 1974 Estimation of linkage disequilibrium in randomly mating populations. Heredity 33: 229–239.[Medline]
HUDSON, R. R., 2001 Two-locus sampling distributions and their application. Genetics 159: 1805–1817.
HUDSON, R. R., 2002 Generating samples under a Wright-Fisher neutral model of genetic variation. Bioinformatics 18: 337–338.
HUDSON, R. R., and N. L. KAPLAN, 1985 Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics 111: 147–164.
HUDSON, R. R., M. KREITMAN and M. AGUADE, 1987 A test of neutral molecular evolution based on nucleotide data. Genetics 116: 153–159.
HUTTLEY, G. A., M. W. SMITH, M. CARRINGTON and S. J. O'BRIEN, 1999 A scan for linkage disequilibrium across the human genome. Genetics 152: 1711–1722.
INGVARSSON, P. K., 2005 Nucleotide polymorphism and linkage disequilibrium within and among natural populations of European aspen (Populus tremula L., Salicaceae). Genetics 169: 945–953.
JUNG, M., A. CHING, D. BHATTRAMAKKI, M. DOLAN, S. TINGEY et al., 2004 Linkage disequilibrium and sequence diversity in a 500-kbp region around the adh1 locus in elite maize germplasm. Theor. Appl. Genet. 109: 681–689.[CrossRef][Medline]
KIM, Y., and R. NIELSEN, 2004 Linkage disequilibrium as a signature of selective sweeps. Genetics 167: 1513–1524.
KIMBER, C., 2000 Origins of domesticated sorghum and its early diffusion to India and China, pp. 3–98 in Sorghum, edited by C. W. SMITH and R. A. FREDERIKSEN. John Wiley & Sons, New York.
KUITTINEN, H., and M. AGUADE, 2000 Nucleotide variation at the CHALCONE ISOMERASE locus in Arabidopsis thaliana. Genetics 155: 863–872.
KUMAR, S., K. TAMURA, I. B. JAKOBSEN and M. NEI, 2001 MEGA2: Molecular Evolutionary Genetics Analysis software. Bioinformatics 17: 1244–1245.
LI, N., and M. STEPHENS, 2003 Modeling linkage disequilibrium and identifying recombination hotspots using single-nucleotide polymorphism data. Genetics 165: 2213–2233.
LIN, J. Z., P. L. MORRELL and M. T. CLEGG, 2002 The influence of linkage and inbreeding on patterns of nucleotide sequence diversity at duplicate alcohol dehydrogenase loci in wild barley (Hordeum vulgare ssp. spontaneum). Genetics 162: 2007–2015.
MULLET, J. E., R. R. KLEIN and P. E. KLEIN, 2002 Sorghum bicolor—an important species for comparative grass genomics and a source of beneficial genes for agriculture. Curr. Opin. Plant Biol. 5: 118–121.[CrossRef][Medline]
NEI, M., 1987 Molecular Evolutionary Genetics. Columbia University Press, New York.
NORDBORG, M., 2000 Linkage disequilibrium, gene trees and selfing: an ancestral recombination graph with partial self-fertilization. Genetics 154: 923–929.
NORDBORG, M., and P. DONNELLY, 1997 The coalescent process with selfing. Genetics 146: 1185–1195.[Abstract]
NORDBORG, M., and S. TAVARE, 2002 Linkage disequilibrium: what history has to tell us. Trends Genet. 18: 83–90.[CrossRef][Medline]
NORDBORG, M., J. O. BOREVITZ, J. BERGELSON, C. C. BERRY, J. CHORY et al., 2002 The extent of linkage disequilibrium in Arabidopsis thaliana. Nat. Genet. 30: 190–193.[CrossRef][Medline]
NORDBORG, M., T. T. HU, Y. ISHINO, J. JHAVERI, C. TOOMAJIAN et al., 2005 The pattern of polymorphism in Arabidopsis thaliana. PLoS Biol. 3: e196.[CrossRef][Medline]
PALAISA, K., M. MORGANTE, S. TINGEY and A. RAFALSKI, 2004 Long-range patterns of diversity and linkage disequilibrium surrounding the maize Y1 gene are indicative of an asymmetric selective sweep. Proc. Natl. Acad. Sci. USA 101: 9885–9890.
PRITCHARD, J. K., and M. PRZEWORSKI, 2001 Linkage disequilibrium in humans: models and data. Am. J. Hum. Genet. 69: 1–14.[CrossRef][Medline]
PRZEWORSKI, M., and J. D. WALL, 2001 Why is there so little intragenic linkage disequilibrium in humans? Genet Res. 77: 143–151.[CrossRef][Medline]
REMINGTON, D. L., J. M. THORNSBERRY, Y. MATSUOKA, L. M. WILSON, S. R. WHITT et al., 2001 Structure of linkage disequilibrium and phenotypic associations in the maize genome. Proc. Natl. Acad. Sci. USA 98: 11479–11484.
ROONEY, W. L., and C. W. SMITH, 2000 Techniques for developing new cultivars, pp. 329–347 in Sorghum, edited by C. W. SMITH and R. A. FREDERIKSEN. John Wiley & Sons, New York.
ROZAS, J., and J. C. SANCHEZ-DELBARRIO, 2003 DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19: 2496–2497.
SEMON, M., R. NIELSEN, M. JONES and S. MCCOUCH, 2004 The population structure of African cultivated rice Oryza glaberrima (Steud.): evidence for elevated levels of LD caused by admixture with O. sativa and ecological adaptation. Genetics 169: 1639–1647.
SHEPARD, K. A., and M. D. PURUGGANAN, 2003 Molecular population genetics of the Arabidopsis CLAVATA2 region: the genomic scale of variation and selection in a selfing species. Genetics 163: 1083–1095.
SWIGONOVA, Z., J. LAI, J. MA, W. RAMAKRISHNA, V. LLACA et al., 2004 Close split of sorghum and maize genome progenitors. Genome Res. 14: 1916–1923.
TAJIMA, F., 1989 Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123: 585–595.
TENAILLON, M. I., M. C. SAWKINS, A. D. LONG, R. L. GAUT, J. F. DOEBLEY et al., 2001 Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp. mays L.). Proc. Natl. Acad. Sci. USA 98: 9161–9166.
TENAILLON, M. I., M. C. SAWKINS, L. K. ANDERSON, S. M. STACK, J. DOEBLEY et al., 2002 Patterns of diversity and recombination along chromosome 1 of maize (Zea mays ssp. mays L.). Genetics 162: 1401–1413.
THORNSBERRY, J. M., M. M. GOODMAN, J. DOEBLEY, S. KRESOVICH, D. NIELSEN et al., 2001 Dwarf8 polymorphisms associate with variation in flowering time. Nat. Genet. 28(3): 286–289.[CrossRef][Medline]
WAKELEY, J., and S. LESSARD, 2003 Theory of the effects of population structure and sampling on patterns of linkage disequilibrium applied to genomic data from humans. Genetics 164: 1043–1053.
WALL, J. D., 2000 A comparison of estimators of the population recombination rate. Mol. Biol. Evol. 17: 156–163.
WALL, J. D., P. ANDOLFATTO and M. PRZEWORSKI, 2002 Testing models of selection and demography in Drosophila simulans. Genetics 162: 203–216.
WHITE, G. M., M. T. HAMBLIN and S. KRESOVICH, 2004 Molecular evolution of the phytochrome gene family in sorghum: changing rates of synonymous and replacement evolution. Mol. Biol. Evol. 21: 716–723.
YAUK, C. L., P. R. BOIS and A. J. JEFFREYS, 2003 High-resolution sperm typing of meiotic recombination in the mouse MHC Ebeta gene. EMBO J. 22: 1389–1397.[CrossRef][Medline]
ZHANG, K., P. CALABRESE, M. NORDBORG and F. SUN, 2002 Haplotype block structure and its applications to association studies: power and study designs. Am. J. Hum. Genet. 71: 1386–1394.[CrossRef][Medline]
Communicating editor: J. BERGELSONThis article has been cited by other articles:
 |
|
 |
|
  A. M. Casa, G. Pressoir, P. J. Brown, S. E. Mitchell, W. L. Rooney, M. R. Tuinstra, C. D. Franks, and S. Kresovich Community Resources and Strategies for Association Mapping in Sorghum Crop Sci., January 16, 2008; 48(1): 30 - 40. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J.-B. Veyrieras, L. Camus-Kulandaivelu, B. Gouesnard, D. Manicacci, and A. Charcosset Bridging Genomics and Genetic Diversity: Linkage Disequilibrium Structure and Association Mapping in Maize and Other Cereals Crop Sci., December 18, 2007; 47(Supplement_3): S-60 - S-71. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 K. A. Mather, A. L. Caicedo, N. R. Polato, K. M. Olsen, S. McCouch, and M. D. Purugganan The Extent of Linkage Disequilibrium in Rice (Oryza sativa L.) Genetics, December 1, 2007; 177(4): 2223 - 2232. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. T. Hamblin, M. G. Salas Fernandez, M. R. Tuinstra, W. L. Rooney, and S. Kresovich Sequence Variation at Candidate Loci in the Starch Metabolism Pathway in Sorghum: Prospects for Linkage Disequilibrium Mapping Crop Sci., July 16, 2007; 47(S2): S-125 - S-134. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. M. Casa, S. E. Mitchell, J. D. Jensen, M. T. Hamblin, A. H. Paterson, C. F. Aquadro, and S. Kresovich Evidence for a Selective Sweep on Chromosome 1 of Cultivated Sorghum Crop Sci., November 1, 2006; 46(Supplement_1): S-27 - S-40. [Abstract] [Full Text] [PDF]
|
|
|
|
|
