CRISPR/Cas9-mediated мultiplexed multi-allelic mutagenesis of genes located on A, B and R subgenomes of hexaploid triticale.

KEY MESSAGE: The first-time generation of hexaploid triticale plants harbouring variable panels of novel mutations in gene families involved in starch biosynthesis has been achieved by the subgenome-independent multiplexed CRISPR/Cas9-mediated editing.

Comparative Characterization of Pseudoroegneria libanotica and Pseudoroegneria tauri Based on Their Repeatome Peculiarities.

Pseudoroegneria species play an important role among Triticeae grasses, as they are the putative donors of the St genome in many polyploid species. Satellite repeats are widely used as a reliable tool for tracking evolutionary changes because they are distributed throughout the genomes of plants. The aim of our work is to perform a comparative characterization of the repeatomes of the closely related species Ps. libanotica and Ps. tauri, and Ps. spicata was also included in the analysis. The overall repeatome structures of Ps. libanotica, Ps. tauri, and Ps. spicata were similar, with some individual peculiarities observed in the abundance of the SIRE (Ty1/Copia) retrotransposons, Mutator and Harbinger transposons, and satellites. Nine new satellite repeats that have been identified from the whole-genome sequences of Ps. spicata and Ps. tauri, as well as the CL244 repeat that was previously found in Aegilops crassa, were localized to the chromosomes of Ps. libanotica and Ps. tauri. Four satellite repeats (CL69, CL101, CL119, CL244) demonstrated terminal and/or distal localization, while six repeats (CL82, CL89, CL168, CL185, CL192, CL207) were pericentromeric. Based on the obtained results, it can be assumed that Ps. libanotica and Ps. tauri are closely related species, although they have individual peculiarities in their repeatome structures and patterns of satellite repeat localization on chromosomes. The evolutionary fate of the identified satellite repeats and their related sequences, as well as their distribution on the chromosomes of Triticeae species, are discussed. The newly developed St genome chromosome markers developed in the present research can be useful in population studies of Ps. libanotica and Ps. tauri; auto- and allopolyploids that contain the St genome, such as Thinopyrum, Elymus, Kengyilia, and Roegneria; and wide hybrids between wheat and related wild species.

Evaluation of phosphate rock as the only source of phosphorus for the growth of tall and semi-dwarf durum wheat and rye plants using digital phenotyping.

Background: Phosphorus nutrition is important for obtaining high yields of crop plants. However, wheat plants are known to be almost incapable of taking up phosphorus from insoluble phosphate sources, and reduced height genes are supposed to decrease this ability further. Methods: We performed a pot experiment using Triticum durum Desf. tall spring variety LD222, its near-isogenic semidwarf line carrying Rht17 (Reduced height 17) gene, and winter rye (Secale cereale L.) variety Chulpan. The individual plants were grown in quartz sand. The phosphorus was provided either as phosphate rock powder mixed with sand, or as monopotassium phosphate solution (normal nutrition control) or was not supplemented at all (no-phosphorus control). Other nutrients were provided in soluble form. During experiment the plants were assessed using the TraitFinder (Phenospex Ltd., Heerlen, Netherlands) digital phenotyping system for a standard set of parameters. Double scan with 90 degrees turns of pots around vertical axis vs. single scan were compared for accuracy of phenotyping. Results: The phenotyping showed that at least 20 days of growth after seedling emergence were necessary to get stable differences between genotypes. After this initial period, phenotyping confirmed poor ability of wheat to grow on substrate with phosphate rock as the only source of phosphorus compared to rye; however, Rht17 did not cause an additional reduction in growth parameters other than plant height under this variant of substrate. The agreement between digital phenotyping and conventionally measured traits was at previously reported level for grasses (R2 = 0.85 and 0.88 for digital biomass and 3D leaf area vs. conventionally measured biomass and leaf area, single scan). Among vegetation indices, only the normalized differential vegetation index (NDVI) and the green leaf index (GLI) showed significant correlations with manually measured traits, including the percentage of dead leaves area. The double scan improved phenotyping accuracy, but not substantially.

Cas9-targeted Nanopore sequencing rapidly elucidates the transposition preferences and DNA methylation profiles of mobile elements in plants.

Transposable element insertions (TEIs) are an important source of genomic innovation by contributing to plant adaptation, speciation, and the production of new varieties. The often large, complex plant genomes make identifying TEIs from short reads difficult and expensive. Moreover, rare somatic insertions that reflect mobilome dynamics are difficult to track using short reads. To address these challenges, we combined Cas9-targeted Nanopore sequencing (CANS) with the novel pipeline NanoCasTE to trace both genetically inherited and somatic TEIs in plants. We performed CANS of the EVADÉ (EVD) retrotransposon in wild-type Arabidopsis thaliana and rapidly obtained up to 40× sequence coverage. Analysis of hemizygous T-DNA insertion sites and genetically inherited insertions of the EVD transposon in the ddm1 (decrease in DNA methylation 1) genome uncovered the crucial role of DNA methylation in shaping EVD insertion preference. We also investigated somatic transposition events of the ONSEN transposon family, finding that genes that are downregulated during heat stress are preferentially targeted by ONSENs. Finally, we detected hypomethylation of novel somatic insertions for two ONSENs. CANS and NanoCasTE are effective tools for detecting TEIs and exploring mobilome organization in plants in response to stress and in different genetic backgrounds, as well as screening T-DNA insertion mutants and transgenic plants.

Erratum: Aegilops crassa Boiss. repeatome characterized using low-coverage NGS as a source of new FISH markers: application in phylogenetic studies of the Triticeae.

[This corrects the article DOI: 10.3389/fpls.2022.980764.].

Allelic Variation of Glu-A1 and Glu-B1 Genes in Winter Durum Wheat and Its Effect on Quality Parameters.

Winter durum wheat is a relatively young crop that is highly adaptable due to its winter type of growth habit. The priority of breeding and genetic improvement of winter durum wheat is to improve grain quality and pasta quality, largely determined by the glutenin storage proteins. In the present study, a collection of 76 accessions of winter durum wheat from P.P. Lukyanenko National Grain Centre was studied. The allelic state of high-molecular-weight glutenin genes, Glu-A1 and Glu-B1, using PCR markers and SDS-PAGE was identified and grain and pasta quality traits were assessed in a two-year field experiment. The positive effect of the Glu-A1a allele and a negative effect of Glu-A1c on the gluten index were shown. It was found that Glu-B1al and Glu-B1f have a positive effect on the quality and quantity of protein and gluten, while the Glu-A1c + Glu-B1al genotypes were closest to the high-quality category in protein-associated quality traits.

Permanent Spreading of 1RS.1AL and 1RS.1BL Translocations in Modern Wheat Breeding.

Wheat-rye translocations 1RS.1BL and 1RS.1AL are used in bread wheat breeding worldwide because a short arm of rye chromosome 1 (1RS) when introgressed into the wheat genome confers resistance to diseases, pests and better performance under drought-stress conditions. However, in durum wheat genotypes, these translocations occur only in experimental lines, although their advantages could enhance the potential of this crop. P.P. Lukyanenko National Grain Centre (NGC) has successfully developed commercially competitive cultivars of bread and durum wheat demanded by many agricultural producers in the South of Russia for decades. Here, 94 accessions of bread and 343 accessions of durum wheat, representing lines and cultivars from collection, competitive variety trials and breeding nursery developed at NGC were screened for 1RS using PCR markers and genomic in situ hybridization. The 1RS.1BL and 1RS.1AL translocations were detected in 38 and 6 bread wheat accessions, respectively. None of the durum wheat accessions showed translocation, despite the fact that some of them had 1RS.1BL donors in their pedigree. The absence of translocations in the studied durum wheat germplasm can be caused by the negative selection of 1RS carriers at different stages of the breeding process due to low quality and difficulties in transferring rye chromatin through wheat gametes.

The launch of satellite: DNA repeats as a cytogenetic tool in discovering the chromosomal universe of wild Triticeae.

Fluorescence in situ hybridization is a powerful tool that enables plant researchers to perform systematic, evolutionary, and population studies of wheat wild relatives as well as to characterize alien introgression into the wheat genome. This retrospective review reflects on progress made in the development of methods for creating new chromosomal markers since the launch of this cytogenetic satellite instrument to the present day. DNA probes based on satellite repeats have been widely used for chromosome analysis, especially for "classical" wheat probes (pSc119.2 and Afa family) and "universal" repeats (45S rDNA, 5S rDNA, and microsatellites). The rapid development of new-generation sequencing and bioinformatical tools, and the application of oligo- and multioligonucleotides has resulted in an explosion in the discovery of new genome- and chromosome-specific chromosome markers. Owing to modern technologies, new chromosomal markers are appearing at an unprecedented velocity. The present review describes the specifics of localization when employing commonly used vs. newly developed probes for chromosomes in J, E, V, St, Y, and P genomes and their diploid and polyploid carriers Agropyron, Dasypyrum, Thinopyrum, Pseudoroegneria, Elymus, Roegneria, and Kengyilia. Particular attention is paid to the specificity of probes, which determines their applicability for the detection of alien introgression to enhance the genetic diversity of wheat through wide hybridization. The information from the reviewed articles is summarized into the TRepeT database, which may be useful for studying the cytogenetics of Triticeae. The review describes the trends in the development of technology used in establishing chromosomal markers that can be used for prediction and foresight in the field of molecular biology and in methods of cytogenetic analysis.

GISH painting of the Y chromosomes suggests advanced phases of sex chromosome evolution in three dioecious Cannabaceae species (Humulus lupulus, H. japonicus, and Cannabis sativa).

In plants, dioecy is relatively rare, and it involves sex chromosome systems that often developed independently over time. These characteristics make dioecious plants an attractive model to study sex chromosome evolution. To clarify the patterns of plant sex chromosome evolution, studies should be performed on a wide range of dioecious species. It is interesting to study the sex chromosomes in related species that evolved during a long period of independent sex chromosome evolution. The Cannabaceae family includes three dioecious species with heteromorphic sex chromosomes. Cannabis sativa and Humulus lupulus use the XX/XY chromosome system, whereas Humulus japonicus contains multiple sex chromosomes (XX/XY1Y2). To better understand sex chromosome evolution and the level of genomic divergence of these three related species, we undertook self-GISH and comparative GISH analyses. The self-GISH allowed visualization of the Y chromosomes of C. sativa, H. lupulus, and H. japonicus. The self-GISH signal was distributed along the entire Y chromosome, excluding the pseudo-autosomal region (PAR). Our results indicate that the male-specific region of the Y chromosome (MSY) spans the overwhelming majority of the Y chromosomes of all three species studied. The self-GISH results reveal the accumulation of repetitive DNA sequences in the Y chromosomes of all three species studied. This sequences presented in autosomes and/or chromosome X at a lower copy number than in Y. In comparative GISH experiments where the probe DNA of one species was applied to another species, a weak signal was exclusively detected on 45S rDNA sites, indicating a high level of genomic differentiation of the species used in this study. We demonstrate small PAR size and opposing large MSY and its positions on Y chromosomes. We also found that these genomes are highly differentiated. Furthermore, the data obtained in this study indicate a long period of independent and advanced sex chromosome evolution. Our study provides a valuable basis for future genomic studies of sex and suggests that the Cannabaceae family offers a promising model to study sex chromosome evolution.

CRISPR/Cas9-induced modification of the conservative promoter region of VRN-A1 alters the heading time of hexaploid bread wheat.

In cereals, the vernalization-related gene network plays an important role in regulating the transition from the vegetative to the reproductive phase to ensure optimal reproduction in a temperate climate. In hexaploid bread wheat (Triticum aestivum L.), the spring growth habit is associated with the presence of at least one dominant locus of VERNALIZATION 1 gene (VRN-1), which usually differs from recessive alleles due to mutations in the regulatory sequences of the promoter or/and the first intron. VRN-1 gene is a key regulator of floral initiation; various combinations of dominant and recessive alleles, especially VRN-A1 homeologs, determine the differences in the timing of wheat heading/flowering. In the present study, we attempt to expand the types of VRN-A1 alleles using CRISPR/Cas9 targeted modification of the promoter sequence. Several mono- and biallelic changes were achieved within the 125-117 bp upstream sequence of the start codon of the recessive vrn-A1 gene in plants of semi-winter cv. 'Chinese Spring'. New mutations stably inherited in subsequent progenies and transgene-free homozygous plants carrying novel VRN-A1 variants were generated. Minor changes in the promoter sequence, such as 1-4 nucleotide insertions/deletions, had no effect on the heading time of plants, whereas the CRISPR/Cas9-mediated 8 bp deletion between -125 and -117 bp of the vrn-A1 promoter shortened the time of head emergence by up to 2-3 days. Such a growth habit was consistently observed in homozygous mutant plants under nonvernalized cultivation using different long day regimes (16, 18, or 22 h), whereas the cold treatment (from two weeks and more) completely leveled the effect of the 8 bp deletion. Importantly, comparison with wild-type plants showed that the implemented alteration has no negative effects on main yield characteristics. Our results demonstrate the potential to manipulate the heading time of wheat through targeted editing of the VRN-A1 gene promoter sequence on an otherwise unchanged genetic background.

Aegilops crassa Boiss. repeatome characterized using low-coverage NGS as a source of new FISH markers: Application in phylogenetic studies of the Triticeae.

Aegilops crassa Boiss. is polyploid grass species that grows in the eastern part of the Fertile Crescent, Afghanistan, and Middle Asia. It consists of tetraploid (4x) and hexaploid (6x) cytotypes (2n = 4x = 28, D1D (Abdolmalaki et al., 2019) XcrXcr and 2n = 6x = 42, D1D (Abdolmalaki et al., 2019) XcrXcrD2D (Adams and Wendel, 2005), respectively) that are similar morphologically. Although many Aegilops species were used in wheat breeding, the genetic potential of Ae. crassa has not yet been exploited due to its uncertain origin and significant genome modifications. Tetraploid Ae. crassa is thought to be the oldest polyploid Aegilops species, the subgenomes of which still retain some features of its ancient diploid progenitors. The D1 and D2 subgenomes of Ae. crassa were contributed by Aegilops tauschii (2n = 2x = 14, DD), while the Xcr subgenome donor is still unknown. Owing to its ancient origin, Ae. crassa can serve as model for studying genome evolution. Despite this, Ae. crassa is poorly studied genetically and no genome sequences were available for this species. We performed low-coverage genome sequencing of 4x and 6x cytotypes of Ae. crassa, and four Ae. tauschii accessions belonging to different subspecies; diploid wheatgrass Thinopyrum bessarabicum (Jb genome), which is phylogenetically close to D (sub)genome species, was taken as an outgroup. Subsequent data analysis using the pipeline RepeatExplorer2 allowed us to characterize the repeatomes of these species and identify several satellite sequences. Some of these sequences are novel, while others are found to be homologous to already known satellite sequences of Triticeae species. The copy number of satellite repeats in genomes of different species and their subgenome (D1 or Xcr) affinity in Ae. crassa were assessed by means of comparative bioinformatic analysis combined with quantitative PCR (qPCR). Fluorescence in situ hybridization (FISH) was performed to map newly identified satellite repeats on chromosomes of common wheat, Triticum aestivum, 4x and 6x Ae. crassa, Ae. tauschii, and Th. bessarabicum. The new FISH markers can be used in phylogenetic analyses of the Triticeae for chromosome identification and the assessment of their subgenome affinities and for evaluation of genome/chromosome constitution of wide hybrids or polyploid species.

Compensatory Effect of the ScGrf3-2R Gene in Semi-Dwarf Spring Triticale (x Triticosecale Wittmack).

The dwarfness in many triticale cultivars is provided by the dominant Ddw1 (Dominant dwarf 1) allele found in rye. However, along with conferring semi-dwarf phenotype to improve resistance to lodging, this gene also reduces grain size and weight and delays heading and flowering. Grf (Growth-regulating factors) genes are plant-specific transcription factors that regulate plant growth, including stem growth, in terms of length and thickness, and leaf and fruit size. In this work, we partially sequenced the rye gene ScGrf3 on chromosome 2R homologous to the wheat Grf3 gene, and found multiple polymorphisms in intron 3 and exon 4 complying with two alternative alleles (haplotypes ScGrf3-2Ra and ScGrf3-2Rb). For the identification of these, we developed a codominant PCR marker. Using a new marker, we studied the effect of ScGrf3-2R alleles in combination with the Ddw1 dwarf gene on economically valuable traits in F4 and F5 recombinant lines of spring triticale from the hybrid combination Valentin 90 x Dublet, grown in the Non-Chernozem zone for 2 years. Allele ScGrf3-2Ra was associated with greater thousand-grain weight, higher spike productivity, and earlier heading and flowering, which makes ScGrf3-2R a perspective compensator for negative effects of Ddw1 on these traits and increases prospects for its involvement in breeding semi-dwarf cultivars of triticale.

The Association of Grain Yield and Agronomical Traits with Genes of Plant Height, Photoperiod Sensitivity and Plastid Glutamine Synthetase in Winter Bread Wheat (Triticum aestivum L.) Collection.

The reduction in plant height caused by mutations in Rht-B1 or Rht-D1 (Reduced height-1) genes in combination with day-length-independent early flowering associated with the Ppd-D1 (Photoperiod-D1) gene were the main factors of the drastic yield increase in bread wheat in the 1960s. Increasing nitrogen use efficiency as well as maintaining high yields under conditions of global climate change are the modern goals of wheat breeding. The glutamine synthetase (GS) enzyme plays a key role in ammonium assimilation in plants. In previous studies, the TaGS2-A1 gene, coding the plastid isoform of GS, was shown to be connected with nitrogen use efficiency in wheat. Using the polymerase chain reaction (PCR) markers, the association of yield and agronomical traits with haplotypes of Rht-B1, Rht-D1, Ppd-D1 and TaGS2-A1 genes was studied in a diverse collection of winter bread wheat cultivars grown in Krasnodar (Russia). In the three-year experiment, semidwarfism and photoperiod insensitivity were confirmed to be highly favorable for the grain yield. The TaGS2-A1b haplotype had a tendency for increased grain yield and lodging resistance, but mainly in plants not possessing the 'green revolution' alleles. Thus, TaGS2-A1b may have potential in breeding wheat cultivars with alternative dwarfing genes or tall cultivars, which may be optimal for growing under certain environments.

Integrating Genetic and Chromosome Maps of Allium cepa: From Markers Visualization to Genome Assembly Verification.

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombinations. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Moreover, it is impossible to assign contigs lacking DNA markers. Here, we report the use of Tyr-FISH to determine the position of the short DNA sequence of markers and non-mapped unique copy sequence on the physical chromosomes of a large-genome onion (Allium cepa L.). In order to minimize potential background masking of the target signal, we improved our earlier developed pipeline for probe design. A total of 23 markers were located on physical chromosomes 2 and 6. The order of markers was corrected by the integration of genetic, pseudochromosome maps and cytogenetic maps. Additionally, the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. The results can be used for the improvement of pseudochromosome 2 and 6 assembly. The present study aims to demonstrate the value of the in situ visualization of DNA sequences in chromosome-scaffold genome assembly.

Developing Oligo Probes for Chromosomes Identification in Hemp (Cannabis sativa L.).

Hemp (Cannabis sativa L., 2n = 20) is a valuable crop that is successfully used as a food, technical and medicinal crop. It is a dioecious plant with an XX\XY sex determination system. Some chromosomes of C. sativa have almost the same lengths and centromeric indexes. Cytogenetic markers help to distinguish similar plant chromosomes, including sex chromosomes, which is important for the breeding process. Two repeats (CS-1 and CS-237) were used to develop labeled oligo-probes for rapid and low-cost oligo-FISH. These oligos can be recommended for use as cytological markers to distinguish sex chromosomes (X and Y) and somatic chromosome pairs 3, 6, and 8 by rapid oligo-FISH in a short time.

Study and Physical Mapping of the Species-Specific Tandem Repeat CS-237 Linked with 45S Ribosomal DNA Intergenic Spacer in Cannabis sativa L.

Hemp (Cannabis sativa L.) is a valuable crop and model plant for studying sex chromosomes. The scientific interest in the plant has led to its whole genome sequencing and the determination of its cytogenetic characteristics. A range of cytogenetic markers (subtelomeric repeat CS-1, 5S rDNA, and 45S rDNA) has been mapped onto hemp's chromosomes by fluorescent in situ hybridization (FISH). In this study, another cytogenetic marker (the tandem repeat CS-237, with a 237 bp monomer) was found, studied, and localized on chromosomes by FISH. The signal distribution and karyotyping revealed that the CS-237 probe was localized in chromosome 6 with one hybridization site and in chromosome 8 with two hybridization sites, one of which colocalizes with the 45S rDNA probe (with which a nucleolus organizer region, NOR, was detected). A BLAST analysis of the genomic data and PCR experiments showed that the modified CS-237 monomers (delCS-237, 208 bp in size) were present in the intergenic spacers (IGSs) of hemp 45S rDNA monomers. Such a feature was firstly observed in Cannabaceae species. However, IGS-linked DNA repeats were found in several plant species of other families (Fabaceae, Solanaceae, and Asteraceae). This phenomenon is discussed in this article. The example of CS-237 may be useful for further studying the phenomenon as well as for the physical mapping of hemp chromosomes.

qPCR as a Selective Tool for Cytogenetics.

qPCR is widely used in quantitative studies of plant genomes and transcriptomes. In this article, this method is considered as an auxiliary step in the preparation and selection of markers for FISH analysis. Several cases from the authors' research on populations of the same species were reviewed, and a comparison of the closely related species, as well as the adaptation of the markers, based on satellite tandem repeats (TRs) using quantitative qPCR data was conducted. In the selected cases, TRs with contrast abundance were identified in the cases of the Dasypyrum, Thinopyrum and Aegilops species, and the transfer of TRs between the wheat and related species was demonstrated. TRs with intraspecific copy number variation were revealed in Thinopyrum ponticum and wheat-wheatgrass partial amphidiploids, and the TR showing predominant hybridization to the sea buckthorn Y chromosome was identified. Additionally, problems such as the absence of a reference gene for qPCR, and low-efficiency and self-complementary primers, were illustrated. In the cases considered here, the qPCR results clearly show high correlation with the subsequent results of the FISH analysis, which confirms the value of this method for cytogenetic studies.

Transposons Hidden in Arabidopsis thaliana Genome Assembly Gaps and Mobilization of Non-Autonomous LTR Retrotransposons Unravelled by Nanotei Pipeline.

Long-read data is a great tool to discover new active transposable elements (TEs). However, no ready-to-use tools were available to gather this information from low coverage ONT datasets. Here, we developed a novel pipeline, nanotei, that allows detection of TE-contained structural variants, including individual TE transpositions. We exploited this pipeline to identify TE insertion in the Arabidopsis thaliana genome. Using nanotei, we identified tens of TE copies, including ones for the well-characterized ONSEN retrotransposon family that were hidden in genome assembly gaps. The results demonstrate that some TEs are inaccessible for analysis with the current A. thaliana (TAIR10.1) genome assembly. We further explored the mobilome of the ddm1 mutant with elevated TE activity. Nanotei captured all TEs previously known to be active in ddm1 and also identified transposition of non-autonomous TEs. Of them, one non-autonomous TE derived from (AT5TE33540) belongs to TR-GAG retrotransposons with a single open reading frame (ORF) encoding the GAG protein. These results provide the first direct evidence that TR-GAGs and other non-autonomous LTR retrotransposons can transpose in the plant genome, albeit in the absence of most of the encoded proteins. In summary, nanotei is a useful tool to detect active TEs and their insertions in plant genomes using low-coverage data from Nanopore genome sequencing.

The Development of New Species-Specific Molecular Markers Based on 5S rDNA in Elaeagnus L. Species.

The Elaeagnus L. species are trees and bushes that mainly grow in temperate zones of Western Europe; Minor, Central, and Southeast Asia; the Far East; and North America. Some species are used as fruit or ornamental plants and have economic value. Problems with the identification of species in the Elaeagnus genus by molecular genetical methods arise in the study of populations, systematics, breeding, and other areas of plant science and practice. Recently, the polymorphism of 5S ribosomal DNA non-transcribed spacers (5S rDNA NTSs) in Elaeagnaceae Adans. has been described. The results were used in our study as a basis for development of new species-specific molecular markers for some members of the Elaeagnus genus. The author's method was applied for finding regions that were potentially applicable for species-specific primer design. As a result, some species-specific molecular markers were developed for Elaeagnus angustifolia L., E. commutata Bernh., E. pungens Thunb., and E. multiflora Thunb. These markers were tested in a range of samples and showed the presence of amplified fragments in lanes of the marked species only. Samples of other species showed no amplifications. Thus, the developed markers may be useful for the species identification of the studied Elaeagnus plants in botanical, dendrological, and genetic research (especially in a leafless period of year), as well as in breeding and hybridization experiments.