Ancient Origin of Two 5S rDNA Families Dominating in the Genus Rosa and Their Behavior in the Canina-Type Meiosis

First insight into the genomes of the Pulmonaria officinalis group (Boraginaceae) provided by repeatome analysis and comparative karyotyping

BACKGROUND

The genus Pulmonaria (Boraginaceae) represents a taxonomically complex group of species in which morphological similarity contrasts with striking karyological variation. The presence of different numbers of chromosomes in the diploid state suggests multiple hybridization/polyploidization events followed by chromosome rearrangements (dysploidy). Unfortunately, the phylogenetic relationships and evolution of the genome, have not yet been elucidated. Our study focused on the P. officinalis group, the most widespread species complex, which includes two morphologically similar species that differ in chromosome number, i.e. P. obscura (2n = 14) and P. officinalis (2n = 16). Ornamental cultivars, morphologically similar to P. officinalis (garden escapes), whose origin is unclear, were also studied. Here, we present a pilot study on genome size and repeatome dynamics of these closely related species in order to gain new information on their genome and chromosome structure.

RESULTS

Flow cytometry confirmed a significant difference in genome size between P. obscura and P. officinalis, corresponding to the number of chromosomes. Genome-wide repeatome analysis performed on genome skimming data showed that retrotransposons were the most abundant repeat type, with a higher proportion of Ty3/Gypsy elements, mainly represented by the Tekay lineage. Comparative analysis revealed no species-specific retrotransposons or striking differences in their copy number between the species. A new set of chromosome-specific cytogenetic markers, represented by satellite DNAs, showed that the chromosome structure in P. officinalis was more variable compared to that of P. obscura. Comparative karyotyping supported the hybrid origin of putative hybrids with 2n = 15 collected from a mixed population of both species and outlined the origin of ornamental garden escapes, presumably derived from the P. officinalis complex.

CONCLUSIONS

Large-scale genome size analysis and repeatome characterization of the two morphologically similar species of the P. officinalis group improved our knowledge of the genome dynamics and differences in the karyotype structure. A new set of chromosome-specific cytogenetic landmarks was identified and used to reveal the origin of putative hybrids and ornamental cultivars morphologically similar to P. officinalis.

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.

Epigenetic histone H3 phosphorylation marks discriminate between univalent- and bivalent-forming chromosomes during canina asymmetrical meiosis

BACKGROUND AND AIMS

Dogroses (Rosa sect. Caninae) are mostly pentaploid, bearing 2n = 5x = 35 chromosomes in somatic cells. They evolved a unique form of asymmetrical meiosis characterized by two types of chromosomes: (1) chromosomes forming bivalents and distributed in the normal sexual way; and (2) chromosomes occurring as univalents and transferred by a female gamete only. In the mature pollen of pentaploid species, seven bivalent-derived chromosomes are transmitted to offspring, and 21 unpaired univalent chromosomes are eliminated during microsporogenesis. To discriminate between bivalent- and univalent-forming chromosomes, we studied histone H3 phosphorylation patterns regulating meiotic chromosome condensation and segregation.

METHODS

We analysed histone modification patterns during male canina meiosis in two representative dogrose species, 5x Rosa canina and 5x Rosa rubiginosa, by immunohistochemical and molecular cytogenetics approaches. Immunostaining of meiotic cells included α-tubulin, histone H3 phosphorylation (H3S10p, H3S28p and H3T3p) and methylation (H3K4me3 and H3K27me3) marks. In addition, fluorescent in situ hybridization was carried out with an 18S rDNA probe.

KEY RESULTS

In the first meiotic division, univalent chromosomes underwent equational division into chromatids, while homologues in bivalents were segregated as regular dyads. In diakinesis, bivalent chromosomes displayed strong H3 phosphorylation signals in proximal regions, spreading to the rest of the chromosome. In contrast, in univalents, the H3 phosphorylation signals were weaker, occurring mostly outside proximal regions largely overlapping with the H3K4me3 signals. Reduced phosphorylation was associated with relative under-condensation of the univalent chromosomes, particularly at early diakinesis.

CONCLUSIONS

We hypothesize that the absence of pairing and/or recombination in univalent chromosomes negatively affects the histone H3 phosphorylation of their chromatin and perhaps the loading of meiotic-specific cohesins. This apparently destabilizes cohesion of sister chromatids, leading to their premature split in the first meiotic division.

Untangling the hedge: Genetic diversity in clonally and sexually transmitted genomes of European wild roses, Rosa L

While European wild roses are abundant and widely distributed, their morphological taxonomy is complicated and ambiguous. In particular, the polyploid Rosa section Caninae (dogroses) is characterised by its unusual meiosis, causing simultaneous clonal and sexual transmission of sub-genomes. This hemisexual reproduction, which often co-occurs with vegetative reproduction, defies the standard definition of species boundaries. We analysed seven highly polymorphic microsatellite loci, scored for over 2 600 Rosa samples of differing ploidy, collected across Europe within three independent research projects. Based on their morphology, these samples had been identified as belonging to 21 dogrose and five other native rose species. We quantified the degree of clonality within species and at individual sampling sites. We then compared the genetic structure within our data to current rose morpho-systematics and searched for hemisexually co-inherited sets of alleles at individual loci. We found considerably fewer copies of identical multi-locus genotypes in dogroses than in roses with regular meiosis, with some variation recorded among species. While clonality showed no detectable geographic pattern, some genotypes appeared to be more widespread. Microsatellite data confirmed the current classification of subsections, but they did not support most of the generally accepted dogrose microspecies. Under canina meiosis, we found co-inherited sets of alleles as expected, but could not distinguish between sexually and clonally inherited sub-genomes, with only some of the detected allele combinations being lineage-specific.

Autotetraploid Origin of Chinese Cherry Revealed by Chromosomal Karyotype and In Situ Hybridization of Seedling Progenies

Polyploidy is considered a driving force in plant evolution and diversification. Chinese cherry [Cerasus pseudocerasus (Lindl.) G.Don], an economically important fruit crop native to China, has evolved at the tetraploid level, with a few pentaploid and hexaploid populations. However, its auto- or allo-polyploid origin remains unclear. To address this issue, we analyzed the ploidy levels and rDNA chromosomal distribution in self- and open-pollinated seedling progenies of tetraploid and hexaploid Chinese cherry. Genomic in situ hybridization (GISH) analysis was conducted to reveal the genomic relationships between Chinese cherry and diploid relatives from the genus Cerasus. Both self- and open-pollinated progenies of tetraploid Chinese cherry exhibited tetraploids, pentaploids, and hexaploids, with tetraploids being the most predominant. In the seedling progenies of hexaploid Chinese cherry, the majority of hexaploids and a few pentaploids were observed. A small number of aneuploids were also observed in the seedling progenies. Chromosome 1, characterized by distinct length characteristics, could be considered the representative chromosome of Chinese cherry. The basic Chinese cherry genome carried two 5S rDNA signals with similar intensity, and polyploids had the expected multiples of this copy number. The 5S rDNA sites were located at the per-centromeric regions of the short arm on chromosomes 4 and 5. Three 45S rDNA sites were detected on chr. 3, 4 and 7 in the haploid complement of Chinese cherry. Tetraploids exhibited 12 signals, while pentaploids and hexaploids showed fewer numbers than expected multiples. Based on the GISH signals, Chinese cherry demonstrated relatively close relationships with C. campanulata and C. conradinae, while being distantly related to another fruiting cherry, C. avium. In combination with the above results, our findings suggested that Chinese cherry likely originated from autotetraploidy.

Reproduction Modes and Conservation Implications in Three Polyploid Sorbus Stenoendemics in Eastern Slovakia (Central Europe)

The remarkable species diversity of the genus Sorbus is a result of polyploidization and frequent hybridization between interacting species of different cytotypes. Moreover, hybridization is possible between several parental taxa. Gametophytic apomixis, which is common among polyploid Sorbus taxa, indicates the role of clonal reproduction in the evolutionary stabilization of hybridogeneous genotypes. The precise determination of the origin of seeds and their quantitative evaluation may elucidate inter-cytotype interactions, the potential role of mixed-cytotype populations in evolutionary success, and the long-term survival of some hybrid species. We investigated the reproduction modes of selected species of Sorbus in mixed-cytotype populations in eastern Slovakia, Central Europe. We determined the pollen quality, seed production rate, and the ploidy level of mature trees, as well as the origin of the embryo and endosperm in seeds of the stenoendemics S. amici-petri, S. dolomiticola, and S. hornadensis. The tetraploids S. amici-petri and S. hornadensis are characterized by regular and highly stainable pollen grains and reproduce predominantly via pseudogamous apomixis. In contrast, triploid S. dolomiticola usually has oval, heterogenous, and weakly stainable pollen grains, suggesting male meiotic irregularities. Although seeds originate via pseudogamous apomixis in S. dolomiticola as well, the ploidy level of sperm cells participating in the fertilization of central cells is usually determined by co-occurring species of different cytotypes. This suggests that maintaining mating partners is necessary for the long-term survival of a triploid species. We documented rare B_III hybrids and the residual sexuality in tetraploids. The distribution of seeds of meiotic and apomeiotic origins in S. amici-petri shows bimodal characteristics; however, genotypes with predominantly sexual seed types are rare. Reproduction modes documented in polyploid stenoendemics of Sorbus and inferred microevolutionary intercytotype relationships highlight the mixed-cytotype populations as the source of biodiversity in apomictic plant complexes. We suggest that conservation efforts should focus on maintaining the species and cytotypic diversity of Sorbus populations, especially when it comes to the conservation of triploid species.

Analysis of 5S rDNA Genomic Organization Through the RepeatExplorer2 Pipeline: A Simplified Protocol

The ribosomal RNA genes (rDNA) are universal genome components with a housekeeping function, given the crucial role of ribosomal RNA in the synthesis of ribosomes and thus for life-on-Earth. Therefore, their genomic organization is of considerable interest for biologists, in general. Ribosomal RNA genes have also been largely used to establish phylogenetic relationships, and to identify allopolyploid or homoploid hybridization.Here, we demonstrate how high-throughput sequencing data, through graph clustering implemented in RepeatExplorer2 pipeline ( https://repeatexplorer-elixir.cerit-sc.cz/galaxy/ ), can be helpful to decipher the genomic organization of 5S rRNA genes. We show that the linear shapes of cluster graphs are reminiscent to the linked organization of 5S and 35S rDNA (L-type arrangement) while the circular graphs correspond to their separate arrangement (S-type). We further present a simplified protocol based on the paper by (Garcia et al., Front Plant Sci 11:41, 2020) about the use of graph clustering of 5S rDNA homoeologs (S-type) to identify hybridization events in the species history. We found that the graph complexity (i.e., graph circularity in this case) is related to ploidy and genome complexity, with diploids typically showing circular-shaped graphs while allopolyploids and other interspecific hybrids display more complex graphs, with usually two or more interconnected loops representing intergenic spacers. When a three-genomic comparative clustering analysis from a given hybrid (homoploid/allopolyploid) and its putative progenitor species (diploids) is performed, it is possible to identify the corresponding homoeologous 5S rRNA gene families, and to elucidate the contribution of each putative parental genome to the 5S rDNA pool of the hybrid. Thus, the analysis of 5S rDNA cluster graphs by RepeatExplorer, together with information coming from other sources (e.g., morphology, cytogenetics) is a complementary approach for the determination of allopolyploid or homoploid hybridization and even ancient introgression events.

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.

Uneven Levels of 5S and 45S rDNA Site Number and Loci Variations across Wild Chrysanthemum Accessions

Ribosomal DNA (rDNA) is an excellent cytogenetic marker owing to its tandem arrangement and high copy numbers. However, comparative studies have focused more on the number of rDNA site variations within the Chrysanthemum genus, and studies on the types of rDNA sites with the same experimental procedures at the species levels are lacking. To further explore the number and types of rDNA site variations, we combined related data to draw ideograms of the rDNA sites of Chrysanthemum accessions using oligonucleotide fluorescence in situ hybridization (Oligo-FISH). Latent variations (such as polymorphisms of 45S rDNA sites and co-localized 5S-45S rDNA) also occurred among the investigated accessions. Meanwhile, a significant correlation was observed between the number of 5S rDNA sites and chromosome number. Additionally, the clumped and concentrated geographical distribution of different ploidy Chrysanthemum accessions may significantly promote the karyotype evolution. Based on the results above, we identified the formation mechanism of rDNA variations. Furthermore, these findings may provide a reliable method to examine the sites and number of rDNA variations among Chrysanthemum and its related accessions and allow researchers to further understand the evolutionary and phylogenetic relationships of the Chrysanthemum genus.

Complex Scenarios of Reticulation, Polyploidization, and Species Diversity within Annual Pansies of Subsect. Bracteolatae (Viola Sect. Melanium, Violaceae) in Italy: Insights from 5S-IGS High-Throughput Sequencing and Plastid DNA Variation

Viola sect. Melanium, the so-called pansy, is an allopolyploid morphologically well-defined lineage of ca. 110 perennial and annual species in the northern hemisphere, characterized by markedly complex genomic configurations. Five annual pansies occur in Italy, four of which are morphologically very similar and belong to the informal 'V. tricolor species complex': V. arvensis (2n = 34), V. hymettia (2n = 16), V. kitaibeliana (2n = 16), and V. tricolor (2n = 26). Their field recognition is difficult and reflects a long-debated taxonomy often resulting in doubtful records in field inventories and across European herbaria. The current lack of comprehensive intra- and interspecific comparative studies and a relative scarcity of appropriate genetic markers coupled with unambiguous cytological descriptions are hindering clear taxa circumscription and phylogenetic inferences within this group. In this work, we tested DNA sequence variation of three highly variable plastid markers and High-Throughput Sequencing (HTS) of the nuclear ribosomal 5S-IGS region in an attempt to decipher species identity within the V. tricolor species complex and to obtain an insight on their genome organization and evolution. Our results document the close relationships within this species group, a reliable molecular resolution for V. tricolor, and the common ancestry of V. arvensis and the poorly differentiated V. kitaibeliana and V. hymettia. Evidence of an important inter-population geographical divergence was recorded in V. tricolor and V. arvensis, pointing at the existence of different eco-cytotypes within these entities. Overall diversity patterns and the occurrence of two to four differently diverging 5S-IGS lineages are discussed in the light of the acknowledged taxonomy and genomic evolutive trajectories of sect. Melanium.

5S Ribosomal DNA of Genus Solanum: Molecular Organization, Evolution, and Taxonomy

The Solanum genus, being one of the largest among high plants, is distributed worldwide and comprises about 1,200 species. The genus includes numerous agronomically important species such as Solanum tuberosum (potato), Solanum lycopersicum (tomato), and Solanum melongena (eggplant) as well as medical and ornamental plants. The huge Solanum genus is a convenient model for research in the field of molecular evolution and structural and functional genomics. Clear knowledge of evolutionary relationships in the Solanum genus is required to increase the effectiveness of breeding programs, but the phylogeny of the genus is still not fully understood. The rapidly evolving intergenic spacer region (IGS) of 5S rDNA has been successfully used for inferring interspecific relationships in several groups of angiosperms. Here, combining cloning and sequencing with bioinformatic analysis of genomic data available in the SRA database, we evaluate the molecular organization and diversity of IGS for 184 accessions, representing 137 species of the Solanum genus. It was found that the main mechanisms of IGS molecular evolution was step-wise accumulation of single base substitution or short indels, and that long indels and multiple base substitutions, which arose repeatedly during evolution, were mostly not conserved and eliminated. The reason for this negative selection seems to be association between indels/multiple base substitutions and pseudogenization of 5S rDNA. Comparison of IGS sequences allowed us to reconstruct the phylogeny of the Solanum genus. The obtained dendrograms are mainly congruent with published data: same major and minor clades were found. However, relationships between these clades and position of some species (S. cochoae, S. clivorum, S. macrocarpon, and S. spirale) were different from those of previous results and require further clarification. Our results show that 5S IGS represents a convenient molecular marker for phylogenetic studies on the Solanum genus. In particular, the simultaneous presence of several structural variants of rDNA in the genome enables the detection of reticular evolution, especially in the largest and economically most important sect. Petota. The origin of several polyploid species should be reconsidered.

5S-IGS rDNA in wind-pollinated trees (Fagus L.) encapsulates 55 million years of reticulate evolution and hybrid origins of modern species

Standard models of plant speciation assume strictly dichotomous genealogies in which a species, the ancestor, is replaced by two offspring species. The reality in wind-pollinated trees with long evolutionary histories is more complex: species evolve from other species through isolation when genetic drift exceeds gene flow; lineage mixing can give rise to new species (hybrid taxa such as nothospecies and allopolyploids). The multi-copy, potentially multi-locus 5S rDNA is one of few gene regions conserving signal from dichotomous and reticulate evolutionary processes down to the level of intra-genomic recombination. Therefore, it can provide unique insights into the dynamic speciation processes of lineages that diversified tens of millions of years ago. Here, we provide the first high-throughput sequencing (HTS) of the 5S intergenic spacers (5S-IGS) for a lineage of wind-pollinated subtropical to temperate trees, the Fagus crenata - F. sylvatica s.l. lineage, and its distant relative F. japonica. The observed 4963 unique 5S-IGS variants reflect a complex history of hybrid origins, lineage sorting, mixing via secondary gene flow, and intra-genomic competition between two or more paralogous-homoeologous 5S rDNA lineages. We show that modern species are genetic mosaics and represent a striking case of ongoing reticulate evolution during the past 55 million years.

Power and Weakness of Repetition - Evaluating the Phylogenetic Signal From Repeatomes in the Family Rosaceae With Two Case Studies From Genera Prone to Polyploidy and Hybridization (Rosa and Fragaria)

Plant genomes consist, to a considerable extent, of non-coding repetitive DNA. Several studies showed that phylogenetic signals can be extracted from such repeatome data by using among-species dissimilarities from the RepeatExplorer2 pipeline as distance measures. Here, we advanced this approach by adjusting the read input for comparative clustering indirectly proportional to genome size and by summarizing all clusters into a main distance matrix subjected to Neighbor Joining algorithms and Principal Coordinate Analyses. Thus, our multivariate statistical method works as a "repeatomic fingerprint," and we proved its power and limitations by exemplarily applying it to the family Rosaceae at intrafamilial and, in the genera Fragaria and Rosa, at the intrageneric level. Since both taxa are prone to hybridization events, we wanted to show whether repeatome data are suitable to unravel the origin of natural and synthetic hybrids. In addition, we compared the results based on complete repeatomes with those from ribosomal DNA clusters only, because they represent one of the most widely used barcoding markers. Our results demonstrated that repeatome data contained a clear phylogenetic signal supporting the current subfamilial classification within Rosaceae. Accordingly, the well-accepted major evolutionary lineages within Fragaria were distinguished, and hybrids showed intermediate positions between parental species in data sets retrieved from both complete repeatomes and rDNA clusters. Within the taxonomically more complicated and particularly frequently hybridizing genus Rosa, we detected rather weak phylogenetic signals but surprisingly found a geographic pattern at a population scale. In sum, our method revealed promising results at larger taxonomic scales as well as within taxa with manageable levels of reticulation, but success remained rather taxon specific. Since repeatomes can be technically easy and comparably inexpensively retrieved even from samples of rather poor DNA quality, our phylogenomic method serves as a valuable alternative when high-quality genomes are unavailable, for example, in the case of old museum specimens.

Variation in Ribosomal DNA in the Genus Trifolium (Fabaceae)

The genus Trifolium L. is characterized by basic chromosome numbers 8, 7, 6, and 5. We conducted a genus-wide study of ribosomal DNA (rDNA) structure variability in diploids and polyploids to gain insight into evolutionary history. We used fluorescent in situ hybridization to newly investigate rDNA variation by number and position in 30 Trifolium species. Evolutionary history among species was examined using 85 available sequences of internal transcribed spacer 1 (ITS1) of 35S rDNA. In diploid species with ancestral basic chromosome number (x = 8), one pair of 5S and 26S rDNA in separate or adjacent positions on a pair of chromosomes was prevalent. Genomes of species with reduced basic chromosome numbers were characterized by increased number of signals determined on one pair of chromosomes or all chromosomes. Increased number of signals was observed also in diploids Trifolium alpestre and Trifolium microcephalum and in polyploids. Sequence alignment revealed ITS1 sequences with mostly single nucleotide polymorphisms, and ITS1 diversity was greater in diploids with reduced basic chromosome numbers compared to diploids with ancestral basic chromosome number (x = 8) and polyploids. Our results suggest the presence of one 5S rDNA site and one 26S rDNA site as an ancestral state.