Hybridization-facilitated genome merger and repeated chromosome fusion after 8 million years

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

Celebrating Mendel, McClintock, and Darlington: On end-to-end chromosome fusions and nested chromosome fusions

The evolution of eukaryotic genomes is accompanied by fluctuations in chromosome number, reflecting cycles of chromosome number increase (polyploidy and centric fissions) and decrease (chromosome fusions). Although all chromosome fusions result from DNA recombination between two or more nonhomologous chromosomes, several mechanisms of descending dysploidy are exploited by eukaryotes to reduce their chromosome number. Genome sequencing and comparative genomics have accelerated the identification of inter-genome chromosome collinearity and gross chromosomal rearrangements and have shown that end-to-end chromosome fusions (EEFs) and nested chromosome fusions (NCFs) may have played a more important role in the evolution of eukaryotic karyotypes than previously thought. The present review aims to summarize the limited knowledge on the origin, frequency, and evolutionary implications of EEF and NCF events in eukaryotes and especially in land plants. The interactions between nonhomologous chromosomes in interphase nuclei and chromosome (mis)pairing during meiosis are examined for their potential importance in the origin of EEFs and NCFs. The remaining open questions that need to be addressed are discussed.

Molecular and Cytogenetic Analysis of rDNA Evolution in Crepis Sensu Lato

Although Crepis was the first model plant group in which chromosomal changes were considered to play an important role in speciation, their chromosome structure and evolution have been barely investigated using molecular cytogenetic methods. The aim of the study was to provide a better understanding of the patterns and directions of Crepis chromosome evolution, using comparative analyses of rDNA loci number and localisation. The chromosome base number and chromosomal organisation of 5S and 35S rDNA loci were analysed in the phylogenetic background for 39 species of Crepis, which represent the evolutionary lineages of Crepis sensu stricto and Lagoseris, including Lapsana communis. The phylogenetic relationships among all the species were inferred from nrITS and newly obtained 5S rDNA NTS sequences. Despite high variations in rDNA loci chromosomal organisation, most species had a chromosome with both rDNA loci within the same (usually short) chromosomal arm. The comparative analyses revealed several independent rDNA loci number gains and loci repositioning that accompanied diversification and speciation in Crepis. Some of the changes in rDNA loci patterns were reconstructed for the same evolutionary lineages as descending dysploidy.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

Genome structure and evolution in the cruciferous tribe Thlaspideae (Brassicaceae)

Whole-genome duplications (WGDs) and chromosome rearrangements (CRs) play the key role in driving the diversification and evolution of plant lineages. Although the direct link between WGDs and plant diversification is well documented, relatively few studies focus on the evolutionary significance of CRs. The cruciferous tribe Thlaspideae represents an ideal model system to address the role of large-scale chromosome alterations in genome evolution, as most Thlaspideae species share the same diploid chromosome number (2n = 2x = 14). Here we constructed the genome structure in 12 Thlaspideae species, including field pennycress (Thlaspi arvense) and garlic mustard (Alliaria petiolata). We detected and precisely characterized genus- and species-specific CRs, mostly pericentric inversions, as the main genome-diversifying drivers in the tribe. We reconstructed the structure of seven chromosomes of an ancestral Thlaspideae genome, identified evolutionary stable chromosomes versus chromosomes prone to CRs, estimated the rate of CRs, and uncovered an allohexaploid origin of garlic mustard from diploid taxa closely related to A. petiolata and Parlatoria cakiloidea. Furthermore, we performed detailed bioinformatic analysis of the Thlaspideae repeatomes, and identified repetitive elements applicable as unique species- and genus-specific barcodes and chromosome landmarks. This study deepens our general understanding of the evolutionary role of CRs, particularly pericentric inversions, in plant genome diversification, and provides a robust base for follow-up whole-genome sequencing efforts.

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation

Plants’ adaptations to and divergence in arid deserts have long fascinated scientists and the general public. Here, we present a genomic analysis of two congeneric desert plant species that clarifies their evolutionary history and shows that their common ancestor arose from a hybrid polyploidization, which provided genomic foundations for their survival in deserts. The whole-genome duplication was followed by translocation-based rearrangements of the ancestral chromosomes. Rapid evolution of genes in these reshuffled chromosomes contributed greatly to the divergences of the two species in desert microhabitats during which gene flow was continuous. Our results provide insights into plant adaptation in the arid deserts and highlight the significance of polyploidy-driven chromosomal structural variations in species divergence.

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

Linked by Ancestral Bonds: Multiple Whole-Genome Duplications and Reticulate Evolution in a Brassicaceae Tribe

Pervasive hybridization and whole-genome duplications (WGDs) influenced genome evolution in several eukaryotic lineages. Although frequent and recurrent hybridizations may result in reticulate phylogenies, the evolutionary events underlying these reticulations, including detailed structure of the ancestral diploid and polyploid genomes, were only rarely reconstructed. Here, we elucidate the complex genomic history of a monophyletic clade from the mustard family (Brassicaceae), showing contentious relationships to the early-diverging clades of this model plant family. Genome evolution in the crucifer tribe Biscutelleae (∼60 species, 5 genera) was dominated by pervasive hybridizations and subsequent genome duplications. Diversification of an ancestral diploid genome into several divergent but crossable genomes was followed by hybridizations between these genomes. Whereas a single genus (Megadenia) remained diploid, the four remaining genera originated by allopolyploidy (Biscutella, Lunaria, Ricotia) or autopolyploidy (Heldreichia). The contentious relationships among the Biscutelleae genera, and between the tribe and other early diverged crucifer lineages, are best explained by close genomic relatedness among the recurrently hybridizing ancestral genomes. By using complementary cytogenomics and phylogenomics approaches, we demonstrate that the origin of a monophyletic plant clade can be more complex than a parsimonious assumption of a single WGD spurring postpolyploid cladogenesis. Instead, recurrent hybridization among the same and/or closely related parental genomes may phylogenetically interlink diploid and polyploid genomes despite the incidence of multiple independent WGDs. Our results provide new insights into evolution of early-diverging Brassicaceae lineages and elucidate challenges in resolving the contentious relationships within and between land plant lineages with pervasive hybridization and WGDs.

A chromosome-scale reference genome of Lobularia maritima, an ornamental plant with high stress tolerance

Lobularia maritima (L.) Desv. is an ornamental plant cultivated across the world. It belongs to the family Brassicaceae and can tolerate dry, poor and contaminated habitats. Here, we present a chromosome-scale, high-quality genome assembly of L. maritima based on integrated approaches combining Illumina short reads and Hi-C chromosome conformation data. The genome was assembled into 12 pseudochromosomes with a 197.70 Mb length, and it includes 25,813 protein-coding genes. Approximately 41.94% of the genome consists of repetitive sequences, with abundant long terminal repeat transposable elements. Comparative genomic analysis confirmed that L. maritima underwent a species-specific whole-genome duplication (WGD) event ~22.99 million years ago. We identified ~1900 species-specific genes, 25 expanded gene families, and 50 positively selected genes in L. maritima. Functional annotations of these genes indicated that they are mainly related to stress tolerance. These results provide new insights into the stress tolerance of L. maritima, and this genomic resource will be valuable for further genetic improvement of this important ornamental plant.

The chromosome-level genome sequence and karyotypic evolution of Megadenia pygmaea (Brassicaceae)

Karyotypic changes in chromosome number and structure are drivers in the divergent evolution of diverse plant species and lineages. This study aimed to reveal the origins of the unique karyotype (2n = 12) and phylogenetic relationships of the genus Megadenia (Brassicaceae). A high-quality chromosome-scale genome was assembled for Megadenia pygmaea using Nanopore long reads and high-throughput chromosome conformation capture (Hi-C). The assembled genome is 215.2 Mb and is anchored on six pseudochromosomes. We annotated a total of 25,607 high-confidence protein-coding genes and corroborated the phylogenetic affinity of Megadenia with the Brassicaceae expanded lineage II, containing numerous agricultural crops. We dated the divergence of Megadenia from its closest relatives to 27.04 (19.11-36.60) million years ago. A reconstruction of the chromosomal composition of the species was performed based on the de novo assembled genome and comparative chromosome painting analysis. The karyotype structure of M. pygmaea is very similar to the previously inferred proto-Calepineae karyotype (PCK; n = 7) of the lineage II. However, an end-to-end translocation between two ancestral chromosomes reduced the chromosome number from n = 7 to n = 6 in Megadenia. Our reference genome provides fundamental information for karyotypic evolution and evolutionary study of this genus.

Multiple and independent rearrangements revealed by comparative cytogenetic mapping in the dysploid Leptostachyus group (Phaseolus L., Leguminosae)

Polyploidy and dysploidy have been reported as the main events in karyotype evolution of plants. In the genus Phaseolus L. (2n = 22), a small monophyletic group of three species, the Leptostachyus group, presents a dysploid karyotype with 2n = 20. It was shown in Phaseolus leptostachyus that the dysploidy was caused by a nested chromosome fusion (NCF) accompanied by several translocations, suggesting a high rate of karyotype evolution in the group. To verify if this karyotype restructuring was a single event or occurred progressively during the evolution of this group, we analysed P. macvaughii, sister to Phaseolus micranthus + P. leptostachyus. Twenty-four genomic clones of P. vulgaris previously mapped on P. leptostachyus, in addition to the 5S and 35S rDNA probes, were used for fluorescence in situ hybridization. Only a single rearrangement was common to the two species: the nested chromosome fusion (NCF) involving chromosomes 10 and 11. The translocation of chromosome 2 is not the same found in P. leptostachyus, and pericentric inversions in chromosomed 3 and 4 were exclusive of P. macvaughii. The other rearrangements observed in P. leptostachyus were not shared with this species, suggesting that they occurred after the separation of these lineages. The presence of private rearrangements indicates a progressive accumulation of karyotype changes in the Leptostachyus group instead of an instant genome-wide repatterning.

Phylogeny of Euclidieae (Brassicaceae) based on plastome and nuclear ribosomal DNA data

Euclidieae, a morphologically diverse tribe in the family Brassicaceae (Cruciferae), consists of 29 genera and more than 150 species distributed mainly in Asia. Prior phylogenetic analyses on Euclidieae are inadequate. In this study, sequence data from the plastid genome and nuclear ribosomal DNA of 72 species in 27 genera of Euclidieae were used to infer the inter- and intra-generic relationships within. The well-resolved and strongly supported plastome phylogenies revealed that Euclidieae could be divided into five clades. Both Cymatocarpus and Neotorularia are polyphyletic in nuclear and plastome phylogenies. Besides, the conflicts of systematic positions of three species of Braya and two species of Solms-laubachia s.l. indicated that hybridization and or introgression might have happened during the evolutionary history of the tribe. Results from divergence-time analyses suggested an early Miocene origin of Euclidieae, and it probably originated from the Central Asia, Pamir Plateau and West Himalaya. In addition, multiple ndh genes loss and pseudogenization were detected in eight species based on comparative genomic study.

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution

Plant genome evolution can be very complex and challenging to describe, even within a genus. Mechanisms that underlie genome variation are complex and can include whole-genome duplications, gene duplication and/or loss, and, importantly, multiple chromosomal rearrangements. Lupins (Lupinus) diverged from other legumes approximately 60 mya. In contrast to New World lupins, Old World lupins show high variability not only for chromosome numbers (2n = 32–52), but also for the basic chromosome number (x = 5–9, 13) and genome size. The evolutionary basis that underlies the karyotype evolution in lupins remains unknown, as it has so far been impossible to identify individual chromosomes. To shed light on chromosome changes and evolution, we used comparative chromosome mapping among 11 Old World lupins, with Lupinus angustifolius as the reference species. We applied set of L. angustifolius-derived bacterial artificial chromosome clones for fluorescence in situ hybridization. We demonstrate that chromosome variations in the species analyzed might have arisen from multiple changes in chromosome structure and number. We hypothesize about lupin karyotype evolution through polyploidy and subsequent aneuploidy. Additionally, we have established a cytogenomic map of L. angustifolius along with chromosome markers that can be used for related species to further improve comparative studies of crops and wild lupins.