Genome size variation at constant chromosome number is not correlated with repetitive DNA dynamism in Anacyclus (Asteraceae)

Expansions and contractions of repetitive DNA elements reveal contrasting evolutionary responses to the polyploid genome shock hypothesis in Brachypodium model grasses

Brachypodium grass species have been selected as model plants for functional genomics of grass crops, and to elucidate the origins of allopolyploidy and perenniality in monocots, due to their small genome sizes and feasibility of cultivation. However, genome sizes differ greatly between diploid or polyploid Brachypodium lineages. We have used genome skimming sequencing data to uncover the composition, abundance, and phylogenetic value of repetitive elements in 44 representatives of the major Brachypodium lineages and cytotypes. We also aimed to test the possible mechanisms and consequences of the "polyploid genome shock hypothesis" (PGSH) under three different evolutionary scenarios of variation in repeats and genome sizes of Brachypodium allopolyploids. Our data indicated that the proportion of the genome covered by the repeatome in the Brachypodium species showed a 3.3-fold difference between the highest content of B. mexicanum-4x (67.97%) and the lowest of B. stacei-2x (20.77%), and that changes in the sizes of their genomes were a consequence of gains or losses in their repeat elements. LTR-Retand and Tekay retrotransposons were the most frequent repeat elements in the Brachypodium genomes, while Ogre retrotransposons were found exclusively in B. mexicanum. The repeatome phylogenetic network showed a high topological congruence with plastome and nuclear rDNA and transcriptome trees, differentiating the ancestral outcore lineages from the recently evolved core-perennial lineages. The 5S rDNA graph topologies had a strong match with the ploidy levels and nature of the subgenomes of the Brachypodium polyploids. The core-perennial B. sylvaticum presents a large repeatome and characteristics of a potential post-polyploid diploidized origin. Our study evidenced that expansions and contractions in the repeatome were responsible for the three contrasting responses to the PGSH. The exacerbated genome expansion of the ancestral allotetraploid B. mexicanum was a consequence of chromosome-wide proliferation of TEs and not of WGD, the additive repeatome pattern of young allotetraploid B. hybridum of stabilized post-WGD genome evolution, and the genomecontraction of recent core-perennials polyploids (B. pinnatum, B. phoenicoides) of repeat losses through recombination of these highly hybridizing lineages. Our analyses have contributed to unraveling the evolution of the repeatome and the genome size variation in model Brachypodium grasses.

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae)

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Unravelling plant diversification: Intraspecific genetic differentiation in hybridizing Anacyclus species in the western Mediterranean Basin

PREMISE

The interfertile species Anacyclus clavatus, A. homogamos, and A. valentinus represent a plant complex coexisting in large anthropic areas of the western Mediterranean Basin with phenotypically mixed populations exhibiting a great floral variation. The goal of this study was to estimate the genetic identity of each species, to infer the role of hybridization in the observed phenotypic diversity, and to explore the effect of climate on the geographic distribution of species and genetic clusters.

METHODS

We used eight nuclear microsatellites to genotype 585 individuals from 31 populations of three Anacyclus species for population genetic analyses by using clustering algorithms based on Bayesian models and ordination methods. In addition, we used ecological niche models and niche overlap analyses for both the species and genetic clusters. We used an expanded data set, including 721 individuals from 129 populations for ecological niche models of the genetic clusters.

RESULTS

We found a clear correspondence between species and genetic clusters, except for A. clavatus that included up to three genetic clusters. We detected individuals with admixed genetic ancestry in A. clavatus and in mixed populations. Ecological niche models predicted similar distributions for species and genetic clusters. For the two specific genetic clusters of A. clavatus, ecological niche models predicted remarkably different areas.

CONCLUSIONS

Gene flow between Anacyclus species likely explains phenotypic diversity in contact areas. In addition, we suggest that introgression could be involved in the origin of one of the two A. clavatus genetic clusters, which also showed ecological differentiation.

Phylogenomics and genome size evolution in Amomum s. s. (Zingiberaceae): Comparison of traditional and modern sequencing methods

BACKGROUND AND AIMS

A targeted enrichment NGS approach was used to construct the phylogeny of Amomum Roxb. (Zingiberaceae). Phylogenies based on hundreds of nuclear genes, the whole plastome and the rDNA cistron were compared with an ITS-based phylogeny. Trends in genome size (GS) evolution were examined, chromosomes were counted and the geographical distribution of phylogenetic lineages was evaluated.

METHODS

In total, 92 accessions of 54 species were analysed. ITS was obtained for 79 accessions, 37 accessions were processed with Hyb-Seq and sequences from 449 nuclear genes, the whole cpDNA, and the rDNA cistron were analysed using concatenation, coalescence and supertree approaches. The evolution of absolute GS was analysed in a phylogenetic and geographical context. The chromosome numbers of 12 accessions were counted.

KEY RESULTS

Four groups were recognised in all datasets though their mutual relationships differ among datasets. While group A (A. subulatum and A. petaloideum) is basal to the remaining groups in the nuclear gene phylogeny, in the cpDNA topology it is sister to group B (A. repoeense and related species) and, in the ITS topology, it is sister to group D (the Elettariopsis lineage). The former Elettariopsis makes a monophyletic group. There is an increasing trend in GS during evolution. The largest GS values were found in group D in two tetraploid taxa, A. cinnamomeum and A. aff. biphyllum (both 2n = 96 chromosomes). The rest varied in GS (2C = 3.54-8.78 pg) with a constant chromosome number 2n = 48. There is a weak connection between phylogeny, GS and geography in Amomum.

CONCLUSIONS

Amomum consists of four groups, and the former Elettariopsis is monophyletic. Species in this group have the largest GS. Two polyploids were found and GS greatly varied in the rest of Amomum.

The relationship between transposable elements and ecological niches in the Greater Cape Floristic Region: A study on the genus Pteronia (Asteraceae)

Non-coding repetitive DNA (repeatome) is an active part of the nuclear genome, involved in its structure, evolution and function. It is dominated by transposable elements (TEs) and satellite DNA and is prone to the most rapid changes over time. The TEs activity presumably causes the global genome reorganization and may play an adaptive or regulatory role in response to environmental challenges. This assumption is applied here for the first time to plants from the Cape Floristic hotspot to determine whether changes in repetitive DNA are related to responses to a harsh, but extremely species-rich environment. The genus Pteronia (Asteraceae) serves as a suitable model group because it shows considerable variation in genome size at the diploid level and has high and nearly equal levels of endemism in the two main Cape biomes, Fynbos and Succulent Karoo. First, we constructed a phylogeny based on multiple low-copy genes that served as a phylogenetic framework for detecting quantitative and qualitative changes in the repeatome. Second, we performed a comparative analysis of the environments of two groups of Pteronia differing in their TEs bursts. Our results suggest that the environmental transition from the Succulent Karoo to the Fynbos is accompanied by TEs burst, which is likely also driving phylogenetic divergence. We thus hypothesize that analysis of rapidly evolving repeatome could serve as an important proxy for determining the molecular basis of lineage divergence in rapidly radiating groups.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Evolutionary Dynamics of the Repeatome Explains Contrasting Differences in Genome Sizes and Hybrid and Polyploid Origins of Grass Loliinae Lineages

The repeatome is composed of diverse families of repetitive DNA that keep signatures on the historical events that shaped the evolution of their hosting species. The cold seasonal Loliinae subtribe includes worldwide distributed taxa, some of which are the most important forage and lawn species (fescues and ray-grasses). The Loliinae are prone to hybridization and polyploidization. It has been observed a striking two-fold difference in genome size between the broad-leaved (BL) and fine-leaved (FL) Loliinae diploids and a general trend of genome reduction of some high polyploids. We have used genome skimming data to uncover the composition, abundance, and potential phylogenetic signal of repetitive elements across 47 representatives of the main Loliinae lineages. Independent and comparative analyses of repetitive sequences and of 5S rDNA loci were performed for all taxa under study and for four evolutionary Loliinae groups [Loliinae, Broad-leaved (BL), Fine-leaved (FL), and Schedonorus lineages]. Our data showed that the proportion of the genome covered by the repeatome in the Loliinae species was relatively high (average ∼ 51.8%), ranging from high percentages in some diploids (68.7%) to low percentages in some high-polyploids (30.7%), and that changes in their genome sizes were likely caused by gains or losses in their repeat elements. Ty3-gypsy Retand and Ty1-copia Angela retrotransposons were the most frequent repeat families in the Loliinae although the relatively more conservative Angela repeats presented the highest correlation of repeat content with genome size variation and the highest phylogenetic signal of the whole repeatome. By contrast, Athila retrotransposons presented evidence of recent proliferations almost exclusively in the Lolium clade. The repeatome evolutionary networks showed an overall topological congruence with the nuclear 35S rDNA phylogeny and a geographic-based structure for some lineages. The evolution of the Loliinae repeatome suggests a plausible scenario of recurrent allopolyploidizations followed by diploidizations that generated the large genome sizes of BL diploids as well as large genomic rearrangements in highly hybridogenous lineages that caused massive repeatome and genome contractions in the Schedonorus and Aulaxyper polyploids. Our study has contributed to disentangling the impact of the repeatome dynamics on the genome diversification and evolution of the Loliinae grasses.

Characterisation of LTR-Retrotransposons of Stevia rebaudiana and Their Use for the Analysis of Genetic Variability

Stevia rebaudiana is one of the most important crops belonging to the Asteraceae family. Stevia is cultivated all over the world as it represents a valid natural alternative to artificial sweeteners thanks to its leaves, which produce steviol glycosides that have high sweetening power and reduced caloric value. In this work, the stevia genome sequence was used to isolate and characterise full-length long-terminal repeat retrotransposons (LTR-REs), which account for more than half of the genome. The Gypsy retrotransposons were twice as abundant as the Copia ones. A disproportionate abundance of elements belonging to the Chromovirus/Tekay lineage was observed among the Gypsy elements. Only the SIRE and Angela lineages represented significant portions of the genome among the Copia elements. The dynamics with which LTR-REs colonised the stevia genome were also estimated; all isolated full-length elements turned out to be relatively young, with a proliferation peak around 1–2 million years ago. However, a different analysis conducted by comparing sequences encoding retrotranscriptase showed the occurrence of an older period in which there was a lot of LTR-RE proliferation. Finally, a group of isolated full-length elements belonging to the lineage Angela was used to analyse the genetic variability in 25 accessions of S. rebaudiana using the Inter-Retrotransposon Amplified Polymorphism (IRAP) protocol. The obtained fingerprints highlighted a high degree of genetic variability and were used to study the genomic structures of the different accessions. It was hypothesised that there are four ancestral subpopulations at the root of the analysed accessions, which all turned out to be admixed. Overall, these data may be useful for genome sequence annotations and for evaluating genetic variability in this species, which may be useful in stevia breeding.

Comparative Analysis of Transposable Elements and the Identification of Candidate Centromeric Elements in the Prunus Subgenus Cerasus and Its Relatives

The subgenus Cerasus and its relatives include many crucial economic drupe fruits and ornamental plants. Repetitive elements make up a large part of complex genomes, and some of them play an important role in gene regulation that can affect phenotypic variation. However, the variation in their genomes remains poorly understood. This work conducted a comprehensive repetitive sequence identification across the draft genomes of eight taxa of the genus Prunus, including four of the Prunus subgenus Cerasus (Prunus pseudocerasus, P. avium, P. yedoensis and P. × yedoensis) as well as congeneric species (Prunus salicina, P. armeniaca, P. dulcis and P. persica). Annotation results showed high proportions of transposable elements in their genomes, ranging from 52.28% (P. armeniaca) to 61.86% (P. pseudocerasus). The most notable differences in the contents of long terminal repeat retrotransposons (LTR-RTs) and tandem repeats (TRs) were confirmed with de novo identification based on the structure of each genome, which significantly contributed to their genome size variation, especially in P. avium and P.salicina. Sequence comparisons showed many similar LTR-RTs closely related to their phylogenetic relationships, and a highly similar monomer unit of the TR sequence was conserved among species. Additionally, the predicted centromere-associated sequence was located in centromeric regions with FISH in the 12 taxa of Prunus. It presented significantly different signal intensities, even within the diverse interindividual phenotypes for Prunus tomentosa. This study provides insight into the LTR-RT and TR variation within Prunus and increases our knowledge about its role in genome evolution.

Genomes shed light on the evolution of Begonia, a mega-diverse genus

Clarifying the evolutionary processes underlying species diversification and adaptation is a key focus of evolutionary biology. Begonia (Begoniaceae) is one of the most species-rich angiosperm genera with c. 2000 species, most of which are shade-adapted. Here, we present chromosome-scale genome assemblies for four species of Begonia (B. loranthoides, B. masoniana, B. darthvaderiana and B. peltatifolia), and whole genome shotgun data for an additional 74 Begonia representatives to investigate lineage evolution and shade adaptation of the genus. The four genome assemblies range in size from 331.75 Mb (B. peltatifolia) to 799.83 Mb (B. masoniana), and harbor 22 059-23 444 protein-coding genes. Synteny analysis revealed a lineage-specific whole-genome duplication (WGD) that occurred just before the diversification of Begonia. Functional enrichment of gene families retained after WGD highlights the significance of modified carbohydrate metabolism and photosynthesis possibly linked to shade adaptation in the genus, which is further supported by expansions of gene families involved in light perception and harvesting. Phylogenomic reconstructions and genomics studies indicate that genomic introgression has also played a role in the evolution of Begonia. Overall, this study provides valuable genomic resources for Begonia and suggests potential drivers underlying the diversity and adaptive evolution of this mega-diverse clade.

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.

Morphological and Genome-Wide Evidence of Homoploid Hybridisation in Urospermum (Asteraceae)

The genus Urospermum is distributed in the Mediterranean region and Macaronesia, and has been introduced to other extra-Mediterranean regions. Although the two species constituting the genus, U. dalechampii and U. picroides, are frequently found together, hybrids have so far only been reported once, from Morocco. However, we found certain individuals in Catalonia, whose intermediate morphology suggested a potential hybrid origin. In this study, we applied morphological and molecular methods to investigate the origin of those individuals. Intermediate features at phenotype, karyological, cytogenetic, and genomic levels were identified in morphologically intermediate individuals, supporting their homoploid hybrid origin. Chloroplast sequence data suggest that U. dalechampii is the maternal progenitor of the hybrid. Together with the intermediate traits displayed, the lack of fertile seeds suggests that hybrids are probably F1. Future monitoring studies will be, nonetheless, needed to evaluate the extent of hybridisation and its potential impact on the biology of the genus.

ModEst: Accurate estimation of genome size from next generation sequencing data

Accurate estimates of genome sizes are important parameters for both theoretical and practical biodiversity genomics. Here we present a fast, easy-to-implement and accurate method to estimate genome size from the number of bases sequenced and the mean sequencing depth. To estimate the latter, we take advantage of the fact that an accurate estimation of the Poisson distribution parameter lambda is possible from truncated data, restricted to the part of the sequencing depth distribution representing the true underlying distribution. With simulations we show that reasonable genome size estimates can be gained even from low-coverage (10×), highly discontinuous genome drafts. Comparison of estimates from a wide range of taxa and sequencing strategies with flow cytometry estimates of the same individuals showed a very good fit and suggested that both methods yield comparable, interchangeable results.

Genome Size Doubling Arises From the Differential Repetitive DNA Dynamics in the Genus Heloniopsis (Melanthiaceae)

Plant genomes are highly diverse in size and repetitive DNA composition. In the absence of polyploidy, the dynamics of repetitive elements, which make up the bulk of the genome in many species, are the main drivers underpinning changes in genome size and the overall evolution of the genomic landscape. The advent of high-throughput sequencing technologies has enabled investigation of genome evolutionary dynamics beyond model plants to provide exciting new insights in species across the biodiversity of life. Here we analyze the evolution of repetitive DNA in two closely related species of Heloniopsis (Melanthiaceae), which despite having the same chromosome number differ nearly twofold in genome size [i.e., H. umbellata (1C = 4,680 Mb), and H. koreana (1C = 2,480 Mb)]. Low-coverage genome skimming and the RepeatExplorer2 pipeline were used to identify the main repeat families responsible for the significant differences in genome sizes. Patterns of repeat evolution were found to correlate with genome size with the main classes of transposable elements identified being twice as abundant in the larger genome of H. umbellata compared with H. koreana. In addition, among the satellite DNA families recovered, a single shared satellite (HeloSAT) was shown to have contributed significantly to the genome expansion of H. umbellata. Evolutionary changes in repetitive DNA composition and genome size indicate that the differences in genome size between these species have been underpinned by the activity of several distinct repeat lineages.

Descending Dysploidy and Bidirectional Changes in Genome Size Accompanied Crepis (Asteraceae) Evolution

The evolution of the karyotype and genome size was examined in species of Crepis sensu lato. The phylogenetic relationships, inferred from the plastid and nrITS DNA sequences, were used as a framework to infer the patterns of karyotype evolution. Five different base chromosome numbers (x = 3, 4, 5, 6, and 11) were observed. A phylogenetic analysis of the evolution of the chromosome numbers allowed the inference of x = 6 as the ancestral state and the descending dysploidy as the major direction of the chromosome base number evolution. The derived base chromosome numbers (x = 5, 4, and 3) were found to have originated independently and recurrently in the different lineages of the genus. A few independent events of increases in karyotype asymmetry were inferred to have accompanied the karyotype evolution in Crepis. The genome sizes of 33 Crepis species differed seven-fold and the ancestral genome size was reconstructed to be 1C = 3.44 pg. Both decreases and increases in the genome size were inferred to have occurred within and between the lineages. The data suggest that, in addition to dysploidy, the amplification/elimination of various repetitive DNAs was likely involved in the genome and taxa differentiation in the genus.

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights

BACKGROUND AND AIMS

Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species.

METHODS

We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples.

KEY RESULTS

We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small.

CONCLUSIONS

We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae)

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Genome downsizing after polyploidy: mechanisms, rates and selection pressures

An analysis of over 10 000 plant genome sizes (GSs) indicates that most species have smaller genomes than expected given the incidence of polyploidy in their ancestries, suggesting selection for genome downsizing. However, comparing ancestral GS with the incidence of ancestral polyploidy suggests that the rate of DNA loss following polyploidy is likely to have been very low (4-70 Mb/million years, 4-482 bp/generation). This poses a problem. How might such small DNA losses be visible to selection, overcome the power of genetic drift and drive genome downsizing? Here we explore that problem, focussing on the role that double-strand break (DSB) repair pathways (non-homologous end joining and homologous recombination) may have played. We also explore two hypotheses that could explain how selection might favour genome downsizing following polyploidy: to reduce (i) nitrogen (N) and phosphate (P) costs associated with nucleic acid synthesis in the nucleus and the transcriptome and (ii) the impact of scaling effects of GS on cell size, which influences CO₂ uptake and water loss. We explore the hypothesis that losses of DNA must be fastest in early polyploid generations. Alternatively, if DNA loss is a more continuous process over evolutionary time, then we propose it is a byproduct of selection elsewhere, such as limiting the damaging activity of repetitive DNA. If so, then the impact of GS on photosynthesis, water use efficiency and/or nutrient costs at the nucleus level may be emergent properties, which have advantages, but not ones that could have been selected for over generational timescales.

Comparative phylogenetics of repetitive elements in a diverse order of flowering plants (Brassicales)

Genome sizes of plants have long piqued the interest of researchers due to the vast differences among organisms. However, the mechanisms that drive size differences have yet to be fully understood. Two important contributing factors to genome size are expansions of repetitive elements, such as transposable elements (TEs), and whole-genome duplications (WGD). Although studies have found correlations between genome size and both TE abundance and polyploidy, these studies typically test for these patterns within a genus or species. The plant order Brassicales provides an excellent system to further test if genome size evolution patterns are consistent across larger time scales, as there are numerous WGDs. This order is also home to one of the smallest plant genomes, Arabidopsis thaliana-chosen as the model plant system for this reason-as well as to species with very large genomes. With new methods that allow for TE characterization from low-coverage genome shotgun data and 71 taxa across the Brassicales, we confirm the correlation between genome size and TE content, however, we are unable to reconstruct phylogenetic relationships and do not detect any shift in TE abundance associated with WGD.

Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids

Although the evolutionary drivers of genome size change are known, the general patterns and mechanisms of plant genome size evolution are yet to be established. Here we aim to assess the relative importance of proliferation of repetitive DNA, chromosomal variation (including polyploidy), and the type of endoreplication for genome size evolution of the Pleurothallidinae, the most species-rich orchid lineage. Phylogenetic relationships between 341 Pleurothallidinae representatives were refined using a target enrichment hybrid capture combined with high-throughput sequencing approach. Genome size and the type of endoreplication were assessed using flow cytometry supplemented with karyological analysis and low-coverage Illumina sequencing for repeatome analysis on a subset of samples. Data were analyzed using phylogeny-based models. Genome size diversity (0.2-5.1 Gbp) was mostly independent of profound chromosome count variation (2n = 12-90) but tightly linked with the overall content of repetitive DNA elements. Species with partial endoreplication (PE) had significantly greater genome sizes, and genomic repeat content was tightly correlated with the size of the non-endoreplicated part of the genome. In PE species, repetitive DNA is preferentially accumulated in the non-endoreplicated parts of their genomes. Our results demonstrate that proliferation of repetitive DNA elements and PE together shape the patterns of genome size diversity in orchids.

Primeras medidas del tamaño del genoma en Carduncellusy los géneros afines Femeniasia y Phonus (Asteraceae, Cardueae), con datos para 21 táxones

El tamaño del genoma de 18 especies del género Carduncellus, dos especies de los géneros relacionados, Phonus y el género monotípico Femeniasia (F. balearica) ha sido medido por primera vez mediante citometría de flujo. Los niveles de ploidía se asignaron utilizando datos de tamaño del genoma junto con los recuentos de cromosomas previamente reportados. Se construyó un marco filogenético para visualizar la distribución de las características citogenéticas de los táxones. Los resultados confirmaron tres niveles de ploidía (2x, 4x y 6x), con un predominio de los táxones diploides. Los valores de 2C oscilaron entre 3,24 pg en Carduncellus calvus y 11,16 pg en C. eriocephalus, mientras que el tamaño del genoma monoploide (1Cx) osciló entre 1,29 pg en C. duvauxii (4x) y 2,30 pg en Phonus rhiphaeus (2x). La media de los valores 1Cx para los tetraploides fue menor que para los diploides. Los valores de tamaño del genoma de Carduncellus, Femeniasia y Phonus fueron más elevados que los de Carthamus dentro del mismo nivel de ploidía. Este resultado concuerda con una tendencia frecuentemente observada en plantas en la que los táxones con ciclos de vida largos presentan tamaños del genoma más elevados que los táxones relacionados que poseen ciclos de vida cortos.

Horizontally Acquired nrDNAs Persist in Low Amounts in Host Hordeum Genomes and Evolve Independently of Native nrDNA

Nuclear ribosomal DNA (nrDNA) has displayed extraordinary dynamics during the evolution of plant species. However, the patterns and evolutionary significance of nrDNA array expansion or contraction are still relatively unknown. Moreover, only little is known of the fate of minority nrDNA copies acquired between species via horizontal transfer. The barley genus Hordeum (Poaceae) represents a good model for such a study, as species of section Stenostachys acquired nrDNA via horizontal transfer from at least five different panicoid genera, causing long-term co-existence of native (Hordeum-like) and non-native (panicoid) nrDNAs. Using quantitative PCR, we investigated copy number variation (CNV) of nrDNA in the diploid representatives of the genus Hordeum. We estimated the copy number of the foreign, as well as of the native ITS types (ribotypes), and followed the pattern of their CNV in relation to the genus' phylogeny, species' genomes size and the number of nrDNA loci. For the native ribotype, we encountered an almost 19-fold variation in the mean copy number among the taxa analysed, ranging from 1689 copies (per 2C content) in H. patagonicum subsp. mustersii to 31342 copies in H. murinum subsp. glaucum. The copy numbers did not correlate with any of the genus' phylogeny, the species' genome size or the number of nrDNA loci. The CNV was high within the recognised groups (up to 13.2 × in the American I-genome species) as well as between accessions of the same species (up to 4×). Foreign ribotypes represent only a small fraction of the total number of nrDNA copies. Their copy numbers ranged from single units to tens and rarely hundreds of copies. They amounted, on average, to between 0.1% (Setaria ribotype) and 1.9% (Euclasta ribotype) of total nrDNA. None of the foreign ribotypes showed significant differences with respect to phylogenetic groups recognised within the sect. Stenostachys. Overall, no correlation was found between copy numbers of native and foreign nrDNAs suggesting the sequestration and independent evolution of native and non-native nrDNA arrays. Therefore, foreign nrDNA in Hordeum likely poses a dead-end by-product of horizontal gene transfer events.

A phylogenetic analysis of Momordica (Cucurbitaceae) in India based on karyo-morphology, nuclear DNA content and rDNA ITS1-5.8S-ITS2 sequences

The infrageneric delimitation of Momordica, a medicinally important genus of Cucurbitaceae, is ill-defined until date. Momordica chromosomes are extremely small and are difficult to stain and visualize because of the dense cytoplasmic background. We have conducted karyomorphometric analysis by EMA method in five Indian Momordica species, and the nuclear genome sizes were estimated by flow cytometry for the first time. The somatic chromosome numbers ranged from 2n = 18 to 56 in the species. We have resolved previously disputed chromosome numbers in M. cymbalaria and M. dioica as 2n = 18 (lowest) and 2n = 56, respectively. Chromosome counts in the other species were re-confirmed as 2n = 22 in M. charantia, 2n = 28 in M. cochinchinensis and 2n = 56 in M. subangulata. The largest genome size was recorded in M. cymbalaria (3.74 pg 2C^-1), while the smallest size (0.72 pg 2C^-1) was detected in M. charantia var. charantia. The nuclear genome sizes were analysed in comparison to chromosome numbers and total chromosome lengths of the species. Karyomorphometric indices showed comparable symmetric karyotypes in the species except in M. cymbalaria having tendency towards asymmetry. The UPGMA phenogram and principle component analysis based on nuclear DNA contents and karyomorphometric parameters demonstrated interspecies differences, intraspecific distinction within M. charantia varieties and highlighted distinction of M. cymbalaria. This study was further supported by the rDNA ITS sequence-based phylogenetic analysis which revealed the monophyletic origin of the Indian members of Momordica and clarified the intraspecies relationship among the studied members. As a whole, the study brought out new insights on species diversification within the genus Momordica in India and would benefit further studies on biosystematics and plant breeding programmes.

Hybridization and cryptic speciation in the Iberian endemic plant genus Phalacrocarpum (Asteraceae-Anthemideae)

Understanding the role and impact of reticulation in phylogenetic inquiry has improved with extended use of high throughput sequencing data. Yet, due to the dynamism of genomes over evolutionary time, disentangling old hybridization events remains a serious challenge. Phalacrocarpum (DC.) Willk. is one of the 27 Iberian endemic plant genera, currently considered monotypic but including three subspecies. Its uncertain phylogenetic relationships within tribe Anthemideae (Asteraceae) point to an Early Miocene divergence from its sister group, and its persistent taxonomic instability has been proposed to be due to hybridization. We aim at understanding the evolutionary history of this genus using SNPs called from a genotyping-by-sequencing (GBS) analysis, Sanger sequences-from three plastid DNA regions (psbJ-petA, petB-petD, trnH-psbA) and the nuclear ribosomal ITS regions (cloned)-as well as leaf morphometric multivariate analysis. SNP data and Sanger sequences strongly support the unforeseen existence of a cryptic species in the eastern populations of P. oppositifolium subsp. anomalum. Broad molecular and morphometric patterns of variation found in conflictive populations from the Sanabria Valley region convincingly identify a recent previously undocumented hybrid zone. By contrast, evidence is less conclusive on relationships between subspecies hoffmannseggii, oppositifolium and a second conflictive group distributed along the Galician-Portuguese border (Orense massifs). Although genetic clustering analysis of SNP data suggests that the former subspecies was the maternal progenitor in hybridization events that gave rise to the other two groups, we found considerable uniqueness of ITS ribotypes and plastid haplotypes in them. This result, in the context of Pleistocene climatically-driven range shifts in NW Iberian Peninsula, can be due to periods of isolation, genetic bottlenecks and drift superimposed on old hybridization events. Our study confirms the idea that unravelling old hybridization events may be compromised by the suite of evolutionary processes accumulated subsequently, particularly in areas with a history of climatic instability.

Characterization and Dynamics of Repeatomes in Closely Related Species of Hieracium (Asteraceae) and Their Synthetic and Apomictic Hybrids

The repetitive content of the plant genome (repeatome) often represents its largest fraction and is frequently correlated with its size. Transposable elements (TEs), the main component of the repeatome, are an important driver in the genome diversification due to their fast-evolving nature. Hybridization and polyploidization events are hypothesized to induce massive bursts of TEs resulting, among other effects, in an increase of copy number and genome size. Little is known about the repeatome dynamics following hybridization and polyploidization in plants that reproduce by apomixis (asexual reproduction via seeds). To address this, we analyzed the repeatomes of two diploid parental species, Hieracium intybaceum and H. prenanthoides (sexual), their diploid F1 synthetic and their natural triploid hybrids (H. pallidiflorum and H. picroides, apomictic). Using low-coverage next-generation sequencing (NGS) and a graph-based clustering approach, we detected high overall similarity across all major repeatome categories between the parental species, despite their large phylogenetic distance. Medium and highly abundant repetitive elements comprise ∼70% of Hieracium genomes; most prevalent were Ty3/Gypsy chromovirus Tekay and Ty1/Copia Maximus-SIRE elements. No TE bursts were detected, neither in synthetic nor in natural hybrids, as TE abundance generally followed theoretical expectations based on parental genome dosage. Slight over- and under-representation of TE cluster abundances reflected individual differences in genome size. However, in comparative analyses, apomicts displayed an overabundance of pararetrovirus clusters not observed in synthetic hybrids. Substantial deviations were detected in rDNAs and satellite repeats, but these patterns were sample specific. rDNA and satellite repeats (three of them were newly developed as cytogenetic markers) were localized on chromosomes by fluorescence in situ hybridization (FISH). In a few cases, low-abundant repeats (5S rDNA and certain satellites) showed some discrepancy between NGS data and FISH results, which is due partly to the bias of low-coverage sequencing and partly to low amounts of the satellite repeats or their sequence divergence. Overall, satellite DNA (including rDNA) was markedly affected by hybridization, but independent of the ploidy or reproductive mode of the progeny, whereas bursts of TEs did not play an important role in the evolutionary history of Hieracium.