Optimized In Vitro Restriction Digestion Protocol for Preparing Maize and Barley ddRAD-Seq Libraries

: In recent years, high-throughput sequencing methods have become increasingly popular in molecular biology laboratories, mainly due to the relatively low cost of small, benchtop platforms, the simplicity of library preparation, and the low price per unit of information. Sequencing huge and complex genomes, such as cereal genomes, remains challenging and may not always be necessary. Therefore, several techniques have been developed to sequence a reduced representation of the genome. The most ﬂexible and widely used of these is ddRAD-Seq, which uses a pair of restriction enzymes to generate a pool of DNA fragments. The aim of this study was to validate in vitro the efﬁcacy of different combinations of restriction enzymes for ddRAD-Seq library construction in barley and maize. Eleven pairs of restriction enzymes were selected and tested to determine the concentrations of fragments with the expected length range and to select suitable pairs for sampling the genomes of these two cereals using ddRAD-Seq. For the selected pairs, i.e., PstI—MspI and HindIII—FspBI for barley and maize, respectively, libraries were prepared for NGS sequencing on Illumina MiSeq. Sequencing conﬁrmed the suitability of the selected enzymes to perform ddRAD-Seq in different genotypes. The results presented can be used for extensive research on these important cereal species.


Introduction
In the nearly 70 years since James Watson and Francis Crick discovered the structure of DNA [1], several different technologies have been developed to study the sequence of nucleotides found in DNA molecules.The first attempts at DNA sequencing date back to the 1950s, while the late 20th century witnessed intensive advances in sequencing methods, leading to the establishment of three generations of sequencing methods.However, whole genome sequencing of every specimen in the population of interest is still time-consuming, costly, and often unnecessary for many applications.The focus has, therefore, been on developing methods for reduced representation genome sequencing (RRS) of genomes.This approach has been extremely successful and has revolutionized research in evolutionary and conservation genetics and molecular ecology over the past decade [2].These methods allow reproducible pools of sequenced DNA fragments to be obtained from consecutive samples analyzed and provide information on single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs).Depending on the approach, reproducibility has been achieved by targeted sequencing of PCR products, library enrichment by probe hybridization or restriction enzyme (RE) fragmentation [3][4][5].
Technological advances over the past 15 years have made sequencing a routine and widely available analysis in molecular biology, and a milestone has been the invention of parallel, high-throughput next-generation sequencing (NGS) technology.The beginning of the 21st century exhibited a swift advancement in sequencing methods.Second-generation sequencing technologies, such as Illumina, have facilitated the assembly of over 200 plant Agronomy 2023, 13, 2956 2 of 15 genomes.Genomic development has enabled comprehensive analysis of plant genomes and examination of the genetic foundation of agronomic characteristics [6].Genome sequencing provides an opportunity to gain insight into the process of domestication, whereby human selection has resulted in the cultivation of plants that are more suitable for agriculture than for survival in natural environments.The comparison of genome sequences among plant species forms the basis for determining evolutionary relationships.The genome of Oryza sativa, commonly known as rice, was the first crop to undergo sequencing.This decision was made due to its importance as a major crop with a genome that is relatively small in size.The sequencing of rice subsequently served as a model for the sequencing of other cereal and grass genomes.Brachypodium distachyon, a model grass genome, was similarly sequenced and is relevant to the sequencing of the wheat genome [7].
Restriction enzyme-based methods use one or two restriction enzymes to cleave genomic DNA in the initial step of library preparation [8].One of these methods is doubledigest restriction site-associated DNA sequencing (ddRAD-Seq), which uses two different restriction enzymes (RE) to cut the DNA [9].The first RE is a rare cutter, i.e., an RE with a rare restriction recognition site of 6-8 bp in length, which is used to generate long fragments.The second is a common cutter, i.e., an RE with a much more common four bp restriction site, whose role is to reduce the length of the resulting fragments.This eliminates random shearing and end repair of genomic DNA and minimizes the number of fragments bound by the same restriction sites at both ends.In addition, the method uses precise size selection of genomic fragments, providing greater control over the proportion of fragments represented in the resulting library [9].As a result, ddRAD-Seq allows the generation of a library consisting of only a subset of hundreds to hundreds of thousands of DNA fragments generated after the double restriction cut.In addition, precise and reproducible size selection increases the ability to sequence the same regions from multiple samples and reduces the risk of sequencing adjacent regions.This is reflected in the minimization of the required sequencing depth for optimal coverage.The percentage of missing loci is also reduced compared to other methods, such as restriction site-associated DNA sequencing (RADseq) [9][10][11].The ddRAD-Seq technology can be successfully applied to the analysis of both species with known and unknown genome sequences [11,12].The reference genome sequence is an advantage right from the start, i.e., it offers the possibility to perform in silico restriction analysis to select REs on this basis [13].For species with an unknown genome sequence, it is important to optimize the protocol, especially regarding the selection of REs [14].Therefore, ddRAD-Seq is highly customizable in terms of the total number of loci, based on the RE selection and the range of fragment sizes selected.It allows for the production of a large number of fragments, thus enabling the identification of variants at tens of thousands of orthologous loci [15,16].Consequently, this method is used for genetic mapping, genetic diversity detection, population genomics, and molecular evolution studies [17][18][19][20][21]. Another application of ddRAD-Seq can be the detection of methylation levels using methylation-sensitive restriction enzymes, as described by Marconi et al. [22] and Dimond and Roberts [23].
The aim of this study was to experimentally confirm the efficacy of restriction enzymes used in constructing ddRAD-Seq libraries to evaluate genetic diversity in maize (Zea mays L.) and barley (Hordeum vulgare L.).In parallel, the results obtained by in silico restriction analysis were verified.Barley and maize are among the most important crops in terms of acreage and protein sources for both humans and livestock.As these species are of interest to scientists and breeders, it is crucial to efficiently break down their huge genomes, i.e., about 5.3 Gbp for barley and about 2.4 Gbp for maize, into fragments that are fractionally accurate for case-or gene-specific ddRAD-Seq library construction.Presenting the efficiency of different RE combinations is crucial, especially for low-budget laboratories that cannot afford an additional assay before implementing the ddRAD-Seq methodology.Based on this, the optimal set of enzymes can be selected depending on the expected coverage and depth of sequencing.

Plant Materials and DNA Extraction
The material for the study were barley (Hordeum vulgare sp.vulgare L.) accessions from the National Centre for Plant Genetic Resources, that is, the Polish Gene Bank and maize (Zea mays L.) inbred lines obtained from DSc El żbieta Kocha ńska-Czembor (Plant Breeding and Acclimatization Institute-NRI, Poland) (Table 1).The maize genotypes utilized in constructing the ddRad-Seq library were outlined in Czembor et al., 2019 [24].Detailed information on all barley accessions can be found in the EGISET database (https://wyszukiwarka.ihar.edu.pl/pl,accessed date: 27 October 2023).The materials used in the study's first phase included a maize line, 30IDT, with moderate resistance to ear rot caused by Fusarium verticillioides and a barley accession numbered 41268.This landrace was collected during an expedition from south-eastern Poland and introduced into the long-term storage in 1978.It was found to be resistant to steam rust, net blotch, and scald, and identified as Hordeum vulgare var.nutans.
DNA was isolated from the first leaf using the CTAB protocol followed by treatment with RNase A (10 mg/mL) (Thermo Fisher Scientific, Waltham, MA, USA) at 37 • C for 30 min [25].The isolated genetic material was dissolved in DNAse-free water, quantified using the NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).and qualified by agarose gel electrophoresis analysis.Samples were diluted to 20 ng/µL.

Evaluation of Restriction Enzymes Efficiency
The products of enzymatic digestion with the above combinations of restriction enzymes were analyzed using an Agilent 2100 Bioanalyzer with a High Sensitivity DNA kit.Within the 300-600 bp range, the concentration, molarity, and percentage of the resulting fragment pool were calculated for each sample.

Preparation of Adaptors
Oligonucleotide-forming adaptors were annealed prior to library construction to form the duplex of the P1 adaptor.The oligonucleotides P1.1/P1.2 and P2.1/P2.2 were designed in pairs for each barcode, respectively, to form their duplex.The reactions were carried out according to Table 3. Mixtures of P1 and P2 adaptors were incubated in a thermocycler at 97.5 • C for 2.

Library Construction and Sequencing
In the next step, fragments resulting from PstI-MspI and HindIII-FspBI digestion

Library Construction and Sequencing
In the next step, fragments resulting from PstI-MspI and HindIII-FspBI digestion for barley and maize genomic DNA, respectively, were ligated with adaptors compatible with the restriction site.A ligation reaction mix was prepared according to Table 4.One µL of P1 adaptor was added directly to 20 µL of digested DNA, and then the mixture was added to the ligation mix.Different adapters (P1) were used for each accession (Table S1), and then 9 µL of mix was added.The reaction was incubated at 16 • C for 2 h and then held at 4 • C. A diagram of ddRAD-Seq library preparation can be found in Figure S1.DNA fragments were then purified from the post-reaction components using 1× Ampure XP beads (Beckman Coulter, Pasadena, CA, USA), and the purified product was suspended in 20 µL of water.The purified library was amplified with 5 µM Illumina PCR1, PCR2 primers, and 25 µL NebNext Ultra II Q5 master mix.The reaction was denatured at 98 • C for 30 s.This was followed by 16 cycles of denaturation at 98 • C for 10 s and combined with annealing at 65 • C for 75 s.The next step was a final extension at 65 • C for 5 min and a final incubation at 4 • C. For this purpose, a PCR mix was prepared.
The library was again purified with 1× Ampure XP beads.The purified product was then suspended in 20 µL of water.Size selection was then performed using a 1.5% DF cassette with Marker K Pippin Prep.The prepared libraries were analyzed qualitatively and quantitatively using the Qubit fluorimeter (ThermoFisher Scientific) with the Qubit dsDNA HS Assay Kit (ThermoFisher Scientific) and the Bioanalyzer 2100 automated electrophoresis system with the High Sensitivity DNA Kit (Agilent, Santa Clara, CA, USA).For normalization, the library was diluted to 4nM and then denatured with 0.2N NaOH.After normalization, the libraries were sequenced on MiSeq (Illumina, San Diego, CA, USA) at 2 × 150 bp paired-end using the MiSeq Reagent Kit v3 (600 cycles).Raw reads were qualitatively evaluated using FastQC software v 0.11.9 [26].

Efficacy of Restriction Enzymes
For barley and maize, RE efficacy was evaluated separately.A fragment fraction of 300-700 bp is often used to construct ddRAD-Seq libraries, but it can be modified depending on the case, enzyme, species, etc. [29,30].In general, the size of fragments in the library should not be less than 200 bp to avoid overlapping of paired-end sequences, which can hinder bioinformatic analysis [31].However, to avoid degradation of base quality in Illumina paired-end sequencing, fragments should also not exceed 800 bp [32]; a selection of fragments in the range of 300 to 600 bp was used here.It is worth noting here that the parameters obtained from automated electrophoresis show the aggregate presence of all types of fragments that can be generated during digestion with a pair of REs, i.e., those with the same restriction sites specific for one of the REs at both ends, two different restriction sites, and fragments that were mechanically broken off during DNA handling.In silico analysis, on the other hand, only shows the presence of fragments with different restriction sites at both ends.Such fragments are indeed sequenced later.

Maize
Restriction of maize DNA using seven RE in 11 combinations produced a diverse range of fragment concentrations and molarities.The distribution of fragments obtained by 11 pairs of RE is shown in Figure 2. Based on the results of automated electrophoresis, the concentration, molarity, and percentage of the desired fraction of fragments were determined (Table 5).The highest values of molarity, concentration, and percentage of fragments of appropriate length were observed for a HindIII-FspBI, i.e., above 124,398 pmol/L 48,290 pg/µL and 46%, respectively.Acceptable fragmentation results were also obtained for the PstI-FspBI.Low fragment concentrations in the range of 300-600 bp were observed after using PstI-MseI, EcoRI-MspI, and HindIII-MspI RE pairs.This means that no fragments of this length were obtained.
In silico analysis showed that the most fragments in the expected range of 300-600 bp should be obtained using HindIII-FspBI.This was confirmed by in vitro analysis.In contrast, for RE pairs such as PstI-MseI and HindIII-MspI, significant discrepancies were observed between the wet analysis performed in the laboratory and the computer simulation.Fragmentation efficiency with PstI-MseI was unsatisfactory for the sample tested here, although in silico results indicated a better yield.

Barley
The effectiveness of the different RE combinations was also investigated for barley DNA.Different sizes and amounts of fragments were also observed in the chromatograms resulting from automated on-chip gel electrophoresis (Figure 3).The highest concentration and molarity of fragments in the range of 300-600 bp occurred for EcoRI-FspBI (Table 6), while the lowest efficiency was shown by the EcoRI-MspI pair.The highest proportion of fragments with the required length was obtained in the sample treated with PstI-MspI.Meanwhile, the results of the in silico analysis indicated that the HindIII-MseI RE pair should be the most efficient in digesting the barley genome into 300-600 bp fragments.In silico treatment of the reference genome with PstI-MspI RE was identified as the least effective.

NGS Sequencing
Based on the results of the efficiency analysis of the RE pairs and the literature data, PstI-MspI and HindIII-FspBI for barley and maize, respectively, were selected for further testing.To verify that the selected restriction systems were effective regardless of the genotype tested, the number of DNA samples was increased to 17 for wheat and 24 for barley.From the resulting fragment mixtures after restriction digestion, a total of 41

NGS Sequencing
Based on the results of the efficiency analysis of the RE pairs and the literature data, PstI-MspI and HindIII-FspBI for barley and maize, respectively, were selected for further testing.To verify that the selected restriction systems were effective regardless of the genotype tested, the number of DNA samples was increased to 17 for wheat and 24 for barley.From the resulting fragment mixtures after restriction digestion, a total of 41 ddRAD-Seq libraries were prepared and subjected to NGS sequencing.After sequencing, a total of 17,648,717 raw reads were obtained for maize and 23,876,096 for barley.The Phred score for all reads from the ddRAD-Seq libraries was above 30 for the entire length of the fragments (Figure 4).Bioinformatic analysis and SNP calling revealed that approximately 50,500 fragments with allele frequencies ≥0.05 and at least 5× [28] coverage were obtained for barley.For maize, 10,500 such fragments were obtained.

Discussion
The ddRAD-Seq method can be an essential and simple answer to the needs of scientists and breeders.Due to its high reproducibility and low cost per sample, the method can be an ideal analytical solution for teams with tighter budgets and struggling to analyze large plant genomes [9].The use of a pair of RE during genome sampling, followed by adaptor ligation and magnetic bead purification, generates reproducible genomic targets [33].The ddRAD-Seq method allows the sequencing of the same loci within the list of sequencing samples of complex genomes with a high rate of repeatable fragments and accurate mapping [34].This approach increases the possibility for easy and cost-effective development of molecular markers.

Discussion
The ddRAD-Seq method can be an essential and simple answer to the needs of scientists and breeders.Due to its high reproducibility and low cost per sample, the method can be an ideal analytical solution for teams with tighter budgets and struggling to analyze large plant genomes [9].The use of a pair of RE during genome sampling, followed by adaptor ligation and magnetic bead purification, generates reproducible genomic targets [33].The ddRAD-Seq method allows the sequencing of the same loci within the list of sequencing samples of complex genomes with a high rate of repeatable fragments and accurate mapping [34].This approach increases the possibility for easy and cost-effective development of molecular markers.
Although the original ddRAD-Seq protocol developed by Peterson et al. [9] used EcoRI-MspI, the method needs to be optimized for each species and/or application, especially in terms of RE selection, as genomes vary in size, complexity, and density of restriction sites.The appropriate pair of RE allows efficient sequencing by constructing libraries that contain precise DNA length fragments and cover the appropriate portion of the genome.This optimization should be prepared for the species, DNA fragment range, sequencing platform, and strategy [11].In general, the goal is to generate an appropriate number of sequencable fragments, according to the downstream analysis of interest.
In silico analysis is often used to pre-identify the optimal combination of RE for genome fragmentation in the ddRAD-Seq protocol.Peterson et al. [9] performed an in silico restriction analysis for 17 species and 5 pairs of RE.Maize was one of the species analyzed.The researchers estimated that EcoRI-MspI digestion would yield 30,000 fragments.Yang et al. [33] indicated the potential presence of about 120,000 fragments for the same REs.The in silico analysis presented here predicted more than 165,000 fragments.Similar inconsistencies appeared in the abundance of fragments obtained by PstI-MspI restriction.According to Yang et al. [33], this combination should generate a yield of about 130,000 fragments, while here, it was estimated to yield about 160,000 fragments.Differences in the in silico estimation of the number of fragments are due to the analysis of different maize genome assemblies, variable sequencing parameters, different length ranges of the fragments analyzed, and the use of different bioinformatics tools.Therefore, the results of the in silico analysis should be treated with appropriate caution, and the possibility that the results of the in vitro analysis may not agree with the estimate should be considered.This is illustrated by the comparison of parameters obtained during the automated electrophoresis of maize DNA.Here, the results clearly indicated that the efficiency of EcoRI-MspI digestion was very low for the expected fragment length, and the efficiency of PstI-MspI was even below the detection threshold of the Agilent Bioanalyzer 2100 with the High Sensitivity DNA chip loaded.However, there is evidence of a high agreement between in silico estimated and laboratory-analyzed results [9,35].Based on both types of results, further analysis of the maize lines was performed using HindIII-FspBI and satisfactory NGS results were obtained.Other RE combinations that have been successfully used in reduced representation protocols for maize genome sequencing can also be found in the literature, i.e., SphI-MluCI [36].
Furthermore, in barley, different teams have used different combinations of RE, i.e., PstI-MseI and PstI-MspI, to effectively reduce genome complexity [37,38].Here, the concentration and molarity of these two combinations were similar in the desired length range.However, the percentage of 300-600 bp fragments was higher for PstI-MspI, while the in silico analysis indicated more than double the number of fragments for PstI-MseI.The PstI-FspBI combination, which, according to in silico estimation, should generate a comparable number of fragments to PstI-MspI in in vitro analysis, was characterized by a significantly lower digestion efficiency.Thus, discrepancies between laboratory analysis and computer simulations can be observed, as in the case of maize.High efficiencies were observed for EcoRI-FspBI and HindIII-MseI.However, considering the size of the barley genome (5.3 Gbp), the planned sequencing depth (1.25 million reads/sample) and the aim of the study (diversity analysis of the landrace collection), the most efficient RE combinations were discarded in order to obtain a higher fraction of fragments with at least a five-fold sequence coverage.Finally, the combination PstI-MspI was chosen because of the possibility of comparing the results with those obtained by Milner et al. for 22,626 accessions from the gene bank [37].
As can be seen, the selection of RE for library construction is not a trivial matter.Different teams have used different pairs of RE to analyze even the same species.The efficiency of the restriction digestion in generating the appropriate fragment length is one of the most important conditions for success in constructing a library containing only a representation of the genome.However, the selection of RE for analysis should also be guided by the availability of data from the experimental work by other teams.Using the same RE greatly increases the chances of obtaining compatible results from different experiments, facilitating co-analysis and inference.The lower number of fragments generated by this RE pair also leads to a reduction in sequencing costs.A lower total number of reads per sample is required to achieve adequate coverage of the resulting fragments.
Another factor to consider is the methylation sensitivity of the RE and its consequences.The PstI-MspI RE pair used in the barley analysis consists of both a methylation-sensitive rare-cutting RE and a common-cutting RE.As the results show, this arrangement has certain advantages, such as a combination resulted in greater uniformity in the depth of reads across loci and provided a higher quality of acquired genetic information compared to pairs containing no methylation-sensitive RE as well as those containing only one methylationsensitive RE [39].In addition, the use of a methylation-sensitive RE pair increases the correlation between SNP and gene density.This may be a derivative of the methylation sensitivity itself, which increases the pool of under-methylated regions of the genome in the pool of sequenced fragments and their association with a higher gene density [39,40].However, the methylation sensitivity of the RE also has drawbacks.In maize, all RE combinations containing MspI showed low efficiency in generating fragments in the length range analyzed.MspI is completely blocked by hemi-and fully methylated mCCGG and mCmCGG sequences.Thus, a highly methylated genome may cause a low cutting efficiency and significantly reduce the number of fragments useful for ddRAD-Seq analysis of maize.However, attention to the efficiency of generating the required fragment length remains critical.If an enzyme pair does not generate fragments of the appropriate length, as in the case of PstI-MseI in maize, it does not matter what the longer/shorter fragments contain and how often they occur, as they will be discarded and not sequenced anyway due to the fragment length selection step in the ddRAD-Seq protocol.
The inconsistency of in silico and experimental genome restriction analysis can arise from several factors.In silico tools use algorithms to predict restriction sites by recognizing restriction enzyme sequence motifs.However, the accuracy of these algorithms is limited, particularly when dealing with complex or repetitive DNA sequences [41].The accuracy of in silico restriction analysis depends on the quality of the genome sequence used as input.If the genome sequence is incomplete or inaccurate, it may result in the misidentification or omission of true restriction sites [41].The substantial genetic distance between the reference and the genotype(s) being examined could also contribute to the observed discrepancies [27].This could lead to the existence of polymorphisms in restriction sites, as well as the emergence of new restriction sites [42,43].Moreover, natural adaptations to specific environmental conditions may lead to structural changes.Such changes may also arise from natural and artificial introgression of genome fragments from another species [44,45].Additionally, mobile genetic elements, both in quantity and type, significantly contribute to these modifications of the genome [46].DNA modifications, such as methylation, can modify recognition sequences for restriction enzymes.These modifications can hinder cleavage at forecasted restriction sites, thereby resulting in inconsistencies between in silico and experimental outcomes [27,47].The discrepancy between the in silico predictions and the experimental results may also be attributed to varying experimental conditions, such as temperature, pH, and reaction time, which could impact the activity of the restriction enzyme.Under conditions that are far from optimal for REs activity, they can cleave sequences that are similar but not identical to their canonical recognition sequences.This altered specificity has been termed "star activity".Star sites are associated with a canonical recognition site, but usually differ by one or more bases.This results in non-specific fragments which presence was not considered in the in silico prediction [48].Initial DNA quality is also a crucial determinant for sequencing reduced representative fractions of the genome.The number of fragments obtained, as well as the depth of sequencing, significantly decrease only after significant genomic DNA degradation.This effect is observed alongside a decrease in the number of SNPs [49].Another factor, that could cause discrepancies between laboratory analysis and computer simulations is the presence of non-target DNA, i.e., residuals of microbial DNA.Restriction enzymes may recognize and cleave such sequences and can lead to the generation of unexpected fragments [50].The final outcome of these factors could be an inadequate initial estimation of the number of fragments to be included in the library, resulting in a different number of reads and sequencing depth than expected.In instances with low sequencing depth, heterozygous loci might be miscategorized as sequencing errors during the analysis of results.In the context of research on genetic variation or population structure, this may result in a falsely elevated level of homozygosity [51].This could potentially suggest the presence of a bottleneck phenomenon that did not actually occur.Insufficient fragments of the appropriate length may result in too few reads, leading to errors in downstream applications such as gene mapping or, in the worst case, making the analysis of the sequences obtained completely impossible.In addition, selecting an inappropriate enzyme pair might lead to inadequate or partial digestion of DNA, significantly increasing the likelihood of poorly reproducible results [52].
The implementation of a ddRAD-Seq methodology, in which restriction enzymes are selected to produce a large number of fragments in the range of 300-600 bp, allows the generation of an effective ddRAD-Seq library.By sequencing such libraries, substantial amounts of data are obtained, which then allows for an in-depth analysis of the variability of SNP types.The outcomes derived from the optimization process, which helped identify the optimal ERs for maize and barley, were deployed during the study conducted by Czembor et al. in 2019 [24].A study on the use of specific REs in barley is being prepared for publication.

Conclusions
For optimal sequencing results, the restriction REs used in ddRAD-Seq analysis should be matched to the genome under study.This is a key step in ddRAD-Seq library construction that affects the number of fragments generated for sequencing.The in silico analyses used to assess the number of genomic fragments resulting from enzymatic digestion may differ from the in vitro results of experimental studies.There is no universal pair of restriction REs that can be recommended for ddRAD-Seq library construction for all species to produce an efficient number of fragments.Therefore, to minimize the inconsistencies between in silico and experimental genome restriction analysis results, it is important to use highquality genome sequences, carefully control experimental conditions, and validate in silico predictions with experimental data whenever possible.
5 min and then cooled to 21 • C at a rate of no more than 3 • C per minute, held at 21 • C for 10 min, and stored at 4 • C. A diagram of adaptor ligation is shown in Figure 1.Agronomy 2023, 13, x FOR PEER REVIEW 5 of 15

Figure 1 .
Figure 1.Scheme of double digestion of genomic DNA and ligation of adaptors.

Figure 1 .
Figure 1.Scheme of double digestion of genomic DNA and ligation of adaptors.

Figure 2 .
Figure 2. Results of enzymatic digestion of maize DNA.Automated electrophoresis was performed on an Agilent 2100 Bioanalyzer using the High Sensitivity DNA kit.(A) Agilent Bioanalyzer chromatograms of DNA fragment size (FU-Arbitrary fluorescence units); (B) Agilent 2100 Bioanalyzer gel-like image of DNA after automated electrophoresis.

Figure 2 .
Figure 2. Results of enzymatic digestion of maize DNA.Automated electrophoresis was performed on an Agilent 2100 Bioanalyzer using the High Sensitivity DNA kit.(A) Agilent Bioanalyzer chromatograms of DNA fragment size (FU-Arbitrary fluorescence units); (B) Agilent 2100 Bioanalyzer gel-like image of DNA after automated electrophoresis.

Figure 3 .
Figure 3. Results of enzymatic digestion of barley DNA.Automated on-chip electrophoresis was performed on an Agilent 2100 Bioanalyzer using the High Sensitivity DNA kit.(A) Agilent Bioanalyzer chromatograms of DNA fragment size (FU-Arbitrary fluorescence units); (B) Agilent 2100 Bioanalyzer gel-like image of DNA after automated electrophoresis.

Figure 3 .
Figure 3. Results of enzymatic digestion of barley DNA.Automated on-chip electrophoresis was performed on an Agilent 2100 Bioanalyzer using the High Sensitivity DNA kit.(A) Agilent Bioanalyzer chromatograms of DNA fragment size (FU-Arbitrary fluorescence units); (B) Agilent 2100 Bioanalyzer gel-like image of DNA after automated electrophoresis.

Agronomy 2023 ,
13,  x FOR PEER REVIEW 10 of 15 approximately 50,500 fragments with allele frequencies ≥0.05 and at least 5×[28] coverage were obtained for barley.For maize, 10,500 such fragments were obtained.

Figure 4 .
Figure 4. Quality of reads obtained from sequencing ddRAD-Seq maize (A) and barley libraries (B).

Figure 4 .
Figure 4. Quality of reads obtained from sequencing ddRAD-Seq maize (A) and barley libraries (B).

Table 1 .
List of materials used in the construction of the ddRAD-Seq libraries.

Table 3 .
Composition of the adaptors duplex.

Table 3 .
Composition of the adaptors duplex.

Table 4 .
Composition of the adaptors' ligation reaction.

Table 5 .
Results of in vitro and in silico enzymatic digestion of the maize genome for fragments in the range of 300-600 bp.

Table 6 .
Results of in vitro and in silico enzymatic digestion of the barley genome for fragments in the range of 300-600 bp.