Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium.

The "genomic shock" hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession "Bhyb26." We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species' similarity in transposable element load may account for the subtlety of the observed genome dominance.

Evolutionary Analysis of the YABBY Gene Family in Brassicaceae.

The YABBY gene family is one of the plant transcription factors present in all seed plants. The family members were extensively studied in various plants and shown to play important roles in plant growth and development, such as the polarity establishment in lateral organs, the formation and development of leaves and flowers, and the response to internal plant hormone and external environmental stress signals. In this study, a total of 364 YABBY genes were identified from 37 Brassicaceae genomes, of which 15 were incomplete due to sequence gaps, and nine were imperfect (missing C2C2 zinc-finger or YABBY domain) due to sequence mutations. Phylogenetic analyses resolved these YABBY genes into six compact clades except for a YAB3-like gene identified in Aethionema arabicum. Seventeen Brassicaceae species each contained a complete set of six basic YABBY genes (i.e., 1 FIL, 1 YAB2, 1 YAB3, 1 YAB5, 1 INO and 1 CRC), while 20 others each contained a variable number of YABBY genes (5-25) caused mainly by whole-genome duplication/triplication followed by gene losses, and occasionally by tandem duplications. The fate of duplicate YABBY genes changed considerably according to plant species, as well as to YABBY gene type. These YABBY genes were shown to be syntenically conserved across most of the Brassicaceae species, but their functions might be considerably diverged between species, as well as between paralogous copies, as demonstrated by the promoter and expression analysis of YABBY genes in two Brassica species (B. rapa and B. oleracea). Our study provides valuable insights for understanding the evolutionary story of YABBY genes in Brassicaceae and for further functional characterization of each YABBY gene across the Brassicaceae species.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids.

Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.

A high-contiguity Brassica nigra genome localizes active centromeres and defines the ancestral Brassica genome.

It is only recently, with the advent of long-read sequencing technologies, that we are beginning to uncover previously uncharted regions of complex and inherently recursive plant genomes. To comprehensively study and exploit the genome of the neglected oilseed Brassica nigra, we generated two high-quality nanopore de novo genome assemblies. The N50 contig lengths for the two assemblies were 17.1 Mb (12 contigs), one of the best among 324 sequenced plant genomes, and 0.29 Mb (424 contigs), respectively, reflecting recent improvements in the technology. Comparison with a de novo short-read assembly corroborated genome integrity and quantified sequence-related error rates (0.2%). The contiguity and coverage allowed unprecedented access to low-complexity regions of the genome. Pericentromeric regions and coincidence of hypomethylation enabled localization of active centromeres and identified centromere-associated ALE family retro-elements that appear to have proliferated through relatively recent nested transposition events (<1 Ma). Genomic distances calculated based on synteny relationships were used to define a post-triplication Brassica-specific ancestral genome, and to calculate the extensive rearrangements that define the evolutionary distance separating B. nigra from its diploid relatives.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors.

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

Whole-Genome Resequencing of a Worldwide Collection of Rapeseed Accessions Reveals the Genetic Basis of Ecotype Divergence.

Rapeseed (Brassica napus), an important oilseed crop, has adapted to diverse climate zones and latitudes by forming three main ecotype groups, namely winter, semi-winter, and spring types. However, genetic variations underlying the divergence of these ecotypes are largely unknown. Here, we report the global pattern of genetic polymorphisms in rapeseed determined by resequencing a worldwide collection of 991 germplasm accessions. A total of 5.56 and 5.53 million single-nucleotide polymorphisms (SNPs) as well as 1.86 and 1.92 million InDels were identified by mapping reads to the reference genomes of "Darmor-bzh" and "Tapidor," respectively. We generated a map of allelic drift paths that shows splits and mixtures of the main populations, and revealed an asymmetric evolution of the two subgenomes of B. napus by calculating the genetic diversity and linkage disequilibrium parameters. Selective-sweep analysis revealed genetic changes in genes orthologous to those regulating various aspects of plant development and response to stresses. A genome-wide association study identified SNPs in the promoter regions of FLOWERING LOCUS T and FLOWERING LOCUS C orthologs that corresponded to the different rapeseed ecotype groups. Our study provides important insights into the genomic footprints of rapeseed evolution and flowering-time divergence among three ecotype groups, and will facilitate screening of molecular markers for accelerating rapeseed breeding.

Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family.

Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus.

Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype.

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (Ar ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with Ar , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

Assembly and comparison of two closely related Brassica napus genomes.

As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor-bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor-bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.

Genome-wide identification, evolution and expression analysis of RING finger protein genes in Brassica rapa.

More and more RING finger genes were found to be implicated in various important biological processes. In the present study, a total of 731 RING domains in 715 predicted proteins were identified in Brassica rapa genome (AA, 2n = 20), which were further divided into eight types: RING-H2 (371), RING-HCa (215), RING-HCb (47), RING-v (44), RING-C2 (38), RING-D (10), RING-S/T (5) and RING-G (1). The 715 RING finger proteins were further classified into 51 groups according to the presence of additional domains. 700 RING finger protein genes were mapped to the 10 chromosomes of B. rapa with a range of 47 to 111 genes for each chromosome. 667 RING finger protein genes were expressed in at least one of the six tissues examined, indicating their involvement in various physiological and developmental processes in B. rapa. Hierarchical clustering analysis of RNA-seq data divided them into seven major groups, one of which includes 231 members preferentially expressed in leaf, and constitutes then a panel of gene candidates for studying the genetic and molecular mechanisms of leafy head traits in Brassica crops. Our results lay the foundation for further studies on the classification, evolution and putative functions of RING finger protein genes in Brassica species.

Surviving a Genome Collision: Genomic Signatures of Allopolyploidization in the Recent Crop Species Brassica napus.

Polyploidization has played a major role in crop plant evolution, leading to advantageous traits that have been selected by humans. Here, we describe restructuring patterns in the genome of L., a recent allopolyploid species. Widespread segmental deletions, duplications, and homeologous chromosome exchanges were identified in diverse genome sequences from 32 natural and 20 synthetic accessions, indicating that homeologous exchanges are a major driver of postpolyploidization genome diversification. Breakpoints of genomic rearrangements are rich in microsatellite sequences that are known to interact with the meiotic recombination machinery. In both synthetic and natural , a subgenome bias was observed toward exchanges replacing larger chromosome segments from the C-subgenome by their smaller, homeologous A-subgenome segments, driving postpolyploidization genome size reduction. Selection in natural favored segmental deletions involving genes associated with immunity, reproduction, and adaptation. Deletions affecting mismatch repair system genes, which are assumed to control homeologous recombination, were also found to be under selection. Structural exchanges between homeologous subgenomes appear to be a major source of novel genetic diversity in de novo allopolyploids. Documenting the consequences of genomic collision by genomic resequencing gives insights into the adaptive processes accompanying allopolyploidization.

Recreating Stable Brachypodium hybridum Allotetraploids by Uniting the Divergent Genomes of B. distachyon and B. stacei.

Brachypodium hybridum (2n = 30) is a natural allopolyploid with highly divergent sub-genomes derived from two extant diploid species, B. distachyon (2n = 10) and B. stacei (2n = 20) that differ in chromosome evolution and number. We created synthetic B. hybridum allotetraploids by hybridizing various lines of B. distachyon and B. stacei. The initial amphihaploid F1 interspecific hybrids were obtained at low frequencies when B. distachyon was used as the maternal parent (0.15% or 0.245% depending on the line used) and were sterile. No hybrids were obtained from reciprocal crosses or when autotetraploids of the parental species were crossed. Colchicine treatment was used to double the genome of the F1 amphihaploid lines leading to allotetraploids. The genome-doubled F1 plants produced a few S1 (first selfed generation) seeds after self-pollination. S1 plants from one parental combination (Bd3-1×Bsta5) were fertile and gave rise to further generations whereas those of another parental combination (Bd21×ABR114) were sterile, illustrating the importance of the parental lineages crossed. The synthetic allotetraploids were stable and resembled the natural B. hybridum at the phenotypic, cytogenetic and genomic levels. The successful creation of synthetic B. hybridum offers the possibility to study changes in genome structure and regulation at the earliest stages of allopolyploid formation in comparison with the parental species and natural B. hybridum.

Co-linearity and divergence of the A subgenome of Brassica juncea compared with other Brassica species carrying different A subgenomes.

BACKGROUND: There are three basic Brassica genomes (A, B, and C) and three parallel sets of subgenomes distinguished in the diploid Brassica (i.e.: B. rapa, A(r)A(r); B. nigra, B(ni)B(ni); B. oleracea, C(o)C(o)) and the derived allotetraploid species (i.e.: B. juncea, A(j)A(j)B(j)B(j); B. napus, A(n)A(n)C(n)C(n); B. carinata, B(c)B(c)C(c)C(c)). To understand subgenome differentiation in B. juncea in comparison to other A genome-carrying Brassica species (B. rapa and B. napus), we constructed a dense genetic linkage map of B. juncea, and conducted population genetic analysis on diverse lines of the three A-genome carrying Brassica species using a genotyping-by-sequencing approach (DArT-seq). RESULTS: A dense genetic linkage map of B. juncea was constructed using an F2 population derived from Sichuan Yellow/Purple Mustard. The map included 3329 DArT-seq markers on 18 linkage groups and covered 1579 cM with an average density of two markers per cM. Based on this map and the alignment of the marker sequences with the physical genome of Arabidopsis thaliana, we observed strong co-linearity of the ancestral blocks among the different A subgenomes but also considerable block variation. Comparative analyses at the level of genome sequences of B. rapa and B. napus, and marker sequence anchored on the genetic map of B. juncea, revealed a total of 30 potential inversion events across large segments and 20 potential translocation events among the three A subgenomes. Population genetic analysis on 26 accessions of the three A genome-carrying Brassica species showed that the highest genetic distance were estimated when comparing A(j)-A(n) than between A(n)-A(r) and A(j)-A(r) subgenome pairs. CONCLUSIONS: The development of the dense genetic linkage map of B. juncea with informative DArT-seq marker sequences and availability of the reference sequences of the A(r), and A(n)C(n) genomes allowed us to compare the A subgenome structure of B. juncea (A(j)) . Our results suggest that strong co-linearity exists among the three A Brassica genomes (A(r), A(n) and A(j)) but with apparent subgenomic variation. Population genetic analysis on three A-genome carrying Brassica species support the idea that B. juncea has distinct genomic diversity, and/or evolved from a different A genome progenitor of B. napus.

Identification, evolution, and expression partitioning of miRNAs in allopolyploid Brassica napus.

The recently published genome of Brassica napus offers for the first time the opportunity to gain insights into the genomic organization and the evolution of miRNAs in oilseed rape. In this study, 12 small RNA libraries from two B. napus cultivars (Tapidor and Ningyou7) and their four double-haploid lines were sequenced, employing the newly sequenced B. napus genome, together with genomes of its progenitors Brassica rapa and Brassica oleracea. A total of 645 miRNAs including 280 conserved and 365 novel miRNAs were identified. Comparative analysis revealed a high level of genomic conservation of MIRNAs (75.9%) between the subgenomes of B. napus and its two progenitors' genomes, and MIRNA lost/gain events (133) occurred in B. napus after its speciation. Furthermore, significant partitioning of miRNA expressions between the two subgenomes in B. napus was detected. The data of degradome sequencing, miRNA-mediated cleavage, and expression analyses support specific interactions between miRNAs and their targets in the modulation of diverse physiological processes in roots and leaves, as well as in biosynthesis of, for example, glucosinolates and lipids in oilseed rape. These data provide a first genome-wide view on the origin, evolution, and genomic organization of B. napus MIRNAs.

Construction of Brassica A and C genome-based ordered pan-transcriptomes for use in rapeseed genomic research.

This data article reports the establishment of the first pan-transcriptome resources for the Brassica A and C genomes. These were developed using existing coding DNA sequence (CDS) gene models from the now-published Brassica oleracea TO1000 and Brassica napus Darmor-bzh genome sequence assemblies representing the chromosomes of these species, along with preliminary CDS models from an updated Brassica rapa Chiifu genome sequence assembly. The B. rapa genome sequence scaffolds required splitting and re-ordering to match the expected genome organisation based on a high density SNP linkage map, but the B. oleracea assembly was used unchanged. The resulting B. rapa (A genome) pseudomolecules contained 47,656 ordered CDS models and the B. oleracea (C genome) pseudomolecules contained 54,766 ordered CDS models. Interpolation of B. napus CDS models not already represented by orthologues resulted in 52,790 and 63,308 ordered CDS models in the A and C pan-transcriptomes, an increase of 13,676 overall. Comparison of the organisation of this resource with publicly available genome sequences for B. napus showed excellent consistency for the B. napus Darmor-bzh resource, but more breakdown of collinearity for the B. napus ZS11 resource. CDS datasets comprising the pan-transcriptomes are available with this article (B. rapa) or from public repositories (B. oleracea and B. napus).

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome.

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.

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

BACKGROUND: Brassica 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. RESULTS: We 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. CONCLUSIONS: Differential 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.

Update on the genomics and basic biology of Brachypodium: International Brachypodium Initiative (IBI).

The scientific presentations at the First International Brachypodium Conference (abstracts available at http://www.brachy2013.unimore.it) are evidence of the widespread adoption of Brachypodium distachyon as a model system. Furthermore, the wide range of topics presented (genome evolution, roots, abiotic and biotic stress, comparative genomics, natural diversity, and cell walls) demonstrates that the Brachypodium research community has achieved a critical mass of tools and has transitioned from resource development to addressing biological questions, particularly those unique to grasses.

Meiotic gene evolution: can you teach a new dog new tricks?

Meiosis, the basis of sex, evolved through iterative gene duplications. To understand whether subsequent duplications have further enriched the core meiotic "tool-kit," we investigated the fate of meiotic gene duplicates following whole genome duplication (WGD), a common occurrence in eukaryotes. We show that meiotic genes return to a single copy more rapidly than genome-wide average in angiosperms, one of the lineages in which WGD is most vividly exemplified. The rate at which duplicates are lost decreases through time, a tendency that is also observed genome-wide and may thus prove to be a general trend post-WGD. The sharpest decline is observed for the subset of genes mediating meiotic recombination; however, we found no evidence that the presence of these duplicates is counterselected in two recent polyploid crops selected for fertility. We therefore propose that their loss is passive, highlighting how quickly WGDs are resolved in the absence of selective duplicate retention.