Cytogenomics of Deschampsia P. Beauv. (Poaceae) Species Based on Sequence Analyses and FISH Mapping of CON/COM Satellite DNA Families

Integration of Repeatomic and Cytogenetic Data on Satellite DNA for the Genome Analysis in the Genus Salvia (Lamiaceae)

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.

A Pipeline NanoTRF as a New Tool for De Novo Satellite DNA Identification in the Raw Nanopore Sequencing Reads of Plant Genomes

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.

Integration of Genomic and Cytogenetic Data on Tandem DNAs for Analyzing the Genome Diversity Within the Genus Hedysarum L. (Fabaceae)

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.

Repeatome Analyses and Satellite DNA Chromosome Patterns in Deschampsia sukatschewii, D. cespitosa, and D. antarctica (Poaceae)

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.

Chromosomal Differentiation of Deschampsia (Poaceae) Based on Four Satellite DNA Families

Diverse families of satellite DNA (satDNA) were detected in heterochromatin regions of Deschampsia. This kind of repetitive DNA consists of tandem repeat sequences forming big arrays in genomes, and can contribute to lineages differentiation. The differentiation between types of satDNA is related to their sequence identity, the size and number of monomers forming the array, and their chromosomal location. In this work, four families of satDNA (D2, D3, D12, D13), previously isolated by genomic analysis, were studied on chromosomal preparations of 12 species of Deschampsia (D. airiformis, D. antarctica, D. cespitosa, D. cordillerarum, D. elongata, D. kingii, D. laxa, D. mendocina, D. parvula, D. patula, D. venustula, and Deschampsia sp) and one of Deyeuxia (D. eminens). Despite the number of satDNA loci showing interspecific variation, the general distribution pattern of each satDNA family is maintained. The four satDNA families are AT-rich and associated with DAPI + heterochromatin regions. D2, D3, and D12 have mainly subterminal distribution, while D13 is distributed in intercalary regions. Such conservation of satDNA patterns suggests a not random distribution in genomes, where the variation between species is mainly associated with the array size and the loci number. The presence of satDNA in all species studied suggests a low genetic differentiation of sequences. On the other hand, the variation of the distribution pattern of satDNA has no clear association with phylogeny. This may be related to high differential amplification and contraction of sequences between lineages, as explained by the library model.