<figure><img src="image_1.png" /></figure>

HAL
open science

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

Alessia Russo, Baptiste Mayjonade, Daniel Frei, Giacomo Potente, Roman T Kellenberger, Léa Frachon, Dario Copetti, Bruno Studer, Jürg E Frey, Ueli Grossniklaus, et al.

► To cite this version:

Alessia Russo, Baptiste Mayjonade, Daniel Frei, Giacomo Potente, Roman T Kellenberger, et al. Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families. Frontiers in Plant Science, 2022, 13, 12 p. 10.3389/fpls.2022.883897. [https://doi.org/10.3389/fpls.2022.883897](https://doi.org/10.3389/fpls.2022.883897). [hal-03748027](https://HAL.inrae.fr/hal-03748027)

HAL Id: hal-03748027

https://hal.inrae.fr/hal-03748027

Submitted on 9 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

University of
Zurich UZH

Zurich Open Repository and
Archive

University of Zurich
University Library
Strickhofstrasse 39
CH-8057 Zurich
www.zora.uzh.ch

Year: 2022

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

Russo, Alessia ; Mayjonade, Baptiste ; Frei, Daniel ; Potente, Giacomo ; Kellenberger, Roman T ; Frachon, Léa ; Copetti, Dario ; Studer, Bruno ; Frey, Jürg E ; Grossniklaus, Ueli ; Schlüter, Philipp M

**Abstract:** Long-read DNA sequencing technologies require high molecular weight (HMW) DNA of adequate purity and integrity, which can be difficult to isolate from plant material. Plant leaves usually contain high levels of carbohydrates and secondary metabolites that can impact DNA purity, affecting downstream applications. Several protocols and kits are available for HMW DNA extraction, but they usually require a high amount of input material and often lead to substantial DNA fragmentation, making sequencing suboptimal in terms of read length and data yield. We here describe a protocol for plant HMW DNA extraction from low input material (0.1 g) which is easy to follow and quick (2.5 h). This method successfully enabled us to extract HMW from four species from different families (Orchidaceae, Poaceae, Brassicaceae, Asteraceae). In the case of recalcitrant species, we show that an additional purification step is sufficient to deliver a clean DNA sample. We demonstrate the suitability of our protocol for long-read sequencing on the Oxford Nanopore Technologies PromethION<sup>6</sup> platform, with and without the use of a short fragment depletion kit.

DOI: https://doi.org/10.3389/fpls.2022.883897

Posted at the Zurich Open Repository and Archive, University of Zurich

ZORA URL: https://doi.org/10.5167/uzh-219157

Journal Article

Published Version

<figure><img src="image_4.png" /></figure>

The following work is licensed under a Creative Commons: Attribution 4.0 International (CC BY 4.0) License.

Originally published at:

Russo, Alessia; Mayjonade, Baptiste; Frei, Daniel; Potente, Giacomo; Kellenberger, Roman T; Frachon, Léa; Copetti, Dario; Studer, Bruno; Frey, Jürg E; Grossniklaus, Ueli; Schlüter, Philipp M (2022). Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families. Frontiers in Plant Science, 13:883897.

DOI: https://doi.org/10.3389/fpls.2022.883897
<figure><img src="image_5.png" /></figure>

Check for
updates

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

Alessia Russo<fnref n="1" />,<fnref n="2" />,³*, Baptiste Mayjonade<fnref n="4" />, Daniel Frei⁵, Giacomo Potente³, Roman T. Kellenberger⁶, Léa Frachon³, Dario Copetti⁷, Bruno Studer⁷, Jürg E. Frey⁵, Ueli Grossniklaus<fnref n="1" /> and Philipp M. Schlüter<fnref n="2" />,³*

## OPEN ACCESS

**Edited by:**

Ruslan Kalendar,
University of Helsinki, Finland

**Reviewed by:**

Cássio Van Den Berg,
State University of Feira de Santana,
Brazil

Ken Naito,
National Agriculture and Food
Research Organization (NARO), Japan

*Correspondence:

Alessia Russo
alessia.russo@systbot.uzh.ch

Philipp M. Schlüter
philipp.schlueter@uni-hohenheim.de

**Specialty section:**

This article was submitted to Technical Advances in Plant Science, a section of the journal Frontiers in Plant Science

Received: 25 February 2022

Accepted: 21 April 2022

Published: 19 May 2022

**Citation:**

Russo A, Mayjonade B, Frei D, Potente G, Kellenberger RT, Frachon L, Copetti D, Studer B, Frey JE, Grossniklaus U and Schlüter PM (2022) Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families. *Front. Plant Sci.* 13:883897.
doi: 10.3389/fpls.2022.883897

Long-read DNA sequencing technologies require high molecular weight (HMW) DNA of adequate purity and integrity, which can be difficult to isolate from plant material. Plant leaves usually contain high levels of carbohydrates and secondary metabolites that can impact DNA purity, affecting downstream applications. Several protocols and kits are available for HMW DNA extraction, but they usually require a high amount of input material and often lead to substantial DNA fragmentation, making sequencing suboptimal in terms of read length and data yield. We here describe a protocol for plant HMW DNA extraction from low input material (0.1 g) which is easy to follow and quick (2.5 h). This method successfully enabled us to extract HMW from four species from different families (Orchidaceae, Poaceae, Brassicaceae, Asteraceae). In the case of recalcitrant species, we show that an additional purification step is sufficient to deliver a clean DNA sample. We demonstrate the suitability of our protocol for long-read sequencing on the Oxford Nanopore Technologies PromethION® platform, with and without the use of a short fragment depletion kit.

**Keywords:** DNA extraction, DNA sequencing, nanopore sequencing, Circulomics, plant genome, ONT long read sequencing, PacBio, genome assembly

## INTRODUCTION

Long-read sequencing technologies have reshaped the research landscape of plant biology over the last few years. With the recent increase in sequencing read length, decrease in sequencing cost, and newly developed bioinformatics tools suitable for these technologies, de novo assembly of large and complex plant genomes of non-model species is now feasible (Jiao and Schneeberger, 2017; Kersey, 2019). This offers unprecedented opportunities to investigate genome structure and function, and focus on molecular and evolutionary questions in organisms that were previously inaccessible (Belser et al., 2018). We are now gaining a deeper understanding of genomic diversity,

<footnote n="1">¹ Department of Plant and Microbial Biology and Zurich-Basel Plant Science Centre, University of Zurich, Zurich, Switzerland,</footnote>

<footnote n="2">² Department of Plant Evolutionary Biology, Institute of Biology, University of Hohenheim, Stuttgart, Germany, ³ Department of Systematic and Evolutionary Botany and Zurich-Basel Plant Science Centre, University of Zurich, Zurich, Switzerland,</footnote>

<footnote n="4">⁴ Laboratoire des Interactions Plantes Microbes Environnement (LIPME), INRAE, Toulouse, France, ⁵ Department of Method Development and Analytics, Agroscope, Wädenswil, Switzerland, ⁶ Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom, ⁷ Institute of Agricultural Sciences and Zurich-Basel Plant Science Centre, ETH Zürich, Zurich, Switzerland</footnote>
evolution, and gene function by sequencing more genomes at a higher resolution (Zhang et al., 2017; Chawla et al., 2021). Lately, the possibility to release high quality reference genome assemblies has led to initiatives such as the European Reference Genome Atlas (Formenti et al., 2022) or the Earth BioGenome Project, which “aims to sequence, catalog and characterise the genomes of all of Earth’s eukaryotic biodiversity,” to study evolution and preserve biodiversity (Exposito-Alonso et al., 2020).

Long-read sequencers are able to generate reads of 10 kbp or longer. The recently developed PacBio® HiFi technology from Pacific Biosciences can provide reads up to 25 kbp (the older CLR mode provides longer but less accurate reads)<fnref n="1" />, while Oxford Nanopore Technologies® (ONT®) nanopore sequencing technology provides the longest reads, up to the current record of 4.2 Mbp<fnref n="2" />. Such long reads are able to unambiguously capture complex and repetitive regions in plant genomes, allowing the exploration of genomic regions that were previously inaccessible (Belser et al., 2018; Goerner-Potvin and Bourque, 2018). This, together with their ability to resolve highly heterozygous regions, has enabled the assembly of large plant genomes at the chromosome level (Hu et al., 2019; Hasing et al., 2020; Pu et al., 2020; Niu et al., 2022).

Long reads are also reshaping the way we approach population genetic studies. Structural Variants (SVs) represent a major form of genetic variation, and may contribute to phenotypic variation as much as – or even more so than – single nucleotide polymorphisms (SNPs) (Chawla et al., 2021). However, it is challenging to reliably detect large SVs using short-read sequencing (Saxena et al., 2014). Structural variants can now be captured with long reads, thus enabling the fine-scale characterisation of genomic rearrangements responsible for trait variation in plants (Zhang et al., 2016; Sedlazeck et al., 2018; Heller and Vinron, 2019). Moreover, since PacBio® and ONT® are able to detect chemical modifications on nucleotides, they provide a new method to directly profile patterns of DNA methylation across genomes and allow epigenetic studies (Flusberg et al., 2010; Simpson et al., 2017). However, to fully exploit the potential of long-read sequencing, it is critical to obtain high molecular weight (HMW) DNA of adequate purity and integrity.

Extraction of HMW DNA from plant material can be challenging. First, plant cells have a cell wall composed of polysaccharide polymers, such as cellulose and pectin, as well as glycoproteins and lignin (Zhang et al., 2021), making the cell wall rigid and hard to break. Thus, steps that achieve effective mechanical disruption of the cell wall are necessary. As a chemical defence against herbivores, plants also produce polysaccharides and phenols, which tend to accumulate in leaves and, upon cell lysis, can bind DNA and affect downstream molecular analyses (Katterman and Shattuck, 1983; Varma et al., 2007; Moreira et al., 2011). The presence of polysaccharides has been

shown to inhibit restriction enzyme activity (Pandey et al., 1996). Thus, purification of DNA from plant material requires careful optimisation. Several commercial DNA isolation kits and protocols are available on the market. Many protocols rely on the isolation of nuclei with an osmotic nuclear isolation buffer, and subsequent lysis of nuclear membranes with a detergent to release DNA. These methods are time-consuming, hazardous, and/or require high amounts of input material (Zerpa-Catanho et al., 2021). Furthermore, the resulting genomic DNA (gDNA) is often highly oxidised and, therefore, unsuitable for long-read applications. Kit-based extraction methods are offered by several companies, and are intended to easily remove contaminants; but they are costly and there is a risk of losing DNA during column washes. A previously published protocol (Mayjonade et al., 2016) presented a method to extract plant HMW DNA via a sodium dodecyl sulphate (SDS)-based lysis buffer and magnetic bead-based purification. The described method is easy and quick, taking only 1.5 h to complete DNA isolation from harvested plant material. Moreover, it requires only 0.1 g starting material. Unfortunately, this method did not yield DNA of sufficient quality for long-read sequencing in the plant species used as study organisms by our groups.

Here, we present an improved HMW DNA extraction method based upon the protocol by Mayjonade et al. (2016). We introduced two simple but effective major modifications: the addition of β-mercaptoethanol, which prevents oxidative damage to nucleic acids (Gerstein, 2001) and prevents nuclease activity, and a phenol:chloroform extraction. To demonstrate the robustness of the method, we applied it to plants from diverse families across both monocots and eudicots, including the Mediterranean early spider orchid (Ophrys sphegodes, Orchidaceae), Italian ryegrass (Lolium multiflorum, Poaceae), wild cabbage (Brassica incana, Brassicaceae), and South African beetle daisy (Gorteria diffusa, Asteraceae). We successfully extracted HMW DNA of high purity and integrity from all four species, showing that our protocol can be applied to a broad range of angiosperm species. To demonstrate the suitability of our protocol for long-read technology, we sequenced one sample each from *O. sphegodes* and *L. multiflorum* on the ONT® PromethION® platform. Finally, we assessed the use of Circulomics' Short Read Eliminator Kit and its impact on sequenced read length in these species. We show that short fragment removal can improve the average read length and increase the proportion of ultra-long sequenced fragments [*>100 kb* (Prall et al., 2021)], thus improving sequencing efficiency. A schematic overview of our protocol is provided in **Figure 1** and a step-by-step version is available at dx.doi.org/10.17504/protocols.io.5t7g6rn online.

# METHODS

## Species Used for DNA Extraction

We used plants from four families for DNA extraction, specifically *Ophrys sphegodes* Mill. (Orchidaceae), *Lolium multiflorum* Lam. (Poaceae), *Brassica incana* Ten. (Brassicaceae) and *Gorteria diffusa* Thunb. (Asteraceae). Since bee orchids

<footnote n="1">¹https://www.pacb.com/literature/application-brief-whole-genome-sequencing-for-de novo-assembly-best-practices/</footnote>

<footnote n="2">²https://nanoporetech.com/about-us/news/ncm-announcements-include-single-read-accuracy-991-new-chemistry-and-sequencing</footnote>
## DNA extraction

① 100 mg Leaf tissue

<figure><img src="image_6.png" /></figure>

② + 600 µl fresh lysis buffer

<figure><img src="image_7.png" /></figure>

65 °C

<figure><img src="image_8.png" /></figure>

Vortex
00:05

③ + 4 µl 100mg/ml RNAse A

<figure><img src="image_9.png" /></figure>

30:00

Invert 10 times

④ Contaminant precipitation
+200 µl 5M KAc

<figure><img src="image_10.png" /></figure>

Pipette up and down

<figure><img src="image_11.png" /></figure>

Invert 25 times

## Contaminant precipitation and Phenol:Chloroform purification 2x

⑤ + 800 µl Phenol:Chloroform:IAA*

<figure><img src="image_12.png" /></figure>

⑥ Incubate @ 20rpm

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

10:00

⑦ Centrifuge @ 10,000g

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

10:00

⑧ Transfer supernatant into a new 2ml tube

<figure><img src="image_17.png" /></figure>

## Genomic DNA purification with beads

⑨ Vortex bead solution

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

+1 Volume bead solution

⑩ Incubate @ 10rpm

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

30:00

⑪ **3x** Wait until beads bind to the magnet

<figure><img src="image_22.png" /></figure>

⑫ 3x
Invert 25 times

<figure><img src="image_23.png" /></figure>

⑬ Remove supernatant + 50 µl Tris-HCl

<figure><img src="image_24.png" /></figure>

⑭ Resuspend beads

<figure><img src="image_25.png" /></figure>

Vortex

⑭ **3x** Remove supernatant + 1 ml 70% EtOH

<figure><img src="image_26.png" /></figure>

⑮ Incubate @ 10rpm

<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /></figure>

15:00

<figure><img src="image_29.png" /></figure>

37°C

⑮ Wait until beads bind to the magnet

<figure><img src="image_30.png" /></figure>

⑯ Transfer eluted DNA into a new tube

<figure><img src="image_31.png" /></figure>

⑰ O/N to homogenize and relax

<figure><img src="image_32.png" /></figure>

⑱ Quality Control

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /><figcaption>**FIGURE 1** | Schematic overview of the DNA extraction method. The individual steps are described in detail in the see section “Methods” and in the accompanying online protocol.</figcaption></figure>

(Ophrys) hardly produce 2 g total leaf mass per individual, common DNA extraction procedures that need ≥ 1 g starting material make it difficult to reach a final DNA amount suitable

for long-read sequencing without pooling individuals. This, coupled with the large genome size (ca. 1C = 5 Gbp) and high heterozygosity, would make a genome assembly project
unfeasible. Italian ryegrass is one of the most important forage grasses (Gilliland et al., 2007) and represents a major feed source for livestock farming due to its high digestibility and biomass yield (Wilkins, 1991; Frei et al., 2021). The Brassicaceae family includes diverse plant species widely cultivated for oilseed production and vegetable consumption, including oilseed rape, kale, broccoli and cauliflower. The major challenge of extracting DNA from these plants comes from the high level of secondary metabolites in leaves that interfere with sample purity (Zhao et al., 2020). The leaves of *G. diffusa* are even more challenging and contain milky latex (and other unknown secondary metabolites) in high amounts.

## High Molecular Weight DNA Extraction for *Ophrys sphegodes*

Plant material used for this experiment was collected from a greenhouse-grown *O. sphegodes* individual. Young (2 weeks old) leaves were collected, flash frozen in liquid nitrogen (LN₂), and stored at -80°C until DNA extraction. On the day of the experiment, a fresh SDS lysis buffer was prepared as in Mayjonade et al. (2016), supplemented with β-mercaptoethanol (β-ME): 1% polyvinylpyrrolidone 40 (PVP40), 1% sodium metabisulphite, 0.5 M sodium chloride, 100 mM Tris-HCl (pH 8), 50 mM EDTA (pH 8), 2% β-ME, 1.5% sodium dodecyl sulphate (SDS), in ddH₂O to a final volume of 10 ml of stock solution (see Supplementary Table 1). The lysis buffer was incubated at 65°C for 1 h to ensure total dissolution of reagents in the buffer. Meanwhile, 100 mg frozen leaf tissue was ground with mortar and pestle (pre-cooled at -80°C for > 1 h) in LN₂ until a fine powder was obtained (Step 1, Figure 1). Note that it is not advisable to use more than 100 mg starting material, as this decreases DNA purity (as measured by the A₂₆₀/₂₃₀ absorbance ratio using a NanoDrop® spectrophotometer). Also, since the grinding step is crucial to ensure optimal outcomes in terms of final DNA yield and integrity, we provide some tips. First, it is critical to keep the sample submerged in LN₂. If the LN₂ evaporates, plant material will thaw (as indicated by a colour change – it becomes dark green when it thaws) and the DNA will degrade. Second, a fine, flour-like texture of ground plant material is optimal to ensure maximal DNA yield. Hence, the plant material was first crushed in a mortar by gently pounding it with a pestle, until small pieces (<5 × 5 mm) were obtained. Then, plant pieces were rubbed against the mortar with circular movements of the pestle to obtain a final powder with a flour-like texture. To avoid thawing, LN₂ was added every half minute (or when LN₂ had almost evaporated). Depending on the plant material and the pressure applied, the grinding can take up to 30 min. After grinding, the powder was immediately transferred to a sterile 2 ml plastic tube with a chilled metal spatula and mixed with 600 µl of the pre-warmed (65°C) SDS lysis buffer (Step 2 in Figure 1). The sample was vortexed for 3-5 s and incubated on a thermomixer with gentle agitation (400 rpm, 20 min at 55°C) to inactivate DNases and remove polyphenols that could bind DNA. Afterward, 4 µl of 100 mg/ml DNase-free RNase A (Qiagen, Germantown, MD, United States) were added, and the sample was incubated for 10 additional minutes at 55°C (Step 3, Figure 1). To fully

precipitate proteins and polysaccharides that form complexes with SDS (Otzen, 2011), 200 µl of 5 M potassium acetate (KAc) were added, and the solution was mixed by inverting the tube 25 times (Step 4, Figure 1). Next, the sample was purified via a phenol/chloroform extraction as follows. Under a fume hood, 800 µl of a phenol:chloroform:isoamyl alcohol mixture (25:24:1 v/v, pH 8) was added, and the sample was incubated for 10 min at room temperature (RT) with gentle agitation on a tube rotator at 20 rpm (Steps 5 and 6, Figure 1). The sample was then centrifuged for 10 min at 10,000 × g at RT (Step 7, Figure 1). Afterward, the supernatant was transferred into a new 2 ml tube using a 1,000 µl wide-bore pipette tip, to avoid shearing the DNA (Step 8, Figure 1). A second purification step was then undertaken by the addition of 800 µl chloroform:isoamyl alcohol (24:1 v/v), followed by a second incubation (10 min at RT at 20 rpm) and centrifugation (10 min at 10,000 × g at RT). Finally, the supernatant was transferred to a new 2 ml tube and the final volume was recorded (~ 700–800 µl).

## Carboxyl Magnetic Bead Purification

The supernatant recovered from the previous step was purified with carboxylated magnetic beads (Sera-Mag SpeedBeads™ Carboxyl Magnetic Beads, GE Healthcare 65152105050250, Fisher Scientific). The bead stock solution was prepared as in Schalamun et al. (2019) and added to the sample in a 1:1 ratio (Step 9, Figure 1) to remove shorter fragments. Note that a complete resuspension of beads in the stock solution was crucial for optimal DNA yield. The sample tube was incubated for 30 min on a rotator at 10 rpm at RT, spun down for 1 s in a benchtop microcentrifuge, and placed into a magnetic rack until all beads migrated toward the magnet and the solution became clear (Step 10 and 11, Figure 1). This step can take several minutes, as the viscosity of the solution may slow down the beads' migration. Afterward, an ethanol (EtOH) washing step was carried out as follows: 1 ml of 70% EtOH was added to the tube; then, the sample tube was removed from the magnetic rack (without extended incubation), mixed by inverting it 25 times to resuspend the bead pellet, spun down for 1 s, and placed back into the magnetic rack. When the solution became clear, the supernatant was discarded, and the washing step repeated for a total of three times (Step 12–14, Figure 1). We note that these steps can be quite challenging to perform, as beads tend to aggregate, making it difficult to separate and resuspend them properly. In this case, it helps to gently flick the tube to help the beads separate and to avoid prolonged incubation in EtOH. Thereafter, DNA was eluted by addition of 50 µl of 10 mM Tris-HCl pH 8.5 (using commercial buffer EB from Qiagen) preheated to 50°C to the tube, followed by a last incubation at 37°C for 15 min (Steps 14–16, Figure 1). The warm temperature is intended to favour the elution of DNA from the magnetic beads. Finally, to collect the eluted DNA, the tube was placed back into the magnetic rack until the solution became clear (Step 17, Figure 1). This step lasted ca. 30 min, as the long fragments migrate slowly in such a small volume. Note that the slow separation speed is an indication of successful HMW DNA extraction. The eluted DNA was gently pipetted into a new 2 ml tube with a wide-bore 1000 µl pipette tip (Step 18, Figure 1). The sample was very viscous at this point,
indicating highly concentrated HMW DNA. The tube was left on the bench overnight at RT to allow DNA to homogenise and relax (Step 19, **Figure 1**). The next day, the sample was ready for quality control (Step 20, **Figure 1**). A total of six *O. sphegodes* samples were prepared (named OPH_1-6). The total amount of time taken for DNA extraction and clean-up was approximately 2.5 h.

## High Molecular Weight DNA Extraction for Other Plant Species

To evaluate the efficacy and reproducibility of the protocol, we applied it to plant material from four different families: Poaceae (Italian ryegrass, *Lolium multiflorum*), Brassicaceae (wild cabbage, *Brassica incana*) and Asteraceae (beetle daisy, *Gorteria diffusa*). When extracting HMW DNA from Italian ryegrass leaf material, we noticed that the ground powder easily clumped when in contact with SDS lysis buffer (Step 2, **Figure 1**), thereby impacting the efficiency of the lysis step. To avoid powder clumping, it was critical to vortex the sample tube immediately after transferring the powder into the pre-warmed SDS lysis buffer. A total of three DNA samples were extracted (named RAB_1-3).

Wild cabbage leaves contain high amounts of polysaccharides that negatively impact final DNA purity and yield. To reduce the level of polysaccharides in leaf tissues, the plant was placed in the dark for ca. 18 h before harvesting (final sample named BRI_1).

Milky latex and other unknown secondary metabolites are present in high amounts in *G. diffusa* leaves. Carryover of these substances negatively impacts sample purity, resulting in NanoDrop A260/230 ratios outside the optimal range for long-read sequencing. We noticed that *G. diffusa* plants perish quickly without light, making a prolonged dark treatment infeasible. To address these issues, we reduced the amount of starting material and included an additional purification step with magnetic beads. Starting with 70–80 mg young leaves, the protocol was performed as described before until the final DNA elution in 50 µl of 10 mM Tris–HCl buffer (pH 8.5) (Steps 1–18, **Figure 1**). A total of eight tubes were prepared. DNA from two tubes each were pooled and carried through an additional bead purification step. Briefly, 1 volume of magnetic bead solution was added to each sample, followed by incubation at RT for 10–15 min (Steps 9 and 10, **Figure 1**). The tubes were placed into the magnetic rack until the solution became clear, and three washing steps were performed with 1 ml of 70% EtOH without removing the tubes from the rack (Step 11–12, **Figure 1**). After the last wash, the tubes were spun down for 1 s and placed back on the magnetic rack. The beads were resuspended in 50 µl 10 mM Tris–HCl buffer (pH 8.5) (Step 14–15, **Figure 1**) and incubated at 37°C for 15 min. Afterward, the tubes were placed on a magnetic rack for final DNA elution (Step 16–19, **Figure 1**) (samples named GOR_1-4).

## Quality Control Prior to Sequencing

Genomic DNA was evaluated for purity on a NanoDrop® spectrophotometer ND-1000 (Thermo Fisher Scientific, MA, United States). Absorbance at 230, 260, and 280 nm was measured, and A260/280 and A260/230 ratios were assessed to determine DNA purity. Genomic DNA concentration was

measured via NanoDrop® and confirmed with a Qubit® 3.0 fluorometer (Invitrogen, CA, United States) using the dsDNA BR Assay Kit (Thermo Fisher Scientific, MA, United States, Q32850). Note that, if the DNA is pure, the measurements of DNA concentrations from NanoDrop® and Qubit® should be identical.

DNA integrity was assessed on a TapeStation 4200 system (Agilent, CA, United States) with a Genomic DNA ScreenTape Assay. Here, DNA quality is assessed using a DNA Integrity Number (DIN) that ranges from 1 (highly degraded DNA) to 10 (intact DNA). Fragment length was measured on a Femto Pulse v-1.0.0.32 system (Agilent, CA, United States, Cat N° M5330AA) using the Genomic DNA 165 kb Ladder Fast Separation assay with a separation time of 70 min (Agilent, CA, United States, Cat N° FP-1002-0275). For *G. diffusa*, DNA integrity was determined on a Pippin Pulse™ electrophoresis system (Sage Science, MA, United States, Cat N° PP10200) with program 5–80 kb and the Bio-Rad CHEF 5 kb DNA Size Standard.

## Size Selection With the Circulomics Kit

We tested the impact of short DNA fragment depletion on final sequencing results by applying the Circulomics Short Read Eliminator Kit (Circulomics, MD, United States, SS-100-101-01) on one of the two samples selected for sequencing (OPH_3 and RAB_2). The kit was applied to the *L. multiflorum* sample RAB_2 before library preparation to remove small DNA fragments. According to supplier information, the kit uses size-selective precipitation to reduce the amount of DNA fragments below 25 kbp in length<fnref n="3" />. Potentially, the kit can thus significantly enhance average read length during sequencing. The kit was used according to the manufacturer's recommendations (handbook v2.0, 07/2019). Briefly, 60 µl of Buffer SRE were added to the sample tube (60 µl volume), gently mixed and the tube centrifuged at 10,000 × g for 30 min at RT. After supernatant removal, two washing steps were performed with 200 µl of 70% EtOH and a centrifugation at 10000 × g for 2 min at RT. Finally, 100 µl Qiagen Buffer EB were added and the tube was incubated at RT overnight to ensure efficient DNA elution (sample named RAB_2_Circ). No Circulomics kit was applied to sample OPH_3.

## ONT® Library Preparation

We tested the suitability of our protocol by sequencing HMW DNA from samples OPH_3 and RAB_2_Circ. Sequencing library preparation was carried out following the general guidelines from Oxford Nanopore Technologies® for 1D Genomic DNA sequencing, with modifications proposed by New England Biolabs® (NEB) to ensure high data yield production and long-fragment sequencing. For library preparation, the following reagents were used: Ligation Sequencing Kit SQK-LSK109 (Oxford Nanopore Technologies®), NEBNext® Companion Module for Oxford Nanopore Technologies® Ligation Sequencing (NEB, MA, United States, Cat N° E7180S), and AMPure XP beads (Beckman Coulter Inc., CA, United States). A DNA amount of 1.5 µg was collected from samples OPH_3 and RAB_2_Circ (since OPH_3 was highly concentrated, 13 µl of DNA were diluted in 35 µl of 10 mM Tris–HCl buffer to reach

<footnote n="3">³https://www.circulomics.com/store/Short-Read-Eliminator-Kit-p131401036</footnote>
a total volume of 48 µl) and transferred into a 0.2 ml thin-walled PCR tube. DNA fragments were repaired and end-prepped as follows: 3.5 µl NEBNext® FFPE DNA Repair Buffer, 2 µl NEBNext® FFPE DNA Repair Mix, 3.5 µl NEBNext® Ultra™ II End Prep Reaction Buffer, and 3 µl NEBNext® Ultra™ II End Prep Enzyme Mix were added to each tube. After mixing and spinning down, the samples were incubated at 20°C for 30 min, followed by a second incubation at 65°C for 30 min. The original recommendations from NEB were followed, instead of the ONT® guidelines, as preliminary experiments showed better results under NEB supplier conditions: prolonged incubation time allowed recovery of longer fragments. After incubation, the solution from each tube was transferred to a clean 1.5 ml Eppendorf DNA LoBind® tube (Eppendorf AG, Hamburg, Germany) for clean-up. First, a stock solution of AMPure XP Beads was prepared as in Schalamun et al. (2019), and 60 µl were added to each tube. The samples were then incubated on a HulaMixer™ sample mixer (Thermo Fisher Scientific, MA, United States, 15920D) for 20 min at RT, until the solution was homogenised. Bead clean-up was performed with two washing steps on a magnetic rack, each time pipetting off the supernatant and adding 200 µl of freshly prepared 70% EtOH. The pellet was resuspended in 61 µl nuclease-free water (1 µl was then taken out for quantification) and incubated for 10 min at RT on a HulaMixer™. Tubes were placed on a magnetic rack to collect the final eluate. For adapter ligation and clean-up, 60 µl DNA from the previous step was combined with 25 µl Ligation Buffer LNB, 5 µl Adapter Mix AMX, and 10 µl NEBNext® Quick T4 DNA Ligase (240 µl) in a 1.5 ml Eppendorf DNA LoBind tube, and incubated for 20 min at RT. A second AMPure bead clean-up step was carried out by adding 45 µl of bead solution to each tube, followed by incubation on a HulaMixer™ for 20 min at RT. After pipetting off the supernatant on a magnet rack, the beads were washed twice with 250 µl Long Fragment Buffer LFB. Finally, the supernatant was discarded, the pellet resuspended in 25 µl Elution Buffer EB, and incubated for 10 min at RT to collect the final library.

## Long-Read Sequencing

The ONT® PromethION PTC0031 sequencing platform (Alpha-Beta model, OS Ubuntu 16.06, Intel CPU) was used to sequence samples OPH_3 and RAB_2_Circ. For each sample, 0.8–1 µg of the prepared library was loaded onto a R9.4.1 chemistry PromethION flow cell. Sequencing ran until the flow cell buffer was exhausted (63.10 h for OPH_3, 40.14 h for RAB_2_Circ). MinKNOW v-2.1 was used for data acquisition, real-time analysis and sample tracking. Base calling was performed with Guppy v-3.0.4.

## RESULTS

### DNA Purity and Quantity

We compared our DNA extraction protocol with the original version of Mayjonade et al. (2016) and the protocol of Schalamun et al. (2019) which was developed for recalcitrant plant species. We used *O. sphegodes* leaves and compared the results in terms

of purity on a NanoDrop® device. When using the original protocol (Mayjonade et al., 2016), we could recover only small amounts of DNA (1,090 ng DNA total) and failed to recover pure DNA (Figure 2A). The A260/280 = 1.68 indicated protein contamination, and the low A260/230 value indicated presence of polysaccharides, polyphenols or other organic compounds. Using the protocol of Schalamun et al. (2019), we were not able to recover any DNA at all (Figure 2B). With our improved extraction protocol, all purity criteria (A260/280 = 1.80–2.0 and A260/230 = 2.0–2.2) were met in all six *O. sphegodes* samples (Figure 2C, OPH_1-6 in Table 1). Remarkably, we were able to recover an average DNA amount of 5,710 ng per sample (3,950 to 8,450 ng). The NanoDrop® and Qubit® measurements were comparable (NanoDrop®/Qubit® ratio close to 1; Table 1), indicating a clean DNA sample. Since Qubit® measures only double-stranded DNA, a ratio of 1 further indicates that the DNA molecules measured are double stranded and that other molecules with absorption at 260 nm are absent (O'Neill et al., 2011). When the ratio was below 1 (OPH_2 = 0.87, OPH_5 = 0.71), we considered the Qubit® values to be more accurate.

To evaluate the efficacy and reproducibility of our method, we applied it to other plant species. All three Italian ryegrass samples met the quality criteria, and the DNA amount recovered was on average 2,400 ng, with a NanoDrop®/Qubit® ratio of ~ 1 (RAB_1-3, Table 1). When DNA was extracted from *B. incana* leaves, A260/280 was 1.85 and A260/230 was 2.19. Total DNA amount was 7,300 ng, confirmed by a NanoDrop®/Qubit® = 0.97 (Table 1). In *G. diffusa*, we included an additional purification step with magnetic beads after DNA elution. This extra step allowed all samples to meet the purity criteria necessary for sequencing (GOR_1-4, average DNA amount 5,550 ng per sample, A260/280 = ~1.83, A260/230 = ~2.13, Table 1).

### DNA Integrity and Fragment Lengths

Genomic DNA integrity was assessed by determining the degree of fragmentation of the sample. On a TapeStation, all *O. sphegodes* samples had a DIN value of 9.3 or above (with 10 = highly intact DNA, Figure 2D). Fragment lengths of samples OPH_3, RAB_2 and BRI_1 were measured on a Femto Pulse system. Because the fast separation assay was used, it was not possible to distinguish any fragment lengths above 165 kbp, where only a compression band was visible. OPH_3 showed peaks at 110 and 153 kbp, and RAB_2 clearly displayed a peak at around 165 kbp (Figures 3A,B). BRI_1 showed a peak at around 50 kbp and a more heterogeneous distribution of fragment sizes (Figure 3C). DNA of GOR samples showed a clear band above the size of the largest marker of 68 kbp on a Pippin Pulse gel (Figure 3D).

### Size Selection With Circulomics and PromethION Sequencing

We tested our DNA extraction protocol by sequencing samples OPH_3 and RAB_2_Circ. Each library was injected into one flow cell. Flow cell behaviour was comparable, with 2,791 and 2,684 out of 3,000 active channels for OPH_3 and RAB_2_Circ, respectively. Since the sequencing run was longer for OPH_3, the
<figure><img src="image_35.png" /></figure>

<figure><img src="image_36.png" /></figure>

<figure><img src="image_37.png" /></figure>

<figure><img src="image_38.png" /><figcaption>**FIGURE 2** | Comparison of DNA extraction performance for different extraction methods. (A–C) Output from a NanoDrop spectrophotometer. (A) DNA extraction results with the Mayjonade et al. (2016) protocol, without β-ME and phenol:chloroform purification step. (B) DNA extraction with lysis buffer as described by Schalamun et al. (2019). (C) DNA extracted with the protocol described in this study (sample OPH_3). (D) TapeStation results showing the fragment size distribution and the DNA Integrity Number (DIN), for six *O. sphegodes* samples.</figcaption></figure>

<table><caption>**TABLE 1** | Summary of DNA quality measurements for all extracted samples.</caption><thead><tr><th>Species/Sample</th><th>Qubit conc. [ng/µL]</th><th>NanoDrop conc. [ng/µL]</th><th>Output [µg]</th><th>A<sub>260/280</sub> ratio</th><th>A<sub>260/230</sub> ratio</th><th>NanoDrop/Qubit conc. ratio</th><th>Qubit conc. after Circulomics [ng/µL]</th></tr></thead><tbody><tr><td colspan="8"><em>Ophrys sphegodes</em></td></tr><tr><td>OPH_1</td><td>79.0</td><td>114.4</td><td>3.95</td><td>1.83</td><td>1.97</td><td>1.45</td><td>–</td></tr><tr><td>OPH_2</td><td>115.0</td><td>99.8</td><td>5.75</td><td>1.83</td><td>2.23</td><td>0.87</td><td>–</td></tr><tr><td>OPH_3</td><td>169.0</td><td>164.3</td><td>8.45</td><td>1.82</td><td>2.14</td><td>0.97</td><td>–</td></tr><tr><td>OPH_4</td><td>86.3</td><td>85.9</td><td>4.32</td><td>1.81</td><td>2.15</td><td>0.99</td><td>–</td></tr><tr><td>OPH_5</td><td>152.0</td><td>107.9</td><td>7.60</td><td>1.79</td><td>1.97</td><td>0.71</td><td>–</td></tr><tr><td>OPH_6</td><td>83.8</td><td>86.4</td><td>4.19</td><td>1.79</td><td>1.79</td><td>1.03</td><td>–</td></tr><tr><td colspan="8"><em>Lolium multiflorum</em></td></tr><tr><td>RAB_1</td><td>42.0</td><td>40.0</td><td>2.10</td><td>1.95</td><td>2.06</td><td>0.95</td><td>17.0</td></tr><tr><td>RAB_2</td><td>56.0</td><td>53.0</td><td>2.80</td><td>1.92</td><td>2.15</td><td>0.95</td><td>23.0</td></tr><tr><td>RAB_3</td><td>46.0</td><td>42.0</td><td>2.30</td><td>1.85</td><td>2.30</td><td>0.91</td><td>18.0</td></tr><tr><td colspan="8"><em>Brassica incana</em></td></tr><tr><td>BRI_1</td><td>150.0</td><td>146.0</td><td>7.50</td><td>1.85</td><td>2.19</td><td>0.97</td><td>–</td></tr><tr><td colspan="8"><em>Gorteria diffusa</em><sup>1</sup></td></tr><tr><td>GOR_1</td><td>146.0</td><td>182.6</td><td>7.30</td><td>1.84</td><td>2.07</td><td>1.25</td><td>–</td></tr><tr><td>GOR_2</td><td>93.6</td><td>97.0</td><td>4.68</td><td>1.83</td><td>2.18</td><td>1.04</td><td>–</td></tr><tr><td>GOR_3</td><td>102.0</td><td>108.4</td><td>5.10</td><td>1.82</td><td>2.15</td><td>1.06</td><td>–</td></tr><tr><td>GOR_4</td><td>103.0</td><td>112.8</td><td>5.15</td><td>1.82</td><td>2.12</td><td>1.09</td><td>–</td></tr></tbody></table>

final sequencing yield in terms of Gbp data produced was higher for OPH_3 (66 Gbp versus 50 Gbp for RAB_2_Circ; Table 2). The effect of the Circulomics Kit in depleting short fragments

is evident when comparing the fragment size distribution of the samples before (Figures 3A,B) and after (Figure 4) sequencing. When DNA molecules from samples OPH_3 and RAB_2 were

<footnote n="1">¹For *G. diffusa*, 70–80 mg of input material were used, whereas for the other species, it was ~100 mg.</footnote>
<figure><img src="image_39.png" /></figure>

<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

**FIGURE 3** | Fragment length profiles of extracted HMW DNA. (A–C) Femto Pulse profiles showing integrity and size of gDNA. (A) *O. sphegodes* sample OPH_3 with peaks at 110 and 153 kbp. (B) *L. multiflorum* (sample RAB_2) with a peak at 167 kbp. (C) *B. incana* (sample BRI_1) showing a broader fragment size distribution centred at 50 kbp. (D) Gel pictures for *G. diffusa* DNA extracts (GOR_1-4, from a pool of two initial extractions).

evaluated on a Femto Pulse system, both showed the presence of smaller fragments (2–10 kbp). After short read depletion in sample RAB_2_Circ, the amount of small fragments decreased

drastically, and as a result, the number of sequenced small reads was appreciably lower than in OPH_3 (Figure 4A and Table 2). Because short fragments were largely removed in RAB_2_Circ, more than half (51.84%) of the total reads were longer than 10 kbp, and one third (33.08%) were longer than 25 kbp. The read length statistics in OPH_3 were very different, with only 15.75% of the total reads longer than 10 kbp, and a read length N50 roughly half that of RAB_2_Circ (27,196 vs. 51,861 bp, Table 2). In both runs, we were able to recover long to ultra-long reads (> 100 kbp) (Table 2). The longest reads were 464 kbp for RAB_2_Circ, and 1.7 Mbp for OPH_3.

<table><caption>**TABLE 2** | Nanopore (PromethION) sequencing statistics for sequenced samples.</caption><thead><tr><th>Sample<br>Species<br>Treatment</th><th>OPH_3<br><i>O. sphegodes</i><br>No Circulomics</th><th>RAB_2_Circ<br><i>L. multiflorum</i><br>With Circulomics</th></tr></thead><tbody><tr><td>Run duration (h)</td><td>63.10</td><td>40.14</td></tr><tr><td>Active channels</td><td>2,791</td><td>2,684</td></tr><tr><td>Total reads</td><td>9,131,684</td><td>2,181,501</td></tr><tr><td>Sequencing yield (Gbp)</td><td>66.4</td><td>50.4</td></tr><tr><td>Read length N50 (bp)</td><td>27,196</td><td>51,861</td></tr><tr><td>Mean read length (bp)</td><td>7,270.60</td><td>23,122.70</td></tr><tr><td>Median read length (bp)</td><td>2,157.00</td><td>11,119.00</td></tr><tr><td># Reads ≥ 10 kbp</td><td>1,438,821</td><td>1,130,927</td></tr><tr><td># Reads ≥ 25 kbp</td><td>635,963</td><td>721,668</td></tr><tr><td># Reads ≥ 50 kbp</td><td>261,570</td><td>331,789</td></tr><tr><td># Reads ≥ 100 kbp</td><td>55,592</td><td>59,998</td></tr><tr><td># Reads ≥ 200 kbp</td><td>2,011</td><td>1,204</td></tr><tr><td># Reads ≥ 500 kbp</td><td>39</td><td>0</td></tr><tr><td>% Reads ≥ 10 kbp</td><td>15.75%</td><td>51.84%</td></tr><tr><td>% Reads ≥ 25 kbp</td><td>6.96%</td><td>33.08%</td></tr><tr><td>% Reads ≥ 50 kbp</td><td>2.86%</td><td>15.21%</td></tr><tr><td>% Reads ≥ 100 kbp</td><td>0.61%</td><td>2.75%</td></tr><tr><td>% Reads ≥ 200 kbp</td><td>0.02%</td><td>0.06%</td></tr><tr><td>% Reads ≥ 500 kbp</td><td>0.00%</td><td>0.00%</td></tr><tr><td>Longest read (Mbp)</td><td>1.7</td><td>0.464</td></tr></tbody></table>

## DISCUSSION

Long-read sequencing technologies offer a new array of opportunities to study plant genomes in ways that were not feasible before. Our study presents an improved HMW DNA extraction method suitable for a wide variety of plant species and shows how pure high-quality DNA, together with optimised library preparation parameters and size selection, are key for reaching high-throughput ultra-long reads for sequencing projects. Isolated DNA suitable for long-read sequencing has to meet two important criteria: (i) DNA purity and quantity and (ii) DNA integrity.

## DNA Purity and Quantity

Nucleic acid purification from plant tissue can be challenging. In particular, when chemicals such as ethanol or phenol are involved in a DNA extraction procedure, or when proteins or
<figure><img src="image_43.png" /><figcaption>**B** Normalized histogram of log transformed read lengths</figcaption></figure>

<figure><img src="image_44.png" /><figcaption>**FIGURE 4** | Nanopore sequencing read length distributions. Read length distribution for *O. sphegodes* OPH_3 (blue; without Circulomics kit) and *L. multiflorum* RAB_2_Circ (red; after Circulomics kit), showing (**A**) the entire normalised histogram of log-transformed read lengths and (**B**) a zoom-in into the section showing the longest reads (OPH_3 = 1.7 Mbp; RAB_2_Circ = 464 kbp).</figcaption></figure>

other organic components from plant material are present, final nucleic acid purity can be compromised. As a result, the purity of DNA extracts needs to be measured to avoid sequencing samples of inferior quality. The purity of isolated nucleic acids is commonly determined with a spectrophotometer, measuring three UV absorbance (A) values: absorbance at 260, 280, and 230 nm. Nucleic acids absorb UV light at a wavelength with a peak at 260 nm and hence an absorbance spectrum with a 260 nm peak indicates pure DNA. Protein (specifically the aromatic amino acids tryptophan and tyrosine) and phenols absorb UV light at 280 nm. When proteins and phenols contaminate a nucleic acid sample, the absorbance peak at 280 nm decreases by ca. 10-fold (Koetsier and Cantor, 2019). Organic components such as carbohydrates, buffer salts from DNA extraction (like Tris), EtOH, and EDTA strongly absorb at 230 nm. If one or more of these compounds are present in a nucleic acid extract, a lower 230 nm absorbance is detected. Overall, the ratios of the absorbance values A260/280 and A260/230 allow estimation of

the purity of a nucleic acid sample. Pure dsDNA has absorbance ratios of A260/280 = 1.8 – 2.0 and A260/230 = 2.0 – 2.2; lower values indicate the presence of contaminants, while higher A260/280 values can be indicative of RNA contamination (pure RNA has an A260/280 = 2.0 – 2.2) (Glasel, 1995).

To obtain pure genomic DNA that meets the aforementioned characteristics, we first tried the protocols of both Mayjonade et al. (2016) and Schalamun et al. (2019) on *Ophrys* leaves, and then produced an improved and more robust method. The DNA extraction protocol from Mayjonade et al. (2016) is divided into three main steps: cell membrane disruption with SDS lysis buffer, contaminant precipitation with 5 M KAc, and final purification of gDNA with Sera-Mag SpeedBeads magnetic beads. Schalamun et al. (2019) introduced changes in lysis buffer composition, and used different incubation times at different temperatures. When our *Ophrys* sample was extracted with the protocol of Mayjonade et al. (2016), DNA concentration was low, and the UV absorbance spectrum showed severe contamination. A260/280 and A260/230
are unreliable at DNA concentrations < 20 ng/µl (Koetsier and Cantor, 2019), and in general such weakly concentrated samples are not suitable for long-read sequencing (as per manufacturers' protocols; **Figure 2A**).

To produce a DNA extract of high concentration and purity, we modified the lysis buffer by increasing SDS concentration to 1.5% and adding 2% β-ME. A higher SDS concentration ensures a more effective rupture of cell walls after a first mechanical breakage by grinding under LN₂, while β-ME is a reducing agent that denatures proteins by breaking the disulphide bonds between cysteine residues. Together, those two reagents increased the recovery of DNA from cells. Moreover, it has been reported that high levels of β-ME successfully remove polyphenols (Khanuja et al., 1999) and other organic compounds, such as tannins, from plant tissue (John, 1992; Moreira et al., 2011). A second modification was introduced after protein precipitation with 5M KAc. We added a phenol:chloroform:isoamyl alcohol (25:24:1 v/v) purification to remove other carryover contaminants. The separation of the solution in two phases allowed us to extract the isolated DNA in the aqueous phase, while leaving carryover contaminants from the extraction in the organic phase. Residual phenol was then removed during the bead purification step. These changes proved to be effective in delivering highly purified DNA in all our study species (see absorbance ratios in **Table 1**). All samples met the A₂₆₀/2₈₀ and A₂₆₀/2₃₀ criteria suitable for long-read sequencing, regardless of the plant species used. Remarkably, final DNA amounts were on average ~ 5,000 ng per sample (from Qubit® measurement; ~100 ng/µl) and concentrations ranged from 42 to 169 ng/µl, depending on the species (**Table 1**). Hence, another positive effect of the improvements was the relatively high DNA quantity we were able to recover. ONT® recommends using at least 1 µg DNA per library preparation. Thus, the final DNA amount per sample was enough for several parallel library preparations. Since PacBio® recommends using 5 µg DNA for CLR and CCS sequencing (less if the plant genome of interest is 500 Mbp - 1 Gbp in size), sample pooling would have been necessary for sequencing with this platform.

The same DNA extraction procedure may result in different yields in different plant species, as different plants have different tissue characteristics. If the extracted DNA appears to be suboptimal in terms of purity, we suggest reducing the amount of starting material to 70–80 mg per sample, as we did with *G. diffusa*. In this way, the level of contaminants from plant material that may interfere with SDS during cell lysis is reduced. An additional purification step with magnetic beads can clean the DNA further. As long fragments bind to the magnetic beads and are not washed away, DNA yield is hardly impacted, while residual contaminants are removed. It may often be advisable to use softer and younger rather than older and tougher leaves. In this regard, plant secondary metabolites that negatively impact DNA extraction have been shown to accumulate over time as leaves age (Boege and Marquis, 2005; Moreira et al., 2011), although in other cases higher concentrations of phenolics have been found in young leaves (Barton et al., 2019). It is important to stress that the grinding step is crucial to recover a high DNA yield and that insufficient grinding can reduce final DNA yield.

Finally, heating the elution buffer to 37°C before usage can help to increase the elution efficiency in the final step of the protocol.

## Integrity of DNA Molecules and Effect on Sequencing

In long-read sequencing, the fraction of long reads (10 kbp or longer) depends strongly on the integrity (i.e., degree of fragmentation) of the DNA molecules used for library preparation. Therefore, one of the goals of an HMW DNA extraction method suitable for long-read sequencing is to preserve and maximise long and ultra-long DNA molecules. One critical step is to prevent DNA damage by thawing. When harvesting, plant material should immediately be flash-frozen, stored at -80°C, and transferred to a cool mortar only before immediate use. During grinding, LN₂ should not evaporate completely, as this can thaw plant material and cause DNA damage. Final extracted DNA can be stored at 4°C if actively used. For long-term storage, gDNA is best stored at -20°C and only thawed when used, as repeated cycles of freezing and thawing can degrade DNA molecules. Another source of DNA degradation is nuclease activity. During mechanical tissue rupture and chemical cell lysis, enzymes such as DNases are released along with DNA. For this reason, β-ME was added to the lysis buffer, which disrupts disulphide bonds, thus inactivating DNases (Price et al., 1969). For best results, we recommend to quickly homogenise the ground tissue powder immediately after adding lysis buffer by vortexing the sample tube. To prevent DNA fragmentation, it is also advisable to use wide bore pipette tips (or P1000 tips with cut ends). Such tips have a wider opening for aspirating and dispensing viscous solutions, and thus they can decrease DNA shearing due to pipetting. For the same reason, we advise carefully pipetting the DNA with slow, gentle movements of the pipette plunger. Since elution time plays a role in long fragment recovery, we prolonged the elution time during bead purification and library preparation. This helped recovering a higher fraction of long fragments and an overall higher DNA amount. As shown in **Figure 4B** and **Table 2**, these improvements allowed us to obtain long to ultra-long reads in both sequencing runs.

Finally, we recommend to remove short DNA molecules (<10 kbp) before sequencing. The presence of these “short” fragments does not result in poorer quality sequence, but it does affect a sequencing run’s fragment size distribution (cf. **Figure 4**). Short fragments compete with longer ones for pore occupancy in the flow cell, decreasing the throughput of long reads per sequencing run. Although shorter reads can be discarded bioinformatically, this approach may not be optimal in terms of cost efficiency when the aim of long-read sequencing is to generate data for genome assembly.

## CONCLUSION

Long-read sequencing approaches are starting to dominate as the primary tools for genome sequencing projects. Thanks to continuous improvements in sequencing chemistry and technology, long reads are used in a plethora of applications,
from *de novo* genome assembly to methylome analysis, to field applications. Here, we provide a robust method to extract purified high molecular weight DNA for long-read sequencing. Our DNA extraction method successfully worked on diverse monocot and eudicot plant species, indicating that the method is effective in a wide variety of plant families. For particularly challenging species, we show that adding an extra purification step allows the user to achieve the purity required for long-read sequencing, while maintaining a high integrity of the DNA molecules. A step-by-step version of the protocol is available online at dx.doi.org/10.17504/protocols.io.5t7g6rn.

## DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

## AUTHOR CONTRIBUTIONS

AR and PS: conceptualisation and writing—original draft. AR, BM, GP, and DF: methodology and investigation. RK, LF, and DF: validation. AR: data curation and visualisation. PS and UG:

## REFERENCES

- Barton, K. E., Edwards, K. F., and Koricheva, J. (2019). Shifts in woody plant defence syndromes during leaf development. *Funct. Ecol.* 33, 2095–2104. doi: 10.1111/1365-2435.13435

- Belser, C., Istace, B., Denis, E., Dubarry, M., Baurens, F. C., Falentin, C., et al. (2018). Chromosome-scale assemblies of plant genomes using nanopore long reads and optical maps. *Nat. Plants* 4, 879–887. doi: 10.1038/s41477-018-0289-4

- Boege, K., and Marquis, R. J. (2005). Facing herbivory as you grow up: the ontogeny of resistance in plants. *Trends Ecol. Evol.* 20, 441–448. doi: 10.1016/j.tree.2005.05.001

- Chawla, H. S., Lee, H. T., Gabur, I., Vollrath, P., Tamilselvan-Nattar-Amutha, S., Obermeier, C., et al. (2021). Long-read sequencing reveals widespread intragenic structural variants in a recent allopolyploid crop plant. *Plant Biotechnol. J.* 19, 240–250. doi: 10.1111/pbi.13456

- Exposito-Alonso, M., Drost, H. G., Burbano, H. A., and Weigel, D. (2020). The earth biogenome project: opportunities and challenges for plant genomics and conservation. *Plant J.* 102, 222–229. doi: 10.1111/tpj.14631

- Flusberg, B. A., Webster, D., Lee, J., Travers, K., Olivares, E., Clark, A., et al. (2010). Direct detection of DNA methylation during single-molecule, real-time sequencing. *Nat. Methods* 7, 461–465. doi: 10.1038/nmeth.1459.

- Formenti, G., Theissinger, K., Fernandes, C., Bista, I., Bleidorn, C., Angelica, C., et al. (2022). The era of reference genomes in conservation genomics. *Trends Ecol. Evol.* 37, 197–202. doi: 10.1016/j.tree.2021.11.008

- Frei, D., Veekman, E., Grogg, D., Stoffel-Studer, I., Morishima, A., Shimizu-Inatsugi, R., et al. (2021). Ultralong Oxford Nanopore reads enable the development of a reference-grade perennial ryegrass genome assembly. *Genome Biol. Evol.* 13, 1–6. doi: 10.1093/gbe/evab159

- Gerstein, A. S. (2001). *Molecular Biology Problem Solver: a Laboratory Guide*. Chichester: Wiley-Liss, Inc. doi: 10.1002/0471223905

- Gilliland, T. J., Johnston, J., and Connolly, C. (2007). A review of forage grass and clover seed use in Northern Ireland, UK between 1980 and 2004. *Grass Forage Sci.* 62, 239–254. doi: 10.1111/j.1365-2494.2007.00588.x

- Glasel, J. (1995). Validity of nucleic acid purities monitored by 260/280 absorbance ratios. *BioTechniques* 18, 62–63.

supervision. DC, BS, JF, UG, and PS: resources. BS, JF, UG, and PS: funding acquisition. All authors: writing, review, and editing.

## FUNDING

This work was supported by the University of Zurich, a Ph.D. project funded by the University Research Priority Programme “Evolution in Action” (to UG and PS) and the Swiss National Science Foundation (mobility grant P400PB_194374 to RK).

## ACKNOWLEDGMENTS

We are grateful to Salvatore Cozzolino for providing plant material and to Kelsey Byers and Sergio Ramos for comments on the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.883897/full#supplementary-material

- Goerner-Potvin, P., and Bourque, G. (2018). Computational tools to unmask transposable elements. *Nat. Rev. Genet.* 19, 688–704. doi: 10.1038/s41576-018-0050-x

- Hasing, T., Tang, H., Brym, M., Khazi, F., Huang, T., and Chambers, A. H. (2020). A phased *Vanilla planifolia* genome enables genetic improvement of flavour and production. *Nat. Food* 1, 811–819. doi: 10.1038/s43016-020-00197-2

- Heller, D., and Vingron, M. (2019). SVIM: structural variant identification using mapped long reads. *Bioinformatics* 35, 2907–2915. doi: 10.1093/bioinformatics/btz041

- Hu, L., Xu, Z., Wang, M., Fan, R., Yuan, D., Wu, B., et al. (2019). The chromosome-scale reference genome of black pepper provides insight into piperine biosynthesis. *Nat. Commun.* 10:4702. doi: 10.1038/s41467-019-12607-6

- Jiao, W. B., and Schneeberger, K. (2017). The impact of third generation genomic technologies on plant genome assembly. *Curr. Opin. Plant Biol.* 36, 64–70. doi: 10.1016/j.pbi.2017.02.002

- John, M. E. (1992). An efficient method for isolation of RNA and DNA from plants containing polyphenolics. *Nucleic Acids Res.* 20:2381. doi: 10.1093/nar/20.9.2381

- Katterman, F. R. H., and Shattuck, V. I. (1983). An effective method of dna isolation from the mature leaves of gossypium species that contain large amounts of phenolic terpenoids and tannins. *Prep. Biochem.* 13, 347–359. doi: 10.1080/00327488308068177

- Kersey, P. J. (2019). Plant genome sequences: past, present, future. *Curr. Opin. Plant Biol.* 48, 1–8. doi: 10.1016/j.pbi.2018.11.001

- Khanuja, S., Shasany, A., Darokar, M. P., and Kumar, S. (1999). Rapid isolation of DNA from dry and fresh samples of plants producing large amounts of secondary metabolites and essential oils. *Plant Mol. Biol. Rep.* 17, 74–1–74. doi: 10.1023/A:1007528101452

- Koetsier, G., and Cantor, E. (2019). *A Practical Guide to Analyzing Nucleic Acid Concentration and Purity with Microvolume Spectrophotometers*. New England Biolabs Tech. Note 7/2019. Available online at: https://prd-sccd01-international.neb.com/-/media/nebus/files/application-notes/technote_mvs_analysis_of_nucleic_acid_concentration_and_purity.pdf?rev=c24cea043416420d84fb6bf7b554dbbb (accessed February 24, 2022).
Mayjonade, B., Muños, S., Donnadieu, C., Marande, W., Langlade, N., Callot, C., et al. (2016). Extraction of high-molecular-weight genomic DNA for long-read sequencing of single molecules. *Biotechniques* 61, 203–205. doi: 10.2144/000114460

Moreira, P. A., Oliveira, D. A., Federal, U., Gerais, D. M., Horizonte, B., Estadual, U., et al. (2011). Leaf age affects the quality of DNA extracted from *Dimorphandra Mollis* (*Fabaceae*), a tropical tree species from the Cerrado region of Brazil. *Genet. Mol. Res.* 10, 353–358. doi: 10.4238/vol10-1gmr 1030

Niu, S., Li, J., Bo, W., Yang, W., Zuccolo, A., Giacomello, S., et al. (2022). The Chinese pine genome and methylome unveil key features of conifer evolution. *Cell* 185, 204.e-217.e. doi: 10.1016/j.cell.2021.12.006

O'Neill, M., McPartlin, J., Arthure, K., Riedel, S., and McMillan, N. D. (2011). Comparison of the TLDA with the Nanodrop and the reference Qubit system. *J. Phys. Conf. Ser.* 307:012047. doi: 10.1088/1742-6596/307/1/012047

Otzen, D. (2011). Protein-surfactant interactions: a tale of many states. *Biochim. Biophys. Acta* 1814, 562-591. doi: 10.1016/j.BBAPAP.2011.03.003

Pandey, R. N., Adams, R. P., and Flournoy, L. E. (1996). Inhibition of random amplified polymorphic DNAs (RAPDs) by plant polysaccharides. *Plant Mol. Biol. Rep.* 14, 17-22. doi: 10.1007/BF02671898

Prall, T. M., Neumann, E. K., Karl, J. A., Shortreed, C. G., Baker, D. A., Bussan, H. E., et al. (2021). Consistent ultra-long DNA sequencing with automated slow pipetting. *BMC Genomics* 22:182. doi: 10.1186/S12864-021-07500-W/ FIGURES/4

Price, P. A., Stein, W. H., and Moore, S. (1969). Effect of divalent cations on the reduction and re-formation of the disulfide bonds of deoxyribonuclease. *J. Biol. Chem.* 244, 929-932. doi: 10.1016/s0021-9258(18)91875-2

Pu, X., Li, Z., Tian, Y., Gao, R., Hao, L., Hu, Y., et al. (2020). The honeysuckle genome provides insight into the molecular mechanism of carotenoid metabolism underlying dynamic flower coloration. *New Phytol.* 227, 930-943. doi: 10.1111/nph.16552

Saxena, R. K., Edwards, D., and Varshney, R. K. (2014). Structural variations in plant genomes. *Brief. Funct. Genomics* 13, 296–307. doi: 10.1093/BFGP/EL U016

Schalamun, M., Nagar, R., Kainer, D., Beavan, E., Eccles, D., Rathjen, J. P., et al. (2019). Harnessing the MinION: an example of how to establish long-read sequencing in a laboratory using challenging plant tissue from *Eucalyptus Pauciflora*. *Mol. Ecol. Resour.* 19, 77–89.

Sedlazeck, F. J., Rescheneder, P., Smolka, M., Fang, H., Nattestad, M., Von Haeseler, A., et al. (2018). Accurate detection of complex structural variations using single-molecule sequencing. *Nat. Methods* 15, 461–468. doi: 10.1038/s41592- 018-0001-7

Simpson, J. T., Workman, R. E., Zuzarte, P. C., David, M., Dursi, L. J., and Timp, W. (2017). Detecting DNA cytosine methylation using nanopore sequencing. *Nat. Methods* 14, 407–410. doi: 10.1038/nmeth.4184

Varma, A., Padh, H., and Shrivastava, N. (2007). Plant genomic DNA isolation: an art or a science. *Biotechnol. J.* 2, 386-392. doi: 10.1002/biot.200600195

Wilkins, P. W. (1991). Breeding perennial ryegrass for agriculture. *Euphytica* 52, 201-214.

Zerpa-Catanho, D., Zhang, X., Song, J., Hernandez, A. G., and Ming, R. (2021). Ultra-long DNA molecule isolation from plant nuclei for ultra-long read genome sequencing. *STAR Protoc.* 2:100343. doi: 10.1016/j.xpro.2021.100343

Zhang, B., Gao, Y., Zhang, L., and Zhou, Y. (2021). The plant cell wall: biosynthesis, construction, and functions. *J. Integr. Plant Biol.* 63, 251-272. doi: 10.1111/JIPB. 13055

Zhang, G. Q., Liu, K. W., Li, Z., Lohaus, R., Hsiao, Y. Y., Niu, S. C., et al. (2017). The Apostasia genome and the evolution of orchids. *Nature* 549, 379-383. doi: 10.1038/nature23897

Zhang, J., Chen, L. L., Xing, F., Kudrna, D. A., Yao, W., Copetti, D., et al. (2016). Extensive sequence divergence between the reference genomes of two elite indica rice varieties zhenshan 97 and minghui 63. *Proc. Natl. Acad. Sci. U.S.A.* 113, E5163-E5171. doi: 10.1073/pnas.1611012113

Zhao, Y., Yue, Z., Zhong, X., Lei, J., Tao, P., and Li, B. (2020). Distribution of primary and secondary metabolites among the leaf layers of headed cabbage (*Brassica oleracea* var. *capitata*). *Food Chem.* 312:126028. doi: 10.1016/j.foodchem.2019.126028

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Russo, Mayjonade, Frei, Potente, Kellenberger, Frachon, Copetti, Studer, Frey, Grossniklaus and Schlüter. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.