Harry Belcram

Early allopolyploid evolution in the post-neolithic Brassica napus oilseed genome

The genomic origins of rape oilseed Many domesticated plants arose through the meeting of multiple genomes through hybridization and genome doubling, known as polyploidy. Chalhoub et al. sequenced the polyploid genome of Brassica napus, which originated from a recent combination of two distinct genomes approximately 7500 years ago and gave rise to the crops of rape oilseed (canola), kale, and rutabaga. B. napus has undergone multiple events affecting differently sized genetic regions where a gene from one progenitor species has been converted to the copy from a second progenitor species. Some of these gene conversion events appear to have been selected by humans as part of the process of domestication and crop improvement. Science, this issue p. 950 The polyploid genome of oilseed rape exhibits evolution through homologous gene conversion. Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes

Polyploidization has provided much genetic variation for plant adaptive evolution, but the mechanisms by which the molecular evolution of polyploid genomes establishes genetic architecture underlying species differentiation are unclear. Brassica is an ideal model to increase knowledge of polyploid evolution. Here we describe a draft genome sequence of Brassica oleracea, comparing it with that of its sister species B. rapa to reveal numerous chromosome rearrangements and asymmetrical gene loss in duplicated genomic blocks, asymmetrical amplification of transposable elements, differential gene co-retention for specific pathways and variation in gene expression, including alternative splicing, among a large number of paralogous and orthologous genes. Genes related to the production of anticancer phytochemicals and morphological variations illustrate consequences of genome duplication and gene divergence, imparting biochemical and morphological variation to B. oleracea. This study provides insights into Brassica genome evolution and will underpin research into the many important crops in this genus.

Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea

BackgroundBrassica oleracea is a valuable vegetable species that has contributed to human health and nutrition for hundreds of years and comprises multiple distinct cultivar groups with diverse morphological and phytochemical attributes. In addition to this phenotypic wealth, B. oleracea offers unique insights into polyploid evolution, as it results from multiple ancestral polyploidy events and a final Brassiceae-specific triplication event. Further, B. oleracea represents one of the diploid genomes that formed the economically important allopolyploid oilseed, Brassica napus. A deeper understanding of B. oleracea genome architecture provides a foundation for crop improvement strategies throughout the Brassica genus.ResultsWe generate an assembly representing 75% of the predicted B. oleracea genome using a hybrid Illumina/Roche 454 approach. Two dense genetic maps are generated to anchor almost 92% of the assembled scaffolds to nine pseudo-chromosomes. Over 50,000 genes are annotated and 40% of the genome predicted to be repetitive, thus contributing to the increased genome size of B. oleracea compared to its close relative B. rapa. A snapshot of both the leaf transcriptome and methylome allows comparisons to be made across the triplicated sub-genomes, which resulted from the most recent Brassiceae-specific polyploidy event.ConclusionsDifferential expression of the triplicated syntelogs and cytosine methylation levels across the sub-genomes suggest residual marks of the genome dominance that led to the current genome architecture. Although cytosine methylation does not correlate with individual gene dominance, the independent methylation patterns of triplicated copies suggest epigenetic mechanisms play a role in the functional diversification of duplicate genes.

The first meiosis of resynthesized Brassica napus, a genome blender

Polyploidy promotes the restructuring of merged genomes within initial generations of resynthesized Brassica napus, possibly caused by homoeologous recombination at meiosis. However, little is known about the impact of the first confrontation of two genomes at the first meiosis which could lead to genome exchanges in progeny. Here, we assessed the role of the first meiosis in the genome instability of synthetic B. napus. We used three different newly resynthesized B. napus plants and established meiotic pairing frequencies for the A and C genomes. We genotyped the three corresponding progenies in a cross to a natural B. napus on the two homoeologous A1 and C1 chromosomes. Pairing at meiosis in a set of progenies with various rearrangements was scored. Here, we confirmed that the very first meiosis of resynthesized plants of B. napus acts as a genome blender, with many of the meiotic-driven genetic changes transmitted to the progenies, in proportions that depend significantly on the cytoplasm background inherited from the progenitors. We conclude that the first meiosis generates rearrangements on both genomes and promotes subsequent restructuring in further generations. Our study advances the knowledge on the timing of genetic changes and the mechanisms that may bias their transmission.

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.

Types and Rates of Sequence Evolution at the High-Molecular-Weight Glutenin Locus in Hexaploid Wheat and Its Ancestral Genomes

The Glu-1 locus, encoding the high-molecular-weight glutenin protein subunits, controls bread-making quality in hexaploid wheat (Triticum aestivum) and represents a recently evolved region unique to Triticeae genomes. To understand the molecular evolution of this locus region, three orthologous Glu-1 regions from the three subgenomes of a single hexaploid wheat species were sequenced, totaling 729 kb of sequence. Comparing each Glu-1 region with its corresponding homologous region from the D genome of diploid wheat, Aegilops tauschii, and the A and B genomes of tetraploid wheat, Triticum turgidum, revealed that, in addition to the conservation of microsynteny in the genic regions, sequences in the intergenic regions, composed of blocks of nested retroelements, are also generally conserved, although a few nonshared retroelements that differentiate the homologous Glu-1 regions were detected in each pair of the A and D genomes. Analysis of the indel frequency and the rate of nucleotide substitution, which represent the most frequent types of sequence changes in the Glu-1 regions, demonstrated that the two A genomes are significantly more divergent than the two B genomes, further supporting the hypothesis that hexaploid wheat may have more than one tetraploid ancestor.

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.

Duplication and partitioning in evolution and function of homoeologous Q loci governing domestication characters in polyploid wheat

The Q gene encodes an AP2-like transcription factor that played an important role in domestication of polyploid wheat. The chromosome 5A Q alleles (5AQ and 5Aq) have been well studied, but much less is known about the q alleles on wheat homoeologous chromosomes 5B (5Bq) and 5D (5Dq). We investigated the organization, evolution, and function of the Q/q homoeoalleles in hexaploid wheat (Triticum aestivum L.). Q/q gene sequences are highly conserved within and among the A, B, and D genomes of hexaploid wheat, the A and B genomes of tetraploid wheat, and the A, S, and D genomes of the diploid progenitors, but the intergenic regions of the Q/q locus are highly divergent among homoeologous genomes. Duplication of the q gene 5.8 Mya was likely followed by selective loss of one of the copies from the A genome progenitor and the other copy from the B, D, and S genomes. A recent V329-to-I mutation in the A lineage is correlated with the Q phenotype. The 5Bq homoeoalleles became a pseudogene after allotetraploidization. Expression analysis indicated that the homoeoalleles are coregulated in a complex manner. Combined phenotypic and expression analysis indicated that, whereas 5AQ plays a major role in conferring domestication-related traits, 5Dq contributes directly and 5Bq indirectly to suppression of the speltoid phenotype. The evolution of the Q/q loci in polyploid wheat resulted in the hyperfunctionalization of 5AQ, pseudogenization of 5Bq, and subfunctionalization of 5Dq, all contributing to the domestication traits.

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.

New insights into the origin of the B genome of hexaploid wheat: Evolutionary relationships at the SPA genomic region with the S genome of the diploid relative Aegilops speltoides

BackgroundSeveral studies suggested that the diploid ancestor of the B genome of tetraploid and hexaploid wheat species belongs to the Sitopsis section, having Aegilops speltoides (SS, 2n = 14) as the closest identified relative. However molecular relationships based on genomic sequence comparison, including both coding and non-coding DNA, have never been investigated. In an attempt to clarify these relationships, we compared, in this study, sequences of the Storage Protein Activator (SPA) locus region of the S genome of Ae. speltoides (2n = 14) to that of the A, B and D genomes co-resident in the hexaploid wheat species (Triticum aestivum, AABBDD, 2n = 42).ResultsFour BAC clones, spanning the SPA locus of respectively the A, B, D and S genomes, were isolated and sequenced. Orthologous genomic regions were identified as delimited by shared non-transposable elements and non-coding sequences surrounding the SPA gene and correspond to 35 268, 22 739, 43 397 and 53 919 bp for the A, B, D and S genomes, respectively. Sequence length discrepancies within and outside the SPA orthologous regions are the result of non-shared transposable elements (TE) insertions, all of which inserted after the progenitors of the four genomes divergence.ConclusionOn the basis of conserved sequence length as well as identity of the shared non-TE regions and the SPA coding sequence, Ae speltoides appears to be more evolutionary related to the B genome of T. aestivum than the A and D genomes. However, the differential insertions of TEs, none of which are conserved between the two genomes led to the conclusion that the S genome of Ae. speltoides has diverged very early from the progenitor of the B genome which remains to be identified.