BACKGROUND:Aegilops biuncialisis tetraploid grass species with UbUbMbMb genome constitution, distributed in Mediterranean and the Middle East. It carries many valuable traits such as disease resistance, drought tolerance, high micronutrient content, which are in demand of wheat breeding. Transfer of genetic material from Ae. biuncialis is not easy due to substantial modification of Ub and especially Mb genomes. Cytogenetic markers permitting easy and reliable chromosome identification will be helpful for successful manipulation with alien genetic material and introgression of useful traits into wheat. Development of chromosome nomenclature is complicated by the significant karyotype diversity of Ae. biuncialis. RESULTS:We used various combinations of eleven DNA probes for studying intraspecific karyotype divergence of Ae. biuncialis by FISH; among them pTa-566, pTa713 and pSc119.2 probes proved to be most informative for chromosome identification and analysis of karyotype evolution. FISH discriminated three chromosomal groups designated A, B and C and showed that the Ub genome of Ae. biuncialis is less modified relative to the parental compared to Mb genome. Based on the obtained results we suggested that Ae. biuncialis originated via multiple hybridization events and the Ub and Mb genomes of group A, B, and C were contributed by different forms of Ae. umbellulata and Ae. comosa. The Mb genome of groups A and C probably derived from Ae. comosa subsp. comosa, whereas in the B-group - from subsp. heldreichii. Divergence of chromosomal groups of Ae. biuncialis was also accompanied by structural chromosome rearrangements. Using multiple DNA probes, we showed that reciprocal translocation between chromosomes 1MbL and 7MbL followed by pericentric inversion of modified 1Mb occurred in group A. Intraspecific divergence of Ae. biuncialis was also associated with amplification/ elimination/ redistribution of repetitive DNA families. CONCLUSIONS:Our study revealed complex genome structure of wild tetraploid grass Ae. biuncialis which might have occurred via multiple hybridization events resulting in formation of three distinct chromosomal groups. Their divergence was accompanied by different chromosomal rearrangements which resulted in formation of highly distinct karyotypes. By using FISH markers, we evaluated the relations between chromosomes of the three groups and developed genetic nomenclature of Ae. biuncialis chromosomes.

The genus Salvia L. (Lamiaceae) is characterized by complex taxonomy and controversial phylogeny. This genus includes about a thousand species with worldwide distribution and high ecological, structural, functional and morphological diversity. Because of their high content of essential oils, various Salvia plants are widely used in medicine, as well as in the food, perfume, cosmetic, and paint industries; they also are valuable melliferous resources. The present study reviews the taxonomic history of the genus Salvia and the phylogenetic relationships between the taxa within the subgenera Salvia, Sclarea, and Glutinaria. Among the Salvia species, three basic chromosome numbers, x = 7, x = 8, and x = 11, were most common, although other basic chromosome numbers (x = 6–19) were determined, which was probably due to events of dysploidy, aneupoidy, and/or polyploidy occurring during speciation. Recent molecular cytogenetic studies based on Next Generation Sequencing technologies have clarified the chromosomal organization of several Salvia species. The patterns of chromosome distribution of 45S rDNA, 5S rDNA, and satellite DNAs made it possible to assess their intra- and interspecific chromosome diversity. However, further cytogenetic studies are needed to characterize the chromosomes in the genomes of other Salvia species and specify the genomic relationships among them.

Amaranthus L. includes valuable and promising crops of multi-purpose use, having high morphological diversity and complicated taxonomy. Their karyotypes and genomic relationships remain insufficiently studied. For the first time, a comparative repeatome analysis of Amaranthus tricolor L., Amaranthus cruentus L., and Amaranthus hypochondriacus L. was performed based on the high-throughput sequencing data obtained via bioinformatic analyses using the RepeatExplorer2/TAREAN/DANTE_LTR pipelines. Interspecific variations in the abundance of Ty1 Copia and Ty3 Gypsy retroelements, DNA transposons, and ribosomal and satellite DNA (satDNA) were detected. Based on fluorescence in situ hybridization (FISH), chromosome mapping of 45S rDNA, 5S rDNA, and satDNAs AmC9 and AmC70, and unique karyograms of A. tricolor, A. cruentus, Amaranthus paniculatus L., and A. hypochondriacus were constructed. The analysis of the interspecies genome diversity/similarity in DNA repeat contents, sequences of the identified satDNAs, and chromosome distribution patterns of the studied molecular markers indicated that these species might also share a common evolutionary ancestor. However, the genomes of A. cruentus, A. paniculatus, and A. hypochondriacus were more similar compared to A. tricolor, which aligns with the previous phylogenetic data. Our results demonstrate that cytogenomic studies might provide important data on Amaranthus species relationships elucidating taxonomy and evolution of these valuable crops.

Polemonium L. (Polemoniaceae) is a widespread genus native to subarctic and arctic regions of the Northern Hemisphere. The taxonomy and genome relationships within Polemonium are still unclear. We analyzed genomes of three species from each Polemonium caeruleum and Polemonium pulcherrimum complex using bioinformatic analysis by RepeatExplorer2/TAREAN pipelines of next-generation sequencing data. The repeatomes of all studied species were similar in type and number of repeats. Satellite DNAs (satDNAs) demonstrated high sequence identity within the studied species. FISH chromosome mapping of 45S rDNA, 5S rDNA, and two satDNAs Pol_C 33 and Pol_C 46 allowed us to construct the species karyograms and assess the genome diversity within the P. caeruleum complex and P. pulcherrimum complex, and also confirm the taxonomic status of P. kiushianum as an independent species. Our findings demonstrate a close genomic relationship among the species from P. caeruleum and P. pulcherrimum complexes, indicating the presence of a common ancestral genome; additionally, our results provide cytogenetic evidence for the monophyletic origin of these sections and also complex evolutionary history of the genus Polemonium. The developed approach may be a valuable framework for further investigation of the chromosomal organization of karyotypes in other species of the genus Polemonium.

Polemonium caeruleum L. (Polemoniaceae) is a perennial flowering plant native to Eurasia and North America, which is used as a fodder, medicinal, and ornamental plant. Many issues related to the taxonomy and origin of this valuable species still remain unclear. The intraspecific genetic variability of P. caeruleum and chromosomal organization of its genome are insufficiently studied. For the first time, we analyzed NGS genomic data of P. caeruleum using ReapeatExplorer2/TAREAN/DANTE Pipelines. In its repeatome, we identified 66.08% of Class I retrotransposons; 0.57% of Class II transposons; 0.42% of ribosomal DNA; and 0.87% of satellite DNA (six high-confident and three low-confident putative satellite DNAs). FISH chromosome mapping of seven tandem DNAs was carried out in two P. caeruleum varieties and two wild populations. Our results demonstrated the effectiveness of using satDNAs Pol_C 46 and Pol_C 33 in combination with 45S rDNA and 5S rDNA for precise chromosome identification. This approach allowed us to study intraspecific chromosomal variability and detect chromosomal rearrangements in the studied accessions of P. caeruleum, which could be related to the speciation process. These novel molecular markers are important for chromosome studies within Polemonium to clarify its taxonomy and phylogeny, and also, they expand the potential of different breeding programs.

The species of the genus Deschampsia are difficult for identification, and the genus is difficult for taxonomic treatment. The regions of 35S rRNA genes were studied for the species of the genus Deschampsia of different geographical origin with a method of sequencing by Sanger (ITS1-5.8S rRNA gene-ITS2, 14 species) and with a method of a locus-specific next-generation sequencing (NGS) on the Illumina platform (ITS1-5.8S rRNA, 7 species). All species of Deschampsia formed one clade; the species, referred by some authors on the basis of morphological characters to the species D. cespitosa s.l., entered one subclade. Subantarctic species formed a separate subclade and their ribotypes formed their own subnetwork. Avenella flexuosa, earlier referred to Deschampsia, entered the other clade, though this species contains some ribotypes common with some Deschampsia species. Deschampsia pamirica and related mountain species have their own specific ribotype groups. On the network of the ribotypes, one can see that D. cespitosa from Great Britain forms a network with some species, but D. cespitosa from the USA forms its own network. Ribotype analysis of each sample revealed traces of introgression with Deyeuxia/Calamagrostis in D. cespitosa and with A. flexuosa and probable introgression of Northern and subantarctic species.

Amaranthus cruentus L. and Amaranthus hypochondriacus L. are valuable and promising food crops for multi-purpose use that are distributed worldwide in temperate, subtropical, and tropical zones. However, their karyotypes and genomic relationships still remain insufficiently studied. For the first time, a comparative repeatome analysis of A. cruentus and A. hypochondriacus was performed based on the available NGS data; bioinformatic analyses using RepeatExplorer/TAREAN pipelines; and chromosome FISH mapping of 45S rDNA, 5S rDNA, and the most abundant satellite DNAs. In the repeatomes of these species, interspecific variations in the amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were detected. In the repeatomes of both species, shared satDNAs with high sequence similarity were identified. The chromosome distribution patterns of four effective molecular markers, 45S rDNA, 5S rDNA, AmC4, and AmC9, allowed us to identify all chromosome pairs in the species karyotypes, construct unique karyograms of A. cruentus and A. hypochondriacus, and confirm the close relationship between their genomes. These results are important for comparative karyotypic studies within the genus Amaranthus. Our findings demonstrated that cytogenomic analyses might provide important data on genomic relationships within Amaranthus and increase knowledge on genome organization in these valuable crops.

The genus Hedysarum L. (Fabaceae) includes about 200 species of annual and perennial herbs distributed in Asia, Europe, North Africa, and North America. Many species of this genus are valuable medicinal, melliferous, and forage resources. In this review, we consider the taxonomic history of the genus Hedysarum, the chromosomal organization of the species from the sections Hedysarum and Multicaulia, as well as phylogenetic relationships between these sections. According to morphological, genetic, and phylogenetic data, the genus Hedysarum is divided into three main sections: Hedysarum (= syn. Gamotion), Multicaulia, and Stracheya. In species of this genus, two basic chromosome numbers, x = 7 (section Hedysarum) and x = 8 (sections Multicaulia and Stracheya), were determined. The systematic positions of some species within the sections are still uncertain due to their morphological similarities. The patterns of distribution of molecular chromosomal markers (45S rDNA, 5S rDNA, and different satellite DNAs) in karyotypes of various Hedysarum species made it possible to determine their ploidy status and also specify genomic relationships within the sections Hedysarum and Multicaulia. Recent molecular phylogenetic studies clarified significantly the taxonomy and evolutionary development of the genus Hedysarum.

The cosmopolitan genus Hedysarum L. (Fabaceae) is divided into sections Hedysarum, Stracheya, and Multicaulia. This genus includes many valuable medicinal, melliferous, and forage species. The species taxonomy and genome relationships within the sections are still unclear. We examined intra- and interspecific diversity in the section (sect.) Hedysarum based on repeatome analyses using NGS data, bioinformatic technologies, and chromosome FISH mapping of 35S rDNA, 5S rDNA, and the identified satellite DNA families (satDNAs). A comparison of repeatomes of H. alpinum, H. theinum, and H. flavescens revealed differences in their composition. However, similarity in sequences of most satDNAs indicated a close relationship between genomes within sect. Hedysarum. New effective satDNA chromosomal markers were detected, which is important for karyotype analyses within Hedysarum. Intra- and interspecific variability in the chromosomal distribution patterns of the studied markers were revealed, and species karyograms were constructed. These results provided new insight into the karyotype structures and genomic diversity within sect. Hedysarum, clarified the systematic position of H. sachalinense and H. arcticum, and confirmed the distant genomic relationships between species from sections Hedysarum and Multicaulia. Our findings are important for further comparative genome studies within the genus Hedysarum.

Artificial hybrids between cultivated Avena species and wild Avena macrostachya that possess genes for resistance to biotic and abiotic stresses can be important for oat breeding. For the first time, a comprehensive study of genomes of artificial fertile hybrids Avena sativa × Avena macrostachya and their parental species was carried out based on the chromosome FISH mapping of satellite DNA sequences (satDNAs) and also analysis of intragenomic polymorphism in the 18S–ITS1–5.8S rDNA region, using NGS data. Chromosome distribution patterns of marker satDNAs allowed us to identify all chromosomes in the studied karyotypes, determine their subgenomic affiliation, and detect several chromosome rearrangements. Based on the obtained cytogenomic data, we revealed differences between two A. macrostachya subgenomes and demonstrated that only one of them was inherited in the studied octoploid hybrids. Ribotype analyses showed that the second major ribotype of A. macrostachya was species-specific and was not represented in rDNA pools of the octoploids, which could be related to the allopolyploid origin of this species. Our results indicate that the use of marker satDNAs in cytogenomic studies can provide important data on genomic relationships within Avena allopolyploid species and hybrids, and also expand the potential for interspecific crosses for breeding.

A comparative analysis of genomes of six Hedysarum species was carried out using the rapid genomic in situ hybridization (rapidGISH) technique, which makes it possible to identify common highly repetitive DNA sequences, and also shows the patterns of their distribution on chromosomes of the studied species. Dispersed localization of genomic DNA of H. flavescens (sect. Hedysarum) was found on chromosomes of H. alpinum and H. theinum (both belonged to sect. Hedysarum). On chromosomes of H. neglectum from the section Hedysarum, clustered hybridization signals of genomic DNA of H. flavescens were observed. The localization of H. alpinum genomic DNA was found along the all chromosomes of H. flavescens and H. theinum. On chromosomes of H. grandiflorum and H. dahuricum (both from sect. Multicaulia), small clustered and dispersed hybridization signals of genomic DNA of H. flavescens and H. alpinum were found. Thus, in the studied species, common highly repetitive DNA sequences of different types of their chromosome organization were revealed. Our findings indicate the presence of common origin of the species from the sections Hedysarum and Multicaulia with different degree of their relationship.

For the first time, a comparative molecular cytogenetic study of four species of Calendula L. (Asteraceae): C. officinalis, C. stellata, C. tripterocarpa, and C. arvensis was carried out. In each species, chromosome numbers were determined: C. stellata (2n = 2х = 14), C. officinalis (2n = 4x = 32), C. tripterocarpa (2n = 2x = 30), and C. arvensis (2n = 4x = 44), and also specific chromosome localization of 45S and 5S rDNA clusters was revealed with the use of fluorescence in situ hybridization. An additional polymorphic minor 45S rDNA hybridization signal was found in the C. officinalis karyotype, which was located median in the short arm of one chromosome. The patterns of chromosome distribution of the major sites of 45S and 5S rDNA in karyotypes of the studied species confirmed the hybrid origin of C. arvensis (2n = 4x = 44) which could be a result of the introgressive hybridization of two other species: C. stellata (2n = 2x = 14) and C. tripterocarpa (2n = 2x = 30) during speciation.

The taxonomically challenging genus Calendula L. (Asteraceae) includes lots of medicinal species characterized by their high morphological and karyological variability. For the first time, a repeatome analysis of a valuable medicinal plant Calendula officinalis L. was carried out using high-throughput genome DNA sequencing and RepeatExplorer/TAREAN pipelines. The FISH-based visualization of the 45S rDNA, 5S rDNA, and satellite DNAs of C. officinalis was performed on the chromosomes of C. officinalis, C. stellata Cav., C. tripterocarpa Rupr., and C. arvensis L. Three satellite DNAs were demonstrated to be new molecular chromosome markers to study the karyotype structure. Karyograms of the studied species were constructed, their ploidy status was specified, and their relationships were clarified. Our results showed that the C. officinalis karyotype differed from the karyotypes of the other three species, indicating its separate position in the Calendula phylogeny. However, the presence of common repeats revealed in the genomes of all the studied species could be related to their common origin. Our findings demonstrated that C. stellata contributed its genome to allotetraploid C. tripterocarpa, and C. arvensis is an allohexaploid hybrid between C. stellata and C. tripterocarpa. At the same time, further karyotype studies of various Calendula species are required to clarify the pathways of chromosomal reorganization that occurred during speciation.

Within the complicated and controversial taxonomy of cosmopolitan genus Salvia L. (Lamiaceae) are valuable species Salvia officinalis L. and Salvia sclarea L., which are important for the pharmaceutical, ornamental horticulture, food, and perfume industries. Genome organization and chromosome structure of these essential oil species remain insufficiently studied. For the first time, the comparative repeatome analysis of S. officinalis and S. sclarea was performed using the obtained NGS data, RepeatExplorer/TAREAN pipelines and FISH-based chromosome mapping of the revealed satellite DNA families (satDNAs). In repeatomes of these species, LTR retrotransposons made up the majority of their repetitive DNA. Interspecific variations in genome abundance of Class I and Class II transposable elements, ribosomal DNA, and satellite DNA were revealed. Four (S. sclarea) and twelve (S. officinalis) putative satDNAs were identified. Based on patterns of chromosomal distribution of 45S rDNA; 5S rDNA and the revealed satDNAs, karyograms of S. officinalis and S. sclarea were constructed. Promising satDNAs which can be further used as chromosome markers to assess inter- and intraspecific chromosome variability in Salvia karyotypes were determined. The specific localization of homologous satDNA and 45S rDNA on chromosomes of the studied Salvia species confirmed their common origin, which is consistent with previously reported molecular phylogenetic data.

Subpolar and polar ecotypes of Deschampsia sukatschewii (Popl.) Roshev, D. cespitosa (L.) P. Beauv, and D. antarctica E. Desv. are well adapted to stressful environmental conditions, which make them useful model plants for genetic research and breeding. For the first time, the comparative repeatome analyses of subpolar and polar D. sukatschewii, D. cespitosa, and D. antarctica was performed using RepeatExplorer/TAREAN pipelines and FISH-based chromosomal mapping of the identified satellite DNA families (satDNAs). In the studied species, mobile genetic elements of class 1 made up the majority of their repetitive DNA; interspecific variations in the total amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were revealed; 12–18 high confident and 7–9 low confident putative satDNAs were identified. According to BLAST, most D. sukatschewii satDNAs demonstrated sequence similarity with satDNAs of D. antarctica and D. cespitosa indicating their common origin. Chromosomal mapping of 45S rDNA, 5S rDNA, and satDNAs of D. sukatschewii allowed us to construct the species karyograms and detect new molecular chromosome markers important for Deschampsia species. Our findings confirmed that genomes of D. sukatschewii and D. cespitosa were more closely related compared to D. antarctica according to repeatome composition and patterns of satDNA chromosomal distribution.

The phylogeny of members of the family Linaceae DC. ex Perleb has not been adequately studied. In particular, data on the phylogenetic relationship between Linum stelleroides Planch. and other representatives of the blue-flowered flax are very controversial. In the present work, to clarify this issue, we obtained DNA sequences of three nuclear loci (IGS and ITS1 + 5.8S rDNA + ITS2 of the 35S rRNA gene and the 5S rRNA gene) and eight chloroplast loci (rbcL, the trnL–trnF intergenic spacer, matK, the 3′ trnK intron, ndhF, trnG, the psbA–trnH intergenic spacer, and rpl16) of 10 Linum L. species (L. stelleroides, L. hirsutum, L. perenne, L. leonii, L. lewisii, L. narbonense, L. decumbens, L. grandiflorum, L. bienne (syn. L. angustifolium), and L. usitatissimum) using high-throughput sequencing data. The phylogenetic analysis showed that L. stelleroides forms a basal branch in the blue-flowered flax clade. Previously found inconsistencies in the position of L. stelleroides and some other species in the Linaceae phylogenetic tree resulted from the erroneous species identification of some of the studied plant samples.

High-copy tandemly organized repeats (TRs), or satellite DNA, is an important but still enigmatic component of eukaryotic genomes. TRs comprise arrays of multi-copy and highly similar tandem repeats, which makes the elucidation of TRs a very challenging task. Oxford Nanopore sequencing data provide a valuable source of information on TR organization at the single molecule level. However, bioinformatics tools for de novo identification of TRs in raw Nanopore data have not been reported so far. We developed NanoTRF, a new python pipeline for TR repeat identification, characterization and consensus monomer sequence assembly. This new pipeline requires only a raw Nanopore read file from low-depth (<1×) genome sequencing. The program generates an informative html report and figures on TR genome abundance, monomer sequence and monomer length. In addition, NanoTRF performs annotation of transposable elements (TEs) sequences within or near satDNA arrays, and the information can be used to elucidate how TR–TE co-evolve in the genome. Moreover, we validated by FISH that the NanoTRF report is useful for the evaluation of TR chromosome organization—clustered or dispersed. Our findings showed that NanoTRF is a robust method for the de novo identification of satellite repeats in raw Nanopore data without prior read assembly. The obtained sequences can be used in many downstream analyses including genome assembly assistance and gap estimation, chromosome mapping and cytogenetic marker development.

The section Multicaulia is the largest clade in the genus Hedysarum L. (Fabaceae). Representatives of the sect. Multicaulia are valuable plants used for medicinal and fodder purposes. The taxonomy and phylogeny of the sect. Multicaulia are still ambiguous. To clarify the species relationships within sect. Multicaulia, we, for the first time, explored repeatomes of H. grandiflorum Pall., H. zundukii Peschkova, and H. dahuricum Turcz. using next-generation sequencing technologies and a subsequent bioinformatic analysis by RepeatExplorer/TAREAN pipelines. The comparative repeatome analysis showed that mobile elements made up 20-24% (Class I) and about 2-2.5% (Class II) of their repetitive DNAs. The amount of ribosomal DNA varied from 1 to 2.6%, and the content of satellite DNA ranged from 2.7 to 5.1%. For each species, five high confident putative tandem DNA repeats and 5-10 low confident putative DNA repeats were identified. According to BLAST, these repeats demonstrated high sequence similarity within the studied species. FISH-based mapping of 35S rDNA, 5S rDNA, and satDNAs made it possible to detect new effective molecular chromosome markers for Hedysarum species and construct the species karyograms. Comparison of the patterns of satDNA localization on chromosomes of the studied species allowed us to assess genome diversity within the sect. Multicaulia. In all studied species, we revealed intra- and interspecific variabilities in patterns of the chromosomal distribution of molecular chromosome markers. In H. gmelinii Ledeb. and H. setigerum Turcz. ex Fisch. et Meyer, similar subgenomes were detected, which confirmed the polyploid status of their genomes. Our findings demonstrated a close genomic relationship among six studied species indicating their common origin and confirmed the taxonomic status of H. setigerum as a subspecies of H. gmelinii as well as the validity of combining the sect. Multicaulia and Subacaulia into one sect. Multicaulia.

Polemonium caeruleum L. (Polemoniaceae) is a valuable medicinal herb with a wide spectrum of biological activities. Under natural conditions, the productivity of this species is rather low. In this study, colchicine-induced tetraploid plants (2n = 4x = 36) of P. caeruleum were obtained, and for the first time, their morphological and cytogenetic characterization was performed. In the tetraploid plants, raw material productivity and also the content of triterpene saponins were significantly higher than in the control diploids. The analysis of chromosome behavior at meiosis and FISH chromosome mapping of 45S and 5S rDNA generally demonstrated stability of both genomes in the tetraploid plants. Based on chromosome morphology and distribution patterns of the studied molecular cytogenetic markers, all chromosome pairs in karyotypes were identified, and chromosome karyograms and idiograms of P. caeruleum were constructed. The revealed specific microdiagnostic characteristics of P. caeruleum (strongly sinuous cells and anomocytic stomata of the leaf epidermis, and also glandular hairs along the veins) could be useful for raw material identification. In the obtained tetraploids, the predominance of large stomata on the lower leaf epidermis was determined. The studied tetraploids can be used in various breeding programs to obtain high-quality pharmaceutical raw materials of P. caeruleum.