Gates Open Res Gates Open Research 2572-4754 F1000 Research Limited London, UK 10.12688/gatesopenres.12856.2 Research Article Articles Insights into population structure of East African sweetpotato cultivars from hybrid assembly of chloroplast genomes

[version 2; peer review: 2 approved, 3 approved with reservations]

 Zhou Chenxi Conceptualization Data Curation Formal Analysis Methodology Software Visualization Writing – Original Draft Preparation 1 Duarte Tania Data Curation Writing – Review & Editing https://orcid.org/0000-0002-3485-742X 1 Silvestre Rocio Data Curation Resources Writing – Review & Editing 2 Rossel Genoveva Data Curation Resources Writing – Review & Editing 2 Mwanga Robert O. M. Funding Acquisition Project Administration Writing – Review & Editing https://orcid.org/0000-0002-3166-0015 3 Khan Awais Funding Acquisition Project Administration Writing – Review & Editing 2 4 George Andrew W. Supervision Writing – Review & Editing https://orcid.org/0000-0001-6916-0684 5 Fei Zhangjun Data Curation Funding Acquisition Project Administration Writing – Review & Editing 6 Yencho G. Craig Funding Acquisition Project Administration Writing – Review & Editing 7 Ellis David Project Administration Resources Writing – Review & Editing 2 Coin Lachlan J. M. Conceptualization Funding Acquisition Project Administration Supervision Writing – Review & Editing https://orcid.org/0000-0002-4300-455X a 1 Institute for Molecular Bioscience, University of Queensland, St Lucia, Brisbane, QLD, 4072, Australia International Potato Center, P.O. Box 1558, Lima 12, Peru International Potato Center, P.O. Box 22274, Kampala, Uganda Plant Pathology and Plant-Microbe Biology Section, School of Integrative Plant Science, Cornell University, Geneva, NY, 14456, USA Data61, CSIRO, Ecosciences Precinct, Brisbane, QLD, 4102, Australia Boyce Thompson Institute, Cornell University, Ithaca, NY, 14853, USA Department of Horticulture, North Carolina State University, Raleigh, North Carolina, 27695, USA l.coin@imb.uq.edu.au

No competing interests were disclosed.

 21 7 2020 2018 2 41 17 7 2020 Copyright: © 2020 Zhou C et al. 2020 This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background: The chloroplast (cp) genome is an important resource for studying plant diversity and phylogeny. Assembly of the cp genomes from next-generation sequencing data is complicated by the presence of two large inverted repeats contained in the cp DNA.

Methods: We constructed a complete circular cp genome assembly for the hexaploid sweetpotato using extremely low coverage (<1×) Oxford Nanopore whole-genome sequencing (WGS) data coupled with Illumina sequencing data for polishing.

Results: The sweetpotato cp genome of 161,274 bp contains 152 genes, of which there are 96 protein coding genes, 8 rRNA genes and 48 tRNA genes. Using the cp genome assembly as a reference, we constructed complete cp genome assemblies for a further 17 sweetpotato cultivars from East Africa and an I. triloba line using Illumina WGS data. Analysis of the sweetpotato cp genomes demonstrated the presence of two distinct subpopulations in East Africa. Phylogenetic analysis of the cp genomes of the species from the Convolvulaceae Ipomoea section Batatas revealed that the most closely related diploid wild species of the hexaploid sweetpotato is I. trifida.

Conclusions: Nanopore long reads are helpful in construction of cp genome assemblies, especially in solving the two long inverted repeats. We are generally able to extract cp sequences from WGS data of sufficiently high coverage for assembly of cp genomes. The cp genomes can be used to investigate the population structure and the phylogenetic relationship for the sweetpotato.

 chloroplast sweetpotato genome assembly Oxford Nanopore sequencing Illumina sequencing phylogenetic analysis Convolvulaceae Ipomoea Gates Foundation OPP1052983 The project was financially supported by Gates Foundation (OPP1052983). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Amendments from Version 1

A number of misspellings and incorrect claims were revised following the reviewers’ comments. Some misleading and confusing sentences were corrected in the new version. Figure 1d has been slightly modified to remove red dots indicating SNPs. The figure legend has been modified accordingly.

 Introduction

The chloroplast (cp) genome has been widely used to study the phylogeography, molecular systematics and the population genetics for plants 1, 2 . The chloroplast DNA (cpDNA) usually displays uniparental inheritance and represents a relatively high degree of conservation in genome structure and gene content 2 . There are over 800 complete cp sequences available for a wide variety of plants from National Center for Biotechnology Information (NCBI) repository ranging in size from 107 to 218 Kb 3 . The cp genomes usually contain 110–130 protein encoding genes (PEGs), about 30 transfer RNA (tRNA) genes and four ribosomal RNA (rRNA) genes, primarily participating in the process of photosynthesis 3, 4 . The cpDNA typically forms a circular quadripartite structure with two inverted repeats (IRs), IRA and IRB, separated by one large single-copy section (LSC) and one small single-copy section (SSC) 5 .

The first cpDNA was sequenced from tobacco ( Nicotiana tabacum) using the bacterial artificial chromosome (BAC) sequencing method in 1986 6 . The two IRs were cloned separately in order to distinguish between them. A plethora of cpDNA had since been sequenced with similar methods 7– 9 . Besides BAC sequencing, an alternative strategy used to sequence cpDNA is whole-cp-genome amplification by rolling-circle amplification (RCA) technology 10– 12 . However, both approaches require complicated library preparation.

The development of next-generation sequencing (NGS) technologies such as Illumina and Roche 454 facilitate faster and cheaper methods to sequence cp genomes 13– 15 . The output of the NGS technologies is short reads of size up to a few hundred base pairs. It is difficult to assemble cp genome with short reads only, especially because of the two large IRs of tens of kilobase pairs. In order to solve this problem, a reference cp genome, normally from a related species, is usually used to anchor the contigs assembled from the short reads 4, 16 . The long reads generated from the third-generation sequencing (TGS) technologies, such as the single-molecule real-time (SMRT) PacBio sequencing and Oxford Nanopore sequencing, can also be used to anchor the contigs and solve the repetitive regions. It is even possible to assemble cp genomes directly from long reads 17 . However, as the sequencing error rate of the long reads from the TGS is typically higher than 10%, it is important to introduce an error correction step to guarantee an accurate genome assembly 18 . The high-quality NGS short reads can be integrated for error correction to improve accuracy 19, 20 .

The aforementioned methods to construct cp genomes from NGS or TGS data assume pure cpDNA were sequenced. More precisely, the cpDNA were isolated from the nuclear DNAs and other organelle DNAs before sequencing 4, 13– 16 . However, whole-genome sequencing (WGS) data generated from NGS or TGS technologies always contains cp sequences at various levels determined by the tissue type and library preparation. Normally we are able to gain enough coverage of cp genome for assembly even from low coverage WGS data. There have since been several studies describing assembly of cp genomes from WGS data 21– 27 . Extraction of cp sequences from the WGS data plays a key role in these methods. The most straightforward idea is to use a reference cp genome. The cp sequences could be extracted by examining the mapping results of the WGS data to the reference cp genome 21, 22 . An alternative strategy relies upon the fact that there are many more copies of the cpDNA than the nuclear DNA and that from other organelles. The entire WGS data is assembled to construct contigs. Contigs that represent significantly higher coverages are treated as cp contigs 23– 25 . NOVOPlasty adopted a seed-and-extend paradigm, where the seed could be a cp read sequence, a conserved gene or a cp genome from a related species 26 . The start and the end of a given seed sequence are iteratively extended with reads that are overlapped with the seed until the circular genome is formed. Izan et al. proposed a K-mer frequency-based selection of cpDNA sequences from WGS data, which was integrated into a reference free cp genome assembler for non-model species 27 .

Sweetpotato ( Ipomoea batatas) ranks among the ten most important food crops worldwide 28 . The total annual production is more than 100 million metric tonnes grown on about 8.6 million hectares around the world in year 2016 29 . Understanding the sweetpotato genomes is of significant importance to achieve the full potential of the sweetpotato 30 . Sweetpotato is a hexaploid (2n=6x=90) with genome size estimated to be between 2,200 to 3,000 Mb 28 . Due to the complex genome structure, the availability of sweetpotato genomic resources is lacking. Under these circumstances, the cp genome provides researchers with an easy and efficient way to study sweetpotato 4, 16, 31, 32 . A number of cp genomes from the genus Ipomoea have been sequenced 4, 16, 33, 34 . Most of them are diploid wild relatives of the sweetpotato. The genome size is around 161 Kb, and the structure represents a standard quadripartite circular with a LSC of 87 Kb, a SSC of 12 Kb and two IRs of 31 Kb 4 . The cp genomes were mainly used to perform phylogenetic analyses 4, 16, 34 .

In the present study, we constructed a complete cp genome assembly for the hexaploid sweetpotato cultivar Tanzania 35 using long reads produced by the Oxford Nanopore sequencing technology. Despite the <1× genome coverage, we obtained approximately 270× data coverage for the cp genome. Illumina sequencing data was integrated to improve the accuracy of the genome assembly. Using the Tanzania cp genome assembly as a reference, we constructed 19 cp genomes for a further 17 sweetpotato cultivars (including a duplicate for one cultivar) and an I. triloba line from paired-end whole genome Illumina sequence data. The assembled sweetpotato cp genomes were combined to perform phylogenetic analysis to investigate the population structure of 18 East African sweetpotato cultivars. Putting together the assembled cp genomes and nine publicly available cp genomes of the sweetpotato and its wild relatives, we performed a phylogenetic analysis to investigate the phylogenetic relationship for species in Convolvulaceae Ipomoea section Batatas.

 Results Extraction of cp genome sequence from whole genome sequencing data

We generated high-coverage (60×) 150 bp paired-end Illumina WGS data, and low-coverage (<1×) Oxford Nanopore WGS data on a single cultivar, referred to as Tanzania 35 (Methods). The cultivar Tanzania was used as one of the parents to develop an F1 outcrossing mapping population (B×T) in the Genomic Tools for Sweetpotato (GT4SP) Improvement Project 30 . Approximately 162,000 Nanopore reads and 1.46 billion Illumina reads were generated (Supplementary Table 1). A total of 6,710 Nanopore reads were identified for cp genome by mapping to 30 publicly available cp genomes of the species from the Convolvulaceae Ipomoea family 4, 16, 33, 36 (Methods, Supplementary Table 2). The total size is ~43.9 Mb, which represents ~270× data coverage for the cp genome. The longest read is ~30 Kb, and the average size is ~6.5 Kb (Supplementary Figure 1). We identified approximately 45 million Illumina reads for cp genome by mapping to the publicly available cp genomes summing to ~6.2 Gb, which were used for error correction for Nanopore reads and the genome assembly. The other parent for the B×T F1 outcrossing mapping population, Beauregard, was subject to whole genome sequencing at 60× coverage (Methods). A total of approximately 1.3 billion 150 bp Illumina reads were generated summing to ~164 Gb, of which approximately 52 million reads were identified as cp sequences with a total size of ~7.2 Gb (Supplementary Table 1). We performed Illumina WGS at 30× coverage for a further 16 sweetpotato cultivars—Wagabolige and New Kawogo 35 , Ejumula and SPK004 37 , NASPOT 1 and NASPOT 5 38 , NASPOT 7 and NASPOT 10 O 39 , NK259L and NASPOT 11 40 , Huarmeyano, Dimbuka-Bukulula and NASPOT 5/58 41 , Resisto 42 , Magabali 43 and Mugande 44 . These cultivars were used as the parental genotypes in the Mwanga Diversity Panel (MDP) which is an 8×8 diallele diversity mating panel constructed by the GT4SP project for genomic selection of the sweetpotato. While the great majority of these sweetpotato cultivars were from East African countries including Uganda and Kenya, Resisto was from USA and Huarmeyano was from Peru (Supplementary Table 3). We have duplicate samples for the cultivar NASPOT 10 O—one was from the screen-house while the other one was from the field. These two NASPOT 10 O samples were analysed separately in this research (Methods). On average, a total of approximately 75 million 251 bp reads were generated for each sample. The number and the total size of the cp reads extracted from the whole genome sequence data, on average, are ~4.4 million and ~1 Gb respectively for each sample (Supplementary Table 1).

We performed Illumina whole genome sequencing at 50× coverage for the I. triloba line, NCNSP-0323 30 (Methods). The raw whole genome sequence data consists of approximately 196 million 150bp reads summing to ~29 Gb. We extracted approximately 13 million reads for the cp genome from the raw sequence data summing to ~2 Gb (Supplementary Table 1).

 Cp genome assembly for the sweetpotato cultivar Tanzania

We combined the Nanopore long reads with Illumina short reads to construct a cp genome assembly for the sweetpotato cultivar Tanzania (Methods). After trimming off the low-quality bases, approximately 2.2 Gb Illumina sequence data remained which was used for error correction for the Nanopore reads with Nanocorr (Supplementary Table 1). A total of 70 low quality Nanopore reads were removed after error correction and the total size reduced to approximately 43.2 Mb ( Figure 1a), which was used to construct a draft genome assembly using Canu. The resulting genome assembly of approximately 218 Kb consists of three contigs of size 46 Kb, 39 Kb and 132 Kb, respectively. Compared to the published sweetpotato cp genome, the assembly is split at the boundaries of the two IRs ( Figure 1b). Utilizing the overlap information between the contigs, the AMOS minimus combined the three contigs and generated a single contig of ~183 Kb ( Figure 1c) (Methods). The contig contains a ~20 Kb redundancy at the ends which was removed after circularization ( Figure 1d). The circularized contig is ~161 Kb, and is highly collinear with the reference cp genome assembly ( Figure 1d). Application of Pilon further identified and corrected 42 single-nucleotide polymorphisms (SNPs) and small indels. To follow the paradigm of the published cp genomes, we restructured the genome assembly so that it starts from the LSC (Methods). The final genome assembly consists of a single circular contig of 161,274 bp ( Figure 1e).

Assembly of the Tanzania chloroplast (cp) genome.

( a) Dot plot of the Nanopore read length versus the alignment identity to reference assembly. The read alignment identity is defined as I = M/L, where M is the total number of base pairs of the exact match and L is the size of the alignment span on the reference genome. The reference genome is the 30 cp genomes downloaded from the NCBI (Supplementary Table 2). The alignment was performed with BWA MEM 45 . The alignment identities were calculated from the Cigar string. The purple and yellow represents before and after error correction with Illumina reads using Nanocorr 20 , respectively. ( b) Dot plot of the reference cp genome versus the contigs produced by Canu 17 . ( c) Dot plot of the reference cp genome versus the contigs produced by AMOS minimus 46 after merging Canu contigs. ( d) Dot plot of the reference cp genome versus the contigs produced by AMOS minimus after circularization. ( e) Dot plot of the reference cp genome versus the final cp genome assembly which was polished with Illumina reads using Pilon 19 and fixed the start at the LSC. For ( b– e), the cp genome assembly of the I. trifida was used as the reference (accession number REM 753, Genbank accession number KF242496) 16 . The green bars on the x-axis indicate positions of the two IRs.

 Cp genome assembly for the other 17 sweetpotato cultivars and the <italic toggle="yes">I. triloba</italic> line NCNSP-0323

The cp sequence data was subjected to quality control before assembled with SPAdes (Methods). After trimming off the low-quality regions, the total sizes of the sequence data of the 19 samples range from approximately 267 Mb to 2.67 Gb (Supplementary Table 1). The contigs generated from SPAdes for the 19 samples vary in numbers and sizes: the minimum number of contigs is 76 for the cultivar NASPOT 7, while the maximum number is 197 for the cultivar Beauregard; and the total sizes of the genome assemblies range from ~169 Kb (cultivars Ejumula and NASPOT 7) to ~229 Kb (cultivar NK259L) (Supplementary Table 4). The SPAdes contigs were then mapped to the Tanzania cp genome assembly for anchoring (Methods). The resulting genome assemblies for the 19 samples are very similar. The largest and the smallest genome assembly is 161,509bp and 161,198bp, derived from the cultivar NASPOT 5 and Beauregard, respectively (Supplementary Table 4).

 Molecular structure and gene content of the sweetpotato cp genome

The gene annotation of the cp genome assembly of the sweetpotato cultivar Tanzania was generated with the web tool DOGMA and further refined with MUSCLE (Methods). The circular plot of the gene annotation is depicted in Figure 2. The sweetpotato cp genome represents a common circular structure with two IRs (IRA and IRB) separating one LSC and one SSC2. The size of the IRA, IRB, LSC and SSC is 30,874, 30,835, 87,489 and 12,076 bp, respectively. The overall GC content of the sweetpotato cp genome is 37.54%. The GC contents in different regions are highly variable. The two IRs represent significantly higher GC content than the single-copy regions: for the LSC and SSC, the GC content is 36.14% and 32.20%, respectively, whereas for the two IRs, the GC content is 40.57%. This is mainly caused by the high GC content ribosomal RNA genes in IR regions, including rrn16, rrn23, rrn4.5 and rrn5 ( Figure 2). We identified 152 genes in the cp genome of which there are 96 protein encoding genes (PEGs), eight rRNA genes and 48 tRNA genes. Table 1 shows a full list of the functional genes. As we can see, the genes can be divided into 16 functional systems. The number of single-copy and double-copy genes is 71 and 11, respectively, and there is one triple-copy gene ( rps12). The results are highly similar to what has been reported for the cultivar Xushu 18 cp genome 4 ; the only difference is that the psbZ gene is not found in the cultivar Xushu 18 cpDNA while the ihbA gene is not found in the cultivar Tanzania cpDNA. It should be noted that the double-copy gene ycf1 was not reported for the cultivar Xushu 18 cp genome4, but this was actually a miss-annotation.

The chloroplast genome of the sweetpotato cultivar Tanzania.

The preliminary annotations were produced by DOGMA 48 . MUSCLE 49 was used to refine the annotations. The plot was generated with OGDRAW 50 .

 List of annotated genes.

The functional systems were adopted from the OGDRAW 50 . Bracketed superscripts represent number of copies.

Functional system Number Gene list
Photosystem I 7 psaA, psaB, psaC, psaI, psaJ, ycf3, ycf4
Photosystem II 15 psbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ
Cytochrome b/f complex 6 petA, petB, petD, petG, petL, petN
ATP synthase 6 atpA, atpB, atpE, atpF, atpH, atpI
NADH dehydrogenase 13 ndhA, ndhB [2], ndhC, ndhD, ndhE, ndhF, ndhG, ndhH [2] , ndhI, ndhJ, ndhK
RubisCO large subunit 1 rbcL
C-type cytochrome synthesis 1 ccsA
RNA polymerase 4 rpoA, rpoB, rpoC1, rpoC2
Ribosomal proteins (LSU) 9 rpl2, rpl14, rpl16, rpl20, rpl22, rpl23, rpl32, rpl33, rpl36
Ribosomal proteins (SSU) 16 rps2, rps3, rps4, rps7 [2], rps8, rps11, rps12 [3] , rps14, rps15 [2] , rps16, rps18, rps19
Maturase K 1 matK
Acetyl-CoA carboxylase carboxyltransferase 1 accD
Clp protease proteolytic subunit 1 clpP
Chloroplast envelope membrane protein 1 cemA
ORFs 6 orf188 [2], orf42 [2], orf56 [2]
Hypothetical chloroplast RF 8 ycf1 [2], ycf15 [2], ycf2 [2], ycf68 [2]
Phylogenetic analysis of the sweetpotato cp genome

We performed a phylogenetic analysis for the Convolvulaceae Ipomoea section Batatas on the basis of the 19 cp genomes of the sweetpotato ( I. batatas) and the cp genome of the I. triloba line NCNSP-0323 assembled in this research, coupled with nine publicly available cp genomes, of which, four of them are for sweetpotato and two of them are for I. trifida and the other three are for I. cordatotriloba, I. splendor-sylvae and I. setosa, respectively 4, 16 (Supplementary Table 2). The resulted phylogenetic tree is depicted in Figure 3. The 18 sweetpotato cultivars used as the parental genotypes for mapping populations in the GT4SP project represent two distinct clades, consisting of 12 and six cultivars, respectively. Here, the length of any branch in a clade is no greater than 2×10 -4 substitutions per bp. The detailed phylogenetic relationship of the 18 sweetpotato cultivars is shown in Figure 4. As we can see, the distance between the two clades is approximately 5×10 -4 substitutions per bp. In the larger clades, the cultivar Tanzania represents a relatively larger distance (2×10 -4 substitutions ber bp) compared to the other cultivars. The population structure discovered here is similar to the one revealed by using simple sequence repeat primers by David et al. with the exception of the classification of the sweetpotato cultivars NK259L, Resisto and Mugande 47 (Supplementary Table 3). For the publicly available sweetpotato cp genomes, PI 561258 and Xushu 18 are closely related to the larger clade, while PI 518474 and PI 508520 have a closer relationship with the smaller clade ( Figure 3). The diploid wild relative of the hexaploid sweetpotato, I. trifida (REM 753), displays a significantly closer relationship to the I. batatas compared to the other species in the Convolvulaceae Ipomoea section Batatas. The other I. trifida accession PI 618966, however, represents a much larger diversity to the I. batatas and shows a close relationship to the I. triloba line NCNSP-0323 assembled in this research. Interestingly, the accession PI 618966 was originally identified as I. triloba and was recently reidentified as I. trifida by the GRIN National Genetic Resources Program. Among the other three species in the Convolvulaceae Ipomoea section Batatas, the I. cordatotriloba (REM 317) is closely related to the I. triloba (NCNSP-0323) and I. trifida (PI 618966) and therefore displays much closer relationship to the I. batatas compared to the I. splendor-sylvae (REM 763) and I. setosa (REM 68).

A phylogenetic tree of the Convolvulaceae <italic toggle="yes">Ipomoea</italic> section <italic toggle="yes">Batatas</italic> on the basis of chloroplast genomes.

The numbers on the branches are bootstrap support values. The branches shorter than 2×10 -4 substitutions per bp were collapsed resulting two clades consisting of 12 and 6 sweetpotato cultivars represented by a big and small solid circle respectively in the plot. The plot was generated with iTOL 52 .

 A phylogenetic tree of the East African sweetpotato cultivars used in the GT4SP project on the basis of chloroplast genomes.

This is a fine-scale representation of the two clades in Figure 3. The numbers on the branches are branch lengths given in terms of substitutions per bp.

 Discussion

The sweetpotato cp genome contains two ~31 Kb IRs which is very difficult for short-read de novo assemblers. There have been a few studies exploring the possibility to perform de novo assembly of organelle genomes with long reads especially with SMRT PacBio sequencing reads 21, 26, 51 . In this study, we constructed a complete sweetpotato cp genome assembly using the long reads generated from Oxford Nanopore sequencing. Nanopore reads proved to be extremely powerful in assembling the cp genome, especially in solving long repetitive regions. The sweetpotato cp genome contains two ~31 Kb IRs, which is very difficult for short-read de novo assemblers. With the overlapping information from long reads; however, the problem can be easily resolved. Canu 17 provides a useful tool set for assembling Nanopore reads, which was used in this research. It is worth noting that although the average depth of coverage of the whole sweetpotato genome is less than 1×, we obtained enough coverage of the cp genome for assembly.

Although long reads are powerful in solving complex genome structures, the error-prone nature of the raw reads necessitates an extra error-correction step. Illumina reads have been widely used to assist long read error-correction 19, 20 . The Illumina read-based correction could be performed either on the raw long reads before assembling 20 or on the draft genome assembly constructed from the raw reads 19 . In the current study, we did both. Before assembling with Canu, the Nanopore reads were corrected with Illumina reads using Nanocorr 20 . After assembling, the draft genome assembly was polished with Illumina reads using Pilon 19 (Methods). With several pipelines examined, we found that to perform error correction both before and after assembling is the best practice to construct the sweetpotato cp genome.

Assembling the cp genome from the short Illumina reads is challenging owing to the two large IRs. Since the structure of the cp genome is generally stable, reference genomes from the closely related species are usually used to perform reference-based assembling 4, 16 . In this study, we used the genome assembly constructed from the Nanopore reads as reference to assemble cp genomes for a further 19 cp genomes including 17 sweetpotato cultivars (including a duplicate for one sample) and the I. triloba line NCNSP-0323. SPAdes 53 was used as the de novo short-read assembler. The contigs generated by SPAdes were fragmented as expected. Among the 19 genome assemblies, the minimum number of contigs was 76. As the two IRs are highly homologous, there was generally only one copy of repetitive regions being assembled. In order to solve this problem, for reference-based scaffolding, we reused some single-copy contigs from the two IR regions to construct complete cp genome assemblies.

The molecular structure and gene content of the cpDNA are relatively conserved in land plants 2 . Many cpDNAs form a circular quadripartite structure with two IRs separated by one large and one small single-copy section 2, 5 . All 20 cp genome assemblies constructed in this research represent this common structure. The size of the two IRs of the sweetpotato cpDNA is approximately 31 Kb each, and is much larger than the other plants such as potato 10 , rice 54 , wheat 55 , and maize 56 , of which the IRs are usually smaller than 26 Kb. This is highly likely due to gene losses in these species. By comparing the gene annotation of the sweetpotato cpDNA in this study ( Figure 2) to the potato cpDNA 10 , we can see that, in the potato cpDNA, the boundary region of the IRA and SSC harbors a deletion of approximately 6Kb involved in the genes, ycf1, rsp15 and ndhH. Meanwhile, these three genes are presented in the symmetric boundary region of the IRB and SSC, which explains why the size of the IRs of the potato cpDNA is approximately 6 Kb smaller than the sweetpotato cpDNA.

The cpDNA usually has uniparental inheritance and undergoes low rates of substitution and recombination, which makes it well suited for phylogenetic analysis. The cp genome has been widely used to perform phylogenetic or comparative analysis in previous studies 2, 10, 16 . In this research, we used the complete cp genome assemblies to study the phylogenetic relationship of the 18 sweetpotato potato cultivars used as the parental genotypes for mapping populations in the GT4SP project, as well as the species from the Convolvulaceae Ipomoea section Batatas. The sweetpotato genotypes from the GT4SP project were classified into two distinct clusters, which guarantees the diversities of mapping populations derived from them. The phylogenetic analysis clearly revealed that the I. trifida is the most closely related diploid wild relatives to the hexaploid sweetpotato, I. batatas, which is consistent with conclusions from the previous studies 32, 57 .

Almost all whole genome sequencing data contains cp sequences, from which we are usually able to obtain cp genome sequences of enough data coverage for de novo assembly. As we can see, all the cp genome assemblies described in this research were constructed using whole genome sequencing data. Given that the cp genome is an important resource for studying plant genomes and whole genome data has gradually become indispensable in modern genome projects, it will be a good practice to construct the cp genome assembly to gain a first insight into the plant genome we are trying to understand before moving to the complex nuclear genome.

 Methods Genome sequencing of the MDP parental genotypes

The 16 sweetpotato cultivars used as the parental genotypes for the MDP diversity panel were subjected to whole genome sequencing. These sweetpotato cultivars were collected from Uganda, Kenya, USA and Peru, and included Wagabolige, New Kawogo, Ejumula, SPK004, NASPOT 1, NASPOT 5, NASPOT 7, NASPOT 10 O, NK259L, NASPOT 11, Huarmeyano, Dimbuka-Bukulula, NASPOT 5/58, Resisto, Magabali and Mugande (Supplementary Table 3). Leaf tissue was ground to a fine powder using the FastPrep-24 TM 5G tissue homogenizer (MP Biomedicals, Santa Ana, California) and DNA extracted from the leaf tissues following published protocols with modifications 58, 59 . Briefly, tissue was suspended in pre-warmed (65°C) CTAB buffer (200mM Tris-CL, 50mM EDTA, 2M NaCl, 2% CTAB and 3% β-mercapto-ethanol), mixed and heated at 65°C for 45 min prior to extraction with chloroform:isoamyl alcohol (24:1) and precipitated with sodium acetate and ethanol. Paired-end genomic libraries were prepared using the Illumina’s Genomic DNA Sample Preparation kit and sequenced on the Illumina HiSeq 2500 system with paired-end mode and read length of 251 bp (Illumina, San Diego, CA).

 10x Genomics’ Chromium sequencing of the sweetpotato cultivar Tanzania and Beauregard

The genomic DNA of Tanzania and Beauregard were extracted using the method cetyltrimethyl ammonium bromide and purified with 1× Agencourt AMPure XP beads (Beckman Coulter), according to manufacturer’s instructions. Before the library preparation, 1.5 µg purified gDNA was size selected using the BluePippin instrument (Sage Science) with the 0.75% Agarose Dye free, Marker U1 High-pass 30–40 kb vs3 protocol followed by a purification step with 0.4× AMPure XP beads. The library preparations for these two samples were done following the Chromium TM Genome Reagent Kits user guide (CG00022, Rev C). In summary, 10 ng of sample DNA was used to generate Gel Bead-In-Emulsions (GEM) in the Chromium TM Controller (10× Genomics) followed by isothermal incubation, post GEM incubation cleanup and quality control (QC). Libraries were constructed with end-repair and A-tailing, adaptor ligation, post ligation cleanup using SPRIselect Reagent (Beckman Coulter, USA), sample index PCR, post PCR cleanup, and QC. We modified the protocol by increasing the number of PCR cycles to nice and adding 105 µl SPRIselect reagent for the Post Sample Index PCR Cleanup, which resulted in the recovery of shorter fragments than it was expected. The libraries were sequenced using the HiSeq X Ten platform (Illumina, San Diego, CA).

 Oxford Nanopore sequencing of the sweetpotato cultivar Tanzania

Before the MinION library preparation, 5.7 µg Tanzania pure DNA was size selected (start selection size: 8Kb) with the same protocol used in 10x Genomics’ Chromium sequencing. The size selected gDNA was purified with 1× AMPure XP beads. The resulting 950 ng of Tanzania gDNA was used in MinION sequencing library preparation with the SQK-LSK108 1D ligation Sequencing kit (May 2017 version). We modified the protocol as follows: 30 min incubation each end-repair step and adapter ligation; 10 min incubation at RT in the end-repair purification step; 0.7× AMPure XP beads used after adapters ligation and ELB buffer (Oxford Nanopore Technologies) warmed up at 50°C previously to use and incubation of the eluted solution at 50°C. A library of 348 ng was loaded into a FLO-MIN106 (R.9.4 version) flowcell used in a MK1B MinION. We run the 1D protocol in the MinKnow software (version 1.5.18) and we basecalled the raw data using Albacore (version 1.1.0).

 Cp genome sequence extraction

WGS data were aligned to 30 publicly available cp genome assemblies of the species from the Ipomoea family 4, 16, 33, 36 (Supplementary Table 2) to extract cp genome reads, using BWA MEM 45 (version 0.7.15). We used the option ‘-x ont2d’ for Nanopore reads, and default options for Illumina reads. For each Nanopore read, the alignment records with at least 500 bp sequence aligned were selected to calculate the total length of the alignment. A Nanopore read was considered as a cp sequence if at least 1 Kb and 80% of the read aligned. A similar strategy was employed for Illumina reads extraction. Both of the two reads of a read pair were required to be aligned. The minimum size of the alignment block was set to 100 bp.

 Cp genome assembly from Nanopore data

We used Nanocorr 20 (version 0.01) to perform error correction for Nanopore reads using the Illumina reads. In order to guarantee the quality of Illumina reads, Trimmomatic 60 (version 0.36) was used to remove the low quality regions. We imposed the quality score of each base pair to be no less than 20 and the length of the reads no less than 100. The corrected Nanopore reads were then used to construct a draft genome assembly with Canu 17 (version 1.5). As the resulting draft genome assembly contained more than one contig, AMOS minimus 46 (version 3.1.0) was used to remove the redundancy and concatenate contigs using the overlap information. The AMOS minimus was also used to circularize the contig. We aligned the Illumina reads to the circularized contig and corrected the SNPs and small indels with Pilon 19 (version 1.22). In order to follow the paradigms of the published cp genomes, we aligned the genome assembly to the published cp genomes with MUMMER 61 (version 3.23) to find homology regions, and let the genome assembly start from the LSC.

 Cp genome assembly from Illumina Hiseq data

The low quality regions of the extracted cp sequences were removed with Trimmomatic 60 (version 0.36). The minimum quality score of each base pair was set to 20 and the minimum length of the reads was set to 100. SPAdes 53 (version 3.10.1) was used to construct contigs from Illumina reads. We excluded the repeat resolve module from SPAdes and used the contigs before repeat resolution as it consistently missed one of the two IRs. The resulting genome assembly contains tens to hundreds of contigs. The size of the contigs ranged from several hundred base pairs to tens of kilobase pairs. Since we know the structure of cp genome is generally stable, the syntenic relationship was used for scaffolding. We mapped the SPAdes contigs to the genome assembly resulting from the Nanopore reads using BWA MEM 45 . The alignments were used to order the contigs. The overlap information between the neighbouring contigs was used to concatenate them.

 Cp genome annotation

We used the web tool Dual Organellar GenoMe Annotator (DOGMA) 48 to generate the preliminarily gene annotations. For each particular gene, we used MUSCLE 49 (version 3.8.31) to align the genuine protein sequences of the gene gained from the NCBI GenBank to the genome assembly to decide the exact boundary positions. The web tool Organellar Genome DRAW (OGDRAW) 50 was used to generate the circular annotation plot of the genome assembly. The hypothetical cp open-reading frame ycf1 was not identified by DOGMA initially. It was added to the annotation on the basis of the MUSCLE alignment results.

 Phylogenetic analysis

Phylogenetic analysis was performed on the 18 sweetpotato cultivars used as the parental genotypes for constructions of mapping populations in GT4SP project as well as the Convolvulaceae Ipomoea section Batatas including the cp genome assemblies constructed in this research and nine publicly available cp genome assemblies. MAFFT 62 (version 7.310) was employed to perform the multiple sequence alignment (MSA) for cp genomes. The phylogenetic structure was constructed with PhyML 63 (version 3.1). Branch certainty was evaluated with 1000 replications of bootstrap resampling. The phylogenetic tree depicted in this research was constructed with the web tool iTOL (version 4) 52 .

 Data availability Underlying data

Nanopore and Illumina reads and the cp genome assemblies are deposited at NCBI BioProject repository, accession number PRJNA438020: http://identifiers.org/bioproject/PRJNA438020.

 Extended data

Supplementary Figure 1. Size distribution of the Nanopore sequencing data of the total DNA. https://doi.org/10.26188/12652034.v2 64

Supplementary Table 1. Statistics of the chloroplast (cp) sequencing data. https://doi.org/10.26188/12652067.v1 65

Supplementary Table 2. List of the 30 publicly available Ipomoea chloroplast (cp) genomes in the NCBI repository. https://doi.org/10.26188/12652079.v1 66

Supplementary Table 3. Description of the parental genotypes of the Mwanga Diversity Panel (MDP). https://doi.org/10.26188/12652085.v1 67

Supplementary Table 4. Statistics of the chloroplast (cp) genome assemblies of the 18 sweetpotato cultivars and the I. triloba line NCNSP-0323. https://doi.org/10.26188/12652094.v1 68

Acknowledgements

We would like to thank Rick Miller for his valuable suggestions on the phylogenetic analysis.

 Martin W Rujan T Richly E : Evolutionary analysis of Arabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus. Proc Natl Acad Sci U S A. 2002;99(19):12246–12251. 12218172 10.1073/pnas.182432999 129430 Shaw J Lickey EB Schilling EE : Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III. Am J Bot. 2007;94(3):275–288. 21636401 10.3732/ajb.94.3.275 Daniell H Lin CS Yu M : Chloroplast genomes: diversity, evolution, and applications in genetic engineering. Genome Biol. 2016;17(1):134. 27339192 10.1186/s13059-016-1004-2 4918201 Yan L Lai X Li X : Analyses of the complete genome and gene expression of chloroplast of sweet potato [ Ipomoea batata]. PLoS One. 2015;10(4):e0124083. 25874767 10.1371/journal.pone.0124083 4398329 Sandelius AS Aronsson H : The chloroplast: interactions with the environment.(Springer Science & Business Media).2008;13. 10.1007/978-3-540-68696-5 Shinozaki K Ohme M Tanaka M : The complete nucleotide sequence of the tobacco chloroplast genome: its gene organization and expression. EMBO J. 1986;5(9):2043–2049. 16453699 10.1002/j.1460-2075.1986.tb04464.x 1167080 Bausher MG Singh ND Lee SB : The complete chloroplast genome sequence of Citrus sinensis (L.) Osbeck var 'Ridge Pineapple': organization and phylogenetic relationships to other angiosperms. BMC Plant Biol. 2006;6:21. 17010212 10.1186/1471-2229-6-21 1599732 Samson N Bausher MG Lee SB : The complete nucleotide sequence of the coffee ( Coffea arabica L.) chloroplast genome: organization and implications for biotechnology and phylogenetic relationships amongst angiosperms. Plant Biotechnol J. 2007;5(2):339–353. 17309688 10.1111/j.1467-7652.2007.00245.x 3473179 Saski C Lee SB Fjellheim S : Complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera, and comparative analyses with other grass genomes. Theor Appl Genet. 2007;115(4):591. 27520091 10.1007/s00122-007-0595-0 Daniell H Lee SB Grevich J : Complete chloroplast genome sequences of Solanum bulbocastanum, Solanum lycopersicum and comparative analyses with other Solanaceae genomes. Theor Appl Genet. 2006;112(8):1503–18. 16575560 10.1007/s00122-006-0254-x Daniell H Wurdack KJ Kanagaraj A : The complete nucleotide sequence of the cassava ( Manihot esculenta) chloroplast genome and the evolution of atpF in Malpighiales: RNA editing and multiple losses of a group II intron. Theor Appl Genet. 2008;116(5):723–37. 18214421 10.1007/s00122-007-0706-y 2587239 Saski C Lee SB Daniell H : Complete chloroplast genome sequence of Gycine max and comparative analyses with other legume genomes. Plant Mol Biol. 2005;59(2):309–322. 16247559 10.1007/s11103-005-8882-0 Moore MJ Dhingra A Soltis PS : Rapid and accurate pyrosequencing of angiosperm plastid genomes. BMC Plant Biol. 2006;6:17. 16934154 10.1186/1471-2229-6-17 1564139 Wu J Liu B Cheng F : Sequencing of chloroplast genome using whole cellular DNA and solexa sequencing technology. Front Plant Sci. 2012;3:243. 23162558 10.3389/fpls.2012.00243 3492724 Cronn R Liston A Parks M : Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology. Nucleic Acids Res. 2008;36(19):e122. 18753151 10.1093/nar/gkn502 2577356 Eserman LA Tiley GP Jarret RL : Phylogenetics and diversification of morning glories (tribe Ipomoeeae, Convolvulaceae) based on whole plastome sequences. Am J Bot. 2014;101(1):92–103. 24375828 10.3732/ajb.1300207 Koren S Walenz BP Berlin K : Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017;27(5):722–736. 28298431 10.1101/gr.215087.116 5411767 Chin CS Alexander DH Marks P : Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10(6):563–569. 23644548 10.1038/nmeth.2474 Walker BJ Abeel T Shea T : Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014;9(11):e112963. 25409509 10.1371/journal.pone.0112963 4237348 Goodwin S Gurtowski J Ethe-Sayers S : Oxford Nanopore sequencing, hybrid error correction, and de novo assembly of a eukaryotic genome. Genome Res. 2015;25(11):1750–1756. 26447147 10.1101/gr.191395.115 4617970 Soorni A Haak D Zaitlin D : Organelle_PBA, a pipeline for assembling chloroplast and mitochondrial genomes from PacBio DNA sequencing data. BMC Genomics. 2017;18(1):49. 28061749 10.1186/s12864-016-3412-9 5219736 Wang W Messing J : High-throughput sequencing of three Lemnoideae (duckweeds) chloroplast genomes from total DNA. PLoS One. 2011;6(9):e24670. 21931804 10.1371/journal.pone.0024670 3170387 Kim K Lee SC Lee J : Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species. Sci Rep. 2015;5: 15655. 26506948 10.1038/srep15655 4623524 Antipov D Hartwick N Shen M : plasmidspades: assembling plasmids from whole genome sequencing data. bioRxiv. 2016;048942. 10.1101/048942 Zhang T Zhang X Hu S : An efficient procedure for plant organellar genome assembly, based on whole genome data from the 454 GS FLX sequencing platform. Plant Methods. 2011;7:38. 22126655 10.1186/1746-4811-7-38 3248859 Dierckxsens N Mardulyn P Smits G : NOVOPlasty: de novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017;45(4):e18. 28204566 10.1093/nar/gkw955 5389512 Izan S Esselink D Visser RGF : De Novo Assembly of Complete Chloroplast Genomes from Non-model Species Based on a K-mer Frequency-Based Selection of Chloroplast Reads from Total DNA Sequences. Front Plant Sci. 2017;8:1271. 28824658 10.3389/fpls.2017.01271 5539191 Loebenstein G Thottappilly G : The sweetpotato. (Springer Science & Business Media).2009. 10.1007/978-1-4020-9475-0 Faostat F : Agriculture organization of the united nations statistics division (2014). [Review date: April 2015].2016. Reference Source Mwanga ROM Andrade MI Carey EE : Sweetpotato ( Ipomoea batatas L.). (Springer International Publishing, Cham).2017;181–218. 10.1007/978-3-319-59819-2_6 Roullier C Rossel G Tay D : Combining chloroplast and nuclear microsatellites to investigate origin and dispersal of New World sweet potato landraces. Mol Ecol. 2011;20(19):3963–3977. 21880085 10.1111/j.1365-294X.2011.05229.x Huang J Sun M : Genetic diversity and relationships of sweetpotato and its wild relatives in ipomoea series batatas (convolvulaceae) as revealed by inter-simple sequence repeat (ISSR) and restriction analysis of chloroplast DNA. Theor Appl Genet. 2000;100(7):1050–1060. 10.1007/s001220051386 Hoshino A Jayakumar V Nitasaka E : Genome sequence and analysis of the Japanese morning glory Ipomoea nil. Nat Commun. 2016;7: 13295. 27824041 10.1038/ncomms13295 5105172 Muñoz-Rodríguez P Carruthers T Wood JR : Reconciling conflicting phylogenies in the origin of sweet potato and dispersal to Polynesia. Curr Biol. 2018;28(8):1246–1256.e12. 29657119 10.1016/j.cub.2018.03.020 Mwanga ROM Odongo B p’Obwoya CO : Release of five sweetpotato cultivars in uganda. HortScience. 2001;36(2):385–386. 10.21273/HORTSCI.36.2.385 McNeal JR Kuehl JV Boore JL : Complete plastid genome sequences suggest strong selection for retention of photosynthetic genes in the parasitic plant genus Cuscuta. BMC Plant Biol. 2007;7:57. 17956636 10.1186/1471-2229-7-57 2216012 Mwanga ROM Odongo B Niringiye C : Release of two orange-fleshed sweetpotato cultivars, ‘spk004’ (‘kakamega’) and ‘ejumula’, in uganda. HortScience. 2007;42(7):1728–1730. 10.21273/HORTSCI.42.7.1728 Mwanga ROM Odongo B Turyamureeba G : Release of six sweetpotato cultivars (‘naspot 1’ to ‘naspot 6’) in uganda. HortScience. 2003;38(3):475–476. 10.21273/HORTSCI.38.3.475 Mwanga ROM Odongo B Niringiye C : ‘naspot 7’, ‘naspot 8’, ‘naspot 9 o’, ‘naspot 10 o’, and ‘dimbuka-bukulula’ sweetpotato. HortScience. 2009;44(3):828–832. Reference Source Mwanga ROM Niringiye C Alajo A : ‘naspot 11’, a sweetpotato cultivar bred by a participatory plant-breeding approach in uganda. HortScience. 2011;46(2):317–321. 10.21273/HORTSCI.46.2.317 Mwanga ROM Kyalo G Ssemakula GN : ‘naspot 12 o’ and ‘naspot 13 o’ sweetpotato. HortScience. 2016;51(3):291–295. 10.21273/HORTSCI.51.3.291 Jones A Dukes PD Schalk JM : ‘resisto’ sweet potato. HortScience. 1983;18(2):251–252. Reference Source Gruneberg WJ Ma D Mwanga ROM : Advances in sweetpotato breeding from 1992 to 2012.(CABI International).2015. 10.1079/9781780644202.0003 Abidin PE van Eeuwijk FA Stam P : Adaptation and stability analysis of sweet potato varieties for low-input systems in uganda. Plant Breed. 2005;124(5):491–497. 10.1111/j.1439-0523.2005.01149.x Li H Durbin R : Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–1760. 19451168 10.1093/bioinformatics/btp324 2705234 Sommer DD Delcher AL Salzberg SL : Minimus: a fast, lightweight genome assembler. BMC Bioinformatics. 2007;8:64. 17324286 10.1186/1471-2105-8-64 1821043 David MC Diaz FC Mwanga ROM : Gene pool subdivision of east african sweetpotato parental material. Crop Science. 2018;58(6):2302–2314. 10.2135/cropsci2017.11.0695 Wyman SK Jansen RK Boore JL : Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004;20(17):3252–3255. 15180927 10.1093/bioinformatics/bth352 Edgar RC : MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–1797. 15034147 10.1093/nar/gkh340 390337 Lohse M Drechsel O Bock R : OrganellarGenomeDRAW (OGDRAW): a tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr Genet. 2007;52(5–6):267–274. 17957369 10.1007/s00294-007-0161-y Li Q Li Y Song J : High-accuracy de novo assembly and SNP detection of chloroplast genomes using a SMRT circular consensus sequencing strategy. New Phytol. 2014;204(4):1041–1049. 25103547 10.1111/nph.12966 Letunic I Bork P : Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–W245. 27095192 10.1093/nar/gkw290 4987883 Bankevich A Nurk S Antipov D : SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–477. 22506599 10.1089/cmb.2012.0021 3342519 Tang J Xia H Cao M : A comparison of rice chloroplast genomes. Plant Physiol. 2004;135(1):412–420. 15122023 10.1104/pp.103.031245 429394 Gogniashvili M Naskidashvili P Bedoshvili D : Complete chloroplast DNA sequences of zanduri wheat ( Triticum spp.). Genet Resour Crop Evol. 2015;62(8):1269–1277. 10.1007/s10722-015-0230-x Strittmauer G Kössel H : Cotranscription and processing of 23S, 4.5S and 5S rRNA in chloroplasts from Zea mays. Nucleic Acids Res. 1984;12(20):7633–7647. 6093045 10.1093/nar/12.20.7633 320190 Srisuwan S Sihachakr D Siljak-Yakovlev S : The origin and evolution of sweet potato ( Ipomoea batatas Lam.) and its wild relatives through the cytogenetic approaches. Plant Sci. 2006;171(3):424–433. 22980213 10.1016/j.plantsci.2006.05.007 Dellaporta SL Wood J Hicks JB : A plant DNA minipreparation: version ii. Plant Mol Biol Report. 1983;1(4):19–21. 10.1007/BF02712670 Mace ES Buhariwalla KK Buhariwalla HK : A high-throughput DNA extraction protocol for tropical molecular breeding programs. Plant Mol Biol Report. 2003;21(4):459–460. 10.1007/BF02772596 Bolger AM Lohse M Usadel B : Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–2120. 24695404 10.1093/bioinformatics/btu170 4103590 Kurtz S Phillippy A Delcher AL : Versatile and open software for comparing large genomes. Genome Biol. 2004;5(2):R12. 14759262 10.1186/gb-2004-5-2-r12 395750 Katoh K Standley DM : MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–780. 23329690 10.1093/molbev/mst010 3603318 Guindon S Dufayard JF Lefort V : New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–321. 20525638 10.1093/sysbio/syq010 Zhou C : Nanopore read length distribution. University of Melbourne. Figure.2020. http://www.doi.org/10.26188/12652034.v2 Zhou C : Statistics of the chloroplast genome sequencing data. University of Melbourne. Dataset.2020. http://www.doi.org/10.26188/12652067.v1 Zhou C : List of the 30 publicly available Ipomoea chloroplast genomes in the NCBI repository. University of Melbourne. Dataset.2020. http://www.doi.org/10.26188/12652079.v1 Zhou C : Description of the parental genotypes of the Mwanga Diversity Panel (MDP). University of Melbourne. Dataset.2020. http://www.doi.org/10.26188/12652085.v1 Zhou C : Statistics of the chloroplast genome assemblies of the 18 sweetpotato cultivars and the I. triloba line NCNSP-0323. University of Melbourne. Dataset.2020. http://www.doi.org/10.26188/12652094.v1 10.21956/gatesopenres.14357.r29627 Reviewer response for version 2 Li Zongyun 1 Referee Institute of Integrative Plant Biology, School of Life Science, Jiangsu Normal University, Xuzhou, Jiangsu 221116, China; Jiangsu Key Laboratory of Phylogenomics and Comparative Genomics, Jiangsu Normal University, Jiangsu, China

Competing interests: No competing interests were disclosed.

 7 10 2020 Copyright: © 2020 Li Z 2020 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

Summary of the key results:

The authors sequenced and assembled 19 cp genomes of sweetpotato cultivars and one wild species using Oxford Nanopore and Illumina sequencing. With the published data, the authors constructed a phylogeny tree and proposed that I. trifida is the most closely related diploid species of the sweetpotato. In addition, 18 sweetpotato cp genomes demonstrated the presence of the two distinct subpopulations in East Africa.

Overall evaluation:

This manuscript provided more information for sweetpotato genome resources which is very important for us to learn more about the genetic diversity, origin, and evolution of sweetpotato. According to the contents of results and discussion, some comments were listed as follows:

Park et al. (2018 1 ) and Sun et al. (2019 2 ) also published some Ipomoea cp genomes. It would be better to cite their works.

There are more genes annotated in this article than in other studies (Eserman et al., 2014 3 ; Yan et al., 2015 4 ; Park et al., 2018 1 ; Sun et al., 2019 2 ). Can the authors discuss the reasons? Which genes were not annotated in previous studies?

I thought one of the highlights of this paper was the assembly of 18 sweetpotato cp genomes. The authors demonstrated the presence of the two distinct subpopulations in East Africa using these cp genomes. However, no other more detailed analysis and discussion about these 18 cp genomes. Would it be better to add more detailed sequences analysis? For example, the factors which impact the genome size, some specific loci between the two subpopulation, etc.

  

Overall, this manuscript is valuable for indexing.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

References : The Complete Chloroplast Genomes of Six Ipomoea Species and Indel Marker Development for the Discrimination of Authentic Pharbitidis Semen (Seeds of I. nil or I. purpurea). Front Plant Sci .2018;9: 10.3389/fpls.2018.00965 965 30026751 10.3389/fpls.2018.00965 : A systematic comparison of eight new plastome sequences from Ipomoea L. PeerJ .2019;7: 10.7717/peerj.6563 e6563 30881765 10.7717/peerj.6563 : Phylogenetics and diversification of morning glories (tribe Ipomoeeae, Convolvulaceae) based on whole plastome sequences. Am J Bot .2014;101(1) : 10.3732/ajb.1300207 92-103 24375828 10.3732/ajb.1300207 : Analyses of the complete genome and gene expression of chloroplast of sweet potato [Ipomoea batata]. PLoS One .2015;10(4) : 10.1371/journal.pone.0124083 e0124083 25874767 10.1371/journal.pone.0124083 10.21956/gatesopenres.14357.r29625 Reviewer response for version 2 Li Qiang 1 Referee https://orcid.org/0000-0003-2923-0136 Xuzhou Institute of Agricultural Sciences in Jiangsu Xuhuai District, Key Laboratory for Biology and Genetic Breeding of Sweetpotato (Xuzhou), Ministry of Agriculture, Sweet Potato Research Institute (SPRI), Chinese Academy of Agricultural Sciences, Xuzhou, China

Competing interests: No competing interests were disclosed.

 7 10 2020 Copyright: © 2020 Li Q 2020 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve-with-reservations

Overall evaluation:

I think that the manuscript is acceptable. It gives us a new method to insight into the population structure of plants, especially sweetpotato with complicated genome. However, I have a number of suggestions to the article.

 

Materials:

 

Just the author’s response to Prof. Aureliano Bombarely, "the primary focus of this study was to investigate the population structure of the East African sweetpotato cultivars used in the GT4SP project", 16 cultivars for MDP and 2 cultivars (Tanzania and Beauregard) for constructing F1 mapping population are suitable for this research. The I. triloba was selected. You can give us reasons why you chose this wild species, not I. trifida, or other wild relatives.

The authors assembled cp genomes 2 times of cultivar NASPOT 10 O - one was from the screen-house while the other one was from the field. What was the aim? You can give us more detail about the difference or consistency between them.

There are total 20 samples, including 19 samples from 18 sweetpotato cultivars (2 samples of the cultivar NASPOT 10 O), and 1 sample of wild relative. I always confuse the number in the article, cultivars or samples.

Methods:

 

"In the present study, we constructed a complete cp genome assembly for the hexaploid sweetpotato cultivar Tanzania using long reads produced by the Oxford Nanopore sequencing technology."  Is this a new method? Is it the first report? It is an efficient tool to assemble complicated genomes in my mind. It is an important part of this article. I suggest the method should be reflected in the title.

 Give more detail in the discussion part about this new tool compared to other methods.

Results:

 

"The sweetpotato cp genome of 161,274 bp contains 152 genes, of which there are 96 protein coding genes, 8 rRNA genes and 48 tRNAgenes..."  this is cp genome of Tanzania. There are a little difference among other cultivars.

Suggest to compare the phylogenetic tree by using cp genome data and nuclear genome data, and to validate the method.

Others:

 

"The only difference is that the psbZ gene is not found in the cultivar Xushu 18 cpDNA while the ihbA gene is not found in the cultivar Tanzania cpDNA. It should be noted that the double-copy gene ycf1 was not reported for the cultivar Xushu 18 cp genome4, but this was actually a miss-annotation." The difference between Xushu 18 and Tanzania are psbZ and ihbA, why the ycf1 was actually a miss-annotation. You should give more information about it.

Just a suggestion from Dr. Yang, the important references should be added.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

 10.21956/gatesopenres.14357.r29622 Reviewer response for version 2 Monden Yuki 1 Referee https://orcid.org/0000-0001-9041-1118 Graduate School of Environmental and Life Science, Okayama University, Okayama, Japan

Competing interests: No competing interests were disclosed.

 5 10 2020 Copyright: © 2020 Monden Y 2020 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

I think this paper is well-written, and the authors have properly revised the manuscript by referring to the comments of the reviewer. This paper performed a complete circular cp genome assembly using NGS and TGS technologies in the hexaploid sweetpotato. Phylogenetic analysis using the cp genomes revealed that there are two distinct clusters of sweetpotato in East Africa and I. trifida is the most closely related diploid wild species to the I. batatas hexaploid sweetpotato. The results of this paper provide insights into the genetic relationships and the population structure of the species from the Convolvulaceae Ipomoea section Batatas. Besides, despite the complexity of the cp genomes by the presence of two large inverted repeats, this research demonstrates the possibility of building the cp genomes using extremely low coverage (<1x) Oxford Nanopore WGS data combined with Illumina short reads. Other comments are shown below.

 

Table 1:

The copy number of rps12 gene should be three, but the bracketed superscript of this gene is two. Please make sure.

Phylogenetic analysis using nuclear genomes

Is it possible to compare the results of phylogenetic analyses based on the cp genomes and the nuclear genomes using the same materials? I think such comparative analysis should provide new insight into evolutionary dynamics on cp and nuclear genomes of Ipomoea species. Do you have any plans for such work?

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

Plant genetics

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

 10.21956/gatesopenres.13938.r26743 Reviewer response for version 1 Bombarely Aureliano 1 Referee https://orcid.org/0000-0001-6257-8914 Department of Horticulture, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA

Competing interests: No competing interests were disclosed.

 20 11 2018 Copyright: © 2018 Bombarely A 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve-with-reservations

Summary of the key results

The manuscript titled “Insights into population structure of East African sweetpotato cultivars from hybrid assembly of chloroplast genomes” presents the chloroplast genome assembly and annotation of the sweet potato cultivar Tanzania and its comparison with seventeen other cultivars and the different species Ipomoea triloba. The final genome assembly of the Tanzania cultivar was 161,274 bp. No major findings were reported except for the deletion of the gene psbZ in the cultivar Xushu 18 and ihbA gene in the Tanzania. The phylogenetic analysis of these cultivars and nine publicly available Ipomoea chloroplast pointed that the diploid I. trifida species is close related to the hexaploid sweet potato ( I. batatas) than other Ipomoea species.

Overall evaluation

This manuscript presents the assembly and analysis of the chloroplast genome of the sweet potato cultivar Tanzania. The I. batatas chloroplast genome was already published in 2015 by Yan et. al. (PLoS One 10:4) so the novelty of the results presented in this manuscript are limited from the point of new of a “new” chloroplast genome. The use of ON for the sequencing of the Tanzania cultivar and the addition of the resequencing data of seventeen other cultivars potentially could add some interesting findings. Nevertheless, the analysis that the authors performed failed in the development of attractive results. Personally, I would propose several analysis that it may help to increase the impact of the manuscript:

Improved phylogenetic analysis partitioning the alignments per gene (or in bins) and using a Bayesian framework (e.g. BEAST). The phylogenetic analysis could include the dating of the divergency time of the different taxa.

Nuclear gene mining. One of the most interesting questions about the polyploids is the origin of those. The use of the resequencing data could potentially derived in the mining of nuclear copy single gene that could help to elucidate a different evolutionary trajectory for these accessions. Additionally, it is interesting the result in which some genes are missing. Maybe they have been transferred to the nuclear genome. It could be interesting that hypothesis.

Positive selection. Each of the chloroplast nuclear genes could be tested for positive selection using PAML and the Ks/Kn ratio.

  

In terms of the manuscript organization and writing, I found confusing some parts. For example, the material and methods are not aligned with the results presented in the manuscript. For example, the section “Extraction of cp genome sequence from whole genome sequencing data” describe the chloroplast data mining from ON and Illumina for the Tanzania accession and then for the Beauregard accession, but the material and methods also describe the use of 10X Genomics Chromium that I am not sure where it comes from. Do the authors used 10X Genomics also? Probably for the genome assembly, the comparison of the Canu assembly with Organelle_PBA (Soorni et al. 2017) could be interesting, to see if the authors obtain only one contig representing the whole chloroplast.

Overall, I think that the manuscript is okay, but there are some space for improvement in the structure of the manuscript and in the results that are presented.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

Coin Lachlan

Competing interests: Not applicable

 12 7 2020

We thank the reviewer for their thoughtful review of  our manuscript. 

Reviewers comment 1. This manuscript presents the assembly and analysis of the chloroplast genome of the sweet potato cultivar Tanzania. The I. batatas chloroplast genome was already published in 2015 by Yan et. al. ( PLoS One 10:4) so the novelty of the results presented in this manuscript are limited from the point of new of a “new” chloroplast genome. The use of ON for the sequencing of the Tanzania cultivar and the addition of the resequencing data of seventeen other cultivars potentially could add some interesting findings. Nevertheless, the analysis that the authors performed failed in the development of attractive results. Personally, I would propose several analysis that it may help to increase the impact of the manuscript:

Improved phylogenetic analysis partitioning the alignments per gene (or in bins) and using a Bayesian framework (e.g. BEAST). The phylogenetic analysis could include the dating of the divergency time of the different taxa.

Nuclear gene mining. One of the most interesting questions about the polyploids is the origin of those. The use of the resequencing data could potentially derived in the mining of nuclear copy single gene that could help to elucidate a different evolutionary trajectory for these accessions. Additionally, it is interesting the result in which some genes are missing. Maybe they have been transferred to the nuclear genome. It could be interesting that hypothesis.

Positive selection. Each of the chloroplast nuclear genes could be tested for positive selection using PAML and the Ks/Kn ratio.

 Response: The primary focus of this study was to investigate the population structure of the East African sweetpotato cultivars used in the GT4SP project. We strongly agree that these suggested analyses will largely increase the impact of the manuscript. However, it is difficult to integrate these suggested analyses in the current study within the constraints of the data. We will include these suggestions in the future directions of the project.  

Reviewers Comment:  In terms of the manuscript organization and writing, I found confusing some parts. For example, the material and methods are not aligned with the results presented in the manuscript. For example, the section “Extraction of cp genome sequence from whole genome sequencing data” describe the chloroplast data mining from ON and Illumina for the Tanzania accession and then for the Beauregard accession, but the material and methods also describe the use of 10X Genomics Chromium that I am not sure where it comes from. Do the authors used 10X Genomics also? Probably for the genome assembly, the comparison of the Canu assembly with Organelle_PBA (Soorni et al. 2017) could be interesting, to see if the authors obtain only one contig representing the whole chloroplast.

Response: We have amended the mansucript to adddress these concerns.   10X Genomics was indeed used to perform the whole genome sequencing for the sweetpotato accessions Tanzania and Beauregard. However, the linked reads information was not utilized in construction of the cp genome assemblies for them. Instead, the sequence data was simply used as paired-end reads to create contigs. The contigs were then used to construct whole cp genome assembly with a cp reference genome.

 10.21956/gatesopenres.13938.r26651 Reviewer response for version 1 Yang Jun 1 Referee https://orcid.org/0000-0002-0371-8814 Yan Mengxiao 2 Co-referee State Key Laboratory of Plant Molecular Genetics, CAS Center for Excellence in Molecular Plant Sciences, Shanghai Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China Shanghai Key Laboratory of Plant Functional Genomics and Resources, Shanghai Chenshan Plant Science Research Center (SCPSRC), Shanghai Chenshan Botanical Garden, Chinese Academy of Sciences, Shanghai, China

Competing interests: No competing interests were disclosed.

 13 9 2018 Copyright: © 2018 Yang J and Yan M 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve-with-reservations

Zhou et al. sequenced 16 sweet potato cultivars in GT4SP project supported by Bill & Melinda Gates Foundation. Here authors presented only partial data about chloroplast (cp) genome assemblies.

I found this study will be a good supplementary work of previous publication in Current Biology 1 . Due to the lack of awareness of this publication, the claims in the manuscript are incorrect and need to be revised.

In this case, the finding of this study about two distinct cp subpopulations in East Africa cultivars is reasonable.

Other comments

“In recent years, the development of next-generation sequencing (NGS) technologies such as Illumina and Roche 454 facilitate faster and cheaper methods to sequence cp genomes 13–15.”

To my knowledge, the Roche 454 already left the market.

“By examining the mapping results of the WGS data, we are able to extract cp sequences 21,22.”

We? Who are we?

“Sweetpotato is a hexaploid (2n=6x=90) with genome size estimated to be between 2,200 to 3,000 Mb 28.”

How about the C-values?

“Due to the complex genome structure, the availability of sweetpotato genomic resources is lacking.”

We do have a published genome, right?

“A number of cp genomes from the Ipomoea family have been sequenced 16,33.”

Dose Ipomoea family mean genus Ipomoea? Or genus Ipomoea Series Batatas?

“Most of them are diploid wild relatives of the sweetpotato. To the best of our knowledge, to date, four cp genomes have been completely sequenced for the hexaploid sweetpotato 4,16; the genome size is around 161 Kb, and the structure represents a standard quadripartite circular with a LSC of 87 Kb, a SSC of 12 Kb and two IRs of 31 Kb 4. The cp genomes were mainly used to perform phylogenetic analyses 4,16.”

The mentioned Current Biology paper has provided hundreds of cp genome sequences of sweet potato and its wild relatives.

“The circularized contig is ~161 Kb, and is highly collinear with the published sweetpotato cp genome assembly (Figure 1d).”

In Figure 1d, it is an I. trifida cp genome, not a published sweet potato cp genome.

“The sweetpotato cp genome represents a common circular structure with two IRs (IRA and IRB) separating one LSC and one SSC2.”

Where does the 2 in SSC2 come from? Convert into right format if it is a citation.

“The red dots represent SNPs between the two cp genomes. The green bars on the x-axis indicate positions of the two IRs”

No red dots there, only black dots.

“It should be noted that the doublecopy gene ycf1 was not reported for the cultivar Xushu 18 cp genome4”

Convert into right format if it is a citation.

“Interestingly, the accession PI 618966 was originally identified as I. triloba and was recently reidentified as I. trifida by the GRIN National Genetic Resources Program.”

The identification of PI 618966 needs to be checked carefully. All individuals of I. trifida formed a monophyletic clade closely related to I. batatas according to Current Biology paper. As the progenitor of sweet potato, it's quite strange that I. trifida is much closer to other species in Series Batatas than I. batatas.

Figure 3 & 4

It will be much clear to add the tip labels rather than collapsed clades on the tree. Figure 4 will be no more informative in this case.

If the tree is not that complicated, it is not suggested to collapse the two clades. Since information about the relationship between within-clade sample and out-clade sample is not visible when one collapse clade. This information will not be illustrated in Figure 4. Clades can be labeled in different colors if one wants to highlight the clades. 

Furthermore, it is not clear to me which place each sample nested on in Figure 4.

“In this study, we used the genome assembly constructed from the Nanopore reads as reference to assemble cp genomes for a further 19 cp genomes including…”

Misleading sentence, authors do rely on published cp genome rather than de novo Nanopore assembly.

“In order to solve this problem, for reference-based scaffolding, we reused some single-copy contigs from the two IR regions to construct complete cp genome assemblies.”

In which cultivar(s), did author investigate the influence on the tree structure?

I agree the population structure of East African sweet potato cultivars is important for GT4SP project. Also obviously, the data organization and visualization could be largely improved to meet the indexing standards.

Is the work clearly and accurately presented and does it cite the current literature?

Partly

If applicable, is the statistical analysis and its interpretation appropriate?

Yes

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

Plant genetics

We confirm that we have read this submission and believe that we have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however we have significant reservations, as outlined above.

References : Reconciling Conflicting Phylogenies in the Origin of Sweet Potato and Dispersal to Polynesia. Curr Biol .2018;28(8) : 10.1016/j.cub.2018.03.020 1246-1256.e12 29657119 10.1016/j.cub.2018.03.020 Coin Lachlan

Competing interests: Not applicable

 12 7 2020

We thank the reviewer for their thoughtful review . 

Comment 1. I found this study will be a good supplementary work of previous publication in Current Biology (Mu, Pablo, et al., 2018). Due to the lack of awareness of this publication, the claims in the manuscript are incorrect and need to be revised.

Response: This publication was cited in the revised version. The incorrect claim was revised (see reviewer’s comment 7).

Comment 2. “In recent years, the development of next-generation sequencing (NGS) technologies such as Illumina and Roche 454 facilitate faster and cheaper methods to sequence cp genomes 13–15.”

To my knowledge, the Roche 454 already left the market.

Response: In the revised version, “In recent years” was deleted to make it more precise.

Comment 3. “By examining the mapping results of the WGS data, we are able to extract cp sequences 21,22.”

We? Who are we?

Response: In the revised version, this sentence was rewritten to “The cp sequences could be extracted by examining the mapping results of the WGS data to the reference cp genome 21, 22.”

Comment 4. “Sweetpotato is a hexaploid (2n=6x=90) with genome size estimated to be between 2,200 to 3,000 Mb 28.”

How about the C-values?

Response: The nuclear genome size is not the key point of this study. The C-value was not investigated.

Comment 5. “Due to the complex genome structure, the availability of sweetpotato genomic resources is lacking.”

We do have a published genome, right?

Response: Even though there is a sweetpotato reference genome published recently (Yang, Jun, et al., 2017), we think the availability of the sweetpotato genome resources is still lacking.  

Comment 6. “A number of cp genomes from the Ipomoea family have been sequenced 16,33.”

Dose Ipomoea family mean genus Ipomoea? Or genus Ipomoea Series Batatas?

Response: “ Ipomoea family” means “genus Ipomoea”. This was made clear in the revised version.

Comment 7. “Most of them are diploid wild relatives of the sweetpotato. To the best of our knowledge, to date, four cp genomes have been completely sequenced for the hexaploid sweetpotato 4, 16; the genome size is around 161 Kb, and the structure represents a standard quadripartite circular with a LSC of 87 Kb, a SSC of 12 Kb and two IRs of 31 Kb 4. The cp genomes were mainly used to perform phylogenetic analyses 4, 16. ”

The mentioned Current Biology paper has provided hundreds of cp genome sequences of sweet potato and its wild relatives.

Response: The claim “To the best of our knowledge, to date, four cp genomes have been completely sequenced for the hexaploid sweetpotato 4, 16;” was removed in the revised version.

Comment 8. “The circularized contig is ~161 Kb, and is highly collinear with the published sweetpotato cp genome assembly (Figure 1d).”

In Figure 1d, it is an I. trifida cp genome, not a published sweet potato cp genome.

Response: This sentence was corrected as “The circularized contig is ~161 Kb, and is highly collinear with the reference cp genome assembly (Figure 1d).”

Comment 9. “The sweetpotato cp genome represents a common circular structure with two IRs (IRA and IRB) separating one LSC and one SSC2.”

Where does the 2 in SSC2 come from? Convert into right format if it is a citation.

Response: This was a citation. It was corrected in the revised version.

Comment 10. “The red dots represent SNPs between the two cp genomes. The green bars on the x-axis indicate positions of the two IRs”

No red dots there, only black dots.

Response: The red dots are hard to see due to the resolution of the image. Since the SNPs are not important in this genome assembly section and are further discussed in the phylogenetic analysis section, the SNPs (red dots) were removed from the figure 1d in the revised version.

Comment 11. “It should be noted that the doublecopy gene ycf1 was not reported for the cultivar Xushu 18 cp genome4”

Convert into right format if it is a citation.

Response: This was a citation. It was corrected in the revised version.

Comment 12. “Interestingly, the accession PI 618966 was originally identified as I. triloba and was recently reidentified as I. trifida by the GRIN National Genetic Resources Program.”

The identification of PI 618966 needs to be checked carefully. All individuals of I. trifida formed a monophyletic clade closely related to I. batatas according to Current Biology paper. As the progenitor of sweet potato, it's quite strange that I. trifida is much closer to other species in Series Batatas than I. batatas.

Response: We agree with the reviewer that the identification of PI 618966 needs to be checked carefully. According to the phylogenetic structure identified in this study (Fig. 3), PI 618966 has a closer phylogenetic relationship to I. triloba instead of I. trifida. However, it was recently reidentified as I. trifida by the GRIN National Genetic Resources Program. A further study is required to fully clarify this.

Comment 13. Figure 3 & 4

It will be much clear to add the tip labels rather than collapsed clades on the tree. Figure 4 will be no more informative in this case.

If the tree is not that complicated, it is not suggested to collapse the two clades. Since information about the relationship between within-clade sample and out-clade sample is not visible when one collapse clade. This information will not be illustrated in Figure 4. Clades can be labeled in different colors if one wants to highlight the clades.

Furthermore, it is not clear to me which place each sample nested on in Figure 4.

Response: The two sweetpotato clades showed in Figure 3 were collapsed since the phylogenetic distances are too small. It will be impossible to see the detail phylogenetic structure of the two East African sweetpotato subpopulations if incorporate Figure 4 into Figure 3.

Comment 14. “In this study, we used the genome assembly constructed from the Nanopore reads as reference to assemble cp genomes for a further 19 cp genomes including…”

Misleading sentence, authors do rely on published cp genome rather than  de novo Nanopore assembly.

Response: The cp genome assembly of the sweetpotato cultivar Tanzania was constructed from the Nanopore reads coupled with a published cp reference genome. The cp genome assembly of the sweetpotato cultivar Tanzania from the Nanopore reads was then used as reference to constructed cp genome assemblies for a further 19 cp genomes.

Comment 15. “In order to solve this problem, for reference-based scaffolding, we reused some single-copy contigs from the two IR regions to construct complete cp genome assemblies.”

In which cultivar(s), did author investigate the influence on the tree structure?

Response: Single-copy contigs were reused in construction of cp genome assembly for all cultivars. The genome assembler SPAdes collapsed the contigs from the two IR regions as they are almost identical. In order to construct the whole cp genome, the contigs from two IR regions need to be reused. This has no influence on the tree structure.