Sexual reproduction relies on meiotic chromosome pairing to form bivalents, a process that is complicated in polyploids owing to the presence of multiple subgenomes1. Uneven ploidy mostly results in sterility due to unbalanced chromosome pairing and segregation during meiosis. However, pentaploid dogroses (Rosa sect. Caninae; 2n = 5x = 35) achieve stable sexual reproduction through a unique mechanism: 14 chromosomes form bivalents and are transmitted biparentally, while the remaining 21 chromosomes are maternally inherited as univalents2,3. Despite being studied for over a century, the role of centromeres in this process has remained unclear. Here we analyse haplotype-resolved chromosome-level genome assemblies for three pentaploid dogroses. Subgenome phasing revealed a bivalent-forming subgenome with two highly homozygous chromosome sets and three divergent subgenomes lacking homologous partners, therefore explaining their meiotic behaviour. Comparative analyses of chromosome synteny, phylogenetic relationships and centromere composition indicate that the subgenomes originated from two divergent clades of the genus Rosa. Pollen genome analysis shows that subgenomes from different evolutionary origins form bivalents, supporting multiple origins of dogroses and highlighting variation in subgenome contributions. We reveal that bivalent-forming centromeres are enriched with ATHILA retrotransposons, contrasting with larger tandem-repeat-based centromeres mainly found in univalents. This centromere structural bimodality possibly contributes to univalent drive during female meiosis. Our findings provide insights into the unique reproductive strategies of dogroses, advancing our understanding of genome evolution, centromere diversity and meiotic mechanisms in organisms with asymmetrical inheritance systems.

• MeSH
• Centromere * genetics   metabolism   MeSH
• Chromosomes, Plant genetics   MeSH
• Phylogeny MeSH
• Genome, Plant genetics   MeSH
• Haplotypes genetics   MeSH
• Meiosis * genetics   MeSH
• Polyploidy * MeSH
• Pollen genetics   cytology   MeSH
• Retroelements genetics   MeSH
• Synteny genetics   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Retroelements MeSH

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.

• Keywords
• Pulmonaria, Comparative karyotyping, Genome size, Repeatome, Satellite DNA,
• MeSH
• Chromosomes, Plant * genetics   MeSH
• Genome Size MeSH
• Phylogeny MeSH
• Genome, Plant * MeSH
• Karyotype MeSH
• Karyotyping * MeSH
• Pulmonaria genetics   MeSH
• Publication type
• Journal Article MeSH
• Comparative Study MeSH

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.

• Keywords
• Rosa genus, 18S ribosomal DNA, canina meiosis, dogroses, euchromatin and heterochromatin, fluorescence in situ hybridization, histone H3 phosphorylation, immunostaining, non-disjunction,
• MeSH
• Chromosomes MeSH
• Epigenesis, Genetic MeSH
• Phosphorylation MeSH
• Histones * genetics   MeSH
• In Situ Hybridization, Fluorescence MeSH
• Meiosis * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Histones * MeSH

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.

• Keywords
• 5S ribosomal RNA genes, Allopolyploid hybridization, Genomic analysis, Graph clustering, Homoploid hybridization, Introgression, Repeatome,
• MeSH
• Phylogeny MeSH
• Genomics * MeSH
• Genes, rRNA MeSH
• DNA, Ribosomal genetics   MeSH
• RNA, Ribosomal, 5S * genetics   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA, Ribosomal MeSH
• RNA, Ribosomal, 5S * MeSH

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.

• Keywords
• Caninae, Fragaria, Rosaceae, graph-based clustering, high-throughput sequencing, phylogenetics, repeatome, repetitive DNA,
• Publication type
• Journal Article MeSH

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.

• Keywords
• 26S rDNA, 5S rDNA, clover, fluorescent in situ hybridization, genome structure, nucleotide polymorphism,
• Publication type
• Journal Article MeSH