Comparative Genome Analysis Reveals Divergent Genome Size Evolution in a Carnivorous Plant Genus

Molecular Phylogenomics Reveals the Deep Evolutionary History of Carnivory across Land Plants

Plastid molecular phylogenies that broadly sampled angiosperm lineages imply that carnivorous plants evolved at least 11 times independently in 13 families and 6 orders. Within and between these clades, the different prey capture strategies involving flypaper and pitfall structures arose in parallel with the subsequent evolution of snap traps and suction bladders. Attempts to discern the deep ontological history of carnivorous structures using multigene phylogenies have provided a plastid-level picture of sister relationships at the family level. Here, we present a molecular phylogeny of the angiosperms based on nuclear target sequence capture data (Angiosperms-353 probe set), assembled by the Kew Plant Trees of Life initiative, which aims to complete the tree of life for plants. This phylogeny encompasses all carnivorous and protocarnivorous families, although certain genera such as Philcoxia (Plantaginaceae) are excluded. This study offers a novel nuclear gene-based overview of relationships within and between carnivorous families and genera. Consistent with previous broadly sampled studies, we found that most carnivorous families are not affiliated with any single family. Instead, they emerge as sister groups to large clades comprising multiple non-carnivorous families. Additionally, we explore recent genomic studies across various carnivorous clades that examine the evolution of the carnivorous syndrome in relation to whole-genome duplication, subgenome dominance, small-scale gene duplication, and convergent evolution. Furthermore, we discuss insights into genome size evolution through the lens of carnivorous plant genomes.

Disentangling Crocus Series Verni and Its Polyploids

Spring crocuses, the eleven species within Crocus series Verni (Iridaceae), consist of di- and tetraploid cytotypes. Among them is a group of polyploids from southeastern Europe with yet-unclear taxonomic affiliation. Crocuses are generally characterized by complex dysploid chromosome number changes, preventing a clear correlation between these numbers and ploidy levels. To reconstruct the evolutionary history of series Verni and particularly its polyploid lineages associated with C. heuffelianus, we used an approach combining phylogenetic analyses of two chloroplast regions, 14 nuclear single-copy genes plus rDNA spacers, genome-wide genotyping-by-sequencing (GBS) data, and morphometry with ploidy estimations through genome size measurements, analysis of genomic heterozygosity frequencies and co-ancestry, and chromosome number counts. Chromosome numbers varied widely in diploids with 2n = 8, 10, 12, 14, 16, and 28 and tetraploid species or cytotypes with 2n = 16, 18, 20, and 22 chromosomes. Crocus longiflorus, the diploid with the highest chromosome number, possesses the smallest genome (2C = 3.21 pg), while the largest diploid genomes are in a range of 2C = 7-8 pg. Tetraploid genomes have 2C values between 10.88 pg and 12.84 pg. Heterozygosity distribution correlates strongly with genome size classes and allows discernment of di- and tetraploid cytotypes. Our phylogenetic analyses showed that polyploids in the C. heuffelianus group are allotetraploids derived from multiple and partly reciprocal crosses involving different genotypes of diploid C. heuffelianus (2n = 10) and C. vernus (2n = 8). Dysploid karyotype changes after polyploidization resulted in the tetraploid cytotypes with 20 and 22 chromosomes. The multi-data approach we used here for series Verni, combining evidence from nuclear and chloroplast phylogenies, genome sizes, chromosome numbers, and genomic heterozygosity for ploidy estimations, provides a way to disentangle the evolution of plant taxa with complex karyotype changes that can be used for the analysis of other groups within Crocus and beyond. Comparing these results with morphometric analysis results in characters that can discern the different taxa currently subsumed under C. heuffelianus.

Transposons and Non-coding Regions Drive the Intrafamily Differences of Genome Size in Insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein–Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

•

The most comprehensive variation pattern of insects genome size (GS) to date

•

GS evolution of the Insecta was reflected by the adaptive Ornstein–Uhlenbeck model

•

Ancestral state of insect GS was estimated to be ∼1 Gb

•

Intrafamily GS variations were driven by recent transpositions and non-coding regions

The complete mitochondrial genome of carnivorous Genlisea tuberosa (Lentibulariaceae): Structure and evolutionary aspects

Sequenced genomic data for carnivorous plants are scarce, especially regarding the mitogenomes (MTs) and further studies are crucial to obtain a better understanding of the topic. In this study, we sequenced and characterized the mitochondrial genome of the tuberous carnivorous plant Genlisea tuberosa, being the first of its genus to be sequenced. The genome comprises 729,765 bp, encoding 80 identified genes of which 36 are protein-coding, 40 tRNA, four rRNA genes, and three pseudogenes. An intronic region from the cox1 gene was identified that encodes an endonuclease enzyme that is present in the other sequenced species of Lentibulariaceae. Chloroplast genes (pseudogene and complete) inserted in the MT genome were identified, showing possible horizontal transfer between organelles. In addition, 50 pairs of long repeats from 94 to 274 bp are present, possibly playing an important role in the maintenance of the MT genome. Phylogenetic analysis carried out with 34 coding mitochondrial genes corroborated the positioning of the species listed here within the family. The molecular dynamism in the mitogenome (e.g. the loss or pseudogenization of genes, insertion of foreign genes, the long repeats as well as accumulated mutations) may be reflections of the carnivorous lifestyle where a significant part of cellular energy was shifted for the adaptation of leaves into traps molding the mitochondrial DNA. The sequence and annotation of G. tuberosa's MT will be useful for further studies and serve as a model for evolutionary and taxonomic clarifications of the group as well as improving our comprehension of MT evolution.

A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

Background

Parasite evolution has been conceptualized as a process of genetic loss and simplification. Contrary to this model, there is evidence of expansion and conservation of gene families related to essential functions of parasitism in some parasite genomes, reminiscent of widespread mosaic evolution—where subregions of a genome have different rates of evolutionary change. We found evidence of mosaic genome evolution in the cnidarian Myxobolus honghuensis, a myxozoan parasite of fish, with extremely simple morphology.

Results

We compared M. honghuensis with other myxozoans and free-living cnidarians, and determined that it has a relatively larger myxozoan genome (206 Mb), which is less reduced and less compact due to gene retention, large introns, transposon insertion, but not polyploidy. Relative to other metazoans, the M. honghuensis genome is depleted of neural genes and has only the simplest animal immune components. Conversely, it has relatively more genes involved in stress resistance, tissue invasion, energy metabolism, and cellular processes compared to other myxozoans and free-living cnidarians. We postulate that the expansion of these gene families is the result of evolutionary adaptations to endoparasitism. M. honghuensis retains genes found in free-living Cnidaria, including a reduced nervous system, myogenic components, ANTP class Homeobox genes, and components of the Wnt and Hedgehog pathways.

Conclusions

Our analyses suggest that the M. honghuensis genome evolved as a mosaic of conservative, divergent, depleted, and enhanced genes and pathways. These findings illustrate that myxozoans are not as genetically simple as previously regarded, and the evolution of some myxozoans is driven by both genomic streamlining and expansion.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01249-8.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Recent ecophysiological, biochemical and evolutional insights into plant carnivory

BACKGROUND

Carnivorous plants are an ecological group of approx. 810 vascular species which capture and digest animal prey, absorb prey-derived nutrients and utilize them to enhance their growth and development. Extant carnivorous plants have evolved in at least ten independent lineages, and their adaptive traits represent an example of structural and functional convergence. Plant carnivory is a result of complex adaptations to mostly nutrient-poor, wet and sunny habitats when the benefits of carnivory exceed the costs. With a boost in interest and extensive research in recent years, many aspects of these adaptations have been clarified (at least partly), but many remain unknown.

SCOPE

We provide some of the most recent insights into substantial ecophysiological, biochemical and evolutional particulars of plant carnivory from the functional viewpoint. We focus on those processes and traits in carnivorous plants associated with their ecological characterization, mineral nutrition, cost-benefit relationships, functioning of digestive enzymes and regulation of the hunting cycle in traps. We elucidate mechanisms by which uptake of prey-derived nutrients leads to stimulation of photosynthesis and root nutrient uptake.

CONCLUSIONS

Utilization of prey-derived mineral (mainly N and P) and organic nutrients is highly beneficial for plants and increases the photosynthetic rate in leaves as a prerequisite for faster plant growth. Whole-genome and tandem gene duplications brought gene material for diversification into carnivorous functions and enabled recruitment of defence-related genes. Possible mechanisms for the evolution of digestive enzymes are summarized, and a comprehensive picture on the biochemistry and regulation of prey decomposition and prey-derived nutrient uptake is provided.

Plant Pan-Genomics Comes of Age

A pan-genome is the nonredundant collection of genes and/or DNA sequences in a species. Numerous studies have shown that plant pan-genomes are typically much larger than the genome of any individual and that a sizable fraction of the genes in any individual are present in only some genomes. The construction and interpretation of plant pan-genomes are challenging due to the large size and repetitive content of plant genomes. Most pan-genomes are largely focused on nontransposable element protein coding genes because they are more easily analyzed and defined than noncoding and repetitive sequences. Nevertheless, noncoding and repetitive DNA play important roles in determining the phenotype and genome evolution. Fortunately, it is now feasible to make multiple high-quality genomes that can be used to construct high-resolution pan-genomes that capture all the variation. However, assembling, displaying, and interacting with such high-resolution pan-genomes will require the development of new tools.

Regarding the F-word: The effects of data filtering on inferred genotype-environment associations

Genotype-environment association (GEA) methods have become part of the standard landscape genomics toolkit, yet, we know little about how to best filter genotype-by-sequencing data to provide robust inferences for environmental adaptation. In many cases, default filtering thresholds for minor allele frequency and missing data are applied regardless of sample size, having unknown impacts on the results, negatively affecting management strategies. Here, we investigate the effects of filtering on GEA results and the potential implications for assessment of adaptation to environment. We use empirical and simulated data sets derived from two widespread tree species to assess the effects of filtering on GEA outputs. Critically, we find that the level of filtering of missing data and minor allele frequency affect the identification of true positives. Even slight adjustments to these thresholds can change the rate of true positive detection. Using conservative thresholds for missing data and minor allele frequency substantially reduces the size of the data set, lessening the power to detect adaptive variants (i.e., simulated true positives) with strong and weak strengths of selection. Regardless, strength of selection was a good predictor for GEA detection, but even some SNPs under strong selection went undetected. False positive rates varied depending on the species and GEA method, and filtering significantly impacted the predictions of adaptive capacity in downstream analyses. We make several recommendations regarding filtering for GEA methods. Ultimately, there is no filtering panacea, but some choices are better than others, depending on the study system, availability of genomic resources, and desired objectives.

An Intron of Invertebrate Microphthalmia Transcription Factor Gene Is Evolved from a Longer Ancestral Sequence

Introns are highly variable in number and size. Sequence simulation is an effective method to elucidate intron evolution patterns. Previously, we have reported that introns are more likely to evolve through mutation-and-deletion (MD) rather than through mutation-and-insertion (MI). In the present study, we further studied evolution models by allowing insertion in the MD model and by allowing deletion in the MI model at various frequencies. It was found that all deletion-biased models with proper parameter settings could generate sequences with attributes matchable to 16 invertebrate introns from the microphthalmia transcription factor gene, whereas all insertion-biased models with any parameter settings failed to generate such sequences. We conclude that the examined invertebrate introns may have evolved from a longer ancestral sequence in a deletion-biased pattern. The constructed models are useful for studying the evolution of introns from other genes and/or from other taxonomic groups. (C++ scripts of all deletion- and insertion-biased models are available upon request.)

The tight genome size of ants: diversity and evolution under ancestral state reconstruction and base composition

The mechanisms and processes driving change and variation in the genome size (GS) are not well known, and only a small set of ant species has been studied. Ants are an ecologically successful insect group present in most distinct ecosystems worldwide. Considering their wide distribution and ecological plasticity in different environmental contexts, we aimed to expand GS estimation within Formicidae to examine distribution patterns and variation in GS and base composition and to reconstruct the ancestral state of this character in an attempt to elucidate the generalized pattern of genomic expansions. Genome size estimates were generated for 99 ant species, including new GS estimates for 91 species of ants, and the mean GS of Formicidae was found to be 0.38 pg. The AT/GC ratio was 62.40/37.60. The phylogenetic reconstruction suggested an ancestral GS of 0.38 pg according to the Bayesian inference/Markov chain Monte Carlo method and 0.37 pg according to maximum likelihood and parsimony methods; significant differences in GS were observed between the subfamilies sampled. Our results suggest that the evolution of GS in Formicidae occurred through loss and accumulation of non-coding regions, mainly transposable elements, and occasionally by whole genome duplication. However, further studies are needed to verify whether these changes in DNA content are related to colonization processes, as suggested at the intraspecific level.

Chromosome-scale genome assembly for the duckweed Spirodela intermedia, integrating cytogenetic maps, PacBio and Oxford Nanopore libraries

Duckweeds are small, free-floating, morphologically highly reduced organisms belonging to the monocot order Alismatales. They display the most rapid growth among flowering plants, vary ~ 14-fold in genome size and comprise five genera. Spirodela is the phylogenetically oldest genus with only two mainly asexually propagating species: S. polyrhiza (2n = 40; 160 Mbp/1C) and S. intermedia (2n = 36; 160 Mbp/1C). This study combined comparative cytogenetics and de novo genome assembly based on PacBio, Illumina and Oxford Nanopore (ON) reads to obtain the first genome reference for S. intermedia and to compare its genomic features with those of the sister species S. polyrhiza. Both species' genomes revealed little more than 20,000 putative protein-coding genes, very low rDNA copy numbers and a low amount of repetitive sequences, mainly Ty3/gypsy retroelements. The detection of a few new small chromosome rearrangements between both Spirodela species refined the karyotype and the chromosomal sequence assignment for S. intermedia.

Evolutionary dynamics of genome size in a radiation of woody plants

PREMISE

Plant genome size ranges widely, providing many opportunities to examine how genome size variation affects plant form and function. We analyzed trends in chromosome number, genome size, and leaf traits for the woody angiosperm clade Viburnum to examine the evolutionary associations, functional implications, and possible drivers of genome size.

METHODS

Chromosome counts and genome size estimates were mapped onto a Viburnum phylogeny to infer the location and frequency of polyploidization events and trends in genome size evolution. Genome size was analyzed with leaf anatomical and physiological data to evaluate the influence of genome size on plant function.

RESULTS

We discovered nine independent polyploidization events, two reductions in base chromosome number, and substantial variation in genome size with a slight trend toward genome size reduction in polyploids. We did not find strong relationships between genome size and the functional and morphological traits that have been highlighted at broader phylogenetic scales.

CONCLUSIONS

Polyploidization events were sometimes associated with rapid radiations, demonstrating that polyploid lineages can be highly successful. Relationships between genome size and plant physiological function observed at broad phylogenetic scales may be largely irrelevant to the evolutionary dynamics of genome size at smaller scales. The view that plants readily tolerate changes in ploidy and genome size, and often do so, appears to apply to Viburnum.

Global analysis of repetitive DNA from unassembled sequence reads using RepeatExplorer2

RepeatExplorer2 is a novel version of a computational pipeline that uses graph-based clustering of next-generation sequencing reads for characterization of repetitive DNA in eukaryotes. The clustering algorithm facilitates repeat identification in any genome by using relatively small quantities of short sequence reads, and additional tools within the pipeline perform automatic annotation and quantification of the identified repeats. The pipeline is integrated into the Galaxy platform, which provides a user-friendly web interface for script execution and documentation of the results. Compared to the original version of the pipeline, RepeatExplorer2 provides automated annotation of transposable elements, identification of tandem repeats and enhanced visualization of analysis results. Here, we present an overview of the RepeatExplorer2 workflow and provide procedures for its application to (i) de novo repeat identification in a single species, (ii) comparative repeat analysis in a set of species, (iii) development of satellite DNA probes for cytogenetic experiments and (iv) identification of centromeric repeats based on ChIP-seq data. Each procedure takes approximately 2 d to complete. RepeatExplorer2 is available at https://repeatexplorer-elixir.cerit-sc.cz .

Is the evolution of carnivory connected with genome size reduction?

PREMISE

As repeatedly shown, the remarkable variation in the genome size of angiosperms can be shaped by extrinsic selective pressures, including nutrient availability. Carnivory has evolved independently in 10 angiosperm clades, but all carnivorous plants share a common affinity to nutrient-poor habitats. As such, carnivory and genome reduction could be responses to the same environmental pressure. Indeed, the smallest genomes among flowering plants are found in the carnivorous family Lentibulariaceae, where a unique mutation in cytochrome c oxidase (COX) is suspected to promote genome miniaturization. Despite these hypotheses, a phylogenetically informed test of genome size and nutrient availability across carnivorous clades has so far been missing.

METHODS

Using linear mixed models, we compared genome sizes of 127 carnivorous plants from 7 diverse angiosperm clades with 1072 of their noncarnivorous relatives. We also tested whether genome size in Lentibulariaceae reflects the presence of the COX mutation.

RESULTS

The genome sizes of carnivorous plants do not differ significantly from those of their noncarnivorous relatives. Based on available data, no significant association between the COX mutation and genome miniaturization could be confirmed, not even when considering polyploidy.

CONCLUSIONS

Carnivory alone does not seem to significantly affect genome size decrease. Plausibly, it might actually counterbalance the effect of nutrient limitation on genome size evolution. The role of the COX mutation in genome miniaturization needs to be evaluated by analysis of a broader data set because current knowledge of its presence across Lentibulariaceae covers less than 10% of the species diversity in this family.

Super-Resolution Microscopy Reveals Diversity of Plant Centromere Architecture

Centromeres are essential for proper chromosome segregation to the daughter cells during mitosis and meiosis. Chromosomes of most eukaryotes studied so far have regional centromeres that form primary constrictions on metaphase chromosomes. These monocentric chromosomes vary from point centromeres to so-called "meta-polycentromeres", with multiple centromere domains in an extended primary constriction, as identified in Pisum and Lathyrus species. However, in various animal and plant lineages centromeres are distributed along almost the entire chromosome length. Therefore, they are called holocentromeres. In holocentric plants, centromere-specific proteins, at which spindle fibers usually attach, are arranged contiguously (line-like), in clusters along the chromosomes or in bands. Here, we summarize findings of ultrastructural investigations using immunolabeling with centromere-specific antibodies and super-resolution microscopy to demonstrate the structural diversity of plant centromeres. A classification of the different centromere types has been suggested based on the distribution of spindle attachment sites. Based on these findings we discuss the possible evolution and advantages of holocentricity, and potential strategies to segregate holocentric chromosomes correctly.

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.

Kmasker plants - a tool for assessing complex sequence space in plant species

Many plant genomes display high levels of repetitive sequences. The assembly of these complex genomes using short high-throughput sequence reads is still a challenging task. Underestimation or disregard of repeat complexity in these datasets can easily misguide downstream analysis. Detection of repetitive regions by k-mer counting methods has proved to be reliable. Easy-to-use applications utilizing k-mer counting are in high demand, especially in the domain of plants. We present Kmasker plants, a tool that uses k-mer count information as an assistant throughout the analytical workflow of genome data that is provided as a command-line and web-based solution. Beside its core competence to screen and mask repetitive sequences, we have integrated features that enable comparative studies between different cultivars or closely related species and methods that estimate target specificity of guide RNAs for application of site-directed mutagenesis using Cas9 endonuclease. In addition, we have set up a web service for Kmasker plants that maintains pre-computed indices for 10 of the economically most important cultivated plants. Source code for Kmasker plants has been made publically available at https://github.com/tschmutzer/kmasker. The web service is accessible at https://kmasker.ipk-gatersleben.de.

Differential Genome Size and Repetitive DNA Evolution in Diploid Species of Melampodium sect. Melampodium (Asteraceae)

Plant genomes vary greatly in composition and size mainly due to the diversity of repetitive DNAs and the inherent propensity for their amplification and removal from the host genome. Most studies addressing repeatome dynamics focus on model organisms, whereas few provide comprehensive investigations across the genomes of related taxa. Herein, we analyze the evolution of repeats of the 13 species in Melampodium sect. Melampodium, representing all but two of its diploid taxa, in a phylogenetic context. The investigated genomes range in size from 0.49 to 2.27 pg/1C (ca. 4.5-fold variation), despite having the same base chromosome number (x = 10) and very strong phylogenetic affinities. Phylogenetic analysis performed in BEAST and ancestral genome size reconstruction revealed mixed patterns of genome size increases and decreases across the group. High-throughput genome skimming and the RepeatExplorer pipeline were utilized to determine the repeat families responsible for the differences in observed genome sizes. Patterns of repeat evolution were found to be highly correlated with phylogenetic position, namely taxonomic series circumscription. Major differences found were in the abundances of the SIRE (Ty1-copia), Athila (Ty3-gypsy), and CACTA (DNA transposon) lineages. Additionally, several satellite DNA families were found to be highly group-specific, although their overall contribution to genome size variation was relatively small. Evolutionary changes in repetitive DNA composition and genome size were complex, with independent patterns of genome up- and downsizing throughout the evolution of the analyzed diploids. A model-based analysis of genome size and repetitive DNA composition revealed evidence for strong phylogenetic signal and differential evolutionary rates of major lineages of repeats in the diploid genomes.

The Terrestrial Carnivorous Plant Utricularia reniformis Sheds Light on Environmental and Life-Form Genome Plasticity

Utricularia belongs to Lentibulariaceae, a widespread family of carnivorous plants that possess ultra-small and highly dynamic nuclear genomes. It has been shown that the Lentibulariaceae genomes have been shaped by transposable elements expansion and loss, and multiple rounds of whole-genome duplications (WGD), making the family a platform for evolutionary and comparative genomics studies. To explore the evolution of Utricularia, we estimated the chromosome number and genome size, as well as sequenced the terrestrial bladderwort Utricularia reniformis (2n = 40, 1C = 317.1-Mpb). Here, we report a high quality 304 Mb draft genome, with a scaffold NG50 of 466-Kb, a BUSCO completeness of 87.8%, and 42,582 predicted genes. Compared to the smaller and aquatic U. gibba genome (101 Mb) that has a 32% repetitive sequence, the U. reniformis genome is highly repetitive (56%). The structural differences between the two genomes are the result of distinct fractionation and rearrangements after WGD, and massive proliferation of LTR-retrotransposons. Moreover, GO enrichment analyses suggest an ongoing gene birth-death-innovation process occurring among the tandem duplicated genes, shaping the evolution of carnivory-associated functions. We also identified unique patterns of developmentally related genes that support the terrestrial life-form and body plan of U. reniformis. Collectively, our results provided additional insights into the evolution of the plastic and specialized Lentibulariaceae genomes.

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.

Small, but surprisingly repetitive genomes: transposon expansion and not polyploidy has driven a doubling in genome size in a metazoan species complex

BACKGROUND

The causes and consequences of genome size variation across Eukaryotes, which spans five orders of magnitude, have been hotly debated since before the advent of genome sequencing. Previous studies have mostly examined variation among larger taxonomic units (e.g., orders, or genera), while comparisons among closely related species are rare. Rotifers of the Brachionus plicatilis species complex exhibit a seven-fold variation in genome size and thus represent a unique opportunity to study such changes on a relatively short evolutionary timescale. Here, we sequenced and analysed the genomes of four species of this complex with nuclear DNA contents spanning 110-422 Mbp. To establish the likely mechanisms of genome size change, we analysed both sequencing read libraries and assemblies for signatures of polyploidy and repetitive element content. We also compared these genomes to that of B. calyciflorus, the closest relative with a sequenced genome (293 Mbp nuclear DNA content).

RESULTS

Despite the very large differences in genome size, we saw no evidence of ploidy level changes across the B. plicatilis complex. However, repetitive element content explained a large portion of genome size variation (at least 54%). The species with the largest genome, B. asplanchnoidis, has a strikingly high 44% repetitive element content, while the smaller B. plicatilis genomes contain between 14 and 25% repetitive elements. According to our analyses, the B. calyciflorus genome contains 39% repetitive elements, which is substantially higher than previously reported (21%), and suggests that high repetitive element load could be widespread in monogonont rotifers.

CONCLUSIONS

Even though the genome sizes of these species are at the low end of the metazoan spectrum, their genomes contain substantial amounts of repetitive elements. Polyploidy does not appear to play a role in genome size variations in these species, and these variations can be mostly explained by changes in repetitive element content. This contradicts the naïve expectation that small genomes are streamlined, or less complex, and that large variations in nuclear DNA content between closely related species are due to polyploidy.

Genome Size Diversity and Its Impact on the Evolution of Land Plants

Genome size is a biodiversity trait that shows staggering diversity across eukaryotes, varying over 64,000-fold. Of all major taxonomic groups, land plants stand out due to their staggering genome size diversity, ranging ca. 2400-fold. As our understanding of the implications and significance of this remarkable genome size diversity in land plants grows, it is becoming increasingly evident that this trait plays not only an important role in shaping the evolution of plant genomes, but also in influencing plant community assemblages at the ecosystem level. Recent advances and improvements in novel sequencing technologies, as well as analytical tools, make it possible to gain critical insights into the genomic and epigenetic mechanisms underpinning genome size changes. In this review we provide an overview of our current understanding of genome size diversity across the different land plant groups, its implications on the biology of the genome and what future directions need to be addressed to fill key knowledge gaps.

Comparative Genomics of an Unusual Biogeographic Disjunction in the Cotton Tribe (Gossypieae) Yields Insights into Genome Downsizing

Long-distance insular dispersal is associated with divergence and speciation because of founder effects and strong genetic drift. The cotton tribe (Gossypieae) has experienced multiple transoceanic dispersals, generating an aggregate geographic range that encompasses much of the tropics and subtropics worldwide. Two genera in the Gossypieae, Kokia and Gossypioides, exhibit a remarkable geographic disjunction, being restricted to the Hawaiian Islands and Madagascar/East Africa, respectively. We assembled and use de novo genome sequences to address questions regarding the divergence of these two genera from each other and from their sister-group, Gossypium. In addition, we explore processes underlying the genome downsizing that characterizes Kokia and Gossypioides relative to other genera in the tribe. Using 13,000 gene orthologs and synonymous substitution rates, we show that the two disjuncts last shared a common ancestor ∼5 Ma, or half as long ago as their divergence from Gossypium. We report relative stasis in the transposable element fraction. In comparison to Gossypium, there is loss of ∼30% of the gene content in the two disjunct genera and a history of genome-wide accumulation of deletions. In both genera, there is a genome-wide bias toward deletions over insertions, and the number of gene losses exceeds the number of gains by ∼2- to 4-fold. The genomic analyses presented here elucidate genomic consequences of the demographic and biogeographic history of these closest relatives of Gossypium, and enhance their value as phylogenetic outgroups.

Evolutionary modes of emergence of short interspersed nuclear element (SINE) families in grasses

Short interspersed nuclear elements (SINEs) are non-autonomous transposable elements which are propagated by retrotransposition and constitute an inherent part of the genome of most eukaryotic species. Knowledge of heterogeneous and highly abundant SINEs is crucial for de novo (or improvement of) annotation of whole genome sequences. We scanned Poaceae genome sequences of six important cereals (Oryza sativa, Triticum aestivum, Hordeum vulgare, Panicum virgatum, Sorghum bicolor, Zea mays) and Brachypodium distachyon to examine the diversity and evolution of SINE populations. We comparatively analyzed the structural features, distribution, evolutionary relation and abundance of 32 SINE families and subfamilies within grasses, comprising 11 052 individual copies. The investigation of activity profiles within the Poaceae provides insights into their species-specific diversification and amplification. We found that Poaceae SINEs (PoaS) fall into two length categories: simple SINEs of up to 180 bp and dimeric SINEs larger than 240 bp. Detailed analysis at the nucleotide level revealed that multimerization of related and unrelated SINE copies is an important evolutionary mechanism of SINE formation. We conclude that PoaS families diversify by massive reshuffling between SINE families, likely caused by insertion of truncated copies, and provide a model for this evolutionary scenario. Twenty-eight of 32 PoaS families and subfamilies show significant conservation, in particular either in the 5' or 3' regions, across Poaceae species and share large sequence stretches with one or more other PoaS families.

TGTT and AACA: two transcriptionally active LTR retrotransposon subfamilies with a specific LTR structure and horizontal transfer in four Rosaceae species

BACKGROUND

Long terminal repeat retrotransposons (LTR-RTs) are major components of plant genomes. Common LTR-RTs contain the palindromic dinucleotide 5'-'TG'-'CA'-3' motif at the ends. Thus, further analyses of non-canonical LTR-RTs with non-palindromic motifs will enhance our understanding of their structures and evolutionary history.

RESULTS

Here, we report two new LTR-RT subfamilies (TGTT and AACA) with atypical dinucleotide ends of 5'-'TG'-'TT'-3', and 5'-'AA'-'CA'-3' in pear, apple, peach and mei. In total, 91 intact LTR-RTs were identified and classified into four TGTT and four AACA families. A structural annotation analysis showed that the four TGTT families, together with AACA1 and AACA2, belong to the Copia-like superfamily, whereas AACA3 and AACA4 appeared to be TRIM elements. The average amplification time frames for the eight families ranged from 0.05 to 2.32 million years. Phylogenetics coupled with sequence analyses revealed that the TGTT1 elements of peach were horizontally transferred from apple. In addition, 32 elements from two TGTT and three AACA families had detectable transcriptional activation, and a qRT-PCR analysis indicated that their expression levels varied dramatically in different species, organs and stress treatments.

CONCLUSIONS

Two novel LTR-RT subfamilies that terminated with non-palindromic dinucleotides at the ends of their LTRs were identified in four Rosaceae species, and a deep analysis showed their recent activity, horizontal transfer and varied transcriptional levels in different species, organs and stress treatments. This work enhances our understanding of the structural variation and evolutionary history of LTR-RTs in plants and also provides a valuable resource for future investigations of LTR-RTs having specific structures in other species.

How Next-Generation Sequencing Has Aided Our Understanding of the Sequence Composition and Origin of B Chromosomes

Accessory, supernumerary, or-most simply-B chromosomes, are found in many eukaryotic karyotypes. These small chromosomes do not follow the usual pattern of segregation, but rather are transmitted in a higher than expected frequency. As increasingly being demonstrated by next-generation sequencing (NGS), their structure comprises fragments of standard (A) chromosomes, although in some plant species, their sequence also includes contributions from organellar genomes. Transcriptomic analyses of various animal and plant species have revealed that, contrary to what used to be the common belief, some of the B chromosome DNA is protein-encoding. This review summarizes the progress in understanding B chromosome biology enabled by the application of next-generation sequencing technology and state-of-the-art bioinformatics. In particular, a contrast is drawn between a direct sequencing approach and a strategy based on a comparative genomics as alternative routes that can be taken towards the identification of B chromosome sequences.

Endogenous sequence patterns predispose the repair modes of CRISPR/Cas9-induced DNA double-stranded breaks in Arabidopsis thaliana

The possibility to predict the outcome of targeted DNA double-stranded break (DSB) repair would be desirable for genome editing. Furthermore the consequences of mis-repair of potentially cell-lethal DSBs and the underlying pathways are not yet fully understood. Here we study the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-induced mutation spectra at three selected endogenous loci in Arabidopsis thaliana by deep sequencing of long amplicon libraries. Notably, we found sequence-dependent genomic features that affected the DNA repair outcome. Deletions of 1-bp to <1000-bp size and/or very short insertions, deletions >1 kbp (all due to NHEJ) and deletions combined with insertions between 5-bp to >100 bp [caused by a synthesis-dependent strand annealing (SDSA)-like mechanism] occurred most frequently at all three loci. The appearance of single-stranded annealing events depends on the presence and distance between repeats flanking the DSB. The frequency and size of insertions is increased if a sequence with high similarity to the target site was available in cis. Most deletions were linked to pre-existing microhomology. Deletion and/or insertion mutations were blunt-end ligated or via de novo generated microhomology. While most mutation types and, to some degree, their predictability are comparable with animal systems, the broad range of deletion mutations seems to be a peculiar feature of the plant A. thaliana.

Rapid amplification of four retrotransposon families promoted speciation and genome size expansion in the genus Panax

Genome duplication and repeat multiplication contribute to genome evolution in plants. Our previous work identified a recent allotetraploidization event and five high-copy LTR retrotransposon (LTR-RT) families PgDel, PgTat, PgAthila, PgTork, and PgOryco in Panax ginseng. Here, using whole-genome sequences, we quantified major repeats in five Panax species and investigated their role in genome evolution. The diploids P. japonicus, P. vietnamensis, and P. notoginseng and the tetraploids P. ginseng and P. quinquefolius were analyzed alongside their relative Aralia elata. These species possess 0.8-4.9 Gb haploid genomes. The PgDel, PgTat, PgAthila, and PgTork LTR-RT superfamilies accounted for 39-52% of the Panax species genomes and 17% of the A. elata genome. PgDel included six subfamily members, each with a distinct genome distribution. In particular, the PgDel1 subfamily occupied 23-35% of the Panax genomes and accounted for much of their genome size variation. PgDel1 occupied 22.6% (0.8 Gb of 3.6 Gb) and 34.5% (1.7 Gb of 4.9 Gb) of the P. ginseng and P. quinquefolius genomes, respectively. Our findings indicate that the P. quinquefolius genome may have expanded due to rapid PgDel1 amplification over the last million years as a result of environmental adaptation following migration from Asia to North America.

Deletion-bias in DNA double-strand break repair differentially contributes to plant genome shrinkage

In order to prevent genome instability, cells need to be protected by a number of repair mechanisms, including DNA double-strand break (DSB) repair. The extent to which DSB repair, biased towards deletions or insertions, contributes to evolutionary diversification of genome size is still under debate. We analyzed mutation spectra in Arabidopsis thaliana and in barley (Hordeum vulgare) by PacBio sequencing of three DSB-targeted loci each, uncovering repair via gene conversion, single strand annealing (SSA) or nonhomologous end-joining (NHEJ). Furthermore, phylogenomic comparisons between A. thaliana and two related species were used to detect naturally occurring deletions during Arabidopsis evolution. Arabidopsis thaliana revealed significantly more and larger deletions after DSB repair than barley, and barley displayed more and larger insertions. Arabidopsis displayed a clear net loss of DNA after DSB repair, mainly via SSA and NHEJ. Barley revealed a very weak net loss of DNA, apparently due to less active break-end resection and easier copying of template sequences into breaks. Comparative phylogenomics revealed several footprints of SSA in the A. thaliana genome. Quantitative assessment of DNA gain and loss through DSB repair processes suggests deletion-biased DSB repair causing ongoing genome shrinking in A. thaliana, whereas genome size in barley remains nearly constant.

Evolution of genome size and genomic GC content in carnivorous holokinetics (Droseraceae)

BACKGROUND AND AIMS

Studies in the carnivorous family Lentibulariaceae in the last years resulted in the discovery of the smallest plant genomes and an unusual pattern of genomic GC content evolution. However, scarcity of genomic data in other carnivorous clades still prevents a generalization of the observed patterns. Here the aim was to fill this gap by mapping genome evolution in the second largest carnivorous family, Droseraceae, where this evolution may be affected by chromosomal holokinetism in Drosera METHODS: The genome size and genomic GC content of 71 Droseraceae species were measured by flow cytometry. A dated phylogeny was constructed, and the evolution of both genomic parameters and their relationship to species climatic niches were tested using phylogeny-based statistics.

KEY RESULTS

The 2C genome size of Droseraceae varied between 488 and 10 927 Mbp, and the GC content ranged between 37·1 and 44·7 %. The genome sizes and genomic GC content of carnivorous and holocentric species did not differ from those of their non-carnivorous and monocentric relatives. The genomic GC content positively correlated with genome size and annual temperature fluctuations. The genome size and chromosome numbers were inversely correlated in the Australian clade of Drosera CONCLUSIONS: Our results indicate that neither carnivory (nutrient scarcity) nor the holokinetism have a prominent effect on size and DNA base composition of Droseraceae genomes. However, the holokinetic drive seems to affect karyotype evolution in one of the major clades of Drosera Our survey confirmed that the evolution of GC content is tightly connected with the evolution of genome size and also with environmental conditions.

The Power of CRISPR-Cas9-Induced Genome Editing to Speed Up Plant Breeding

Genome editing with engineered nucleases enabling site-directed sequence modifications bears a great potential for advanced plant breeding and crop protection. Remarkably, the RNA-guided endonuclease technology (RGEN) based on the clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9) is an extremely powerful and easy tool that revolutionizes both basic research and plant breeding. Here, we review the major technical advances and recent applications of the CRISPR-Cas9 system for manipulation of model and crop plant genomes. We also discuss the future prospects of this technology in molecular plant breeding.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

Recurrent evolution of heat-responsiveness in Brassicaceae COPIA elements

BACKGROUND

The mobilization of transposable elements (TEs) is suppressed by host genome defense mechanisms. Recent studies showed that the cis-regulatory region of Arabidopsis thaliana COPIA78/ONSEN retrotransposons contains heat-responsive elements (HREs), which cause their activation during heat stress. However, it remains unknown whether this is a common and potentially conserved trait and how it has evolved.

RESULTS

We show that ONSEN, COPIA37, TERESTRA, and ROMANIAT5 are the major families of heat-responsive TEs in A. lyrata and A. thaliana. Heat-responsiveness of COPIA families is correlated with the presence of putative high affinity heat shock factor binding HREs within their long terminal repeats in seven Brassicaceae species. The strong HRE of ONSEN is conserved over millions of years and has evolved by duplication of a proto-HRE sequence, which was already present early in the evolution of the Brassicaceae. However, HREs of most families are species-specific, and in Boechera stricta, the ONSEN HRE accumulated mutations and lost heat-responsiveness.

CONCLUSIONS

Gain of HREs does not always provide an ultimate selective advantage for TEs, but may increase the probability of their long-term survival during the co-evolution of hosts and genomic parasites.

Genome Stability and Evolution: Attempting a Holistic View

The reason why the DNA content, chromosome number and shape, and gene content of eukaryotic genomes vary independently remains a matter of speculation. The same is true for the questions of whether there is a general tendency for increase or decrease of genome size and chromosome number and whether genome size and/or chromosome number have an adaptive value and, if so, what this value is. Here we assume that three strategies of genome evolution (shrinkage, expansion, and equilibrium) have developed to find the optimal balance between genomic stability and plasticity. We suggest various modes of DNA double-strand break (DSB) repair in combination with whole-genome duplication (WGD) and dysploid chromosome number alteration to explain the different strategies of genome size and karyotype evolution.

Genome sequencing of Metrosideros polymorpha (Myrtaceae), a dominant species in various habitats in the Hawaiian Islands with remarkable phenotypic variations

Whole genome sequences, which can be provided even for non-model organisms owing to high-throughput sequencers, are valuable in enhancing the understanding of adaptive evolution. Metrosideros polymorpha, a tree species endemic to the Hawaiian Islands, occupies a wide range of ecological habitats and shows remarkable polymorphism in phenotypes among/within populations. The biological functions of genetic variations observed within this species could provide significant insights into the adaptive radiation found in a single species. Here de novo assembled genome sequences of M. polymorpha are presented to reveal basic genomic parameters about this species and to develop our knowledge of ecological divergences. The assembly yielded 304-Mbp genome sequences, half of which were covered by 19 scaffolds with >5 Mbp, and contained 30 K protein-coding genes. Demographic history inferred from the genome-wide heterozygosity indicated that this species experienced a dramatic rise and fall in the effective population size, possibly owing to past geographic or climatic changes in the Hawaiian Islands. This M. polymorpha genome assembly represents a high-quality genome resource useful for future functional analyses of both intra- and interspecies genetic variations or comparative genomics.

Transposable element evolution in the allotetraploid Capsella bursa-pastoris

PREMISE OF THE STUDY

Shifts in ploidy affect the evolutionary dynamics of genomes in a myriad of ways. Population genetic theory predicts that transposable element (TE) proliferation may follow because the genomewide efficacy of selection should be reduced and the increase in gene copies may mask the deleterious effects of TE insertions. Moreover, in allopolyploids, TEs may further accumulate because of hybrid breakdown of TE silencing. However, to date the evidence of TE proliferation following an increase in ploidy is mixed, and the relative importance of relaxed selection vs. silencing breakdown remains unclear.

METHODS

We used high-coverage whole-genome sequence data to evaluate the abundance, genomic distribution, and population frequencies of TEs in the self-fertilizing recent allotetraploid Capsella bursa-pastoris (Brassicaceae). We then compared the C. bursa-pastoris TE profile with that of its two parental diploid species, outcrossing C. grandiflora and self-fertilizing C. orientalis.

KEY RESULTS

We found no evidence that C. bursa-pastoris has experienced a large genomewide proliferation of TEs relative to its parental species. However, when centromeric regions are excluded, we found evidence of significantly higher abundance of retrotransposons in C. bursa-pastoris along the gene-rich chromosome arms compared with C. grandiflora and C. orientalis.

CONCLUSIONS

The lack of a genomewide effect of allopolyploidy on TE abundance, combined with the increases TE abundance in gene-rich regions, suggests that relaxed selection rather than hybrid breakdown of host silencing explains the TE accumulation in C. bursa-pastoris.

Chromosome identification for the carnivorous plant Genlisea margaretae

Genlisea margaretae, subgenus Genlisea, section Recurvatae (184 Mbp/1C), belongs to a plant genus with a 25-fold genome size difference and an extreme genome plasticity. Its 19 chromosome pairs could be distinguished individually by an approach combining optimized probe pooling and consecutive rounds of multicolor fluorescence in situ hybridization (mcFISH) with bacterial artificial chromosomes (BACs) selected for repeat-free inserts. Fifty-one BACs were assigned to 18 chromosome pairs. They provide a tool for future assignment of genomic sequence contigs to distinct chromosomes as well as for identification of homeologous chromosome regions in other species of the carnivorous Lentibulariaceae family, and potentially of chromosome rearrangements, in cases where more than one BAC per chromosome pair was identified.

Centromere and telomere sequence alterations reflect the rapid genome evolution within the carnivorous plant genus Genlisea

Linear chromosomes of eukaryotic organisms invariably possess centromeres and telomeres to ensure proper chromosome segregation during nuclear divisions and to protect the chromosome ends from deterioration and fusion, respectively. While centromeric sequences may differ between species, with arrays of tandemly repeated sequences and retrotransposons being the most abundant sequence types in plant centromeres, telomeric sequences are usually highly conserved among plants and other organisms. The genome size of the carnivorous genus Genlisea (Lentibulariaceae) is highly variable. Here we study evolutionary sequence plasticity of these chromosomal domains at an intrageneric level. We show that Genlisea nigrocaulis (1C = 86 Mbp; 2n = 40) and G. hispidula (1C = 1550 Mbp; 2n = 40) differ as to their DNA composition at centromeres and telomeres. G. nigrocaulis and its close relative G. pygmaea revealed mainly 161 bp tandem repeats, while G. hispidula and its close relative G. subglabra displayed a combination of four retroelements at centromeric positions. G. nigrocaulis and G. pygmaea chromosome ends are characterized by the Arabidopsis-type telomeric repeats (TTTAGGG); G. hispidula and G. subglabra instead revealed two intermingled sequence variants (TTCAGG and TTTCAGG). These differences in centromeric and, surprisingly, also in telomeric DNA sequences, uncovered between groups with on average a > 9-fold genome size difference, emphasize the fast genome evolution within this genus. Such intrageneric evolutionary alteration of telomeric repeats with cytosine in the guanine-rich strand, not yet known for plants, might impact the epigenetic telomere chromatin modification.