Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

A CRISPR-based strategy for targeted sequencing in biodiversity science

Many applications in molecular ecology require the ability to match specific DNA sequences from single- or mixed-species samples with a diagnostic reference library. Widely used methods for DNA barcoding and metabarcoding employ PCR and amplicon sequencing to identify taxa based on target sequences, but the target-specific enrichment capabilities of CRISPR-Cas systems may offer advantages in some applications. We identified 54,837 CRISPR-Cas guide RNAs that may be useful for enriching chloroplast DNA across phylogenetically diverse plant species. We tested a subset of 17 guide RNAs in vitro to enrich plant DNA strands ranging in size from diagnostic DNA barcodes of 1,428 bp to entire chloroplast genomes of 121,284 bp. We used an Oxford Nanopore sequencer to evaluate sequencing success based on both single- and mixed-species samples, which yielded mean chloroplast sequence lengths of 2,530-11,367 bp, depending on the experiment. In comparison to mixed-species experiments, single-species experiments yielded more on-target sequence reads and greater mean pairwise identity between contigs and the plant species' reference genomes. But nevertheless, these mixed-species experiments yielded sufficient data to provide ≥48-fold increase in sequence length and better estimates of relative abundance for a commercially prepared mixture of plant species compared to DNA metabarcoding based on the chloroplast trnL-P6 marker. Prior work developed CRISPR-based enrichment protocols for long-read sequencing and our experiments pioneered its use for plant DNA barcoding and chloroplast assemblies that may have advantages over workflows that require PCR and short-read sequencing. Future work would benefit from continuing to develop in vitro and in silico methods for CRISPR-based analyses of mixed-species samples, especially when the appropriate reference genomes for contig assembly cannot be known a priori.

A DNA Extraction Method for Nondestructive Testing and Evaluation of Cotton Seeds (Gossypium L.)

Kernels of cotton provide lint and linter for textiles, oil and protein for food and feed. Cotton seed is formed following fertilization between an ovule and a pollen grain. The seed coat is maternal in origin, whereas the embryo and attached cotyledonary leaves are hybrids of parental lines. The extraction of genomic DNA from an ungerminated whole, a portion or mixed seeds are prerequisite in genetic and genomic studies of cotton. As far as our knowledge, there is only one method of nondescriptive DNA extraction from ungerminated cotton seeds without affecting the seed germination capability, but it has technical difficulties and requires special equipment. Furthermore, the amount of DNA extracted using the published method is low and, therefore, it is only suitable for routine marker assisted selection studies. In this study, a DNA extraction protocol referred to as the CTAB-LiCl was developed for single whole cotton seed, a portion of cotton seed and bulked cotton seeds. This protocol uses a combination of CTAB and LiCl to lyse cells and deplete RNAs simultaneously. The CTAB-LiCl DNA extraction method was evaluated in ninety-six individuals of six different cotton cultivars along with two genetic standards of cotton, TM-1 (G. hirsutum L.), Pima 3-79 (G. barbadense L.), and several other plant species of different plant genera. Results revealed that this method produced high quality and amounts of DNA as confirmed by spectrophotometry, agarose gel, restriction enzyme digestion, polymerase chain reaction, and library production for next generation sequencing studies of whole genome bisulfite sequencing. It does not require the use of liquid nitrogen, RNase, proteinase K, or beta-mercaptoethanol and can be completed in approximately 2 h. Small tissues of the chalaza ends of ungerminated cotton seeds could be used to obtain high quality and quantity of DNA ranging from 14 to 28 µg without affecting the seeds' germination ability, allowing marker-assisted selection before planting and flowering.

Correlation between antimicrobial resistance, biofilm formation, and virulence determinants in uropathogenic Escherichia coli from Egyptian hospital

BACKGROUND

Uropathogenic Escherichia coli (UPEC) is the main etiological agent behind community-acquired and hospital-acquired urinary tract infections (UTIs), which are among the most prevalent human infections. The management of UPEC infections is becoming increasingly difficult owing to multi-drug resistance, biofilm formation, and the possession of an extensive virulence arsenal. This study aims to characterize UPEC isolates in Tanta, Egypt, with regard to their antimicrobial resistance, phylogenetic profile, biofilm formation, and virulence, as well as the potential associations among these factors.

METHODS

One hundred UPEC isolates were obtained from UTI patients in Tanta, Egypt. Antimicrobial susceptibility was assessed using the Kirby-Bauer method. Extended-spectrum β-lactamases (ESBLs) production was screened using the double disk synergy test and confirmed with PCR. Biofilm formation was evaluated using the microtiter-plate assay and microscopy-based techniques. The phylogenetic groups of the isolates were determined. The hemolytic activity, motility, siderophore production, and serum resistance of the isolates were also evaluated. The clonal relatedness of the isolates was assessed using ERIC-PCR.

RESULTS

Isolates displayed elevated resistance to cephalosporins (90-43%), sulfamethoxazole-trimethoprim (63%), and ciprofloxacin (53%). Ninety percent of the isolates were multidrug-resistant (MDR)/ extensively drug-resistant (XDR) and 67% produced ESBLs. Notably, there was an inverse correlation between biofilm formation and antimicrobial resistance, and 31%, 29%, 32%, and 8% of the isolates were strong, moderate, weak, and non-biofilm producers, respectively. Beta-hemolysis, motility, siderophore production, and serum resistance were detected in 64%, 84%, 65%, and 11% of the isolates, respectively. Siderophore production was correlated to resistance to multiple antibiotics, while hemolysis was more prevalent in susceptible isolates and associated with stronger biofilms. Phylogroups B2 and D predominated, with lower resistance and stronger biofilms in group B2. ERIC-PCR revealed considerable diversity among the isolates.

CONCLUSION

This research highlights the dissemination of resistance in UPEC in Tanta, Egypt. The evident correlation between biofilm and resistance suggests a resistance cost on bacterial cells; and that isolates with lower resistance may rely on biofilms to enhance their survival. This emphasizes the importance of considering biofilm formation ability during the treatment of UPEC infections to avoid therapeutic failure and/or infection recurrence.

An atypical NLR gene confers bacterial wilt susceptibility in Arabidopsis

Quantitative disease resistance (QDR) remains the most prevalent form of plant resistance in crop fields and wild habitats. Genome-wide association studies (GWAS) have proved to be successful in deciphering the quantitative genetic basis of complex traits such as QDR. To unravel the genetics of QDR to the devastating worldwide bacterial pathogen Ralstonia solanacearum, we performed a GWAS by challenging a highly polymorphic local mapping population of Arabidopsis thaliana with four R. solanacearum type III effector (T3E) mutants, identified as key pathogenicity determinants after a first screen on an A. thaliana core collection of 25 accessions. Although most quantitative trait loci (QTLs) were highly specific to the identity of the T3E mutant (ripAC, ripAG, ripAQ, and ripU), we finely mapped a common QTL located on a cluster of nucleotide-binding domain and leucine-rich repeat (NLR) genes that exhibited structural variation. We functionally validated one of these NLRs as a susceptibility factor in response to R. solanacearum, named it Bacterial Wilt Susceptibility 1 (BWS1), and cloned two alleles that conferred contrasting levels of QDR. Further characterization indicated that expression of BWS1 leads to suppression of immunity triggered by different R. solanacearum effectors. In addition, we showed a direct interaction between BWS1 and RipAC T3E, and BWS1 and SUPPRESSOR OF G2 ALLELE OF skp1 (SGT1b), the latter interaction being suppressed by RipAC. Together, our results highlight a putative role for BWS1 as a quantitative susceptibility factor directly targeted by the T3E RipAC, mediating negative regulation of the SGT1-dependent immune response.

Development of a High-Quality/Yield Long-Read Sequencing-Adaptable DNA Extraction Method for Crop Seeds

Genome sequencing is important for discovering critical genes in crops and improving crop breeding efficiency. Generally, fresh, young leaves are used for DNA extraction from plants. However, seeds, the storage form, are more efficient because they do not require cultivation and can be ground at room temperature. Yet, only a few DNA extraction kits or methods suitable for seeds have been developed to date. In this study, we introduced an improved (IMP) Boom method that is relatively low-cost, simple to operate, and yields high-quality DNA that can withstand long-read sequencing. The method successfully extracted approximately 8 µg of DNA per gram of seed weight from soybean seeds at an average concentration of 48.3 ng/µL, approximately 40-fold higher than that extracted from seeds using a common extraction method kit. The A260/280 and A260/230 values of the DNA were 1.90 and 2.43, respectively, which exceeded the respective quality thresholds of 1.8 and 2.0. The DNA also had a DNA integrity number value (indicating the degree of DNA degradation) of 8.1, higher than that obtained using the kit and cetyltrimethylammonium bromide methods. Furthermore, the DNA showed a read length N50 of 20.96 kbp and a maximum read length of 127.8 kbp upon long-read sequencing using the Oxford Nanopore sequencer, with both values being higher than those obtained using the other methods. DNA extracted from seeds using the IMP Boom method showed an increase in the percentage of the nuclear genome with a decrease in the relative ratio of chloroplast DNA. These results suggested that the proposed IMP Boom method can extract high-quality and high-concentration DNA that can be used for long-read sequencing, which cannot be achieved from plant seeds using other conventional DNA extraction methods. The IMP Boom method could also be adapted to crop seeds other than soybeans, such as pea, okra, maize, and sunflower. This improved method is expected to improve the efficiency of various crop-breeding operations, including seed variety determination, testing of genetically modified seeds, and marker-assisted selection.

Research Progress of Nucleic Acid Detection Technology for Genetically Modified Maize

Genetically modified (GM) maize is one of the earliest GM crops to have achieved large-scale commercial cultivation globally, and it is of great significance to excel in the development and implementation of safety policy regarding GM, and in its technical oversight. This article describes the general situation regarding genetically modified maize, including its varieties, applications, relevant laws and regulations, and so on. From a technical point of view, we summarize and critically analyze the existing methods for detecting nucleic acid levels in genetically modified maize. The nucleic acid extraction technology used for maize is explained, and the introduction of traditional detection techniques, which cover variable-temperature and isothermal amplification detection technology and gene chip technology, applications in maize are described. Moreover, new technologies are proposed, with special attention paid to nucleic acid detection methods using sensors. Finally, we review the current limitations and challenges of GM maize nucleic acid testing and share our vision for the future direction of this field.

First whole genome assembly and annotation of a European common bean cultivar using PacBio HiFi and Iso-Seq data

Common bean (Phaseolus vulgaris L.) is the most important grain legume for direct human consumption worldwide. Flageolet bean originates from France and presents typical organoleptic properties, including the remarkable feature of having small pale green colored seeds. Here, we report the whole-genome data, assembly and annotation of the flageolet bean accession 'Flavert'. High molecular weight DNA and RNA were extracted and subjected to long-read sequencing using PacBio Sequel II platform. The genome consisted of 566,238,753 bp assembled in 13 molecules, including 11 chromosomes plus the mitochondrial and chloroplastic genomes. Annotation predicted 29,549 protein coding genes and 6,958 non-coding RNA. This high-quality genome (99.2% BUSCO completeness) represents a valuable data set for further genomic and genetic studies on common bean and more generally on legumes. To our knowledge, this is the first whole-genome sequence of a common bean accession originating from Europe.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Balancing read length and sequencing depth: Optimizing Nanopore long-read sequencing for monocots with an emphasis on the Liliales

Premise

We present approaches used to generate long-read Nanopore sequencing reads for the Liliales and demonstrate how modifications to standard protocols directly impact read length and total output. The goal is to help those interested in generating long-read sequencing data determine which steps may be necessary for optimizing output and results.

Methods

Four species of Calochortus (Liliaceae) were sequenced. Modifications made to sodium dodecyl sulfate (SDS) extractions and cleanup protocols included grinding with a mortar and pestle, using cut or wide-bore tips, chloroform cleaning, bead cleaning, eliminating short fragments, and using highly purified DNA.

Results

Steps taken to maximize read length can decrease overall output. Notably, the number of pores in a flow cell is correlated with the overall output, yet we did not see an association between the pore number and the read length or the number of reads produced.

Discussion

Many factors contribute to the overall success of a Nanopore sequencing run. We showed the direct impact that several modifications to the DNA extraction and cleaning steps have on the total sequencing output, read size, and number of reads generated. We show a tradeoff between read length and the number of reads and, to a lesser extent, the total sequencing output, all of which are important factors for successful de novo genome assembly.

Modified CTAB protocols for high-molecular-weight DNA extractions from ferns

Premise

Efficient protocols for extracting high-molecular-weight (HMW) DNA from ferns facilitate the long-read sequencing of their large and complex genomes. Here, we perform two cetyltrimethylammonium bromide (CTAB)-based protocols to extract HMW DNA and evaluate their applicability in diverse fern taxa for the first time.

Methods and Results

We describe two modified CTAB protocols, with key adjustments to minimize mechanical disruption during lysis to prevent DNA shearing. One of these protocols uses a small amount of fresh tissue but yields a considerable quantity of HMW DNA with high efficiency. The other accommodates a large amount of input tissue, adopts an initial step of nuclei isolation, and thus ensures a high yield in a short period of time. Both methods were proven to be robust and effective in obtaining HMW DNA from diverse fern lineages, including 33 species in 19 families. The DNA extractions mostly had high DNA integrity, with mean sizes larger than 50 kbp, as well as high purity (A260/A230 and A260/A280 > 1.8).

Conclusions

This study provides HMW DNA extraction protocols for ferns in the hope of facilitating further attempts to sequence their genomes, which will bridge our genomic understanding of land plant diversity.

A high quality, high molecular weight DNA extraction method for PacBio HiFi genome sequencing of recalcitrant plants

BACKGROUND

PacBio HiFi sequencing provides highly accurate long-read sequencing datasets which are of great advantage for whole genome sequencing projects. One limitation of the method is the requirement for high quality, high molecular weight input DNA. This can be particularly challenging for plants that frequently contain common and species-specific secondary metabolites, which often interfere with downstream processes. Cape Primroses (genus Streptocarpus), are some of these recalcitrant plants and are selected here as material to develop a high quality, high molecular weight DNA extraction protocol for long read genome sequencing.

RESULTS

We developed a DNA extraction method for PacBio HiFi sequencing for Streptocarpus grandis and Streptocarpus kentaniensis. A CTAB lysis buffer was employed to avoid guanidine, and the traditional chloroform and phenol purification steps were replaced with pre-lysis sample washes. Best cells/nucleus lysis was achieved with 4 h at 58 °C. The obtained high quality and high molecular weight DNAs were tested in PacBio SMRTBell™ library preparations, which resulted in circular consensus sequencing (CCS) reads from 17 to 27 Gb per cell, and a read length N50 from 14 to 17 kbp. To evaluate the quality of the reads for whole genome sequencing, they were assembled with HiFiasm into draft genomes, with N50 = 49 Mb and 23 Mb, and L50 = 10 and 11. The longest contigs were 95 Mb and 57 Mb respectively, showing good contiguity as these are longer than the theoretical chromosome length (genome size/chromosome number) of 78 Mb and 55 Mb, for S. grandis and S. kentaniensis respectively.

CONCLUSIONS

DNA extraction is a critical step towards obtaining a complete genome assembly. Our DNA extraction method here provided the required high quality, high molecular weight DNA for successful standard-input PacBio HiFi library preparation. The contigs from those reads showed a high contiguity, providing a good starting draft assembly towards obtaining a complete genome. The results obtained here were highly promising, and demonstrated that the DNA extraction method developed here is compatible with PacBio HiFi sequencing and suitable for de novo whole genome sequencing projects of plants.

Putative Signals of Generalist Plant Species Adaptation to Local Pollinator Communities and Abiotic Factors

The reproductive success of flowering plants with generalized pollination systems is influenced by interactions with a diverse pollinator community and abiotic factors. However, knowledge about the adaptative potential of plants to complex ecological networks and the underlying genetic mechanisms is still limited. Based on a pool-sequencing approach of 21 natural populations of Brassica incana in Southern Italy, we combined a genome-environmental association analysis with a genome scan for signals of population genomic differentiation to discover genetic variants associated with the ecological variation. We identified genomic regions putatively involved in the adaptation of B. incana to the identity of local pollinator functional categories and pollinator community composition. Interestingly, we observed several shared candidate genes associated with long-tongue bees, soil texture, and temperature variation. We established a genomic map of potential generalist flowering plant local adaptation to complex biotic interactions, and the importance of considering multiple environmental factors to describe the adaptive landscape of plant populations.

Assembling Quality Genomes of Flax Fungal Pathogens from Oxford Nanopore Technologies Data

Flax (Linum usitatissimum L.) is attacked by numerous devastating fungal pathogens, including Colletotrichum lini, Aureobasidium pullulans, and Fusarium verticillioides (Fusarium moniliforme). The effective control of flax diseases follows the paradigm of extensive molecular research on pathogenicity. However, such studies require quality genome sequences of the studied organisms. This article reports on the approaches to assembling a high-quality fungal genome from the Oxford Nanopore Technologies data. We sequenced the genomes of C. lini, A. pullulans, and F. verticillioides (F. moniliforme) and received different volumes of sequencing data: 1.7 Gb, 3.9 Gb, and 11.1 Gb, respectively. To obtain the optimal genome sequences, we studied the effect of input data quality and genome coverage on assembly statistics and tested the performance of different assembling and polishing software. For C. lini, the most contiguous and complete assembly was obtained by the Flye assembler and the Homopolish polisher. The genome coverage had more effect than data quality on assembly statistics, likely due to the relatively low amount of sequencing data obtained for C. lini. The final assembly was 53.4 Mb long and 96.4% complete (according to the glomerellales_odb10 BUSCO dataset), consisted of 42 contigs, and had an N50 of 4.4 Mb. For A. pullulans and F. verticillioides (F. moniliforme), the best assemblies were produced by Canu-Medaka and Canu-Homopolish, respectively. The final assembly of A. pullulans had a length of 29.5 Mb, 99.4% completeness (dothideomycetes_odb10), an N50 of 2.4 Mb and consisted of 32 contigs. F. verticillioides (F. moniliforme) assembly was 44.1 Mb long, 97.8% complete (hypocreales_odb10), consisted of 54 contigs, and had an N50 of 4.4 Mb. The obtained results can serve as a guideline for assembling a de novo genome of a fungus. In addition, our data can be used in genomic studies of fungal pathogens or plant-pathogen interactions and assist in the management of flax diseases.

Capturing marine microbiomes and environmental DNA: A field sampling guide

The expanding interest in marine microbiome and eDNA sequence data has led to a demand for sample collection and preservation standard practices to enable comparative assessments of results across studies and facilitate meta-analyses. We support this effort by providing guidelines based on a review of published methods and field sampling experiences. The major components considered here are environmental and resource considerations, sample processing strategies, sample storage options, and eDNA extraction protocols. It is impossible to provide universal recommendations considering the wide range of eDNA applications; rather, we provide information to design fit-for-purpose protocols. To manage scope, the focus here is on sampling collection and preservation of prokaryotic and microeukaryotic eDNA. Even with a focused view, the practical utility of any approach depends on multiple factors, including habitat type, available resources, and experimental goals. We broadly recommend enacting rigorous decontamination protocols, pilot studies to guide the filtration volume needed to characterize the target(s) of interest and minimize PCR inhibitor collection, and prioritizing sample freezing over (only) the addition of preservation buffer. An annotated list of studies that test these parameters is included for more detailed investigation on specific steps. To illustrate an approach that demonstrates fit-for-purpose methodologies, we provide a protocol for eDNA sampling aboard an oceanographic vessel. These guidelines can aid the decision-making process for scientists interested in sampling and sequencing marine microbiomes and/or eDNA.

Pushing the limits of HiFi assemblies reveals centromere diversity between two Arabidopsis thaliana genomes

Although long-read sequencing can often enable chromosome-level reconstruction of genomes, it is still unclear how one can routinely obtain gapless assemblies. In the model plant Arabidopsis thaliana, other than the reference accession Col-0, all other accessions de novo assembled with long-reads until now have used PacBio continuous long reads (CLR). Although these assemblies sometimes achieved chromosome-arm level contigs, they inevitably broke near the centromeres, excluding megabases of DNA from analysis in pan-genome projects. Since PacBio high-fidelity (HiFi) reads circumvent the high error rate of CLR technologies, albeit at the expense of read length, we compared a CLR assembly of accession Eyach15-2 to HiFi assemblies of the same sample. The use of five different assemblers starting from subsampled data allowed us to evaluate the impact of coverage and read length. We found that centromeres and rDNA clusters are responsible for 71% of contig breaks in the CLR scaffolds, while relatively short stretches of GA/TC repeats are at the core of >85% of the unfilled gaps in our best HiFi assemblies. Since the HiFi technology consistently enabled us to reconstruct gapless centromeres and 5S rDNA clusters, we demonstrate the value of the approach by comparing these previously inaccessible regions of the genome between the Eyach15-2 accession and the reference accession Col-0.