Fast diploidization in close mesopolyploid relatives of Arabidopsis

Chromosomal dynamics in Senna: comparative PLOP-FISH analysis of tandem repeats and flow cytometric nuclear genome size estimations

Introduction

Tandem repeats (TRs) occur abundantly in plant genomes. They play essential roles that affect genome organization and evolution by inducing or generating chromosomal rearrangements such as duplications, deletions, inversions, and translocations. These impact gene expression and chromosome structure and even contribute to the emergence of new species.

Method

We investigated the effects of TRs on speciation in Senna genus by performing a comparative analysis using fluorescence in situ hybridization (FISH) with S. tora-specific TR probes. We examined the chromosomal distribution of these TRs and compared the genome sizes of seven Senna species (estimated using flow cytometry) to better understand their evolutionary relationships.

Results

Two (StoTR03_159 and StoTR04_55) of the nine studied TRs were not detected in any of the seven Senna species, whereas the remaining seven were found in all or some species with patterns that were similar to or contrasted with those of S. tora. Of these studies species, only S. angulata showed significant genome rearrangements and dysploid karyotypes resembling those of S. tora. The genome sizes varied among these species and did not positively correlate with chromosome number. Notably, S. angulata had the fewest chromosomes (2n = 22) but a relatively large genome size.

Discussion

These findings reveal the dynamics of TRs and provide a cytogenetic depiction of chromosomal rearrangements during speciation in Senna. To further elucidate the dynamics of repeat sequences in Senna, future studies must include related species and extensive repeatomic studies, including those on transposable elements.

Interploidy gene flow via a 'pentaploid bridge' and ploidy reduction in Cystopteris fragilis fern complex (Cystopteridaceae: Polypodiales)

KEY MESSAGE

Our results indicate the existence of interploidy gene flow in Cystopteris fragilis, resulting in sexual triploid and diploid gametophytes from pentaploid parents. Similar evolutionary dynamics might operate in other fern complexes and need further investigation. Polyploidization and hybridization are a key evolutionary processes in ferns. Here, we outline an interploidy gene flow pathway operating in the polyploid Cystopteris fragilis complex. The conditions necessary for the existence of this pathway were tested. A total of 365 C. fragilis individuals were collected covering representatives of all three predominant ploidy levels (tetraploid, pentaploid, and hexaploid), cultivated, had their ploidy level estimated by flow cytometry, and their spores collected. The spores, as well as gametophytes and sporophytes established from them, were analysed by flow cytometry. Spore abortion rate was also estimated. In tetraploids, we observed the formation of unreduced (tetraploid) spores (ca 2%). Collected pentaploid individuals indicate ongoing hybridization between ploidy levels. Pentaploids formed up to 52% viable spores, ca 79% of them reduced, i.e. diploid and triploid. Reduced spores formed viable gametophytes, and, in the case of triploids, filial hexaploid sporophytes, showing evidence of sexual reproduction. Some tetraploid sporophytes reproduce apomictically (based on uniform ploidy of their metagenesis up to filial sporophytes). Triploid and diploid gametophytes from pentaploid parents are able to mate among themselves, or with "normal" reduced gametophytes from the sexual tetraploid sporophytes (the dominant ploidy level in the sporophytes in this populations), to produce tetraploid, pentaploid, and hexaploid sporophytes, allowing for geneflow from the pentaploids to both the tetraploid and hexaploid populations. Similar evolutionary dynamics might operate in other fern complexes and need further investigation.

The evolution of the hypotetraploid Catolobus pendulus genome - the poorly known sister species of Capsella

The establishment of Arabidopsis as the most important plant model has also brought other crucifer species into the spotlight of comparative research. While the genus Capsella has become a prominent crucifer model system, its closest relative has been overlooked. The unispecific genus Catolobus is native to temperate Eurasian woodlands, from eastern Europe to the Russian Far East. Here, we analyzed chromosome number, genome structure, intraspecific genetic variation, and habitat suitability of Catolobus pendulus throughout its range. Unexpectedly, all analyzed populations were hypotetraploid (2n = 30, ~330 Mb). Comparative cytogenomic analysis revealed that the Catolobus genome arose by a whole-genome duplication in a diploid genome resembling Ancestral Crucifer Karyotype (ACK, n = 8). In contrast to the much younger Capsella allotetraploid genomes, the presumably autotetraploid Catolobus genome (2n = 32) arose early after the Catolobus/Capsella divergence. Since its origin, the tetraploid Catolobus genome has undergone chromosomal rediploidization, including a reduction in chromosome number from 2n = 32 to 2n = 30. Diploidization occurred through end-to-end chromosome fusion and other chromosomal rearrangements affecting a total of six of 16 ancestral chromosomes. The hypotetraploid Catolobus cytotype expanded toward its present range, accompanied by some longitudinal genetic differentiation. The sister relationship between Catolobus and Capsella allows comparative studies of tetraploid genomes of contrasting ages and different degrees of genome diploidization.

An in-depth investigation of cryptic taxonomic diversity in the rare endemic mustard Draba maguirei

PREMISE

Previously published evidence suggests that Draba maguirei, a mustard endemic to a few localities in the Bear River, Wellsville, and Wasatch Mountains of northern Utah, may represent a cryptic species complex rather than a single species. Conservation concerns prompted an in-depth systematic study of this taxon and its putative relatives.

METHODS

Sampling most known populations of D. maguirei s.l. (D. maguirei var. maguirei and D. maguirei var. burkei), we integrate data from geography, ecology, morphology, cytogenetics and pollen, enzyme electrophoresis, and the phylogenetic analysis of nuclear internal transcribed spacer sequences to explore potential taxonomic diversity in the species complex.

RESULTS

Draba maguirei var. burkei is shown here to be a distinct species (D. burkei) most closely related to D. globosa, rather than to D. maguirei. Within D. maguirei s.s., the northern (high elevation) and southern (low elevation) population clusters are genetically isolated and morphologically distinguishable, leading to the recognition here of the southern taxon as D. maguirei subsp. stonei.

CONCLUSIONS

Our study reveals that plants traditionally assigned to D. maguirei comprise three genetically divergent lineages (D. burkei and two newly recognized subspecies of D. maguirei), each exhibiting a different chromosome number and occupying a discrete portion of the geographic range. Although previously overlooked and underappreciated taxonomically, the three taxa are morphologically recognizable based on the distribution and types of trichomes present on the leaves, stems, and fruit. Our clarification of the diversity and distribution of these taxa provides an improved framework for conservation efforts.

Polyploidy on islands - concerted evolution and gene loss amid chromosomal stasis

BACKGROUND AND AIMS

Polyploidy is an important process that often generates genomic diversity within lineages, but it can also cause changes that result in loss of genomic material. Island lineages, while often polyploid, typically show chromosomal stasis but have not been investigated in detail regarding smaller-scale gene loss. Our aim was to investigate post-polyploidization genome dynamics in a chromosomally stable lineage of Malvaceae endemic to New Zealand.

METHODS

We determined chromosome numbers and used fluorescence in situ hybridization to localize 18S and 5S rDNA. Gene sequencing of 18S rDNA, the internal transcribed spacers (ITS) with intervening 5.8S rDNA, and a low-copy nuclear gene, GBSSI-1, was undertaken to determine if gene loss occurred in the New Zealand lineage following polyploidy.

KEY RESULTS

The chromosome number for all species investigated was 2n = 42, with the first published report for the monotypic Australian genus Asterotrichion. The five species investigated all had two 5S rDNA signals localized interstitially on the long arm of one of the largest chromosome pairs. All species, except Plagianthus regius, had two 18S rDNA signals localized proximally on the short arm of one of the smallest chromosome pairs. Plagianthus regius had two additional 18S rDNA signals on a separate chromosome, giving a total of four. Sequencing of nuclear ribosomal 18S rDNA and the ITS cistron indicated loss of historical ribosomal repeats. Phylogenetic analysis of a low-copy nuclear gene, GBSSI-1, indicated that some lineages maintained three copies of the locus, while others have lost one or two copies.

CONCLUSIONS

Although island endemic lineages show chromosomal stasis, with no additional changes in chromosome number, they may undergo smaller-scale processes of gene loss and concerted evolution ultimately leading to further genome restructuring and downsizing.

The grapevine LysM receptor-like kinase VvLYK5-1 recognizes chitin oligomers through its association with VvLYK1-1

The establishment of defense reactions to protect plants against pathogens requires the recognition of invasion patterns (IPs), mainly detected by plasma membrane-bound pattern recognition receptors (PRRs). Some IPs, also termed elicitors, are used in several biocontrol products that are gradually being developed to reduce the use of chemicals in agriculture. Chitin, the major component of fungal cell walls, as well as its deacetylated derivative, chitosan, are two elicitors known to activate plant defense responses. However, recognition of chitooligosaccharides (COS) in Vitis vinifera is still poorly understood, hampering the improvement and generalization of protection tools for this important crop. In contrast, COS perception in the model plant Arabidopsis thaliana is well described and mainly relies on a tripartite complex formed by the cell surface lysin motif receptor-like kinases (LysM-RLKs) AtLYK1/CERK1, AtLYK4 and AtLYK5, the latter having the strongest affinity for COS. In grapevine, COS perception has for the moment only been demonstrated to rely on two PRRs VvLYK1-1 and VvLYK1-2. Here, we investigated additional players by overexpressing in Arabidopsis the two putative AtLYK5 orthologs from grapevine, VvLYK5-1 and VvLYK5-2. Expression of VvLYK5-1 in the atlyk4/5 double mutant background restored COS sensitivity, such as chitin-induced MAPK activation, defense gene expression, callose deposition and conferred non-host resistance to grapevine downy mildew (Erysiphe necator). Protein-protein interaction studies conducted in planta revealed a chitin oligomer-triggered interaction between VvLYK5-1 and VvLYK1-1. Interestingly, our results also indicate that VvLYK5-1 mediates the perception of chitin but not chitosan oligomers showing a part of its specificity.

Chromosome Painting Using Chromosome-Specific BAC Clones

Chromosome painting (CP) refers to visualization of large chromosome regions, chromosome arms or entire chromosomes via fluorescence in situ hybridization (FISH) of chromosome-specific DNA sequences. For CP in crucifers (Brassicaceae), typically contigs of chromosome-specific bacterial artificial chromosomes (BAC) from Arabidopsis thaliana are applied as painting probes on chromosomes of A. thaliana or other species (comparative chromosome painting, CCP). CP/CCP enables to identify and trace particular chromosome regions and/or chromosomes throughout all mitotic and meiotic stages as well as corresponding interphase chromosome territories. However, extended pachytene chromosomes provide the highest resolution of CP/CCP. Fine-scale chromosome structure, structural chromosome rearrangements (such as inversions, translocations, centromere repositioning), and chromosome breakpoints can be investigated by CP/CCP. BAC DNA probes can be accompanied by other types of DNA probes, such as repetitive DNA, genomic DNA, or synthetic oligonucleotide probes. Here, we describe a robust step-by-step protocol of CP and CCP which proved to be efficient across the family Brassicaceae, but which is also applicable to other angiosperm families.

POInT: Modeling Polyploidy in the Era of Ubiquitous Genomics

Thirteen years ago, we described an evolutionary modeling tool that could resolve the orthology relationships among the homologous genomic regions created by a whole-genome duplication. This tool, which we subsequently named POInT (the Polyploid Orthology Inference Tool), was originally only useful for studying a genome duplication known from bakers' yeast and its relatives. Now, with hundreds of genome sequences that contain the relicts of ancient polyploidy available, POInT can be used to study dozens of different polyploidies, asking both questions about the history of individual events and about the commonalities and differences seen between those events. In this chapter, I give a brief history of the development of POInT as an illustration of the interconnected nature of computational biology research. I then further describe how POInT operates and some of the strengths and drawbacks of its structure. I close with a few examples of discoveries we have made using it.

Genome assembly of an Australian native grass species reveals a recent whole-genome duplication and biased gene retention of genes involved in stress response

BACKGROUND

The adaptive significance of polyploidy has been extensively debated, and chromosome-level genome assemblies of polyploids can provide insight into this. The Australian grass Bothriochloa decipiens belongs to the BCD clade, a group with a complex history of hybridization and polyploid. This is the first genome assembly and annotation of a species that belongs to this fascinating yet complex group.

FINDINGS

Using Illumina short reads, 10X Genomics linked reads, and Hi-C sequencing data, we assembled a highly contiguous genome of B. decipiens, with a total length of 1,218.22 Mb and scaffold N50 of 42.637 Mb. Comparative analysis revealed that the species experienced a relatively recent whole-genome duplication. We clustered the 20 major scaffolds, representing the 20 chromosomes, into the 2 subgenomes of the parental species using unique repeat signatures. We found evidence of biased fractionation and differences in the activity of transposable elements between the subgenomes prior to hybridization. Duplicates were enriched for genes involved in transcription and response to external stimuli, supporting a biased retention of duplicated genes following whole-genome duplication.

CONCLUSIONS

Our results support the hypotheses of a biased retention of duplicated genes following polyploidy and point to differences in repeat activity associated with subgenome dominance. B. decipiens is a widespread species with the ability to establish across many soil types, making it a prime candidate for climate change- resilient ecological restoration of Australian grasslands. This reference genome is a valuable resource for future population genomic research on Australian grasses.

Origin of the Rare Hybrid Genus ×Trisetokoeleria Tzvelev (Poaceae) According to Molecular Phylogenetic Data

In our article, we analyzed new data on the origin of the hybrid genus ×Trisetokoeleria. According to the morphological criteria ×T. jurtzevii is a hybrid between Koeleria asiatica s. l. and Trisetum spicatum, ×T. taimyrica, and originated from Koeleria asiatica s. l. and Trisetum subalpestre, ×T. gorodkowii, a hybrid between Koeleria asiatica and Trisetum ruprechtianum. Later ×T. taimyrica was transferred to Koeleria. Parental taxa are prone to active hybridization themselves, thus, new methods of next-generation sequencing (NGS) were needed to clarify the relationships of these genera. For NGS we used the fragment 18S rDNA (part)-ITS1-5.8S rDNA (totally 441 accessions). We analyzed ITS1-5.8S rDNA-ITS2 region, trnL-trnF and trnK-rps16 from eight samples of the five species, using the Sanger method: ×Trisetokoeleria jurtzevii, ×T. taimyrica, Koeleria asiatica, Sibirotrisetum sibiricum (=Trisetum sibiricum), and Trisetum spicatum. We also studied the pollen fertility of ×Trisetokoeleria and its possible progenitors. Our data partly contradicted previous assumptions, based on morphological grounds, and showed us a picture of developed introgression within and between Koeleria and Trisetum. ×T. jurtzevii, a totally sterile hybrid formed rather recently. We can suppose that ×T. jurtzevii is a hybrid between K. asiatica and some Trisetum s. str. Species, but not T. spicatum. ×T. gorodkowii, a hybrid in the stage of primary stabilization; it has one unique ribotype related to T. spicatum s. l. The second parental species is unrelated to Trisetum ruprechtianum. ×T. taimyrica and is a stabilized hybrid species; it shares major ribotypes with the T. spicatum/T. wrangelense group and has a minor fraction of rDNA related to genus Deyeuxia s. l.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Chromosome Numbers and Genome Sizes of All 36 Duckweed Species (Lemnaceae)

Usually, chromosome sets (karyotypes) and genome sizes are rather stable for distinct species and therefore of diagnostic value for taxonomy. In combination with (cyto)genomics, both features provide essential cues for genome evolution and phylogenetic relationship studies within and between taxa above the species level. We present for the first time a survey on chromosome counts and genome size measurement for one or more accessions from all 36 duckweed species and discuss the evolutionary impact and peculiarities of both parameters in duckweeds.

Genome diploidization associates with cladogenesis, trait disparity, and plastid gene evolution

Angiosperm genome evolution was marked by many clade-specific whole-genome duplication events. The Microlepidieae is one of the monophyletic clades in the mustard family (Brassicaceae) formed after an ancient allotetraploidization. Postpolyploid cladogenesis has resulted in the extant c. 17 genera and 60 species endemic to Australia and New Zealand (10 species). As postpolyploid genome diploidization is a trial-and-error process under natural selection, it may proceed with different intensity and be associated with speciation events. In Microlepidieae, different extents of homoeologous recombination between the two parental subgenomes generated clades marked by slow ("cold") versus fast ("hot") genome diploidization. To gain a deeper understanding of postpolyploid genome evolution in Microlepidieae, we analyzed phylogenetic relationships in this tribe using complete chloroplast sequences, entire 35S rDNA units, and abundant repetitive sequences. The four recovered intra-tribal clades mirror the varied diploidization of Microlepidieae genomes, suggesting that the intrinsic genomic features underlying the extent of diploidization are shared among genera and species within one clade. Nevertheless, even congeneric species may exert considerable morphological disparity (e.g. in fruit shape), whereas some species within different clades experience extensive morphological convergence despite the different pace of their genome diploidization. We showed that faster genome diploidization is positively associated with mean morphological disparity and evolution of chloroplast genes (plastid-nuclear genome coevolution). Higher speciation rates in perennials than in annual species were observed. Altogether, our results confirm the potential of Microlepidieae as a promising subject for the analysis of postpolyploid genome diploidization in Brassicaceae.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

BACKGROUND

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

RESULTS

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

CONCLUSIONS

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

Idahoa and Subularia: Hidden polyploid origins of two enigmatic genera of crucifers

PREMISE

The monotypic Idahoa (I. scapigera) and the bispecific Subularia (S. aquatica and S. monticola) belong to Brassicaceae with unclear phylogenetic relationships and no tribal assignment. To fill this knowledge gap, we investigated these species and their closest relatives by combining cytogenomic and phylogenomic methods.

METHODS

We used whole plastome sequences in maximum likelihood and Bayesian inference analyses. We tested the phylogenetic informativeness of shared genomic repeats. We combined nuclear gene tree reconciliation and comparative chromosome painting (CCP) to examine the occurrence of past whole-genome duplications (WGDs).

RESULTS

The plastid data set corroborated the sister relationship between Idahoa and Subularia within the crucifer Lineage V but failed to resolve consistent topologies using both inference methods. The shared repetitive sequences provided conflicting pwhylogenetic signals. CCP analysis unexpectedly revealed that Idahoa (2n = 16) has a diploidized mesotetraploid genome, whereas two Subularia species (2n = 28 and 30) have diploidized mesoctoploid genomes. Several ancient allopolyploidy events have also been detected in closely related taxa (Chamira circaeoides, Cremolobeae, Eudemeae, and Notothlaspideae).

CONCLUSIONS

Our results suggest that the contentious phylogenetic placement of Idahoa and Subularia is best explained by two WGDs involving one or more shared parental genomes. The newly identified mesopolyploid genomes highlight the challenges of studying plant clades with complex polyploidy histories and provide a better framework for understanding genome evolution in the crucifer family.

Dynamic enhancer transcription associates with reprogramming of immune genes during pattern triggered immunity in Arabidopsis

BACKGROUND

Enhancers are cis-regulatory elements present in eukaryote genomes, which constitute indispensable determinants of gene regulation by governing the spatiotemporal and quantitative expression dynamics of target genes, and are involved in multiple life processes, for instance during development and disease states. The importance of enhancer activity has additionally been highlighted for immune responses in animals and plants; however, the dynamics of enhancer activities and molecular functions in plant innate immunity are largely unknown. Here, we investigated the involvement of distal enhancers in early innate immunity in Arabidopsis thaliana.

RESULTS

A group of putative distal enhancers producing low-abundance transcripts either unidirectionally or bidirectionally are identified. We show that enhancer transcripts are dynamically modulated in plant immunity triggered by microbe-associated molecular patterns and are strongly correlated with open chromatin, low levels of methylated DNA, and increases in RNA polymerase II targeting and acetylated histone marks. Dynamic enhancer transcription is correlated with target early immune gene expression patterns. Cis motifs that are bound by immune-related transcription factors, such as WRKYs and SARD1, are highly enriched within upregulated enhancers. Moreover, a subset of core pattern-induced enhancers are upregulated by multiple patterns from diverse pathogens. The expression dynamics of putative immunity-related enhancers and the importance of WRKY binding motifs for enhancer function were also validated.

CONCLUSIONS

Our study demonstrates the general occurrence of enhancer transcription in plants and provides novel information on the distal regulatory landscape during early plant innate immunity, providing new insights into immune gene regulation and ultimately improving the mechanistic understanding of the plant immune system.

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.

Genomes, repeatomes and interphase chromosome organization in the meadowfoam family (Limnanthaceae, Brassicales)

The meadowfoam family (Limnanthaceae) is one of the smallest and genomically underexplored families of the Brassicales. The Limnanthaceae harbor about seven species in the genus Limnanthes (meadowfoam) and Floerkea proserpinacoides (false mermaidweed), all native to North America. Because all Limnanthes and Floerkea species have only five chromosome pairs, i.e., a chromosome number rare in Brassicales and shared with Arabidopsis thaliana (Arabidopsis), we examined the Limnanthaceae genomes as a potential model system. Using low-coverage whole-genome sequencing data, we reexamined phylogenetic relationships and characterized the repeatomes of Limnanthaceae genomes. Phylogenies based on complete chloroplast and 35S rDNA sequences corroborated the sister relationship between Floerkea and Limnanthes and two major clades in the latter genus. The genome size of Limnanthaceae species ranges from 1.5 to 2.1 Gb, apparently due to the large increase in DNA repeats, which constitute 60-70% of their genomes. Repeatomes are dominated by long terminal repeat retrotransposons, while tandem repeats represent only less than 0.5% of the genomes. The average chromosome size in Limnanthaceae species (340-420 Mb) is more than 10 times larger than in Arabidopsis (32 Mb). A three-dimensional fluorescence in situ hybridization analysis demonstrated that the five chromosome pairs in interphase nuclei of Limnanthes species adopt the Rabl-like configuration.

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.

A chromosome-level genome assembly of an alpine plant Crucihimalaya lasiocarpa provides insights into high-altitude adaptation

It remains largely unknown how plants adapt to high-altitude habitats. Crucihimalaya (Brassicaceae) is an alpine genus occurring in the Qinghai-Tibet Plateau characterized by cold temperatures and strong ultraviolet radiation. Here, we generated a chromosome-level genome for C. lasiocarpa with a total size of 255.8 Mb and a scaffold N50 size of 31.9 Mb. We first examined the karyotype origin of this species and found that the karyotype of five chromosomes resembled the ancestral karyotype of the Brassicaceae family, while the other three showed strong chromosomal structural variations. In combination with the rough genome sequence of another congener (C. himalaica), we found that the significantly expanded gene families and positively selected genes involved in alpine adaptation have occurred since the origin of this genus. Our new findings provide valuable information for the chromosomal karyotype evolution of Brassicaceae and investigations of high-altitude environment adaptation of the genus.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.

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.

A hypomorphic allele of telomerase uncovers the minimal functional length of telomeres in Arabidopsis

Despite the essential requirement of telomeric DNA for genome stability, the length of telomere tracts between species substantially differs, raising the question of the minimal length of telomeric DNA necessary for proper function. Here, we address this question using a hypomorphic allele of the telomerase catalytic subunit, TERT. We show that although this construct partially restored telomerase activity to a tert mutant, telomeres continued to shorten over several generations, ultimately stabilizing at a bimodal size distribution. Telomeres on two chromosome arms were maintained at a length of 1 kb, while the remaining telomeres were maintained at 400 bp. The longest telomeres identified in this background were also significantly longer in wild-type populations, suggesting cis-acting elements on these arms either promote telomerase processivity or recruitment. Genetically disrupting telomerase processivity in this background resulted in immediate lethality. Thus, telomeres of 400 bp are both necessary and sufficient for Arabidopsis viability. As this length is the estimated minimal length for t-loop formation, our data suggest that telomeres long enough to form a t-loop constitute the minimal functional length.

Chromosome Painting Provides Insights Into the Genome Structure and Evolution of Sugarcane

The genus Saccharum is composed of species with high polyploidy and highly varied chromosome numbers, laying a challenge for uncovering its genomic structure and evolution. We developed a chromosome 2 painting (CP2) probe by designing oligonucleotides covering chromosome 2 of Saccharum spontaneum (2n = 8x = 64). Fluorescence in situ hybridization (FISH) using this CP2 probe revealed six types of ploidies from twenty S. spontaneum clones, including 6x, 8x, 10x, 11x, 12x, and 13x clones. The finding of S. spontaneum clones with uneven of ploid suggested that certain S. spontaneum clones come from hybridization. It renews our knowledge that S. spontaneum is derived from autopolyploidization. Combined with a S. spontaneum-specific probe, chromosome 2-derived chromosome or fragments from either S. spontaneum or Saccharum officinarum can be identified in sugarcane modern cultivars. We revealed unexpected high level of interspecific recombination from introgressive S. spontaneum chromosomes (>50.0%) in cultivars ROC22 and ZZ1, indicating frequent chromosome exchange in cultivars. Intriguingly, we observed interspecific recombination recurring among either homoeologous or non-homoeologous chromosomes in sugarcane cultivars. These results demonstrated that chromosome painting FISH is a powerful tool in the genome dissection of sugarcane and provide new insights into the genome structure and evolution of the complex genus Saccharum.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.

Phenotypic diploidization in plant functional traits uncovered by synthetic neopolyploids in Dianthus broteri

Whole-genome duplication and post-polyploidization genome downsizing play key roles in the evolution of land plants; however, the impact of genomic diploidization on functional traits still remains poorly understood. Using Dianthus broteri as a model, we compared the ecophysiological behaviour of colchicine-induced neotetraploids (4xNeo) to diploids (2x) and naturally occurring tetraploids (4xNat). Leaf gas-exchange and chlorophyll fluorescence analyses were performed in order to asses to what extent post-polyploidization evolutionary processes have affected 4xNat. Genomic diploidization and phenotypic novelty were evident. Distinct patterns of variation revealed that post-polyploidization processes altered the phenotypic shifts directly mediated by genome doubling. The photosynthetic phenotype was affected in several ways but the main effect was phenotypic diploidization (i.e. 2x and 4xNat were closer to each other than to 4xNeo). Overall, our results show the potential benefits of considering experimentally synthetized versus naturally established polyploids when exploring the role of polyploidization in promoting functional divergence.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5'-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

Functional divergence of Brassica napus BnaABI1 paralogs in the structurally conserved PP2CA gene subfamily of Brassicaceae

Group A PP2C (PP2CA) genes form a gene subfamily whose members play an important role in regulating many biological processes by dephosphorylation of target proteins. In this study we examined the effects of evolutionary changes responsible for functional divergence of BnaABI1 paralogs in Brassica napus against the background of the conserved PP2CA gene subfamily in Brassicaceae. We performed comprehensive phylogenetic analyses of 192 PP2CA genes in 15 species in combination with protein structure homology modeling. Fundamentally, the number of PP2CA genes remained relatively constant in these taxa, except in the Brassica genus and Camelina sativa. The expansion of this gene subfamily in these species has resulted from whole genome duplication. We demonstrated a high degree of structural conservation of the PP2CA genes, with a few minor variations between the different PP2CA groups. Furthermore, the pattern of conserved sequence motifs in the PP2CA proteins and their secondary and 3D structures revealed strong conservation of the key ion-binding sites. Syntenic analysis of triplicated regions including ABI1 paralogs revealed significant structural rearrangements of the Brassica genomes. The functional and syntenic data clearly show that triplication of BnaABI1 in B. napus has had an impact on its functions, as well as the positions of adjacent genes in the corresponding chromosomal regions. The expression profiling of BnaABI1 genes showed functional divergence, i.e. subfunctionalization, potentially leading to neofunctionalization. These differences in expression are likely due to changes in the promoters of the BnaABI1 paralogs. Our results highlight the complexity of PP2CA gene subfamily evolution in Brassicaceae.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

Revisiting the Karyotypes of Alligators and Caimans (Crocodylia, Alligatoridae) after a Half-Century Delay: Bridging the Gap in the Chromosomal Evolution of Reptiles

Although crocodilians have attracted enormous attention in other research fields, from the cytogenetic point of view, this group remains understudied. Here, we analyzed the karyotypes of eight species formally described from the Alligatoridae family using differential staining, fluorescence in situ hybridization with rDNA and repetitive motifs as a probe, whole chromosome painting (WCP), and comparative genome hybridization. All Caimaninae species have a diploid chromosome number (2n) 42 and karyotypes dominated by acrocentric chromosomes, in contrast to both species of Alligatorinae, which have 2n = 32 and karyotypes that are predominantly metacentric, suggesting fusion/fission rearrangements. Our WCP results supported this scenario by revealing the homeology of the largest metacentric pair present in both Alligator spp. with two smaller pairs of acrocentrics in Caimaninae species. The clusters of 18S rDNA were found on one chromosome pair in all species, except for Paleosuchus spp., which possessed three chromosome pairs bearing these sites. Similarly, comparative genomic hybridization demonstrated an advanced stage of sequence divergence among the caiman genomes, with Paleosuchus standing out as the most divergent. Thus, although Alligatoridae exhibited rather low species diversity and some level of karyotype stasis, their genomic content indicates that they are not as conserved as previously thought. These new data deepen the discussion of cytotaxonomy in this family.

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.

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

Multimerization variants as potential drivers of neofunctionalization

Whole-genome duplications are common during evolution, creating genetic redundancy that can enable cellular innovations. Novel protein-protein interactions provide a route to diversified gene functions, but, at present, there is limited proteome-scale knowledge on the extent to which variability in protein complex formation drives neofunctionalization. Here, we used protein correlation profiling to test for variability in apparent mass among thousands of orthologous proteins isolated from diverse species and cell types. Variants in protein complex size were unexpectedly common, in some cases appearing after relatively recent whole-genome duplications or an allopolyploidy event. In other instances, variants such as those in the carbonic anhydrase orthologous group reflected the neofunctionalization of ancient paralogs that have been preserved in extant species. Our results demonstrate that homo- and heteromer formation have the potential to drive neofunctionalization in diverse classes of enzymes, signaling, and structural proteins.

Chromosome inheritance and meiotic stability in allopolyploid Brassica napus

Homoeologous recombination, aneuploidy, and other genetic changes are common in resynthesized allopolyploid Brassica napus. In contrast, the chromosomes of cultivars have long been considered to be meiotically stable. To gain a better understanding of the underlying mechanisms leading to stabilization in the allopolyploid, the behavior of chromosomes during meiosis can be compared by unambiguous chromosome identification between resynthesized and natural B. napus. Compared with natural B. napus, resynthesized lines show high rates of nonhomologous centromere association, homoeologous recombination leading to translocation, homoeologous chromosome replacement, and association and breakage of 45S rDNA loci. In both natural and resynthesized B. napus, we observed low rates of univalents, A–C bivalents, and early sister chromatid separations. Reciprocal homoeologous chromosome exchanges and double reductions were photographed for the first time in meiotic telophase I. Meiotic errors were non-uniformly distributed across the genome in resynthesized B. napus, and in particular homoeologs sharing synteny along their entire length exhibited multivalents at diakinesis and polysomic inheritance at telophase I. Natural B. napus appeared to resolve meiotic errors mainly by suppressing homoeologous pairing, resolving nonhomologous centromere associations and 45S rDNA associations before diakinesis, and reducing homoeologous cross-overs.

Chromosomal Mapping of Tandem Repeats Revealed Massive Chromosomal Rearrangements and Insights Into Senna tora Dysploidy

Tandem repeats can occupy a large portion of plant genomes and can either cause or result from chromosomal rearrangements, which are important drivers of dysploidy-mediated karyotype evolution and speciation. To understand the contribution of tandem repeats in shaping the extant Senna tora dysploid karyotype, we analyzed the composition and abundance of tandem repeats in the S. tora genome and compared the chromosomal distribution of these repeats between S. tora and a closely related euploid, Senna occidentalis. Using a read clustering algorithm, we identified the major S. tora tandem repeats and visualized their chromosomal distribution by fluorescence in situ hybridization. We identified eight independent repeats covering ~85 Mb or ~12% of the S. tora genome. The unit lengths and copy numbers had ranges of 7-5,833 bp and 325-2.89 × 10⁶, respectively. Three short duplicated sequences were found in the 45S rDNA intergenic spacer, one of which was also detected at an extra-NOR locus. The canonical plant telomeric repeat (TTTAGGG)_n was also detected as very intense signals in numerous pericentromeric and interstitial loci. StoTR05_180, which showed subtelomeric distribution in Senna occidentalis, was predominantly pericentromeric in S. tora. The unusual chromosomal distribution of tandem repeats in S. tora not only enabled easy identification of individual chromosomes but also revealed the massive chromosomal rearrangements that have likely played important roles in shaping its dysploid karyotype.

The Genome Sequence of Alpine Megacarpaea delavayi Identifies Species-Specific Whole-Genome Duplication

Megacarpaea delavayi (Brassicaceae), a plant found the high mountains of southwest China at high altitudes (3000–4800 m), is used as a vegetable or medicine. Here, we report a draft genome for this species. The assembly genome of M. delavayi is 883 Mb, and 61.59% of the genome is composed of repeat sequences. Annotation of the genome identified a total of 41,114 protein-coding genes. We found that M. delavayi experienced an independent whole-genome duplication (WGD), paralleling those independent WGDs in Iberis, Biscutella, and Anastatica in the early Miocene. Phylogenetic analyses based on the single-copy genes confirmed the position of the genus Megacarpaea within the expanded lineage II of the family and resolved its basal divergence to a subclade consisting of Anastatica, Iberis, and Biscutella. Species-specific and fast-evolving genes in M. delavayi are mainly involved in “DNA repair” and “response to UV-B radiation.” These genetic changes may together help this species survive in high-altitude environments. The reference genome reported here provides a valuable resource for studying adaptation of this and other alpine plants to the high-altitude habitats.

Nested whole-genome duplications coincide with diversification and high morphological disparity in Brassicaceae

Angiosperms have become the dominant terrestrial plant group by diversifying for ~145 million years into a broad range of environments. During the course of evolution, numerous morphological innovations arose, often preceded by whole genome duplications (WGD). The mustard family (Brassicaceae), a successful angiosperm clade with ~4000 species, has been diversifying into many evolutionary lineages for more than 30 million years. Here we develop a species inventory, analyze morphological variation, and present a maternal, plastome-based genus-level phylogeny. We show that increased morphological disparity, despite an apparent absence of clade-specific morphological innovations, is found in tribes with WGDs or diversification rate shifts. Both are important processes in Brassicaceae, resulting in an overall high net diversification rate. Character states show frequent and independent gain and loss, and form varying combinations. Therefore, Brassicaceae pave the way to concepts of phylogenetic genome-wide association studies to analyze the evolution of morphological form and function.

Nested whole-genome duplications coincide with diversification and high morphological disparity in Brassicaceae

Angiosperms have become the dominant terrestrial plant group by diversifying for ~145 million years into a broad range of environments. During the course of evolution, numerous morphological innovations arose, often preceded by whole genome duplications (WGD). The mustard family (Brassicaceae), a successful angiosperm clade with ~4000 species, has been diversifying into many evolutionary lineages for more than 30 million years. Here we develop a species inventory, analyze morphological variation, and present a maternal, plastome-based genus-level phylogeny. We show that increased morphological disparity, despite an apparent absence of clade-specific morphological innovations, is found in tribes with WGDs or diversification rate shifts. Both are important processes in Brassicaceae, resulting in an overall high net diversification rate. Character states show frequent and independent gain and loss, and form varying combinations. Therefore, Brassicaceae pave the way to concepts of phylogenetic genome-wide association studies to analyze the evolution of morphological form and function.

Comparative analyses of DNA repeats and identification of a novel Fesreba centromeric element in fescues and ryegrasses

BACKGROUND

Cultivated grasses are an important source of food for domestic animals worldwide. Increased knowledge of their genomes can speed up the development of new cultivars with better quality and greater resistance to biotic and abiotic stresses. The most widely grown grasses are tetraploid ryegrass species (Lolium) and diploid and hexaploid fescue species (Festuca). In this work, we characterized repetitive DNA sequences and their contribution to genome size in five fescue and two ryegrass species as well as one fescue and two ryegrass cultivars.

RESULTS

Partial genome sequences produced by Illumina sequencing technology were used for genome-wide comparative analyses with the RepeatExplorer pipeline. Retrotransposons were the most abundant repeat type in all seven grass species. The Athila element of the Ty3/gypsy family showed the most striking differences in copy number between fescues and ryegrasses. The sequence data enabled the assembly of the long terminal repeat (LTR) element Fesreba, which is highly enriched in centromeric and (peri)centromeric regions in all species. A combination of fluorescence in situ hybridization (FISH) with a probe specific to the Fesreba element and immunostaining with centromeric histone H3 (CENH3) antibody showed their co-localization and indicated a possible role of Fesreba in centromere function.

CONCLUSIONS

Comparative repeatome analyses in a set of fescues and ryegrasses provided new insights into their genome organization and divergence, including the assembly of the LTR element Fesreba. A new LTR element Fesreba was identified and found in abundance in centromeric regions of the fescues and ryegrasses. It may play a role in the function of their centromeres.

Two combinatorial patterns of telomere histone marks in plants with canonical and non-canonical telomere repeats

Telomeres, nucleoprotein structures at the ends of linear eukaryotic chromosomes, are crucial for the maintenance of genome integrity. In most plants, telomeres consist of conserved tandem repeat units comprising the TTTAGGG motif. Recently, non-canonical telomeres were described in several plants and plant taxons, including the carnivorous plant Genlisea hispidula (TTCAGG/TTTCAGG), the genus Cestrum (Solanaceae; TTTTTTAGGG), and plants from the Asparagales order with either a vertebrate-type telomere repeat TTAGGG or Allium genus-specific CTCGGTTATGGG repeat. We analyzed epigenetic modifications of telomeric histones in plants with canonical and non-canonical telomeres, and further in telomeric chromatin captured from leaves of Nicotiana benthamiana transiently transformed by telomere CRISPR-dCas9-eGFP, and of Arabidopsis thaliana stably transformed with TALE_telo C-3×GFP. Two combinatorial patterns of telomeric histone modifications were identified: (i) an Arabidopsis-like pattern (A. thaliana, G. hispidula, Genlisea nigrocaulis, Allium cepa, Narcissus pseudonarcissus, Petunia hybrida, Solanum tuberosum, Solanum lycopersicum) with telomeric histones decorated predominantly by H3K9me2; (ii) a tobacco-like pattern (Nicotiana tabacum, N. benthamiana, C. elegans) with a strong H3K27me3 signal. Our data suggest that epigenetic modifications of plant telomere-associated histones are related neither to the sequence of the telomere motif nor to the lengths of the telomeres. Nor the phylogenetic position of the species plays the role; representatives of the Solanaceae family are included in both groups. As both patterns of histone marks are compatible with fully functional telomeres in respective plants, we conclude that the described specific differences in histone marks are not critical for telomere functions.

Flexible chromosome painting based on multiplex PCR of oligonucleotides and its application for comparative chromosome analyses in Cucumis

Chromosome painting is a powerful technique for chromosome and genome studies. We developed a flexible chromosome painting technique based on multiplex PCR of a synthetic oligonucleotide (oligo) library in cucumber (Cucumis sativus L., 2n = 14). Each oligo in the library was associated with a universal as well as nested specific primers for amplification, which allow the generation of different probes from the same oligo library. We were also able to generate double-stranded labelled oligos, which produced much stronger signals than single-stranded labelled oligos, by amplification using fluorophore-conjugated primer pairs. Oligos covering cucumber chromosome 1 (Chr1) and chromosome 4 (Chr4) consisting of eight segments were synthesized in one library. Different oligo probes generated from the library painted the corresponding chromosomes/segments unambiguously, especially on pachytene chromosomes. This technique was then applied to study the homoeologous relationships among cucumber, C. hystrix and C. melo chromosomes based on cross-species chromosome painting using Chr4 probes. We demonstrated that the probe was feasible to detect interspecies chromosome homoeologous relationships and chromosomal rearrangement events. Based on its advantages and great convenience, we anticipate that this flexible oligo-painting technique has great potential for the studies of the structure, organization, and evolution of chromosomes in any species with a sequenced genome.

Recurrent gene duplication in the angiosperm tribe Delphinieae (Ranunculaceae) inferred from intracellular gene transfer events and heteroplasmic mutations in the plastid matK gene

The study of intracellular gene transfer may allow for the detection of interesting evolutionary processes such as ancient polyploidization. We compared 24 plastid genomes (plastomes) from tribe Delphinieae, one from tribe Nigelleae and one from tribe Ranunculeae, including five newly sequenced genomes. The functional transfers of the plastids rpl32 and rps16 to the nucleus in tribe Delphinieae were identified. Unexpectedly, we discovered multiple divergent copies of the nuclear-encoded plastid rpl32 in the genus Aconitum. Phylogenetic and synonymous substitution rate analyses revealed that the nuclear-encoded plastid rpl32 underwent two major duplication events. These ancient gene duplication events probably occurred via multiple polyploidization events in Aconitum between 11.9 and 24.7 Mya. Furthermore, our sequence rate analysis indicated that the eight plastid-encoded rpl subunits in Aconitum had a significantly accelerated evolutionary rate compared to those in other genera, suggesting that highly divergent paralogs targeted to the plastid may contribute to an elevated rate of evolution in plastid rpl genes. In addition, heteroplasmy of the plastid matK from two Aconitum species suggested the existence of potentially functional plastid maturases in its plastome. Our results provide insight into the evolutionary history of the tribe Delphinieae.

Chromosome-Scale Assembly of Winter Oilseed Rape Brassica napus

Rapeseed (Brassica napus), the second most important oilseed crop globally, originated from an interspecific hybridization between B. rapa and B. oleracea. After this genome collision, B. napus underwent extensive genome restructuring, via homoeologous chromosome exchanges, resulting in widespread segmental deletions and duplications. Illicit pairing among genetically similar homoeologous chromosomes during meiosis is common in recent allopolyploids like B. napus, and post-polyploidization restructuring compounds the difficulties of assembling a complex polyploid plant genome. Specifically, genomic rearrangements between highly similar chromosomes are challenging to detect due to the limitation of sequencing read length and ambiguous alignment of reads. Recent advances in long read sequencing technologies provide promising new opportunities to unravel the genome complexities of B. napus by encompassing breakpoints of genomic rearrangements with high specificity. Moreover, recent evidence revealed ongoing genomic exchanges in natural B. napus, highlighting the need for multiple reference genomes to capture structural variants between accessions. Here we report the first long-read genome assembly of a winter B. napus cultivar. We sequenced the German winter oilseed rape accession 'Express 617' using 54.5x of long reads. Short reads, linked reads, optical map data and high-density genetic maps were used to further correct and scaffold the assembly to form pseudochromosomes. The assembled Express 617 genome provides another valuable resource for Brassica genomics in understanding the genetic consequences of polyploidization, crop domestication, and breeding of recently-formed crop species.

Chromosome-Scale Assembly of Winter Oilseed Rape Brassica napus

Rapeseed (Brassica napus), the second most important oilseed crop globally, originated from an interspecific hybridization between B. rapa and B. oleracea. After this genome collision, B. napus underwent extensive genome restructuring, via homoeologous chromosome exchanges, resulting in widespread segmental deletions and duplications. Illicit pairing among genetically similar homoeologous chromosomes during meiosis is common in recent allopolyploids like B. napus, and post-polyploidization restructuring compounds the difficulties of assembling a complex polyploid plant genome. Specifically, genomic rearrangements between highly similar chromosomes are challenging to detect due to the limitation of sequencing read length and ambiguous alignment of reads. Recent advances in long read sequencing technologies provide promising new opportunities to unravel the genome complexities of B. napus by encompassing breakpoints of genomic rearrangements with high specificity. Moreover, recent evidence revealed ongoing genomic exchanges in natural B. napus, highlighting the need for multiple reference genomes to capture structural variants between accessions. Here we report the first long-read genome assembly of a winter B. napus cultivar. We sequenced the German winter oilseed rape accession 'Express 617' using 54.5x of long reads. Short reads, linked reads, optical map data and high-density genetic maps were used to further correct and scaffold the assembly to form pseudochromosomes. The assembled Express 617 genome provides another valuable resource for Brassica genomics in understanding the genetic consequences of polyploidization, crop domestication, and breeding of recently-formed crop species.

The story of promiscuous crucifers: origin and genome evolution of an invasive species, Cardamine occulta (Brassicaceae), and its relatives

BACKGROUND AND AIMS

Cardamine occulta (Brassicaceae) is an octoploid weedy species (2n = 8x = 64) originated in Eastern Asia. It has been introduced to other continents including Europe and considered to be an invasive species. Despite its wide distribution, the polyploid origin of C. occulta remained unexplored. The feasibility of comparative chromosome painting (CCP) in crucifers allowed us to elucidate the origin and genome evolution in Cardamine species. We aimed to investigate the genome structure of C. occulta in comparison with its tetraploid (2n = 4x = 32, C. kokaiensis and C. scutata) and octoploid (2n = 8x = 64, C. dentipetala) relatives.

METHODS

Genomic in situ hybridization (GISH) and large-scale CCP were applied to uncover the parental genomes and chromosome composition of the investigated Cardamine species.

KEY RESULTS

All investigated species descended from a common ancestral Cardamine genome (n = 8), structurally resembling the Ancestral Crucifer Karyotype (n = 8), but differentiated by a translocation between chromosomes AK6 and AK8. Allotetraploid C. scutata originated by hybridization between two diploid species, C. parviflora and C. amara (2n = 2x = 16). By contrast, C. kokaiensis has an autotetraploid origin from a parental genome related to C. parviflora. Interestingly, octoploid C. occulta probably originated through hybridization between the tetraploids C. scutata and C. kokaiensis. The octoploid genome of C. dentipetala probably originated from C. scutata via autopolyploidization. Except for five species-specific centromere repositionings and one pericentric inversion post-dating the polyploidization events, the parental subgenomes remained stable in the tetra- and octoploids.

CONCLUSIONS

Comparative genome structure, origin and evolutionary history was reconstructed in C. occulta and related species. For the first time, whole-genome cytogenomic maps were established for octoploid plants. Post-polyploid evolution in Asian Cardamine polyploids has not been associated with descending dysploidy and intergenomic rearrangements. The combination of different parental (sub)genomes adapted to distinct habitats provides an evolutionary advantage to newly formed polyploids by occupying new ecological niches.

Genome invasion by a hypomethylated satellite repeat in Australian crucifer Ballantinia antipoda

Repetitive sequences are ubiquitous components of all eukaryotic genomes. They contribute to genome evolution and the regulation of gene transcription. However, the uncontrolled activity of repetitive sequences can negatively affect genome functions and stability. Therefore, repetitive DNAs are embedded in a highly repressive heterochromatic environment in plant cell nuclei. Here, we analyzed the sequence, composition and the epigenetic makeup of peculiar non-pericentromeric heterochromatic segments in the genome of the Australian crucifer Ballantinia antipoda. By the combination of high throughput sequencing, graph-based clustering and cytogenetics, we found that the heterochromatic segments consist of a mixture of unique sequences and an A-T-rich 174 bp satellite repeat (BaSAT1). BaSAT1 occupies about 10% of the B. antipoda nuclear genome in >250 000 copies. Unlike many other highly repetitive sequences, BaSAT1 repeats are hypomethylated; this contrasts with the normal patterns of DNA methylation in the B. antipoda genome. Detailed analysis of several copies revealed that these non-methylated BaSAT1 repeats were also devoid of heterochromatic histone H3K9me2 methylation. However, the factors decisive for the methylation status of BaSAT1 repeats remain currently unknown. In summary, we show that even highly repetitive sequences can exist as hypomethylated in the plant nuclear genome.

Polyploid evolution: The ultimate way to grasp the nettle

Polyploidy is one of the major forces of plant evolution and widespread mixed-ploidy species offer an opportunity to evaluate its significance. We therefore selected the cosmopolitan species Urtica dioica (stinging nettle), examined its cytogeography and pattern of absolute genome size, and assessed correlations with bioclimatic and ecogeographic data (latitude, longitude, elevation). We evaluated variation in ploidy level using an extensive dataset of 7012 samples from 1317 populations covering most of the species' distribution area. The widespread tetraploid cytotype (87%) was strongly prevalent over diploids (13%). A subsequent analysis of absolute genome size proved a uniform Cx-value of core U. dioica (except for U. d. subsp. cypria) whereas other closely related species, namely U. bianorii, U. kioviensis and U. simensis, differed significantly. We detected a positive correlation between relative genome size and longitude and latitude in the complete dataset of European populations and a positive correlation between relative genome size and longitude in a reduced dataset of diploid accessions (the complete dataset of diploids excluding U. d. subsp. kurdistanica). In addition, our data indicate an affinity of most diploids to natural and near-natural habitats and that the tetraploid cytotype and a small part of diploids (population from the Po river basin in northern Italy) tend to inhabit synanthropic sites. To sum up, the pattern of ploidy variation revealed by our study is in many aspects unique to the stinging nettle, being most likely first of all driven by the greater ecological plasticity and invasiveness of the tetraploid cytotype.

Genome-wide characterization of ALDH Superfamily in Brassica rapa and enhancement of stress tolerance in heterologous hosts by BrALDH7B2 expression

Aldehyde dehydrogenase (ALDH) carries out oxidation of toxic aldehydes using NAD⁺/NADP⁺ as cofactors. In the present study, we performed a genome-wide identification and expression analysis of genes in the ALDH gene family in Brassica rapa. A total of 23 ALDH genes in the superfamily have been identified according to the classification of ALDH Gene Nomenclature Committee (AGNC). They were distributed unevenly across all 10 chromosomes. All the 23 Brassica rapa ALDH (BrALDH) genes exhibited varied expression patterns during treatments with abiotic stress inducers and hormonal treatments. The relative expression profiles of ALDH genes in B. rapa showed that they are predominantly expressed in leaves and stem suggesting their function in the vegetative tissues. BrALDH7B2 showed a strong response to abiotic stress and hormonal treatments as compared to other ALDH genes; therefore, it was overexpressed in heterologous hosts, E. coli and yeast to study its possible function under abiotic stress conditions. Over-expression of BrALDH7B2 in heterologous systems, E. coli and yeast cells conferred significant tolerance to abiotic stress treatments. Results from this work demonstrate that BrALDH genes are a promising and untapped genetic resource for crop improvement and could be deployed further in the development of drought and salinity tolerance in B. rapa and other economically important crops.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.