|
|
|
|
Vol. 10, Issue 10, 1509-1531, October 2000
LETTER
|
| |
ABSTRACT |
|---|
|
|
|---|
The main objectives of the study reported here were to construct a molecular map of wild emmer wheat, Triticum dicoccoides, to characterize the marker-related anatomy of the genome, and to evaluate segregation and recombination patterns upon crossing T. dicoccoides with its domesticated descendant Triticum durum (cultivar Langdon). The total map length exceeded 3000 cM and possibly covered the entire tetraploid genome (AABB). Clusters of molecular markers were observed on most of the 14 chromosomes. AFLP (amplified fragment length polymorphism) markers manifested a random distribution among homologous groups, but not among genomes and chromosomes. Genetic differentiation between T. dicoccoides and T. durum was attributed mainly to the B genome as revealed by AFLP markers. The segregation-distorted markers were mainly clustered on 4A, 5A, and 5B chromosomes. Homeoalleles, differentially conferring the vigor of gametes, might be responsible for the distortion on 5A and 5B chromosomes. Quasilinkage, deviation from free recombination between markers of nonhomologous chromosomes, was discovered. Massive negative interference was observed in most of the chromosomes (an excess of double crossovers in adjacent intervals relative to the expected rates on the assumption of no interference). The general pattern of distribution of islands of negative interference included near-centromeric location, spanning the centromere, and median/subterminal location.
[An appendix describing the molecular marker loci is available as an online supplement at http://www.genome.org.]
| |
INTRODUCTION |
|---|
|
|
|---|
Wild emmer wheat, Triticum dicoccoides
[Triticum turgidum (L) Thell. ssp. dicoccoides
(Koern) Thell.] with genome AABB, was discovered in Northern Israel by
Aaron Aaronsohn in 1906 (Aaronsohn 1910
). It is the tetraploid,
predominantly self-pollinated, wild progenitor from which modern
tetraploid and hexaploid cultivated wheats were derived (Zohary 1970
).
The distribution center of T. dicoccoides is found in the
catchment area of the upper Jordan Valley in Israel and the vicinity
(Nevo and Beiles 1989
). Wild emmer wheat is an important genetic
resource that could be exploited in breeding for resistance to a broad
range of diseases, pests, and for tolerance to poor soil and climatic
factors (Nevo 1983
, 1989
, 1995
). Among many agriculturally important
characteristics already found in T. dicoccoides is its
resistance to stripe rust (Gerechter-Amitai and Stubbs 1970
; Nevo et
al. 1986
; Fahima et al. 1998
), stem rust (Nevo et al. 1991
), and
powdery mildew (Nevo et al. 1985
). Wild emmer wheat contains rich and
variable genetic resources that will play a major role in future wheat
improvement (Nevo 1983
,1989
, 1995
).
Genetic research on important agronomic and quality traits in this
plant species has lagged far behind other cereals, so at present
detailed genetic information is urgently required as a basis for
effective utilization in the breeding programs of cultivated wheat
including bread and durum wheats. Molecular linkage maps of many plant
species have been obtained recently and utilized in quantitative trait
analysis, gene tagging, genome organization, and evolutionary studies,
as well as in improved selection activities (Paterson et al. 1991
;
Whitkus et al. 1994
; Blanco et al. 1998
). Restriction fragment length
polymorphism (RFLP) markers have been used extensively to construct
genetic maps in many cultivated species (Phillips and Vasil 1994
).
Bread wheat, Triticum aestivum (L.) Thell., has received much
attention and several RFLP-based maps either for groups of homoeologous
chromosomes (Chao et al. 1989
; Devos et al. 1992
, 1993
; Nelson et al.
1995
; Van Deynze et al. 1995
; Gill et al. 1996a
,b
) or for all the
chromosome groups (Liu and Tsunewaki 1991
; Anderson et al. 1992
; Gale
et al. 1995
; Mingeot and Jacquemin 1999
) have been reported. In
contrast, tetraploid durum wheat has received little attention and an
RFLP-based linkage map has been published only recently (Blanco et al.
1998
).
Molecular markers can provide a spectacular improvement in the
efficiency and sophistication of plant breeding. It is now generally
accepted that markers represent the most significant advance in
breeding technology that has occurred in the last few decades and are
currently the most important application of molecular biology to plant
breeding (Langridge and Chalmers 1998
). Compared with RFLPs,
microsatellites are PCR-based, easily handled, cheaper to use, suitable
for automation, highly reproducible, and therefore can be used on a
larger scale. Recently, Röder et al. (1998)
developed a set of
hexaploid wheat microsatellite markers, and constructed a molecular map
consisting of 279 microsatellite loci amplified by 230 primer sets. The
efficiency of these primer sets in analysis of T. dicoccoides
genomic DNA was demonstrated previously both for mapping (Chagué
et al. 1999
; Peng et al. 1999
, 2000
) and population genetics purposes
(Fahima et al. 1998
; Li et al. 2000
). Seventy-nine of these
microsatellite markers were integrated into the above-mentioned
RFLP-linkage map in durum wheat (Korzun et al. 1999
).
The amplified fragment length polymorphism (AFLP) technology is based
on the amplification of selected restriction fragments of a total
genomic digest by PCR, and separation of labeled amplified products by
denaturing polyacrylamide gel electrophoresis. A great advantage of the
AFLP approach is that it allows simultaneous identification of a large
number of amplification products (Van Eck et al. 1995
; Vos et al.
1995
). Compared with RFLP or other PCR-based marker systems, AFLP is
fast, reliable, and cost-effective. It may be a good supplement to
other marker systems in species like wheat that give a low-level
polymorphism (such as RFLP). It may be especially useful for high
density mapping in regions containing genes of interest (Ma and Lapitan
1998
). In wheat, a complete AFLP-linkage map has not yet been seen in
formal publications.
The main objectives of the present study were, therefore, to construct a DNA molecular genetic map in wild emmer wheat, T. dicoccoides, by use firstly of microsatellite markers and then of AFLP and random amplified polymorphic DNA (RAPD) markers to fill the gaps, so as to contribute to the understanding of genome evolution in the wheat progenitor and to accelerate the utilization of this important genetic resource in wheat improvement programs. Our aim was also to use the generated information to characterize the marker-related anatomy of T. dicoccoides genome and the segregation and recombination patterns on crossing with its domesticated descendant, Triticum durum (cultivar Langdon).
| |
RESULTS |
|---|
|
|
|---|
Polymorphism between the Parental Lines, T. dicoccoides and T. durum
Of the 203 microsatellite primer pairs, 187 (92.1%) generated
various levels of polymorphism between the parental lines of the
mapping population. All the 33 AFLP primer combinations could detect
the polymorphism between the two parental lines. They amplified 3593 fragments in total, and each of them amplified 109 fragments on
average, with a variation range of 78-128. Among the amplified AFLP
fragments, 24-43 were polymorphic between the two parental lines for
individual primer combinations with an average polymorphism rate of
30.4%, and range of 23%-40%. For the 11 primer combinations chosen
to genotype the mapping population (Table 1), 1241 AFLP fragments were
amplified in total, of which 408 (32.9%) were polymorphic and 315 of
the 408 (77.2%) were mapped onto specific chromosomes (Appendix 1, available as online supplement at www.genome.org). Of the 437 efficient
RAPD primers, 215 were found to produce polymorphism, with an average
polymorphism per primer of 49.2%. However, per-band polymorphism was
only 11%. In total, 14 RAPD primers detected 39 segregating marker
loci (Table 2).
|
|
Construction of Molecular Genetic Map
When 14 previously mapped microsatellite markers (Röder et al.
1998
; Appendix 1) were used as anchors of the 14 chromosomes of T. dicoccoides (one for each chromosome), 543 marker loci, the
YrH52 stripe-rust resistance gene and Ws
gene-conferring spike glaucousness were assigned to the 14 chromosomes
at a minimum LOD value of 2.5. Among the 544 assigned loci, 428 (78.5%) had LOD scores over 10, 110 (20.2%) had LODs ranging from 3.1 to 10, and only seven LODs (1.3%) ranged from 2.5 to 3.0 (Appendix 1, available online at http://www.genome.org). Therefore, the chromosomal assignment of the loci or the construction of the 14 linkage groups proved highly reliable (Lincoln et al. 1992).
Among the 315 assigned AFLP fragments, 20 bands (10 pairs) amplified by the same primer combinations were closely linked in repulsion phase, and were converted into 10 codominant AFLP markers. Thus, a total of 549 loci, comprising 545 PCR-based marker loci, two RFLP loci and two genes, YrH52 and Ws, were included in the mapping analysis. By means of the three- and multipoint analysis of MAPMAKER 3.0b (using Kosambi mapping function) at a minimum LOD of 3.0, two molecular genetic maps for each of the 14 chromosomes of wild emmer wheat, T. dicoccoides, were constructed (Fig. 1), each using codominant markers and dominant markers in coupling phase. The dominant markers of the H and L maps were based on PCR bands amplified from genomic segments of T. dicoccoides (the parental genotype derived from the Hermon population) and T. durum (Langdon), correspondingly. The H and L maps span a distance of 3169 cM and 3180 cM, respectively. The mean interval lengths of H and L maps were 9.9 cM and 10.7 cM, respectively (Table 3).
|
Distribution of Molecular Markers among the Genomes and Chromosomes
The Entire Set of Markers
Among the 549 loci, 231 (42.2%) were located on A genome and 318 (57.8%) on B genome (Table 3). The total lengths of H and L maps for A genome were 1589 cM and 1607 cM, and those for B genome were slightly shorter, 1580 cM and 1573 cM, respectively. The difference of map length between A and B genome was not significant. The mean interval length of A genome (11.1 cM and 12.8 cM for H and L maps, respectively) was significantly greater than that of B genome (8.9 cM and 9.1 cM for H and L maps, respectively). The marker density of B genome was thus obviously higher than that of A genome. The number of marker loci per chromosome ranged from 17 to 58, with all large-marker-number (>40) chromosomes being of the B genome. Except for group 4 chromosomes, B chromosomes amplified obviously more markers than the corresponding A chromosomes. The variations of the number of markers among genomes and among homoeologous groups were highly significant as revealed by log-linear analysis (Table 3).AFLP Markers
From the 11 primer combinations, 21-36 AFLP fragments were mapped onto specific chromosomes, amounting to a total of 315. Among these, 188 (59.7%) were mapped to B genome, and the other 127 (40.3%) to A genome. The number of AFLP markers obviously varied with chromosomes: There were 13-27 AFLP fragments on A chromosomes, and 12-37 on B chromosomes. Log-linear analysis indicated that the effect of the genome on the distribution of AFLP fragments was highly significant, as was the interaction "genome × homologous group". The effect of homoeologous group was not significant (Table 3). Table 4 also shows that AFLP fragments amplified by different primer combinations had various distribution patterns among genomes and chromosomes. The AFLP fragments derived from P55M53, P55M56, P57M49, and P57M51 primer combinations were mainly (>60%) located on the B genome. For P55M60, P56M50, P56M52, P56M53, P56M60, and P57M53 primer combinations, the number of AFLP fragments mapped on the B genome also exceeded 50% (51.5%-57.1%). Only for P57M52 did the numbers of the AFLP fragments mapped onto the A genome slightly exceed those on the B genome.
|
|
Clustering of Markers
Interval length obviously varied within and among chromosomes with a coefficient of variation (CV) of range 63%-125% in the entire genome. The weighted mean CV of B genome was larger than that of A genome (Table 3). This may reflect, to some extent, marker clustering observed on most of the chromosomes for A and B genomes (Fig. 1). The significance of clustering (on the levels of A, B, and the entire genome for both H and L maps) was tested by comparison of the observed distribution of intervals with different marker numbers with Poisson distribution (Korol et al. 1994
|
|
Segregation Distortion of the Molecular Markers
In the present study, a total of 573 marker loci including 203 microsatellites, 326 AFLPs, 39 RAPDs, two RFLPs and one SCAR (Appendix
1), and two genes (YrH52 and Ws) were scored.
2 test was used to check whether the marker segregation
in F2 fitted the Mendelian model (1:2:1 for
codominant and 3:1 for dominant markers). In total, 34 (5.9%) loci
showed significant (P < 0.05) deviation from the expected
ratio, which is close to the number one would expect to get by chance
in a test including many loci, even if no real distortion exists at
all. However, significant deviation at P < 0.01 was
manifested by 19 (19/571 = 3.3%) out of these 34 loci, i.e., three
times that expected by chance (deviating loci are marked * and ** in
Fig.1 and Appendix 1, for P < 0.05 and
P < 0.01, respectively).
Besides the threefold excess of highly distorted loci over the level
expected by chance, the genomic distribution of these loci indicates
the reality of segregation distortion in our material. Indeed, of the
34 segregation-distorted markers, 31 were assigned to specific
chromosomes. Out of these 31 markers, one was assigned to each of 2A,
2B, 3A, 6A, and 7A chromosomes, two were assigned to 1B, four (12.9%)
to 4A, five (16.1%) to 5B, and 15 (48.4%) to 5A chromosome (Table
6, Fig. 1, Appendix 1). Actually, the segregation-distorted markers clustered in some specific regions on the
5A (three clusters along the chromosome), and one cluster was present
in each of the 5BS and 4AS chromosomes (Fig. 1). Distortion favored
T. durum alleles for 23 segregation-distorted loci, including all the loci on 5A and 5B, P57M52u on 1B, P56M50g on 4A, and P57M52x on
6A. Only for five loci, P55M60r on 1B, P56M52o on 2A, P55M53l on 2B,
P56M53a on 3A, and P56M60w on 4A, was the bias toward the T. dicoccoides alleles (Table 6). Thus, the distorted loci showed a
significant (
2 = 11.57, P < 0.001) bias
toward a deficit of alleles of the pollen parent of the hybrid.
|
Nonrandom Segregation of Markers from Nonhomologous Chromosomes
With independent segregation of loci from nonhomologous chromosomes,
the frequency of parental and nonparental combinations of corresponding
alleles is 50%. In 1953 Michie and Wallace observed a departure from
random segregation of markers on nonhomologous chromosomes in crosses
between different strains of house mouse (Mus musculus). Such
a departure from independent assortment of unlinked genes was termed
"quasi-linkage." This phenomenon was observed both in plants and
animals, especially in interspecific hybrids (cotton, tomato, maize,
Coix, mice, mule, dipteran insect Sciara coprophila,
human) (for review, see Sapre and Deshpande 1987
; Korol et al. 1994
).
High genome coverage by molecular markers makes our F2 mapping population (T. durum × T. dicoccoides) especially suitable for testing quasi-linkage. Several pairs of nonhomologs manifested a highly significant deviation of recombination of their markers from the expected r = 0.5 level. These included either excess or deficit of recombinant genotypes, with the following highest deviating r values: r(1A,2A) = 0.701, r(1A,6A) = 0.778, r(2A,2B) = 0.768, r(2A,7B) = 0.271, r(5A,7A) = 0.723, and r(3B,5B) = 0.726.
However, the high pair-wise significance of observed deviations may be
an artifact caused by multiple comparisons. Indeed, if only one marker
per chromosome is considered, a total number of possible nonhomologous
pairs (and tests) for a genome with n = 14 is
n(n
1)/2 = 91. The situation is even more complicated when
each chromosome is represented by multiple markers. To cope with these
complications in evaluating the genome-wise significance of the
observed numerous manifestations of quasi-linkage, we employed a
computing-intensive permutation test. Its result indicate that quasi-linkage in our case could be declared "significant"
(P < 0.02). Spurred by these results, we conducted a few
other tests of quasi-linkage, with hexaploid wheat, maize, and
Arabidopsis (data obtained from http://wheat.pw.usda.gov,
http://ars-genome.cornell.edu/rice, http://www.agron.missouri.edu,
and http://ukcrop.net/agr). It appeared that quasi-linkage might
indeed be a significant phenomenon in wheat: For hexaploid wheat we
obtained genome-wise significance (P < 0.001). The same
test gave P = 0.06 for maize, P = 0.05 for rice,
and P = 0.08 for Arabidopsis.
It is noteworthy that deviations of r values of linked loci to an
unlinked one were positively correlated as expected if the deviations
from r = 0.5 are not related to misclassification of markers. In such
a case, it would be of special interest to conduct a test that will
include not only the "representative markers" of nonhomologs
showing the highest deviations, but to take into account the average
recombination rates characteristic of entire segments (Fig.
3). Application of permutation tests in
such a formulation seems to provide much more reliable results (Table 7). The presented data show that the
observed phenomenon involves quite large segments of nonhomologous
chromosomes. In some cases, even the mean(r) ± 3
interval
build using individual marker-marker recombination rates did not
include the expected 50%. The permutation tests, conducted for the
marker segments, revealed 16 pairs of segments with significant
(P < 0.05) deviation of mean(r) from 50% (Table 7)
should be corrected for multiple comparisons. Indeed, out of 91 chromosome pairs, five pairs may exceed the significance level by
chance. However, the probability to observe 16 such pairs by chance
(H0 hypothesis) is low (<10
5, by binomial
test). Moreover, eight of the foregoing 16 pairs were significant at
P < 0.01 and three at P < 0.001;
probability to obtain such results when H0 is true is
extremely low (<10
7).
|
|
Negative Crossover Interference
The obtained maximum likelihood (ML) estimates of coefficient of
coincidence (c) indicate that negative crossover interference (i.e., c > 1) is characteristic of our T. durum
× T. dicoccoides hybrid. In fact, it is
manifested in some regions of all chromosomes in both genomes A and B,
and for both map versions, H and L (Table 8). Deviations from the `no interference'
Haldane recombination scheme were highly significant. The maximum (per
chromosome) values of
2 (df = 1) ranged,
correspondingly, from 5.65 (with c = 4.35) for one of the
islands of chromosome 5B to 59.25 (c = 3.51) for one of the
islands on chromosome 6B. Clearly, the deviations from Kosambi
interference in such cases should be even more significant. Each
chromosome manifested a few islands of negative interference, from 1 to
2-3. As a rule, these islands were located in proximal regions, either
spanning the centromere (chromosomes 1A, 3B, 4A, 6B, 7A, 7B), or
excluding it. In the latter cases, the location of the islands could be
to one (chromosomes 1B, 2A, 3A, 4B) or both (5A, 5B) sides of the
centromere. When two or more islands were found per chromosomal arm,
the location of extra islands was either median (3B) or
terminal/subterminal (1B, 5A, 5B, and 7A).
|
The revealed trends were confirmed for the foregoing regions by using
larger intervals from one or both sides of the central point of any
considered interval pair and by variation of the position of the
central point. Correspondence between the results of the two map
versions, H and L, was employed as an additional important test for
existence of a tract of negative interference (Table 8). As in our
previous results with chromosome 1B (Peng et al. 1999
), a tendency for
a higher level of negative crossover interference was found in regions
proximal to or spanning the centromere. Likewise, alternation of
negative interference by segments with strong positive interference
found earlier in chromosome 1B (Peng et al. 1999
) appears to be a
general phenomenon. In particular, for many such regions with positive
interference, the estimated value of c was significantly
smaller than predicted by Kosambi mapping function.
| |
DISCUSSION |
|---|
|
|
|---|
High Polymorphism between the Two Parental Lines
Hexaploid bread wheat shows relatively low levels of polymorphism
for RFLP loci, most likely as a result of its narrow genetic base (Chao
et al. 1989
). Usually <10% of all RFLP loci are polymorphic in an
intraspecific context (Röder et al. 1998
). However, in the
genetically distant cross of hexaploid bread wheat, Opata 85 × W7984, 72% of the RFLP probes (Mingeot and Jacquemin 1999
) and 80% of the microsatellite primers (Röder et al. 1998
)
exhibited polymorphism between the parents. In another relatively
distant cross of hexaploid wheat, T. aestivum × T.
spelta, 64% of the RFLP probes were polymorphic (Messmer et al.
1999
). In a distant cross of tetraploid wheat, Messapia × MG4343,
which is similar to our mapping population, 70.1% RFLP probes (Blanco
et al. 1998
) and 84.4% microsatellites (Korzun et al. 1999
) detected
polymorphism between the parents. Therefore, for distant crosses of
both hexaploid and tetraploid wheats, about 70% RFLP probes and 80%
microsatellite primer sets can detect polymorphism.
The level of polymorphism of microsatellite primer sets in the present mapping population was even higher (>90%). Furthermore, all the 33 AFLP primer combinations screened detected polymorphism between the parents, 30.4% of overall AFLP loci were polymorphic, and 49.2% of the RAPD primers detected polymorphism. This clearly shows that highly significant genetic differentiation occurred between T. durum and T. dicoccoides during their evolutionary divergence caused by domestication, even though they share the same A and B genomes and display no genetic or reproductive obstacles.
In the present study, the two parental lines, T. dicoccoides accession Hermon H52 and T. durum cultivar Langdon, were also highly polymorphic in morphological and agronomic traits, i.e., plant height, growth duration, leaf shape, grain characteristic, spike fragility, stripe-rust resistance, and so on (data not shown). The high level of DNA polymorphism between the two parental lines makes it possible to map these domestication-related and agronomically important traits, and hopefully to map many candidate genes in wheat. We now summarize the results of such an analysis, which will be published elsewhere.
Large Genome Coverage of the Molecular Map
It is estimated that the genome sizes of tetraploid wheat (T. dicoccoides and T. turgidum) and hexaploid wheat (T. aestivum) are ~ 1.2 × 1010 and
1.7 × 1010 bp, respectively (Bennett et al. 1998
). If
the average chromosome length is assumed to be 200 cM (Messmer et al.
1999
), the total map size of hexaploid wheat is estimated as 4200 cM.
The sizes of published maps for the populations Chinese Spring
× T. spelta (Liu and Tsunawaki 1991
), Chinese
Spring × Synthetic (Gale et al. 1995
), Forno × Oberkulmer
(Messmer et al. 1999
) and W7984 × Opata 85 (Faris et al. 2000
) in
hexaploid wheat were 1801, 2575, 2469, and 3700 cM, respectively. These
maps cover 43%, 61%, 59%, and 88% of the entire genome of hexaploid
wheat, respectively.
The total map size of tetraploid wheat (T. dicoccoides
and T. turgidum) could be estimated as 2800 cM if the size
of hexaploid wheat is 4200 cM (Messmer et al. 1999
). In tetraploid
durum wheat, the total length of published RFLP-based genetic map is
1352 cM (Blanco et al. 1998
), and has been extended to 2035 cM after
integrating 79 microsatellite markers (Korzun et al. 1999
). These two
maps cover 48% and 73% of the entire genome of tetraploid wheat,
respectively. In the present study, the total lengths for the H and L
molecular maps in tetraploid T. dicoccoides were 3169 and 3180 cM, respectively (Table 3), therefore covering the entire genome
of T. dicoccoides. The present molecular map (Fig. 1) thus has
the largest relative coverage in wheat, and covers the entire genome.
Another advantage of our map is its relatively large population size
[150 F2 individuals compared to 65-120 in other studies
(Gale et al. 1995
; Blanco et al. 1998
; Röder et al. 1998
)]. This
is important for mapping quantitative trait loci (QTL) that confer the
evolutionarily and agronomically important traits of this species (Peng
et al., in prep.).
In the present map, AFLP markers terminate the molecular maps at one or
both ends of all chromosomes except 1B and 3B (Fig. 1). We assume that
the high coverage of the present molecular maps has resulted from the
integration of AFLP markers derived from PstI restriction
enzyme. AFLP is more effective with wheat genomic DNA if PstI
rather than EcoRI is used to generate DNA fragments. This is
probably because of the high G + C content of PstI
recognition sites relative to EcoRI resulting in preferential targeting of low-copy/gene-rich regions of the genome.
Consequently, use of PstI rather than EcoRI allows
better genome coverage and less clustering of marker loci (Langridge
and Chalmers 1998
). In maize, the total map length also increased
significantly due to the addition of PstI AFLP markers
to the telomeric regions where RFLP markers were represented poorly
(Castiglioni et al. 1999
).
Highly Conserved Order of Microsatellite Loci and Structural Changes of Chromosomes
The chromosomal locations and map orders for most of the
microsatellite loci in our present molecular map of T. dicoccoides (Fig. 1) are the same as those in T. aestivum
(Röder et al. 1998
) and in T. durum (Korzun et al.
1999
). This result indicates that the macro-chromosomal organization is
highly conserved among various species in the genus Triticum.
It also further corroborates that wheat microsatellite markers are
mainly genome-specific, and the microsatellite markers developed in
hexaploid wheat are easily transferred to other species in the genus
(Röder et al. 1998
, Korzun et al. 1999
).
However, some changes of chromosomal location and map order were also
observed for a few microsatellite markers in the present map in
comparison with the maps of Röder et al. (1998)
and Korzun et al.
(1999)
. The inconsistency of map order may be explained by the presence
of additional loci and structural changes of chromosomes in the wheat
genome (Korzun et al. 1999
), and limitations caused by the population
size [70 and 65 recombinant inbred lines in Röder et al. (1998)
and Korzun et al. (1999)
, respectively, as against 150 F2
individuals in the present study]. For example, the microsatellite
locus Xgwm550 was mapped to 1BL by Korzun et al. (1999)
. In contrast,
Xgwm550 was mapped to 1BS by Röder et al. (1998)
and in the
present study (Fig. 1). Thus, the Xgwm550 locus in Korzun et al. (1999)
is definitely different from that in Röder et al. (1998)
and in
the present study (Fig. 1). The Xgwm498 was mapped to 1A and 1B by
Korzun et al. (1999)
and Röder et al. (1998)
, respectively. So
there should be two loci for this microsatellite and these two loci
were mapped to chromosome 1A and 1B, respectively, in our present study
(Fig. 1). In wheat, paracentric and pericentric inversions have been
repeatedly observed on chromosome 4A (Liu et al. 1992
; Devos et al.
1995
). Faris et al. (2000)
observed an inversion within a small segment
at the proximal region of T. dicoccoides 5B map involving
several closely linked markers. Thus, the changes of the map order of
Xgwm397, Xgwn610, Xgwm601a, and Xgwm165a, genetically near the
centromere of chromosome 4A, are possibly due to the inversion of this
chromosome segment, but for this we should assume that the inversion is
present in both T. dicoccoides and T. durum.
Joppa et al. (1995)
found that a high proportion (70%) of T. dicoccoides genotypes had translocations, as was evident in the study of Kawahara and Nevo (1996)
. Compared with the map of Röder et al. (1998)
, the changes of chromosomal location of
Xgwm582b, Xgwm344, and Xgwm526 may be explained by the translocations
of chromosomal segments between 1B versus 2A, 7B versus 7A, and 2B versus 5A, respectively. Nonhomoeologous translocations have been repeatedly reported in hexaploid wheat (Liu et al. 1992
; Devos et al.
1993
) and tetraploid wheat (Blanco et al. 1998
). Translocation may be
one of the forces maintaining the high genetic diversity of T. dicoccoides under the diverse natural environments in Israel (Joppa
et al. 1995
; Kawahara and Nevo 1996
). Translocation frequencies (TF) of
various populations were correlated with environmental variables,
primarily with water availability and humidity, and possibly also with
soil type. In general, TF was higher in peripheral populations in
ecologically heterogeneous frontiers of species distribution than in
the central populations located in the catchment area of the upper
Jordan valley (Joppa et al. 1995
).
Distorted Segregation of Molecular Markers
Distorted segregation of molecular markers has been observed in
mapping populations derived from intra- and interspecific hybrids in
many plants, including potato (Gebhardt et al. 1989
), corn (Gardiner et
al. 1993
), rice (Causse et al. 1994
), common bean (Vallejos et al.
1992
), and barley (Heun et al. 1991
). In wheat, this phenomenon has
also been reported repeatedly (Liu and Tsunewaki 1991
; Devos et al.
1993
; Nelson et al. 1995
; Blanco et al. 1998
; Messmer et al. 1999
). In
the present study, the portion of segregation-distorted marker loci
(5.9%) was relatively low compared with the previous studies in wheat
(Liu and Tsunewaki 1991
; Blanco et al. 1998
; Messmer et al. 1999
). The
possible causes for segregation deviation of molecular markers are
chromosomal rearrangement (Tanksley 1984
) and gametic or zygotic
selection (Nakagahra 1986
). Further indications of causes of deviation
such as presence of lethals, meiotic drive, and chromosomal
rearrangements could be obtained from the analysis of additional
populations with segregation distortion (Blanco et al. 1998
).
It is possible to discern whether the cause of distorted segregation is
gametic competition or zygotic selection by estimating the frequencies
of two alleles for a codominant locus. The obtained data indicate that
the gametes carrying T. durum alleles have stronger vigor and
higher competition ability than those with T. dicoccoides
alleles. Therefore, the cause for distorted segregation on 5A and 5B
would be gametic competition (Table 7). It is noteworthy that T. dicoccoides was a pollen-parent of the F1 hybrid, hence the cytoplasm was provided by T. durum. We could speculate
that the indicated regions with presumably gametic selection carry loci
involved in nuclear-cytoplasmic interaction. Faris et al. (1998)
reported a segregation distortion locus within the homoeologous region
of chromosome 5D in Aegilops tauschii. They also postulated a
possible homeoallele of the distortion factor on 5B causing skewed
segregation of two markers (Faris et al. 2000
). It seems that a common
locus confers, via differentially affecting the vigor of gametes, the
segregation distortion on group 5 chromosomes (5A, 5B, and 5D) in wheat.
Nonrandom Segregation of Nonhomologous Chromosomes
Fifty years ago, a departure from random segregation of markers on
nonhomologous chromosomes observed in crosses between different strains
of house mouse, was explained by a mutual attraction of the segregating
chromosomes of the same origin to migrate to the same pole during
meiosis (the "affinity" hypothesis) (Michie 1953
; Wallace 1953
).
Consequently, gametes with parental combinations of alleles at loci of
nonhomologous chromosomes should appear with a higher frequency than
the recombinant gametes. The departure from independent segregation of
unlinked genes was termed quasi-linkage. Reviewing the data in the
literature, Wallace (1960a)
postulated possible chromosome affinity in
cotton. Even earlier, Malinowsky (1927)
explained, based on the
chromosome affinity hypothesis, a number of cases of abnormal
segregation in hybrids of lettuce, peas, tobacco, beans, and wheat. He
was probably the first to introduce the term "affinity." A
pronounced effect of quasi-linkage was reported in an interspecific
hybrid of Coix (Sapre and Deshpande 1987
). There is evidence
for nonrandom segregation of nonhomologs during the first meiotic
division in man (see Driscoll et al. 1979
). Quasi-linkage has also been
established in tomato (Wallace 1960b
; Zhuchenko et al. 1977
; Korol et
al. 1989
, 1994
).
Besides nonrandom assortment (or affinity), quasi-linkage can also be
explained by recourse to mechanisms of differential viability. Let us
consider the MiMj/mimj heterozygote resulting from a cross between
MiMj/MiMj and
mimj/mimj. Suppose that
parental combinations of whole chromosomes (MiMj
and mimj) and recombinant combinations
(Mimj and miMj) have
differential selective values. Experimental evidence of this type has
been obtained in Drosophila (Dobzhansky et al. 1965
).
Selective differences between parental and recombinant combinations may
be reflected in the differential viability of zygotes, embryos, and
adult individuals. The effect can also be expressed at the gamete stage
in the form of differential fertilizing capacity of spermatozoa,
differential rates of pollen tube growth, etc. (Korol et al. 1989
,
1994
).
The major problem with the old data on quasi-linkage was poor genome
coverage of the morphological markers that could be followed up in one
cross. This strongly limited the possibilities of discriminating among
different explanatory hypotheses. Molecular markers solve this problem,
but at the price of another one: The size of the available segregating
populations is usually very small, restricting the detection power of
the tests, so that only big deviations could be declared significant.
Nonetheless, our data (Table 7), together with the few examples on
molecular marker segregation in hexaploid wheat, rice, maize, and
Arabidopsis indicate that quasi-linkage may be a much more
common phenomenon in plants (especially, cereals) than ever thought
before. One practical aspect is related to a possible effect of
quasi-linkage on the rate of false positive detection in QTL mapping.
The widespread approach of multilocus (composite) interval mapping
(Zeng 1994
; Jansen and Stam 1994
) may be helpful in such cases.
Marker Distribution in the Genome
Clustering of Marker Loci
Marker distribution along the present molecular maps of tetraploid wheat was far from uniform, with clusters of tightly linked loci and regions with low density of markers. Most of the 14 chromosomes had marker clusters in the centromeric regions of the genetic maps (Fig. 1). A statistical test proves significance of the marker clustering on the B genome, but not on the A genome despite clusters on a few A chromosomes (Table 5). This feature has been observed in most mapping studies of wheat (Chao et al. 1989Nonrandom Distribution of AFLP Markers
More than half (59.7%) of the AFLP markers in the present study were mapped to the B genome (Table 4). This difference between A and B genome is highly significant (P < 0.01) and implies nonrandom distribution of AFLP markers among A and B genomes of tetraploid wheat. The same trend has also been observed for microsatellites (Röder et al. 1998Negative Crossover Interference
Positive crossover interference, i.e., a reduced frequency of
adjacent double crossovers compared to that expected with the assumption of independence, is a characteristic of meiotic organisms, with only a very few exceptions (Egel-Mitani et al. 1982
).
Consequently, it is generally assumed that negative crossover
interference is mainly associated with intragenic recombination. Still,
cases are known of higher than expected frequency of double crossovers in adjacent segments of small genetic but large physical length. In
Drosophila melanogaster, within a segment 4 cM long accounting for about 25% of the cytological length of chromosome 3 and spanning the centromere, a significant excess of multiple exchanges has been
found (Sinclair 1975
). Similar results have been obtained in other
Drosophila studies with autosomes (Green 1975
; Dennell and
Keppy 1979
; Korol et al. 1994
), but not with the X-chromosome (Lake
1986
). Significant negative crossover interference was found in barley
(Søgaard 1977
). Denell and Keppy (1979)
suggested that negative
chromosome interference could be a characteristic of all regions
exhibiting a very low density of recombination per unit physical
length. This hypothesis fits the results of our study. We found a
significant excess of double exchanges in segments spanning, or
proximal to the centromere, in nearly all chromosomes of both genomes,
A and B (Table 8). These proximal segments comprise about 50%-70% of
the chromosomes cytological length but only 5%-20% of the genetic
length (Lukaszewski and Curtis 1993
; Gill et al. 1996a
). In some
chromosomes, additional islands of negative interference were found in
median or subterminal regions. An alternation of strong negative and
full positive interference was characteristic of our data.
Interest in the problem of coinciding crossovers is due to the current
large-scale genome mapping efforts and the growing evidence that the
length of genetic maps of some plants tends to increase with the number
of molecular markers employed. The simplest explanation is to assume
that double crossovers can occur, at least in some organisms, at much
smaller distances than generally believed. Our data support this view,
corroborating other results (Gill et al. 1996b
; Takahashi et al. 1997
)
and an apparent noncorrespondence between the total chiasma frequency
and genome length of some species, including cereals (Nilsson et al.
1993
).
The simplest explanation that the size of the mapping population (n = 150) is too small to allow reliable conclusions about such kinds of patterns cannot be considered plausible. How could one then explain the high resemblance of different chromosomes and excellent correspondence between the two map versions (H and L) as well as between dominant and codominant markers? Moreover, in light of the unusual fact of massive negative interference manifested by our data, we made some preliminary estimates of interference using recent mapping data available from public domain websites. To our surprise, we discovered that negative interference seems not to be a rare phenomenon among other plant and animal species (Peng et al., in prep.).
The observed chromosomal distribution of islands of negative
interference (either near-centromeric or median/subterminal) and the
alternating positive/negative interference pattern make it possible to
recruit two explanatory models: One is based on the foregoing
hypothesis of Denell and Keppy (1979)
that negative interference could
be a characteristic of regions with low density of recombination per
unit physical length, and the other is based on the recent findings in
cereal genomics (Gill et al. 1996a
,b
; Faris et al. 2000
; Kunzel et al.
2000
) indicating the existence of gene-rich segments in wheat and
barley chromosomes and higher recombination rates in these regions than
in gene-poor segments. We can assume that the positive-negative
interference doubted by Lukaszewski and Curtis (1993)
may be a real
phenomenon in wheat, if double crossovers occur within the foregoing
islands and recombination within an island reduces the chance of
crossover in adjacent gene-poor segments. One important aspect of our
results is that the foregoing recombination pattern was revealed mainly
by using microsatellite markers, which may not necessarily follow the
island-like distribution of structural genes [half of the RFLP markers
used by Gill et al. (1996b)
were cDNA probes]. If negative crossover
interference is indeed a real phenomenon in wheat, then gene
introgression via homologous recombination from relatively close wild
relatives to cultivated wheat can be considered even more optimistic
than the estimates based on the association between gene distribution and recombination density in cereal genomes.
| |
METHODS |
|---|
|
|
|---|
Plant Materials and Genomic DNA Extraction
A highly stripe-rust resistant T. dicoccoides accession,
Hermon H52 (H52), from the Mt. Hermon population, Israel, and a T. durum cultivar, Langdon (Ldn), were used to develop an F2
mapping population consisting of 150 individuals. Ldn was used as
the female parent and H52 as the male parent in the cross. Young and healthy leaf samples were collected from each of F2
individuals, the two parental lines (H52 and Ldn) and their
F1 hybrid in the greenhouse, frozen in liquid nitrogen, and
stored at
80°C. Genomic DNA was extracted by means of plant
genomic DNA isolation reagent DNAzol ES (Molecular Research Center,
Inc.) with some modifications. The scoring of stripe-rust resistance
and spike glaucousness was conducted in the field using F3
families (Peng et al. 1999
).
Molecular Marker Analysis
Microsatellites
Microsatellite genotyping was conducted following Röder et al. (1998)AFLP
The AFLP technique was described by Vos et al. (1995)RAPD
The PCR reactions were performed in a PTC-100 Programmable Thermal Controller (MJ Research, Inc.) according to the protocol described by Chagué et al. (1999)RFLP
To clarify the chromosomal location of stripe-rust resistance gene YrH52 carried by H52, we also used the RFLP probe Nor to genotype the mapping population (Peng et al. 1999Data Analysis
Genetic Mapping
Multiple loci amplified by a single primer/primer combination/primer pair had a suffix a, b, c, ... added following the regular marker names based on the fragment size (bp). The description of all the scored marker loci in the present study is listed in Appendix 1. The mapping analysis was conducted by MAPMAKER 3.0b (Lander et al. 1987
2-test were adopted. The Poisson distribution was
used to test the marker clustering on A, B, and the entire genome
(StatSoft, Inc. 1996Testing for Negative Interference
Kosambi (1944)Testing for Segregation Anomalies
These included deviation of monohybrid segregations from the expected ratios (3:1 and 1:2:1), and disturbance of independent segregation of markers belonging to nonhomologous chromosomes.
2-test was used to analyze both types of
anomalies. In addition, to correct for multiple comparisons when
nonrandom cosegregation of unlinked markers was tested, we employed a
permutation test by reshuffling marker sets of entire chromosomes
relative to each other. Namely, if some deviation from the expected
free recombination was observed in real data, say r = 0.5 +
instead of 0.5, then the genome-wise significance was evaluated as a
proportion of permutation runs resulting in estimates of r outside the
interval 0.5 ±
. This test was also repeated with entire
segments of nonhomologous chromosomes using as a score the average r
over all pairs of markers from the defined segments.
| |
ACKNOWLEDGMENTS |
|---|
This work was supported by the following grants: EMBOGrant no. ASTF9195 to J.H.P.; the Israel Discount Bank Chair of Evolutionary Biology; the Ancell-Teicher Research Foundation for Molecular Genetics and Evolution; the Israeli Ministry of Science (Grant no. 5757-1-95); the Israel Science Foundation (Grant no. 9048/99), the German-Israeli Project Cooperation (DIP project funded by the BMBF and supported by BMBF's International Bureau at the DLR), and the Graduate School of the University of Haifa, Israel. The authors thank T. Krugman and A. Dahan at the Institute of Evolution, University of Haifa, Israel for assistance in the experiment and for kindly providing the SCAR199 primer pairs, respectively; K. Wendehake at the Institute of Plant Genetics and Crop Plant Research (IPK), Gatersleben, Germany for her technical assistance in microsatellite genotyping; Dr. H.J. van Eck at the Department of Plant Breeding, Wageningen Agricultural University, the Netherlands for his kind help in AFLP analysis; and to Dr. L.R. Joppa, Northern Crop Science Laboratory, Fargo, USA for kindly providing the F1 seeds of our mapping population.
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| |
FOOTNOTES |
|---|
3 Corresponding author.
E-MAIL nevo{at}research.haifa.ac.il; FAX 972-4-8246554.
Article and publication are at www.genome.org/cgi/doi/10.1101/gr.150300.
| |
REFERENCES |
|---|
|
|
|---|