# De novo assembly and annotation of three Leptosphaeria genomes using Oxford Nanopore MinION sequencing

Fabien Dutreux, Corinne da Silva, Léo d'Agata, Arnaud Couloux, Elise Gay, Benjamin Istace, Nicolas Lapalu, Arnaud Lemainque, Juliette Linglin, Benjamin Noel, et al.

► To cite this version:

Fabien Dutreux, Corinne da Silva, Léo d'Agata, Arnaud Couloux, Elise Gay, et al. De novo assembly and annotation of three Leptosphaeria genomes using Oxford Nanopore MinION sequencing. Scientific Data, 2018, 5 (180235), 11 p. 10.1038/sdata.2018.235. hal-02621276

HAL Id: hal-02621276

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

Submitted on 26 May 2020

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC DATA
110110
0111101
11011110
011101101

OPEN

# Data Descriptor: *De novo* assembly and annotation of three *Leptosphaeria* genomes using Oxford Nanopore MinION sequencing

Received: 5 April 2018

Accepted: 13 July 2018

Published: 6 November 2018

Fabien Dutreux<fnref n="1" />,², Corinne Da Silva<fnref n="1" />, Léo d'Agata<fnref n="1" />, Arnaud Couloux<fnref n="1" />, Elise J. Gay²,
Benjamin Istace<fnref n="1" />, Nicolas Lapalu², Arnaud Lemainque<fnref n="1" />, Juliette Linglin², Benjamin Noel<fnref n="1" />,
Patrick Wincker³, Corinne Cruaud<fnref n="1" />, Thierry Rouxel², Marie-Hélène Balesdent² &
Jean-Marc Aury<fnref n="1" />

*Leptosphaeria maculans* and *Leptosphaeria biglobosa* are ascomycete phytopathogens of *Brassica napus* (oilseed rape, canola). Here we report the complete sequence of three *Leptosphaeria* genomes (*L. maculans* JN3, *L. maculans* Nz-T4 and *L. biglobosa* G12-14). Nz-T4 and G12-14 genome assemblies were generated *de novo* and the reference JN3 genome assembly was improved using Oxford Nanopore MinION reads. The new assembly of *L. biglobosa* showed the existence of AT rich regions and pointed to a genome compartmentalization previously unsuspected following Illumina sequencing. Moreover nanopore sequencing allowed us to generate a chromosome-level assembly for the *L. maculans* reference isolate, JN3. The genome annotation was supported by integrating conserved proteins and RNA sequencing from Leptosphaeria-infected samples. The newly produced high-quality assemblies and annotations of those three *Leptosphaeria* genomes will allow further studies, notably focused on the tripartite interaction between *L. maculans*, *L. biglobosa* and oilseed rape. The discovery of as yet unknown effectors will notably allow progress in *B. napus* breeding towards *L. maculans* resistance.

<footnote n="1">¹Genoscope, Institut de Biologie François-Jacob, Commissariat à l'Energie Atomique (CEA), Université Paris-Saclay, F-91057 Evry, France. ²UMR BIOGER, INRA, AgroParisTech, Université Paris-Saclay, Avenue Lucien Brétignières, BP 01, F-78850 Thiverval-Grignon, France. ³Génomique Métabolique, Genoscope, Institut de Biologie François Jacob, Commissariat à l'Energie Atomique (CEA), CNRS, Université d'Evry, Université Paris-Saclay, 91057 Evry, France. Correspondence and requests for materials should be addressed to M.-H.B. (email: marie-helene.balesdent@inra.fr) or J.-M.A. (email: jmaury@genoscope.cns.fr)</footnote>
<table><tr><td>Design Type(s)</td><td>sequence assembly objective ● transcription profiling by high throughput sequencing design ● sequence annotation objective</td></tr><tr><td>Measurement Type(s)</td><td>genome assembly ● sequence annotation</td></tr><tr><td>Technology Type(s)</td><td>whole genome sequencing ● RNA sequencing</td></tr><tr><td>Factor Type(s)</td><td>Organism Strain ● Biospecimen Type ● cultivar ● culture medium ● sampling time</td></tr><tr><td>Sample Characteristic(s)</td><td>Leptosphaeria maculans 'brassicae' group ● Leptosphaeria biglobosa 'brassicae' group ● Leptosphaeria maculans species complex</td></tr></table>

## Background & Summary

Fungi have an important ecological and economic role. Many of them are pathogens of animals or plants, especially *Leptosphaeria maculans* which is an ascomycete phytopathogen of *Brassica napus* (oilseed rape, canola). Two genomes of *Leptosphaeria* were already available, the first genome sequence of this fungus (JN3) was obtained in the 2000's using Sanger sequencing¹ and composed of 41 scaffolds totaling 43.76 Mb (excluding ambiguous bases). In-depth analysis of this genome showed an unusual compartmentalization between AT-rich and GC-equilibrated blocks termed isochores. The location of effector genes in the plastic, AT-rich compartments of the genome was instrumental in establishing the “two-speed fungal genome” paradigm¹. According to this paradigm, phytopathogen genomes encompass a plastic compartment enriched in genes involved in niche adaptation. Similar genome organisation was later on found to be a commonality in numerous filamentous phytopathogens, but its identification was slowed down by the use of short-read sequencing that tended to misrepresent the plastic, Transposable Element (TE)-rich genome in the assemblies. Thus the genome sequencing of *Leptosphaeria biglobosa* (B3.5), a species related to *L. maculans* resulted in a poor-quality assembly with complete absence of a putative dispensable genome². This assembly of B3.5 was composed of 606 scaffolds totaling 29.44 Mb (excluding ambiguous bases).

Generally, when using short-read sequencing strategies, TE complicate the assembly step and it results in an under estimation of the repetitive content by collapsing similar copies of repetitions. In this study, we generated genome assembly and annotation of three *Leptosphaeria* isolates. To overcome the limitation of short reads, we generated long reads using the Oxford Nanopore technology (Data Citation 1–3). First, we resequenced the JN3 isolate with the objective of improving the existing assembly in terms of both contiguity and gaps closure (Fig. 1). And then we *de novo* sequenced and assembled two isolates *L. maculans* Nz-T4 and *L. biglobosa* G12-14 (Fig. 1). For each isolate, we generated long and short reads using respectively the MinION device from Oxford Nanopore and the Illumina sequencing technology (Tables 1, 2, 3). Reconciling manually the data from nanopore sequencing, optical and genetic maps allowed the generation of a chromosome-scale assembly for JN3 composed of 19 chromosomes plus the conditionally dispensable chromosome³, with only four missing telomeric extremities. The three genome assemblies were composed of 33, 288 and 156 scaffolds totaling 45.99 Mb, 43.42 Mb and 34.95 Mb for JN3, Nz-T4 and G12-14, respectively (Table 4). Furthermore, the gene prediction was supported by integrating conserved proteins and RNA sequencing from Leptosphaeria-infected samples (Fig. 2 and Data Citation 4). It resulted in 13,047, 14,026 and 12,678 predicted protein-coding genes for JN3, Nz-T4 and G12-14, respectively (Table 5 and Fig. 3).

## Methods

### Biological material

Two *L. maculans brassicae* isolates, JN3 (v23.1.3) and Nz-T4 were sequenced here. JN3 is the reference sequenced isolate¹ of European origin. Nz-T4 has been used in previous genetic analyses⁴ and originates from New Zealand. The *L. biglobosa brassicae* isolate G12-14, also sequenced here, was isolated in 2014 in France (Thiverval-Grignon) as a single-ascospore isolate from naturally infected oilseed rape stubble. G12-14 was chosen as a new *L. biglobosa* reference isolate since the previously sequenced isolate, B3.5², showed a very low level of aggressiveness when used in cotyledon inoculation tests. In contrast, G12-14 rapidly produced on oilseed rape plants typical and intense symptoms of *L. biglobosa brassicae*. In addition it was considered as a better representative of the current *L. biglobosa* populations infecting oilseed rape in France than B3.5.

For RNA sequencing, additional isolates were used, including *L. maculans* JN2 (v23.1.2), a sister isolate of JN3 with a high level of aggressiveness and the *L. biglobosa brassicae* reference isolate, B3.5². Plant materials infected with *L. maculans* and/or *L. biglobosa* were obtained under various conditions (in the field or under controlled conditions) and with three different *B. napus* varieties. Naturally infected plant material was obtained in Thiverval-Grignon experimental fields as either infected stem base at the
<figure><img src="image_2.png" /><figcaption>**Figure 1.** General description of the assembly workflow. Amelioration of the existing JN3 assembly and *de novo* assembly of Nz-T4 and G12-14 isolates.</figcaption></figure>

<table><caption>**Table 1.** Metrics of raw nanopore datasets.</caption><thead><tr><th></th><th></th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans Nz-T4</th><th>Leptosphaeria biglobosa G12-14</th></tr></thead><tbody><tr><td rowspan="7">All</td><td>Number of reads</td><td>1,736,075 (23 flowcells)</td><td>883,625 (30 flowcells)</td><td>931,376 (19 flowcells)</td></tr><tr><td>Cumulative size</td><td>6,044,436,091</td><td>4,813,632,244</td><td>3,362,956,224</td></tr><tr><td>Estimated coverage</td><td>134X</td><td>107X</td><td>96X</td></tr><tr><td>Average Size (bp)</td><td>3,482</td><td>5,448</td><td>3,611</td></tr><tr><td>Longest read (bp)</td><td>1,612,597</td><td>1,769,119</td><td>949,732</td></tr><tr><td>N50 (bp)</td><td>6,441</td><td>7,197</td><td>7,115</td></tr><tr><td># reads > 10Kb</td><td>72,482</td><td>76,176</td><td>61,856</td></tr><tr><td rowspan="7">1D</td><td>Number of reads</td><td>1,278,253</td><td>443,063</td><td>482,374</td></tr><tr><td>Cumulative size</td><td>4,402,377,873</td><td>2,272,788,773</td><td>1,693,737,429</td></tr><tr><td>Estimated coverage</td><td>98X</td><td>51X</td><td>48X</td></tr><tr><td>Average Size (bp)</td><td>3,444</td><td>5,130</td><td>3,511</td></tr><tr><td>Longest read (bp)</td><td>1,612,597</td><td>1,769,119</td><td>949,732</td></tr><tr><td>N50 (bp)</td><td>6,431</td><td>7,497</td><td>7,250</td></tr><tr><td># reads > 10Kb</td><td>52,999</td><td>42,132</td><td>33,148</td></tr><tr><td rowspan="7">2D</td><td>Number of reads</td><td>457,822</td><td>440,562</td><td>449,002</td></tr><tr><td>Cumulative size</td><td>1,642,058,218</td><td>2,540,843,471</td><td>1,669,218,795</td></tr><tr><td>Estimated coverage</td><td>36X</td><td>56X</td><td>48X</td></tr><tr><td>Average Size (bp)</td><td>3,587</td><td>5,767</td><td>3,718</td></tr><tr><td>Longest read (bp)</td><td>237,404</td><td>301,948</td><td>84,061</td></tr><tr><td>N50 (bp)</td><td>6,469</td><td>7,008</td><td>7,005</td></tr><tr><td># reads > 10Kb</td><td>19,483</td><td>34,044</td><td>28,708</td></tr></tbody></table>
<table><caption>Table 2. Metrics of illumina datasets.</caption><thead><tr><th></th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans Nz-T4</th><th>Leptosphaeria biglobosa G12-14</th></tr></thead><tbody><tr><td>Number of reads</td><td>36,421,608</td><td>19,106,872</td><td>22,108,636</td></tr><tr><td>Cumulative size</td><td>8,089,956,155</td><td>4,399,409,082</td><td>4,198,089,452</td></tr><tr><td>Estimated coverage</td><td>180 X</td><td>98 X</td><td>120 X</td></tr><tr><td>Size (bp)</td><td>2×251 # HiSeq2500</td><td>2×251 # MiSeq</td><td>2×251 # MiSeq</td></tr></tbody></table>

<table><caption>Table 3. Metrics of nanopore datasets used for genome assemblies.</caption><thead><tr><th rowspan="2">status</th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans Nz-T4</th><th>Leptosphaeria biglobosa G12-14</th></tr><tr><th>complete dataset</th><th>2D dataset</th><th>2D dataset</th></tr></thead><tbody><tr><td>Number of reads</td><td>1,736,075</td><td>440,562</td><td>449,002</td></tr><tr><td>Cumulative size</td><td>6,044,436,091</td><td>2,540,843,471</td><td>1,669,218,795</td></tr><tr><td>Estimated coverage</td><td>134X</td><td>56X</td><td>48X</td></tr><tr><td>Average Size (bp)</td><td>3,482</td><td>5,767</td><td>3,718</td></tr><tr><td>Longest read (bp)</td><td>1,612,597</td><td>301,948</td><td>84,061</td></tr><tr><td>N50 (bp)</td><td>6,441</td><td>7,008</td><td>7,005</td></tr><tr><td># reads > 10Kb</td><td>72,482</td><td>34,044</td><td>28,708</td></tr></tbody></table>

<table><caption>Table 4. Metrics of the existing and new assemblies.</caption><thead><tr><th rowspan="2">Reference</th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans Nz-T4</th><th>Leptosphaeria biglobosa B3.5</th><th>Leptosphaeria biglobosa G12-14</th></tr><tr><th>Rouxel et al.<sup>1</sup></th><th>This study</th><th>This study</th><th>Grandaubert et al.<sup>2</sup></th><th>This study</th></tr></thead><tbody><tr><td># sequences</td><td>41</td><td>33</td><td>288</td><td>606</td><td>156</td></tr><tr><td>Cumulative size</td><td>44,892,605</td><td>45,986,477</td><td>43,426,637</td><td>31,788,051</td><td>34,950,111</td></tr><tr><td>N50</td><td>1,769,547</td><td>2,437,616</td><td>383,462</td><td>779,070</td><td>462,395</td></tr><tr><td>N90</td><td>1,020,521</td><td>1,391,278</td><td>64,152</td><td>125,916</td><td>123,107</td></tr><tr><td>L50</td><td>10</td><td>8</td><td>38</td><td>14</td><td>22</td></tr><tr><td>L90</td><td>22</td><td>18</td><td>150</td><td>49</td><td>77</td></tr><tr><td># of N's</td><td>1,128,152 (2.51%)</td><td>489,445 (1.06%)</td><td>0 (0%)</td><td>2,343,201 (7.37%)</td><td>0 (0%)</td></tr><tr><td>% GC</td><td>45.24%</td><td>45.22%</td><td>45.69%</td><td>51.39%</td><td>49.13%</td></tr></tbody></table>

end of the growing season, sampled one week before harvest (variety Darmor-Bzh) or as whole stem residues sampled two weeks after harvest (variety Alpaga). Under controlled conditions, two varieties were used, Darmor-Bzh displaying a high level of field resistance to L. maculans, and Bristol, with a high level of susceptibility. Cotyledons and stems were inoculated as previously described⁵ and inoculated material was sampled at different time points after inoculation. For petiole inoculations, the petiole was cut under the leaf blade. Ten µL of inoculum (10⁷ pycnidiospores mL⁻¹) were placed over the wound. Inoculated plants were incubated in darkness for 48 h at room temperature under saturated humidity, and then transferred to a growth chamber at 16 °C (night) and 24 °C (day), with a 16 h photoperiod. Infected petioles were sampled after 7 and 14 days post infection (dpi).

RNA/DNA extraction

Fresh infected cotyledons or freeze-dried cultures of fungal mycelia were ground using a Retsch MM300 mixer mill in Eppendorf tubes, with one tungsten carbide bead per tube, for 45 s with 30 oscillations per second. Petioles were ground with a pestle and mortar. Stem base and stem residues were ground with a Retsch MM300 mixer mill using Zirconium oxide jars and beads, for 40 s with 30 oscillations per second. For RNA extraction, all grinding material was frozen with liquid nitrogen prior grinding. RNA extraction was then performed using the Trizol reagent (Invitrogen, Cergy Pontoise, France) as previously described⁶. DNA was extracted from freeze-dried material using the plant DNeasy mini-kit or the Wizard® Genomic DNA Purification Kit (Promega, Madison, WI, USA) according to the manufacturers'
<figure><img src="image_3.png" /><figcaption>Figure 2. General description of the gene prediction workflow.</figcaption></figure>

<table><caption>Table 5. Metrics of the existing and new gene predictions.</caption><thead><tr><th rowspan="2">Reference</th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans JN3</th><th>Leptosphaeria maculans Nz-T4</th><th>Leptosphaeria biglobosa B3.5</th><th>Leptosphaeria biglobosa G12-14</th></tr><tr><th>Rouxel et al.<sup>1</sup></th><th>This study</th><th>This study</th><th>Grandaubert et al.<sup>2</sup></th><th>This study</th></tr></thead><tbody><tr><td># genes</td><td>12,611</td><td>13,047</td><td>14,026</td><td>11,390</td><td>12,678</td></tr><tr><td># mono-exonic genes</td><td>2,931</td><td>5,204</td><td>5,898</td><td>2,726</td><td>4,630</td></tr><tr><td>Gene length (avg:med)</td><td>1,592:1,278</td><td>1,652:1,341</td><td>1,507:1,211</td><td>1,501:1,217</td><td>1,679:1,337</td></tr><tr><td># exons per gene (avg:med)</td><td>2.94:2</td><td>2.28:2</td><td>2.20:2</td><td>2.67:2</td><td>2.35:2</td></tr><tr><td>CDS length (avg:med)</td><td>1,392:1,110</td><td>1,177:891</td><td>1,107:795</td><td>1,307:1,065</td><td>1,208:954</td></tr><tr><td># introns</td><td>24,475</td><td>16,700</td><td>16,887</td><td>18,998</td><td>17,076</td></tr><tr><td>introns length (avg:med)</td><td>103:63</td><td>127:57</td><td>98:57</td><td>116:57</td><td>151:56</td></tr><tr><td>Coding fraction</td><td>38.9%</td><td>33.4%</td><td>35.8%</td><td>46.8%</td><td>43.8%</td></tr><tr><td>BUSCO (euk)</td><td>89%</td><td>98%</td><td>96%</td><td>93%</td><td>96%</td></tr><tr><td>BUSCO (fungi)</td><td>84%</td><td>96%</td><td>93%</td><td>90%</td><td>95%</td></tr></tbody></table>

instructions. Input DNA for the nanopore sequencing, based on long DNA molecules, was extracted from agarose plugs of concentrated pycnidiospores (>1.10⁹ spore mL⁻¹) partially digested¹,³. Briefly, fungal pycnidiospores were collected from V8-agar plates³ with sterile distilled water, centrifuged (15 min at 10000 × g) and adjusted to 6.10⁹ spore mL⁻¹. The spore suspension was mixed with an equal volume of 2.5% low melting agarose (Seaplaque GTG* in TSE; Tris 25 mM, pH 7.5; Sorbitol 1 M, EDTA 25 mM) maintained at 40 °C. This suspension was poured in disposable plug molds (BioRad). After cooling, the resulting plugs were treated for 20 h in 0.5 M EDTA, SDS 10%, 1 mg.mL⁻¹ proteinase K (Sigma) at 50 °C under mild stirring, then rinsed three times with 0.5 M EDTA, 50 °C and stored at 4 °C in 0.5 M EDTA. For a part of Nanopore sequencing runs, DNA was also extracted from freeze-dried mycelium ground in liquid nitrogen: cells have been then lysed using CTAB (Cetyltrimethylammonium bromide) with 2% of PVP (Polyvinylpyrrolidone) and 0.1% of 2-mercaptoéthanol and DNA was purified following a phenol-chloroform protocol.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Genome browser database for the three *Leptosphaeria* isolates. The genome browser is available at http://www.genoscope.cns.fr/leptolife and contains repeats (green track), coverage of genomic reads (black wiggle), gene prediction (blue track), RNA contigs (dark green track) and protein homologies (salmon track).</figcaption></figure>

## Illumina PCR-Free library preparation and sequencing

DNA (4.5 to 6 µg) was sonicated to a 100- to 1500-bp size range using a Covaris E210 sonicator (Covaris, Woburn, MA, USA). Fragments were end-repaired using the NEBNext End Repair Module (New England Biolabs, Ipswich, MA, USA) and 3'-adenylated with the NEBNext dA-Tailing Module (New England Biolabs). Illumina adapters were then added using the NEBNext Quick Ligation Module (New England Biolabs) and ligation products were purified with AMPure XP beads (Beckmann Coulter Genomics, Danvers, MA, USA). Libraries were sequenced on an Illumina MiSeq (G12-14 and N2T4 genomes) or a HiSeq 2500 (JN3 genome) instrument (San Diego, CA, USA) using 250 base-length read chemistry in a paired-end mode.

## Illumina cDNA library preparation and sequencing

RNA-Seq library preparations were carried out from 1–2 µg total RNA using the TruSeq Stranded mRNA kit (Illumina, San Diego, CA, USA), which allows mRNA strand orientation (sequence reads occur in the same orientation as anti-sense RNA). Briefly, poly(A)⁺ RNA was selected with oligo(dT) beads, chemically fragmented and converted into single-stranded cDNA using random hexamer priming. Then, the second strand was generated to create double-stranded cDNA. cDNA were then 3'-adenylated, and Illumina adapters were added. Ligation products were PCR-amplified. Each library was sequenced using 101 bp paired end reads chemistry on a HiSeq2000 Illumina sequencer.

## Nanopore 8-kb and 20-kb libraries preparation

MinION sequencing libraries were prepared according to the nanopore Sequencing Kit protocol SQK-MAP006 or the Low Input Expansion Pack protocol for genomic DNA. Briefly, 100 ng to 1.5 µg of genomic DNA was sheared to approximately 8 Kb with g-TUBE (Covaris, Woburn, MA, USA) and
cleaned-up using AMPure XP beads (Beckmann Coulter Genomics). For some libraries, DNA fragments were repaired using the NEBNext FFPE Repair Mix (New England Biolabs). After end-reparation and 3'-adenylation with the NEBNext® Ultra™ II End Repair/dA-Tailing Module (New England Biolabs), sequencing adapters provided by Oxford Nanopore Technologies (Oxford Nanopore Technologies Ltd, UK) were ligated using Blunt/TA Ligase Master Mix (New England Biolabs). Libraries were then enriched using Dynabeads MyOne C1 Streptavidin (ThermoFisher Scientific, Wilmington, DE, USA).

Nanopore 20-Kb libraries were prepared according to the same protocol, with the exception that 250 ng to 2.5 µg of genomic DNA was sheared to approximately 20-kb with g-TUBE (Covaris).

## MinION flow cell preparation and sample loading

Each library was mixed with the fuel mix and the running buffer according to the SQK-MAP006 or the Low Input Expansion Pack protocols. The sequencing mix was then added to the R7.3 flowcell for a 48-hour run. Several loading schedules were tested, but the main one used for the MAP006 libraries has been the following: a first loading with 8 µl of the library, then three reloading after 5, 24 and 29 h of run with respectively 4 µl, 8 µl and 4 µl of the DNA library. Regarding the low-input libraries, the main schedule used was a first loading with 10 µl of the DNA library and a reloading after 24 hours of run time with 10 µl of the DNA library.

## MinION sequencing and base calling

Read event data generated by MinKNOW control software (version 0.50.2.15 to 0.51.3.40) were base-called using the Metrichor software (version 2.34.3 to 2.39.3). The data generated (pores metrics, sequencing, and base-calling data) by MinION software were stored and organized using a Hierarchical Data Format. Three types of reads were obtained: template, complement, and two-directions (2D). The template and complement reads correspond to sequencing of the two DNA strands. Metrichor combines template and complement reads to produce a consensus (2D) sequence⁷. FASTA reads were extracted from MinION Hierarchical Data Format files using poretools⁸. Pass and fail reads were both considered as 1D reads. This dataset is described in Table 1 and (Genomic datasets, Data Citation 5).

## Genetic map

A genetic map was built using 93 progeny of the JN2 × Nz-T4 *in vitro* cross⁹. A total of 150 polymorphic micro- and minisatellite markers, three biological markers (Mat locus, *AvrLm3* and *AvrLm7*) and eight markers designed at the junction between two repeated elements¹⁰ were first used to generate an unsaturated genetic map using Mapmaker 3.0, with a LOD score of 3 and a maximum genetic distance of 20 cM. The 93 progeny isolates and the two parental isolates were also sequenced (2 × 100 bp illumina sequencing; 48 samples per line; Beckman Coulter Genomics). SNP calling between parental isolates was done as described in the study¹¹ following mapping of the reads on the JN3 reference genome. A total of 20,103 reliable SNPs were detected and assigned to one of the two parental isolate in the progeny. These SNPs were converted into 2,104 Bins. From all these markers, a saturated map comprising 20 linkage groups was constructed using MSTmap online.

## Long reads-based de novo assembly of *L. biglobosa* G12-14 and *L. maculans* Nz-T4

R7.3 nanopore 2D reads (Table 3) were assembled into contigs using SmartDeNovo (https://github.com/ruanjue/smartdenovo), with a k-mer size parameter set to 21. Eventually, the consensus of the assembly was polished by aligning the Illumina paired-end 2 × 250 bp reads (Table 2) with BWA mem¹² and using the resulting alignment file as input for Pilon¹³. We iteratively polished twice the assembly.

*L. biglobosa* G12-14 final assembly had a cumulative size of 34.9 Mb and a N50 equal to 462 kb. *L. maculans* Nz-T4 final assembly had a cumulative size of 43.4 Mb and a N50 equal to 383 kb. Further usual metrics regarding those assemblies are provided in Table 4.

## Genome assembly update of *L. maculans* JN3

With 33 scaffolds, a N50 equal to 2.40 Mb and a cumulative size of 44.8 Mb, the JN3 rearranged reference genome was already contiguous enough for our analyses. However, it contained gaps that represented 1.1 Mb in total. In addition, both the optical map² and the high-density genetic map indicated that nine mis-assembled super contigs had to be either splitted or fused. The reference genome was scaffolded and gap-filled using PBJelly¹⁴ with the entire nanopore dataset as input (Table 3). As a result we obtained 33 scaffolds with a cumulative size of 45.8 Mb, the N50 increased to 2.43 Mb and we decreased by half the number of N's (570 kb).

As the R7.3 nanopore reads were prone to errors, particularly in the homopolymeric regions¹⁵, we polished the consensus of the gap closed assembly by aligning 150X of Illumina paired-end 2 × 250 bp reads with BWA mem¹² and giving the resulting alignment file to Pilon¹³. The final assembly had a cumulative size of 45.99 Mb with 1% of unknown bases (Table 4).

## Genome annotation

The gene prediction of the three *Leptosphaeria* genomes was done by combining gene models information computed from transcript and protein alignments (Fig. 2).
First, repeated sequences in the reference genomes were masked using multiple tools. Tandem repeats were identified and masked using Tandem Repeats Finder¹⁶. Low complexity regions were identified and masked using Dust¹⁷. Interspersed repeats and other low complexity sequences were masked using RepeatMasker (http://www.repeatmasker.org). Furthermore 121 known transposable elements (identified in the *L. maculans*-*L. biglobosa* species complex in a previous study²) were also masked using RepeatMasker.

The *Alternaria alternata*, *Cochliobolus heterostrophus*, *Pyrenochaeta* sp., *Pyrenophora tritici-repentis* and *Zymoseptoria tritici* proteomes included in UniProt (version available as of 07 September 2016) were used to detect conserved proteins in *L. maculans* JN3, *L. maculans* Nz-T4 and *L. biglobosa* G12-14. *L. maculans* JN3 and swissprot proteins (release available as of 30 November 2016) were also used as a resource to detect conserved proteins in *L. maculans* Nz-T4 and *L. biglobosa* G12-14.

Those proteins were first aligned to *L. maculans* JN3, *L. maculans* Nz-T4 and *L. biglobosa* G12-14 genomes assemblies using BLAT¹⁸. BLAT results were then filtered as follows: for each protein, the best match (based on the BLAT score) and other matches with a score > = 90% of the score of the best match were retained. Genewise¹⁹ was then used to refine matches and identify exon-intron boundaries and alignment were filtered out if less than 50% of the protein length was aligned.

RNAseq reads from 23 samples were assembled using Velvet²⁰ 1.2.07 and Oases²¹ 0.2.08, using a k-mer size of 63 bp. The results of this assembly for each sample are summarized in (Transcriptomic datasets, Data Citation 5). Reads were mapped back to the contigs with BWA-mem¹² and the consistent paired-end reads were selected. Chimeric contigs were identified and splitted (uncovered regions) based on coverage information from consistent paired-end reads. Moreover, open reading frames (ORF) and domains were searched using respectively TransDecoder²² and CDDsearch²³. We only allowed breaks outside ORF and domains. At the end, the read strand information was used to correctly orient the RNA-Seq contigs. In a second step, RNAseq contigs were aligned to the relevant genome assembly (*L. maculans* JN3, *L. maculans* Nz-T4, *L. biglobosa* G12-14) using BLAT. The best matches (based on BLAT score) for each contig and other matches with a score > = 90% of the score of the best match were selected. Then, Est2genome²⁴ was used to refine the alignments and we kept alignments with an identity percent > = 95% and at least 80% of the RNA-Seq contig that can be aligned (Transcriptomic datasets, Data Citation 5).

Gene predictions from mRNA and proteomes were given as input to the GMOVE combiner²⁵ to build the gene models. This tool considers biological data such as RNA-seq from the organisms of interest and proteome from homologous species. Proteins give information about the frame of the coding sequence (CDS) and RNA-seq helps to find splicing sites with more accurate alignments. This tool is designed for all eukaryotes and *de novo* studies with no need of calibration pre-processing. After gene models were produced by GMOVE, a final filter was applied to remove artefactual untranslated regions (UTRs). UTRs were trimmed if it overlapped CDS in predicted gene models for more than 10 bases. Moreover, UTRs overlapping other UTRs for more than 10 bases were also trimmed to avoid large and potentially false UTRs. Final metrics of the gene prediction are available in Table 5 and Fig. 3 is a screenshot of the genome browser.

## Code availability

All the tools used for this study are cited in the method section and they were used with default settings unless options were specified. The in-house tools to clean illumina data are available from the following website: http://www.genoscope.cns.fr/fastxtend/. The integration of resources in the gene prediction workflow was achieved using GMOVE: http://www.genoscope.cns.fr/gmove/.

## Data Records

The authors declare that all data reported herein are fully and freely available from the date of publication. Genomic data have been deposited in the ENA repository (Genomic datasets, Data Citation 5): *L. maculans* JN3 (Data Citation 1), *L. maculans* Nz-T4 (Data Citation 2) and *L. biglobosa* G12-14 (Data Citation 3). RNA data used for gene prediction have been deposited in the ENA repository (Data Citation 4 and Transcriptomic datasets, Data Citation 5). Moreover, a dedicated website brings together all the information: accession numbers, assembly and annotation files as well as a genome browser (Fig. 3): http://www.genoscope.cns.fr/leptolife.

## Technical Validation

### DNA and RNA sample quality

DNA quality was assessed using 1% agarose gel. RNA integrity was assessed using an Agilent 2100 Bioanalyzer electrophoresis system (Agilent, Santa Clara, CA, USA).

### Illumina libraries

Ready-to-sequence Illumina libraries were quantified by qPCR using the KAPA Library Quantification Kit for Illumina Libraries (KapaBiosystems, Wilmington, MA, USA), and libraries profiles evaluated with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
<table><caption>**Table 6.** Genome compartmentalization of old and new assemblies. *Results from Grandaubert *et al.* AT rich regions were manually curated.</caption><thead><tr><th rowspan="2">Reference</th><th rowspan="2"></th><th>LmJN3<sup>*</sup></th><th>LmJN3</th><th>LmNz-T4</th><th>LbB3.5</th><th>LbG12-14</th></tr><tr><th>Rouxel et al.<sup>1</sup></th><th>This study</th><th>This study</th><th>Grandaubert et al.<sup>2</sup></th><th>This study</th></tr></thead><tbody><tr><td rowspan="7">GC blocks</td><td># blocks</td><td>399</td><td>534</td><td>587</td><td>389</td><td>318</td></tr><tr><td>% assembly</td><td>64</td><td>63.1</td><td>66.6</td><td>95.1</td><td>84.5</td></tr><tr><td>mean size (kb)</td><td>70.4</td><td>54.4</td><td>49.2</td><td>77.7</td><td>92.9</td></tr><tr><td>stdev (kb)</td><td>—</td><td>63.1</td><td>61.9</td><td>226</td><td>157</td></tr><tr><td>min (kb)</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>max (kb)</td><td>500</td><td>431</td><td>349</td><td>1708</td><td>1180</td></tr><tr><td># genes</td><td>—</td><td>12,892</td><td>13,790</td><td>11,361</td><td>12,603</td></tr><tr><td rowspan="7">AT blocks</td><td># blocks</td><td>413</td><td>564</td><td>779</td><td>324</td><td>335</td></tr><tr><td>% assembly</td><td>36</td><td>36.9</td><td>33.4</td><td>4.87</td><td>15.5</td></tr><tr><td>mean size (kb)</td><td>38.6</td><td>30.1</td><td>18.6</td><td>4.8</td><td>16.2</td></tr><tr><td>stdev (kb)</td><td>—</td><td>48.1</td><td>23.5</td><td>10.7</td><td>17.7</td></tr><tr><td>min (kb)</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>max (kb)</td><td>320</td><td>319</td><td>245</td><td>120</td><td>142</td></tr><tr><td># genes</td><td>—</td><td>233</td><td>236</td><td>29</td><td>75</td></tr></tbody></table>

## Illumina reads processing and quality filtering

After the Illumina sequencing, an in-house quality control process was applied to the reads that passed the Illumina quality filters. The first step discards low-quality nucleotides ($Q < 20$) from both ends of the reads. Next, Illumina sequencing adapters and primer sequences were removed from the reads. Then, reads shorter than 30 nucleotides after trimming were discarded. These trimming and removal steps were achieved using in-house-designed software based on the FastX package. The last step identifies and discards read pairs that mapped to the phage phiX genome, using SOAP²⁶ and the phiX reference sequence (GenBank: NC_001422.1). This processing resulted in high-quality data and improvement of the subsequent analyses. This dataset is described in Table 2 and (Genomic datasets, Data Citation 5). A specific filter aiming to remove ribosomal reads was applied to data generated from RNA samples. The reads and their mates that mapped onto a ribosomal sequences database were filtered using SortMeRNA²⁷. It contains different rRNA databases and we used it to split the data into two files: rRNA reads in a file (ribo_clean) and other reads in another file (noribo_clean). This dataset is described in (Transcriptomic datasets, Data Citation 5).

## Impact of the gap closing procedure on the JN3 assembly

The JN3 assembly was gap closed using PBJelly and we detected several duplicated regions of several Kb in this gap-closed assembly that were not present in the previous assembly. We first detected tandemly duplicated genes using Blast and revealed 28 candidate regions totaling 150 kb in length with a median size of 5.5 kb and covering 63 genes. Further investigations of these regions using genomic data coverage revealed that they were artefactual duplicates created by PBJelly. We mapped both Illumina and nanopore reads onto the gap-closed assembly and highlighted genomic regions without illumina and nanopore reads coverage. This method revealed 30 regions with a cumulative size of 185 kb and a median size of 5.1 kb spanning a total of 43 genes. Finally, after a careful expertise using the genome browser (Fig. 3), the genomic regions that presented duplicated genes and no proper read support were masked and overlapping genes were filtered out. In total, we masked 41 regions with a cumulative size of 264 kb and a median size of 5.5 kb that spanned 81 genes (Duplicated genomic regions, Data Citation 5).

## Contamination detection

In order to detect possible contamination in our assemblies, we first ran the metagene²⁸ gene predictor on NzT4 and G12-14 contigs and on JN3 scaffolds. Then, we aligned all predicted proteins against the NR database using blastp. Hits were then filtered by their e-value to only conserve the best hit with an e-value of at least $10^{-10}$. Finally, a contig or scaffold was flagged as non-contaminant if at least 50% of its predicted proteins were assigned as fungi. Following this procedure, we couldn't detect any contaminant sequences in JN3 and G12-14 assemblies. We detected 13 small contigs that were not assigned as fungi in the NzT4 assembly, but further investigations (blat against the JN3 assembly and blastx using the NCBI blast server) revealed that they were not contaminants.
## Comparison with existing assemblies

An AT/GC region segmentation was performed on all assemblies with OcculterCut²⁹ to detect potential genome compartmentalization as described in the following study¹. *L. maculans* JN3 and Nz-T4 isolates showed similar genome architecture with a third of their assembly dispatched into AT blocks of 30.1 kb and 18.6 kb on average, respectively. The number of genes present within the AT-blocks was also steady, and ranged between 233 and 236. In contrast, there was a drastic modification in genome architecture and genome size between *L. biglobosa* B3.5, sequenced using only Illumina, and G12-14 sequenced here. The two assemblies differed by more than 10% of their genome content. This difference corresponded to 3 Mb missing in B3.5 and recovered in G12-14 thanks to the nanopore long reads (Table 6). The difference mostly comprised of AT-blocks similar to those found in *L. maculans* and amounting to 15.5% of the G12-14 assembly vs. only 4.9% in the B3.5 assembly (Table 6). The increase in AT-block content is consistent with the lowering of the average GC content of the *L. biglobosa* genome in the new assembly (Table 4).

## Assessment of gene prediction

First, protein sequences of each genome (*L. maculans* JN3, *L. maculans* Nz-T4 and *L. biglobosa* G12-14) were obtained by translating gene CDS and aligned to protein databases (uniprot_sprot and uniprot_trembl 2016_11 release however *L. maculans* matches were excluded to prevent the introduction of known biases from the previous annotation effort) using blastp. Gene models for which proteins aligned with an e-value < 10⁻¹⁰ or with a CDS size > 200 nt were retained.

Moreover gene completion for the newly generated annotations was evaluated using BUSCO³⁰ v1.1b1 fungal and eukaryotic conserved genes databases. The relevance in the increase in gene number between the “old” and “new” gene prediction reflects a higher accuracy of the “new” annotation as evidenced by careful examination of genes encoding effectors. Effector genes encode for small (less than 300 amino-acids) secreted proteins (SSPs). They are usually poorly predicted *de novo* and may lack transcriptomic support due to expression in certain stages of parasitism only. The data obtained here drastically increased the number of predicted SSPs in both species (1,082 for *L. maculans* vs. 651 in the “old” annotation, 875 for *L. biglobosa* vs. 665 in the previous annotation) and drastically improved the gene models (data not shown). Also this new annotation allowed the identification of 75 genes hosted within the AT-block of the genome of *L. biglobosa*, previously unidentified (Table 6).

## References

1. Rouxel, T. et al. Effector diversification within compartments of the *Leptosphaeria maculans* genome affected by Repeat-Induced Point mutations. *Nature communications* **2**, 202, doi:10.1038/ncomms1189 (2011).

2. Grandaubert, J. et al. Transposable element-assisted evolution and adaptation to host plant within the *Leptosphaeria maculans* *Leptosphaeria biglobosa* species complex of fungal pathogens. *BMC genomics* **15**, 891, doi:10.1186/1471-2164-15-891 (2014).

3. Leclair, S., Ansan-Melayah, D., Rouxel, T. & Balesdent, M. Meiotic behaviour of the minichromosome in the phytopathogenic ascomycete *Leptosphaeria maculans*. *Current genetics* **30**, 541–548 (1996).

4. Parlange, F. et al. *Leptosphaeria maculans* avirulence gene *AvrLm4-7* confers a dual recognition specificity by the *Rlm4* and *Rlm7* resistance genes of oilseed rape, and circumvents *Rlm4*-mediated recognition through a single amino acid change. *Molecular microbiology* **71**, 851–863, doi:10.1111/j.1365-2958.2008.06547.x (2009).

5. Gervais, J. et al. Different waves of effector genes with contrasted genomic location are expressed by *Leptosphaeria maculans* during cotyledon and stem colonization of oilseed rape. *Molecular plant pathology* **18**, 1113–1126, doi:10.1111/mpp.12464 (2017).

6. Fudal, I. et al. Heterochromatin-like regions as ecological niches for avirulence genes in the *Leptosphaeria maculans* genome: map-based cloning of *AvrLm6*. *Mol plant-microbe interact.*: *MPMI* **20**, 459–470, doi:10.1094/MPMI-20-4-0459 (2007).

7. Quick, J., Quinlan, A. R. & Loman, N. J. A reference bacterial genome dataset generated on the MinION portable single-molecule nanopore sequencer. *GigaScience* **3**, 22, doi:10.1186/2047-217X-3-22 (2014).

8. Loman, N. J. & Quinlan, A. R. Poretools: a toolkit for analyzing nanopore sequence data. *Bioinformatics* **30**, 3399–3401, doi:10.1093/bioinformatics/btu555 (2014).

9. Kuhn, M. L. et al. Genetic linkage maps and genomic organization in *Leptosphaeria maculans*. *Eur J Plant Pathol* **114**, 17–31, doi:10.1007/s10658-005-3168-6 (2006).

10. Plissonneau, C. et al. A game of hide and seek between avirulence genes *AvrLm4-7* and *AvrLm3* in *Leptosphaeria maculans*. *New Phytologist* **209**, 1613–1624, doi:10.1111/nph.13736 (2016).

11. Zhong, Z. et al. A small secreted protein in *Zymoseptoria tritici* is responsible for avirulence on wheat cultivars carrying the *Stb6* resistance gene. *New Phytologist* **214**, 619–631, doi:10.1111/nph.14434 (2017).

12. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* **25**, 1754–1760, doi:10.1093/bioinformatics/btp324 (2009).

13. Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. *PloS one* **9**, e112963, doi:10.1371/journal.pone.0112963 (2014).

14. English, A. C. et al. Mind the gap: upgrading genomes with Pacific Biosciences RS long-read sequencing technology. *PloS one* **7**, e47768, doi:10.1371/journal.pone.0047768 (2012).

15. Istace, B. et al. de novo assembly and population genomic survey of natural yeast isolates with the Oxford Nanopore MinION sequencer. *GigaScience* **6**, 1–13, doi:10.1093/gigascience/giw018 (2017).

16. Benson, G. Tandem repeats finder: a program to analyze DNA sequences. *Nucleic acids research* **27**, 573–580 (1999).

17. Morgulis, A., Gertz, E. M., Schaffer, A. A. & Agarwala, R. A fast and symmetric DUST implementation to mask low-complexity DNA sequences. *Journal of Computational Biology* **13**, 1028–1040, doi:10.1089/cmb.2006.13.1028 (2006).

18. Kent, W. J. BLAT--the BLAST-like alignment tool. *Genome Res.* **12**, 656–664, doi:10.1101/gr.229202 (2002).

19. Birney, E., Clamp, M. & Durbin, R. GeneWise and Genomewise. *Genome Res.* **14**, 988–995, doi:10.1101/gr.1865504 (2004).

20. Zerbino, D. R. & Birney, E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. *Genome Res.* **18**, 821–829, doi:10.1101/gr.074492.107 (2008).

21. Schulz, M. H., Zerbino, D. R., Vingron, M. & Birney, E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. *Bioinformatics* **28**, 1086–1092, doi:10.1093/bioinformatics/bts094 (2012).
22. Haas, B. J. et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. *Nature protocols* **8**, 1494–1512, doi:10.1038/nprot.2013.084 (2013).

23. Marchler-Bauer, A. & Bryant, S. H. CD-Search: protein domain annotations on the fly. *Nucleic acids research* **32**, W327–W331, doi:10.1093/nar/gkh454 (2004).

24. Mott, R. EST_GENOME: a program to align spliced DNA sequences to unspliced genomic DNA. *Computer applications in the biosciences: CABIOS* **13**, 477–478 (1997).

25. Dubarry, M, N. B. et al. Gmove a tool for eukaryotic gene predictions using various evidences. *F1000Research*, doi:10.7490/f1000research.1111735.1 (2016).

26. Li, R. et al. SOAP2: an improved ultrafast tool for short read alignment. *Bioinformatics* **25**, 1966–1967, doi:10.1093/bioinformatics/btp336 (2009).

27. Kopylova, E., Noe, L. & Touzet, H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. *Bioinformatics* **28**, 3211–3217, doi:10.1093/bioinformatics/bts611 (2012).

28. Noguchi, H., Park, J. & Takagi, T. MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. *Nucleic acids research* **34**, 5623–5630, doi:10.1093/nar/gkl723 (2006).

29. Testa, A. C., Oliver, R. P. & Hane, J. K. OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes. *Genome biology and evolution* **8**, 2044–2064, doi:10.1093/gbe/evw121 (2016).

30. Simao, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. *Bioinformatics* **31**, 3210–3212, doi:10.1093/bioinformatics/btv351 (2015).

## Data Citations

1. *European Nucleotide Archive* PRJEB24468 (2018).

2. *European Nucleotide Archive* PRJEB24469 (2018).

3. *European Nucleotide Archive* PRJEB24467 (2018).

4. *European Nucleotide Archive* PRJEB21682 (2018).

5. Dutreux, F. et al. figshare https://doi.org/10.6084/m9 figshare.c.4124543 (2018).

## Acknowledgements

The authors are grateful to Oxford Nanopore Technologies Ltd for providing early access to the MinION device through the MinION Access Programme, and we thank the staff of Oxford Nanopore Technology Ltd for technical help. BI, CC, PW, and JMA are part of the MinION Access Programme (MAP) and JMA received travel and accommodation expenses to speak at Oxford Nanopore Technologies conferences. This work was supported by the Genoscope, the Commissariat à l'Energie Atomique et aux Energies Alternatives (CEA) and France Génomique (ANR-10-INBS-09-08). F. Dutreux was funded by a grant from PROMOSOL (Association pour la promotion de la sélection des plantes oléagineuses). The genetic map was obtained through a funding of the French National research agency (ANR project ANR-12-ADAP-0009 'Gandalf').

## Author Contributions

C.C. and A.L. optimized and performed the sequencing. J.L. prepared all RNA samples. M.H.B. supervised the production of RNA and DNA samples and performed the genetic map. F.D., C.D.S., B.I., A.C., N.L., E.J.G., L.A. and J.M.A. performed the bioinformatic analyses. T.R. and M.H.B. conceived the Leptolife project and led the Promosol project. F.D. and J.M.A. wrote the article. P.W. and J.M.A. supervised the study.

## Additional Information

**Competing interests:** The authors declare no competing interests.

**How to cite this article:** Dutreux, F. et al. De novo assembly and annotation of three *Leptosphaeria* genomes using Oxford Nanopore MinION sequencing. *Sci. Data.* 5:180235 doi: 10.1038/sdata.2018.235 (2018).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files made available in this article.

© The Author(s) 2018