Alternating between even and odd ploidy levels switches on and off the recombination control, even near the centromeres.

Meiotic recombination is a key biological process in plant evolution and breeding, as it generates genetic diversity in each generation through the formation of crossovers (COs). However, due to their importance in genome stability, COs are highly regulated in frequency and distribution. We previously demonstrated that this strict regulation of COs can be modified, both in terms of CO frequency and distribution, in allotriploid Brassica hybrids (2n = 3x = 29; AAC) resulting from a cross between Brassica napus (2n = 4x = 38; AACC) and Brassica rapa (2n = 2x = 20; AA). Using the recently updated B. napus genome now including pericentromeres, we demonstrated that COs occur in these cold regions in allotriploids, as close as 375 kb from the centromere. Reverse transcription quantitative PCR (RT-qPCR) of various meiotic genes indicated that Class I COs are likely involved in the increased recombination frequency observed in allotriploids. We also demonstrated that this modified recombination landscape can be maintained via successive generations of allotriploidy (odd ploidy level). This deregulated meiotic behavior reverts to strict regulation in allotetraploid (even ploidy level) progeny in the second generation. Overall, we provide an easy way to manipulate tight recombination control in a polyploid crop.

Epigenomic and structural events preclude recombination in Brassica napus.

Meiotic recombination is a major evolutionary process generating genetic diversity at each generation in sexual organisms. However, this process is highly regulated, with the majority of crossovers lying in the distal chromosomal regions that harbor low DNA methylation levels. Even in these regions, some islands without recombination remain, for which we investigated the underlying causes. Genetic maps were established in two Brassica napus hybrids to detect the presence of such large nonrecombinant islands. The role played by DNA methylation and structural variations in this local absence of recombination was determined by performing bisulfite sequencing and whole genome comparisons. Inferred structural variations were validated using either optical mapping or oligo fluorescence in situ hybridization. Hypermethylated or inverted regions between Brassica genomes were associated with the absence of recombination. Pairwise comparisons of nine B. napus genome assemblies revealed that such inversions occur frequently and may contain key agronomic genes such as resistance to biotic stresses. We conclude that such islands without recombination can have different origins, such as DNA methylation or structural variations in B. napus. It is thus essential to take into account these features in breeding programs as they may hamper the efficient combination of favorable alleles in elite varieties.

Screening diversity and distribution of Copia retrotransposons reveals a specific amplification of BARE1 elements in genomes of the polyploid Hordeum murinum complex.

We explored diversity, distribution and evolutionary dynamics of Ty1-Copia retrotransposons in the genomes of the Hordeum murinum polyploid complex and related taxa. Phylogenetic and fluorescent in situ hybridization (FISH) analyses of reverse transcriptase sequences identified four Copia families in these genomes: the predominant BARE1 (including three groups or subfamilies, A, B and C), and the less represented RIRE1, IKYA and TAR-1. Within the BARE1 family, BARE1-A elements and a subgroup of BARE1-B elements (named B1) have proliferated in the allopolyploid members of the H. murinum complex (H. murinum and H. leporinum), and in their extant diploid progenitor, subsp. glaucum. Moreover, we found a specific amplification of BARE1-B elements within each Hordeum species surveyed. The low occurrence of RIRE1, IKYA and TAR-1 elements in the allopolyploid cytotypes suggests that they are either weakly represented or highly degenerated in their diploid progenitors. The results demonstrate that BARE1-A and BARE1-B1 Copia elements are particularly well represented in the genomes of the H. murinum complex and constitute its genomic hallmark. No BARE1-A and -B1 homologs were detected in the reference barley genome. The similar distribution of RT-Copia probes across chromosomes of diploid, tetraploid and hexaploid taxa of the murinum complex shows no evidence of proliferation following polyploidization.

Amplifying recombination genome-wide and reshaping crossover landscapes in Brassicas.

Meiotic recombination by crossovers (COs) is tightly regulated, limiting its key role in producing genetic diversity. However, while COs are usually restricted in number and not homogenously distributed along chromosomes, we show here how to disrupt these rules in Brassica species by using allotriploid hybrids (AAC, 2n = 3x = 29), resulting from the cross between the allotetraploid rapeseed (B. napus, AACC, 2n = 4x = 38) and one of its diploid progenitors (B. rapa, AA, 2n = 2x = 20). We produced mapping populations from different genotypes of both diploid AA and triploid AAC hybrids, used as female and/or as male. Each population revealed nearly 3,000 COs that we studied with SNP markers well distributed along the A genome (on average 1 SNP per 1.25 Mbp). Compared to the case of diploids, allotriploid hybrids showed 1.7 to 3.4 times more overall COs depending on the sex of meiosis and the genetic background. Most surprisingly, we found that such a rise was always associated with (i) dramatic changes in the shape of recombination landscapes and (ii) a strong decrease of CO interference. Hybrids carrying an additional C genome exhibited COs all along the A chromosomes, even in the vicinity of centromeres that are deprived of COs in diploids as well as in most studied species. Moreover, in male allotriploid hybrids we found that Class I COs are mostly responsible for the changes of CO rates, landscapes and interference. These results offer the opportunity for geneticists and plant breeders to dramatically enhance the generation of diversity in Brassica species by disrupting the linkage drag coming from limits on number and distribution of COs.

Crossover rate between homologous chromosomes and interference are regulated by the addition of specific unpaired chromosomes in Brassica.

Recombination is a major mechanism generating genetic diversity, but the control of the crossover rate remains a key question. In Brassica napus (AACC, 2n = 38), we can increase the homologous recombination between A genomes in AAC hybrids. Hypotheses for this effect include the number of C univalent chromosomes, the ratio between univalents and bivalents and, finally, which of the chromosomes are univalents. To test these hypotheses, we produced AA hybrids with zero, one, three, six or nine additional C chromosomes and four different hybrids carrying 2n = 32 and 2n = 35 chromosomes. The genetic map lengths for each hybrid were established to compare their recombination rates. The rates were 1.4 and 2.7 times higher in the hybrids having C6 or C9 alone than in the control (0C). This enhancement reached 3.1 and 4.1 times in hybrids carrying six and nine C chromosomes, and it was also higher for each pair of hybrids carrying 2n = 32 or 2n = 35 chromosomes, with a dependence on which chromosomes remained as univalents. We have shown, for the first time, that the presence of one chromosome, C9 , affects significantly the recombination rate and reduces crossover interference. This result will have fundamental implications on the regulation of crossover frequency.

Molecular cytogenetic identification of B genome chromosomes linked to blackleg disease resistance in Brassica napus × B. carinata interspecific hybrids.

KEY MESSAGE: Provide evidence that the Brassica B genome chromosome B3 carries blackleg resistance gene, and also the B genome chromosomes were inherited several generations along with B. napus chromosomes. Blackleg disease caused by fungus Leptosphaeria maculans causes significant yield losses in Brassica napus. Brassica carinata possesses excellent resistance to this disease. To introgress blackleg resistance, crosses between B. napus cv. Westar and B. carinata were done. The interspecific-hybrids were backcrossed twice to Westar and self-pollinated three times to produce BC2S3 families. Doubled haploid lines (DH1) were produced from one blackleg resistant family. SSR markers were used to study the association between B genome chromosome(s) and blackleg resistance. The entire B3 chromosome of B. carinata was associated with blackleg resistance in DH1. A second DH population (DH2) was produced from F1s of resistant DH1 lines crossed to blackleg susceptible B. napus cv. Polo where resistance was found to be associated with SSR markers from the middle to bottom of the B3 and top of the B8 chromosomes. The results demonstrated that the B3 chromosome carried gene(s) for blackleg resistance. Genomic in situ hybridization (GISH) and GISH-like analysis of the DH2 lines revealed that susceptible lines, in addition to B. napus chromosomes, possessed one pair of B genome chromosomes (2n = 40), while resistant lines had either one (2n = 40) or two pairs (2n = 42) of B chromosomes. The molecular and GISH data suggested that the B chromosome in the susceptible lines was B7, while it was difficult to confirm the identity of the B chromosomes in the resistant lines. Also, B chromosomes were found to be inherited over several generations along with B. napus chromosomes.

Prevalence of gene expression additivity in genetically stable wheat allohexaploids.

The reprogramming of gene expression appears as the major trend in synthetic and natural allopolyploids where expression of an important proportion of genes was shown to deviate from that of the parents or the average of the parents. In this study, we analyzed gene expression changes in previously reported, highly stable synthetic wheat allohexaploids that combine the D genome of Aegilops tauschii and the AB genome extracted from the natural hexaploid wheat Triticum aestivum. A comprehensive genome-wide analysis of transcriptional changes using the Affymetrix GeneChip Wheat Genome Array was conducted. Prevalence of gene expression additivity was observed where expression does not deviate from the average of the parents for 99.3% of 34,820 expressed transcripts. Moreover, nearly similar expression was observed (for 99.5% of genes) when comparing these synthetic and natural wheat allohexaploids. Such near-complete additivity has never been reported for other allopolyploids and, more interestingly, for other synthetic wheat allohexaploids that differ from the ones studied here by having the natural tetraploid Triticum turgidum as the AB genome progenitor. Our study gave insights into the dynamics of additive gene expression in the highly stable wheat allohexaploids.

Crossovers get a boost in Brassica allotriploid and allotetraploid hybrids.

Meiotic crossovers are necessary to generate balanced gametes and to increase genetic diversity. Even if crossover number is usually constrained, recent results suggest that manipulating karyotype composition could be a new way to increase crossover frequency in plants. In this study, we explored this hypothesis by analyzing the extent of crossover variation in a set of related diploid AA, allotriploid AAC, and allotetraploid AACC Brassica hybrids. We first used cytogenetic methods to describe the meiotic behavior of the different hybrids. We then combined a cytogenetic estimation of class I crossovers in the entire genome by immunolocalization of a key protein, MutL Homolog1, which forms distinct foci on meiotic chromosomes, with genetic analyses to specifically compare crossover rates between one pair of chromosomes in the different hybrids. Our results showed that the number of crossovers in the allotriploid AAC hybrid was higher than in the diploid AA hybrid. Accordingly, the allotetraploid AACC hybrid showed an intermediate behavior. We demonstrated that this increase was related to hybrid karyotype composition (diploid versus allotriploid versus allotetraploid) and that interference was maintained in the AAC hybrids. These results could provide another efficient way to manipulate recombination in traditional breeding and genetic studies.

Immediate unidirectional epigenetic reprogramming of NORs occurs independently of rDNA rearrangements in synthetic and natural forms of a polyploid species Brassica napus.

The dynamics of genome modification that occurred from the initial hybridization event to the stabilization of allopolyploid species remains largely unexplored. Here, we studied inheritance and expression of rDNA loci in the initial generations of Brassica napus allotetraploids (2n = 38, AACC) resynthesized from Brassica oleracea (2n = 18, CC) and B. rapa (2n = 20, AA) and compared the patterns to natural forms. Starting already from F1 generation, there was a strong uniparental silencing of B. oleracea genes. The epigenetic reprogramming was accompanied with immediate condensation of C-genome nucleolar organizer region (NOR) and progressive transgeneration hypermethylation of polymerase I promoters, mainly at CG sites. No such changes were observed in the A-genome NORs. Locus loss and gains affecting mainly non-NOR loci after the first allotetraploid meiosis did not influence established functional status of NORs. Collectively, epigenetic and genetic modifications in synthetic lines resemble events that accompanied formation of natural allopolyploid species.

Diversity and evolution of the Hordeum murinum polyploid complex in Algeria.

Population diversity and evolutionary relationships in the Hordeum murinum L. polyploid complex were explored in contrasted bioclimatic conditions from Algeria. A multidisciplinary approach based on morphological, cytogenetic, and molecular data was conducted on a large population sampling. Distribution of diploids (subsp. glaucum) and tetraploids (subsp. leporinum) revealed a strong correlation with a North-South aridity gradient. Most cytotypes exhibit regular meiosis with variable irregularities in some tetraploid populations. Morphological analyses indicate no differentiation among taxa but high variability correlated with bioclimatic parameters. Two and three different nuclear sequences (gene coding for an unspliced genomic protein kinase domain) were isolated in tetraploid and hexaploid cytotypes, respectively, among which one was identical with that found in the diploid subsp. glaucum. The tetraploids (subsp. leporinum and subsp. murinum) do not exhibit additivity for 5S and 45S rDNA loci comparative with the number observed in the related diploid (subsp. glaucum). The subgenomes in the tetraploid taxa could not be differentiated using genomic in situ hybridization (GISH). Results support an allotetraploid origin for subsp. leporinum and subsp. murinum that derives from the diploid subsp. glaucum and another unidentified diploid parent. The hexaploid (subsp. leporinum) has an allohexaploid origin involving the two genomes present in the allotetraploids and another unidentified third diploid progenitor.

Genome structure affects the rate of autosyndesis and allosyndesis in AABC, BBAC and CCAB Brassica interspecific hybrids.

Gene introgression into allopolyploid crop species from diploid or polyploid ancestors can be accomplished through homologous or homoeologous chromosome pairing during meiosis. We produced trigenomic Brassica interspecific hybrids (genome complements AABC, BBAC and CCAB) from the amphidiploid species Brassica napus (AACC), Brassica juncea (AABB) and Brassica carinata (BBCC) in order to test whether the structure of each genome affects frequencies of homologous and homoeologous (both allosyndetic and autosyndetic) pairing during meiosis. AABC hybrids produced from three genotypes of B. napus were included to assess the genetic control of homoeologous pairing. Multi-colour fluorescent in situ hybridisation was used to quantify homologous pairing (e.g. A-genome bivalents in AABC), allosyndetic associations (e.g. B-C in AABC) and autosyndetic associations (e.g. B-B in AABC) at meiosis. A high percentage of homologous chromosomes formed pairs (97.5-99.3%), although many pairs were also involved in autosyndetic and allosyndetic associations. Allosyndesis was observed most frequently as A-C genome associations (mean 4.0 per cell) and less frequently as A-B genome associations (0.8 per cell) and B-C genome associations (0.3 per cell). Autosyndesis occurred most frequently in the haploid A genome (0.75 A-A per cell) and least frequently in the haploid B genome (0.13 B-B per cell). The frequency of C-C autosyndesis was greater in BBAC hybrids (0.75 per cell) than in any other hybrid. The frequency of A-B, A-C and B-C allosyndesis was affected by the genomic structure of the trigenomic hybrids. Frequency of allosyndesis was also influenced by the genotype of the B. napus paternal parent for the three AABC (B. juncea × B. napus) hybrid types. Homoeologous pairing between the Brassica A, B and C genomes in interspecific hybrids may be influenced by complex interactions between genome structure and allelic composition.

Evolutionary dynamics of transposable elements and satellite DNAs in polyploid Spartina species.

Repeated sequences and polyploidy play a central role in plant genome dynamics. Here, we analyze the evolutionary dynamics of repeats in tetraploid and hexaploid Spartina species that diverged during the last 10 million years within the Chloridoideae, one of the poorest investigated grass lineages. From high-throughput genome sequencing, we annotated Spartina repeats and determined what sequence types account for the genome size variation among species. We examined whether differential genome size evolution correlated with ploidy levels and phylogenetic relationships. We also examined the tempo of repeat sequence dynamics associated with allopatric speciation over the last 3-6 million years between hexaploid species that diverged on the American and European Atlantic coasts and tetraploid species from North and South America. The tetraploid S. spartinae, whose phylogenetic placement has been debated, exhibits a similar repeat content as hexaploid species, suggesting common ancestry. Genome expansion or contraction resulting from repeat dynamics seems to be explained mostly by the contrasting divergence times between species, rather than by genome changes triggered by ploidy level change per se. One 370 bp satellite may be exhibiting 'meiotic drive' and driving chromosome evolution in S. alterniflora. Our results provide crucial insights for investigating the genetic and epigenetic consequences of such differential repeat dynamics on the ecology and distribution of the meso- and neopolyploid Spartina species.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency.

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

Newly synthesized wheat allohexaploids display progenitor-dependent meiotic stability and aneuploidy but structural genomic additivity.

To understand key mechanisms leading to stabilized allopolyploid species, we characterized the meiotic behaviour of wheat allohexaploids in relation to structural and genetic changes. For that purpose, we analysed first generations of synthetic allohexaploids obtained through interspecific hybridization, followed by spontaneous chromosome doubling, between several genotypes of Triticum turgidum and Aegilops tauschii wheat species, donors of AB and D genomes, respectively. As expected for these Ph1 (Pairing homoeologous 1) gene-carrying allopolyploids, chromosome pairing at metaphase I of meiosis essentially occurs between homologous chromosomes. However, the different synthetic allohexaploids exhibited progenitor-dependent meiotic irregularities, such as incomplete homologous pairing, resulting in univalent formation and leading to aneuploidy in the subsequent generation. Stability of the synthetic allohexaploids was shown to depend on the considered genotypes of both AB and D genome progenitors, where few combinations compare to the natural wheat allohexaploid in terms of regularity of meiosis and euploidy. Aneuploidy represents the only structural change observed in these synthetic allohexaploids, as no apparent DNA sequence elimination or rearrangement was observed when analysing euploid plants with molecular markers, developed from expressed sequence tags (ESTs) as well as simple sequence repeat (SSR) and transposable element sequences.

Genome-wide gene expression changes in genetically stable synthetic and natural wheat allohexaploids.

*The present study aims to understand regulation of gene expression in synthetic and natural wheat (Triticum aestivum) allohexaploids, that combines the AB genome of Triticum turgidum and the D genome of Aegilops tauschii; and which we have recently characterized as genetically stable. *We conducted a comprehensive genome-wide analysis of gene expression that allowed characterization of the effect of variability of the D genome progenitor, the intergenerational stability as well as the comparison with natural wheat allohexaploid. We used the Affymetrix GeneChip Wheat Genome Array, on which 55 049 transcripts are represented. *Additive expression was shown to represent the majority of expression regulation in the synthetic allohexaploids, where expression for more than c. 93% of transcripts was equal to the mid-parent value measured from a mixture of parental RNA. This leaves c. 2000 (c. 7%) transcripts, in which expression was nonadditive. No global gene expression bias or dominance towards any of the progenitor genomes was observed whereas high intergenerational stability and low effect of the D genome progenitor variability were revealed. *Our study suggests that gene expression regulation in wheat allohexaploids is established early upon allohexaploidization and highly conserved over generations, as demonstrated by the high similarity of expression with natural wheat allohexaploids.

Quantitative resistance increases the durability of qualitative resistance to Leptosphaeria maculans in Brassica napus.

It has frequently been hypothesized that quantitative resistance increases the durability of qualitative (R-gene mediated) resistance but supporting experimental evidence is rare. To test this hypothesis, near-isogenic lines with/without the R-gene Rlm6 introduced into two Brassica napus cultivars differing in quantitative resistance to Leptosphaeria maculans were used in a 5-yr field experiment. Recurrent selection of natural fungal populations was done annually on each of the four plant genotypes, using crop residues from each genotype to inoculate separately the four series of field trials for five consecutive cropping seasons. Severity of phoma stem canker was measured on each genotype and frequencies of avirulence alleles in L. maculans populations were estimated. Recurrent selection of virulent isolates by Rlm6 in a susceptible background rendered the resistance ineffective by the third cropping season. By contrast, the resistance was still effective after 5 yr of selection by the genotype combining this gene with quantitative resistance. No significant variation in the performance of quantitative resistance alone was noted over the course of the experiment. We conclude that quantitative resistance can increase the durability of Rlm6. We recommend combining quantitative resistance with R-gene mediated resistance to enhance disease control and crop production.

Multicolored Fluorescent In Situ Hybridization to Assess Pairing Configurations at Metaphase I in Brassica Hybrids.

Genetic diversity can be introduced into polyploid crop species through meiotic recombination by exchanges between homologous or homoeologous chromosomes. Fluorescent in situ hybridization (FISH) enables the characterization of these homologous and homoeologous chromosome pairs during meiosis and identification of structural rearrangements during mitosis in metaphase I. In this chapter, we describe a protocol for the multicolored fluorescent labeling of chromosome spreads. This protocol allows the characterization of each A and C homoeologous subgenomes in a polyploid species using a genome-specific BAC combined with specific chromosome labeling BAC sequences.

Genome Size Variation and Comparative Genomics Reveal Intraspecific Diversity in Brassica rapa.

Traditionally, reference genomes in crop species rely on the assembly of one accession, thus occulting most of intraspecific diversity. However, rearrangements, gene duplications, and transposable element content may have a large impact on the genomic structure, which could generate new phenotypic traits. Comparing two Brassica rapa genomes recently sequenced and assembled using long-read technology and optical mapping, we investigated structural variants and repetitive content between the two accessions and genome size variation among a core collection. We explored the structural consequences of the presence of large repeated sequences in B. rapa 'Z1' genome vs. the B. rapa 'Chiifu' genome, using comparative genomics and cytogenetic approaches. First, we showed that large genomic variants on chromosomes A05, A06, A09, and A10 are due to large insertions and inversions when comparing B. rapa 'Z1' and B. rapa 'Chiifu' at the origin of important length differences in some chromosomes. For instance, lengths of 'Z1' and 'Chiifu' A06 chromosomes were estimated in silico to be 55 and 29 Mb, respectively. To validate these observations, we compared using fluorescent in situ hybridization (FISH) the two A06 chromosomes present in an F1 hybrid produced by crossing these two varieties. We confirmed a length difference of 17.6% between the A06 chromosomes of 'Z1' compared to 'Chiifu.' Alternatively, using a copy number variation approach, we were able to quantify the presence of a higher number of rDNA and gypsy elements in 'Z1' genome compared to 'Chiifu' on different chromosomes including A06. Using flow cytometry, the total genome size of 12 Brassica accessions corresponding to a B. rapa available core collection was estimated and revealed a genome size variation of up to 16% between these accessions as well as some shared inversions. This study revealed the contribution of long-read sequencing of new accessions belonging to different cultigroups of B. rapa and highlighted the potential impact of differential insertion of repeat elements and inversions of large genomic regions in genome size intraspecific variability.

Homoeologous recombination within bread wheat to develop novel combinations of HMW-GS genes: transfer of the Glu-A1 locus to chromosome 1D.

In an attempt to improve the bread-making quality within hexaploid wheat by elaborating novel high-molecular weight glutenin subunits (HMW-GS) combinations useful in wheat-breeding programmes, a 1A chromosome fragment carrying the Glu-A1 locus encoding the subunit Ax2*, was translocated to the long arm of chromosome 1D. The partially isohomoeoallelic line, designated RR239, had a meiotic behaviour as regular as cv. Courtot. It was characterised using genomic in situ hybridization and microsatellite markers as well as biochemical and proteomic approaches. The translocated 1D chromosome had an interstitial 1AL segment representing in average 30% of the recombinant arm length that was confirmed by molecular analysis. The genetic length of the removed segment in chromosome 1DL was estimated to be at least 51 cM, and that of the interstitial 1AL translocation to be at least 33 cM. Proteome analysis performed on total endosperm proteins revealed variation in amounts, 8 spots and 1 spot being up- and downregulated, respectively. Quantitative variations in HMW-GS were observed for the Glu-A1 (Ax2*) and Glu-B1 (Bx7 + By8) loci in response to duplication of the Glu-A1 locus.

Dynamics and differential proliferation of transposable elements during the evolution of the B and A genomes of wheat.

Transposable elements (TEs) constitute >80% of the wheat genome but their dynamics and contribution to size variation and evolution of wheat genomes (Triticum and Aegilops species) remain unexplored. In this study, 10 genomic regions have been sequenced from wheat chromosome 3B and used to constitute, along with all publicly available genomic sequences of wheat, 1.98 Mb of sequence (from 13 BAC clones) of the wheat B genome and 3.63 Mb of sequence (from 19 BAC clones) of the wheat A genome. Analysis of TE sequence proportions (as percentages), ratios of complete to truncated copies, and estimation of insertion dates of class I retrotransposons showed that specific types of TEs have undergone waves of differential proliferation in the B and A genomes of wheat. While both genomes show similar rates and relatively ancient proliferation periods for the Athila retrotransposons, the Copia retrotransposons proliferated more recently in the A genome whereas Gypsy retrotransposon proliferation is more recent in the B genome. It was possible to estimate for the first time the proliferation periods of the abundant CACTA class II DNA transposons, relative to that of the three main retrotransposon superfamilies. Proliferation of these TEs started prior to and overlapped with that of the Athila retrotransposons in both genomes. However, they also proliferated during the same periods as Gypsy and Copia retrotransposons in the A genome, but not in the B genome. As estimated from their insertion dates and confirmed by PCR-based tracing analysis, the majority of differential proliferation of TEs in B and A genomes of wheat (87 and 83%, respectively), leading to rapid sequence divergence, occurred prior to the allotetraploidization event that brought them together in Triticum turgidum and Triticum aestivum, <0.5 million years ago. More importantly, the allotetraploidization event appears to have neither enhanced nor repressed retrotranspositions. We discuss the apparent proliferation of TEs as resulting from their insertion, removal, and/or combinations of both evolutionary forces.