Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing
Genome Biology volume 23, Article number: 258 (2022)
Advancing crop genomics requires efficient genetic systems enabled by high-quality personalized genome assemblies. Here, we introduce RagTag, a toolset for automating assembly scaffolding and patching, and we establish chromosome-scale reference genomes for the widely used tomato genotype M82 along with Sweet-100, a new rapid-cycling genotype that we developed to accelerate functional genomics and genome editing in tomato. This work outlines strategies to rapidly expand genetic systems and genomic resources in other plant species.
Recent advances in genome sequencing and editing enable the interrogation and manipulation of crop genomes with unprecedented accuracy. Pan-genomes can capture diverse alleles within crop species but studying their phenotypic consequences is limited by efficient functional genetic systems in relevant and diverse genotypes. Tomato is a model system for dissecting the genetics of crop domestication and quantitative traits. Resequencing hundreds of tomato genomes has uncovered vast genomic diversity [1, 2] but chromosome-scale genomes are only available for a few accessions [3,4,5] and there is a historical discrepancy between the reference genome (Heinz 1706) and genotypes that are commonly used for genetic and molecular experimentation (e.g., cultivars M82, Moneymaker, and Ailsa Craig). Furthermore, large-scale genetics experiments in genotypes with larger fruits are limited by the labor-intensive phenotyping and the requirement for extensive growth facilities to accommodate large plants with long generation times. The ultra-dwarfed variety Micro-tom overcomes some of these limitations , but its highly mutagenized background, severe hormonal and developmental abnormalities, and low fruit quality undermine its value for studying many phenotypes of translational agronomic importance, such as shoot, inflorescence, and fruit development (Fig. 1a and Additional file 1: Fig. S1a-f). As a consequence, the large-fruited cultivar M82 remains a primary reference for genetic, metabolic, and developmental analyses [7, 8] but the lack of a high-quality genome assembly causes reference bias and false signals in genomics analyses.
Genome assemblies are typically built from PacBio High Fidelity (HiFi) and/or Oxford Nanopore long-reads (ONT) . HiFi reads average 15 kbp in length, are highly accurate (~0.1% error), and can produce contiguous draft genome assemblies . However, HiFi-based assemblies often fragment at large and homogenous repeats as well as known sequence-specific coverage dropouts . Built from much longer, though noisier reads with a distinct error profile, ONT-based assemblies can complement HiFi-based assemblies by resolving some larger repeats or compensating for HiFi coverage dropouts . Combining complementary genome assemblies is sometimes referred to as “reconciliation” and several tools have been developed for this task [12, 13]. Even after combining complementary assemblies, modern draft genome assemblies rarely achieve complete chromosome scale. Longer and ultimately chromosome-scale sequences are produced by scaffolding, the process of ordering and orienting genome assembly contigs, and placing gaps between adjacent contigs. Scaffolding is usually achieved by comparing a genome assembly to genome maps encoding the relative distances of genomic markers along chromosomes. Linkage, physical (including optical maps and reference assemblies), and spatial proximity maps (from Chromatin Conformation Capture, or Hi-C data) are popular and effective for scaffolding assemblies. However, because genome maps are noisy and scaffolding methods are fallible, automated scaffolding usually results in incomplete or misassembled scaffolds and researchers often rely on laborious manual curation to correct these shortcomings [14, 15].
Results and discussion
To address the limitations for large-scale genetics experiments in tomato and illustrate how new genomic systems can be rapidly developed, we established the small-fruited tomato cultivar Sweet-100 (S100) as a new system for genome editing and functional genomics. Previously, we used clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9 genome editing to engineer mutations in the paralogous flowering repressor genes SELF PRUNING (SP) and SELF PRUNING 5G (SP5G) in S100 to induce fast flowering and compact growth . We confirmed that null mutations in these genes cause rapid cycling and compact growth without severe developmental abnormalities (Fig. 1a and Additional file 1: Fig. S1a-g). Furthermore, we quantified flowering and ripening time in S100, and found that the first ripe fruits mature in only ~65 days after sowing, allowing up to five generations per year compared to at most three generations for most large-fruited genotypes (Fig. 1a–d). The short generation time and compact growth habit of S100 allows both greenhouse and field growth at high density and reduced input. Notably, S100 yields ripe fruits ~3 weeks earlier than Micro-Tom and produces more seeds per fruit (Fig. 1b and Additional file 1: Fig. S1d,f). Together, these characteristics make S100 a highly efficient system for genetics and a valuable complement to the widely used M82 cultivar. However, a corresponding high-quality reference genome assembly is required to actualize the utility of this new S100 model system.
To help overcome the limitations of current assembly scaffolding approaches and rapidly generate reference genomes for S100 and M82, we developed RagTag, an improved version of RaGOO, to automate scaffolding and improve modern genome assemblies (Fig. 1e). RagTag supersedes our previously published RaGOO scaffolder by implementing improvements to the homology-based correction and scaffolding modules (see Methods section “RagTag overview”)  and by providing two new scaffolding tools called “patch” and “merge”. RagTag “patch” uses one genome assembly to make scaffolding joins and fill gaps in a second genome assembly (Fig. 1f). This is especially useful for genome assembly projects with complementary sequencing technology types, such as HiFi and ONT, as we demonstrate by accurately patching the CHM13 “Telomere-2-Telomere” human reference assembly (Additional files 1 and 2) (see Methods section “Patching a human genome assembly”) . However, while utilizing complementary types of sequencing data is ideal, RagTag “patch” is generic and can be applied to any pair of assemblies. RagTag “merge” is an extension of the CAMSA scaffolder that reconciles multiple candidate scaffolds for a given assembly (Fig. 1g) . This allows users to scaffold an assembly with potentially several map or map-specific technical parameters and biases and then synergistically combine results into a single scaffolding solution. The single scaffolding solution will likely have fewer false joins than the input scaffolds, which reduces the manual labor necessary to correct mis-assemblies that may be present in any of the individual scaffolding solutions. The input scaffolds are encoded within a “scaffold graph” with contigs as nodes and edges weighted by the confidence of the placement of contigs. RagTag analyzes this graph to resolve ambiguous paths, optionally using Hi-C data to re-weight the edges . RagTag can be used for assemblies of varying quality, though benchmarking with several Arabidopsis thaliana Columbia-0 draft genome assemblies indicated that higher quality input assemblies led to more accurate RagTag results and ultimately more contiguous assemblies (Additional file 2: Tables S1, S2 and Additional file 1: Fig. S2-S10).
We used RagTag to produce high-quality chromosome-scale reference genomes for M82 and S100. For each genotype, we assembled HiFi and ONT data independently. After screening the HiFi primary contigs for contaminates and organellar sequences, we used RagTag “patch” to patch the HiFi contigs with the ONT contigs, ultimately increasing the N50 from 20.1 to 40.8 Mbp and 12.6 to 27.8 Mbp in S100 and M82, respectively, without introducing any gaps (Fig. 1h, Additional file 1: Fig. S11, Additional file 2: Table S3). After patching, S100 chromosomes 1 and 5 and M82 chromosome 7 were each represented in a single chromosome-scale contig (Fig. 1i). A comparison to other patching tools showed that RagTag “patch” outperformed DENTIST  and it was competitive with SAMBA  and Quickmerge  (Additional file 2: Table S4). We incorporated a variety of physical and spatial proximity maps into RagTag “merge” to build 12 chromosome-scale pseudomolecules for each assembly (Additional file 2: Table S5). For S100, RagTag synergistically combined three distinct scaffolding solutions (two homology-based and one Hi-C based) to generate a new solution. While two distinct scaffolding solutions were used as input for M82 (homology-based and Hi-C-based), the output was identical to the homology-based input scaffolds. This agreement between the homology-based and the Hi-C-based scaffolds served as a validation of the homology-based scaffolds. Finally, each set of scaffolds was manually validated and corrected with Juicebox Assembly Tools and the assemblies were packaged according to our Pan-Sol specification (https://github.com/pan-sol/pan-sol-spec) . The final S100 and M82 assemblies had a QV score of 56.6 and 53.1, respectively, and re-mapping Hi-C data indicated broad structural accuracy for each assembly (Additional file 2: Table S3, Additional file 1: Fig. S12), demonstrating that RagTag allows fast and accurate generation of chromosome-scale assemblies with little manual intervention.
Using a read mapping approach, we previously reported that S100 and M82 are admixed due to historical breeding, and are thus structurally distinct from the Heinz 1706 reference genome . When comparing all three genomes, we confirmed elevated rates of structural variation across broad chromosomal regions, indicating introgressions from wild relatives (Fig. 1i–k). Multiple chromosomes, such as chromosomes 4, 9, 11, and 12 in S100 and chromosomes 4, 5, and 11 in M82, are nearly entirely introgressed from wild relatives, and within introgressions, we detected several large inversions. The largest inversion, a ~8.6-Mbp inversion observed on chromosome 9 of S100, was recently discovered in the wild tomato S. pimpinellifolium accession LA2093 and was genotyped in 99% of S. pimpinellifolium accessions, reinforcing the contribution of S. pimpinellifolium and other wild tomato species to the S100 and M82 genomes . Such widespread structural variation between these three tomato accessions highlights the need for personalized genomes to mitigate reference bias and false signals in genomics experiments.
Powerful experimental systems for genetics and functional genomics allow routine genetic manipulation. Using the new S100 genome assembly as a foundation, we adapted our tomato transformation and genome editing protocols to genetically modify S100 (see Methods section “Plant transformation”). We obtained transgenic plants in less than 4 months, comparable to previously published protocols (Additional file 1: Fig. S13a, b) [6, 24]. To test the efficiency of CRISPR-Cas9 genome editing in S100, we utilized our new S100 assembly for accurate guide-RNA (gRNA) design and targeted the tomato homolog of Arabidopsis APETALA3 (SlAP3, Solyc04g081000) on the chromosome 4 introgression (Fig. 2a). In Arabidopsis, AP3 activity is essential for petal and stamen development  and we observed the expected abnormal or missing petals and stamens on all seven ap3 CRISPR-Cas9 (ap3CR) first-generation (T0) transgenic plants (Fig. 2b). We identified multiple ap3CR mutant alleles by Sanger sequencing and observed germline transmission in the next generation, demonstrating efficient and robust editing in S100 (Fig. 2c, d and Additional file 1: Fig. S13c-e). We next explored the ability to delete entire gene loci, which is important to mitigate potential confounding genetic compensation responses to mutant allele transcripts . We targeted the floral identity gene ANANTHA (AN)  183 bp upstream and 18 bp downstream of the protein coding sequence (Fig. 2e). From four T0 transgenics we identified a complete 1568 bp gene deletion allele (anCR-1568), which was transmitted to the next generation (Fig. 2f–h and Additional file 1: Fig. S13f). Second-generation individuals that carried the anCR-1568 allele developed cauliflower-like inflorescences due to floral meristem overproliferation that is characteristic of an mutants. (Fig. 2h). Finally, we tested the potential for mutating multigene families and targeted the three floral regulator genes JOINTLESS2 (J2), ENHANCER OF J2 (EJ2), and LONG INFLORESCENCE (LIN), which belong to the MADS-box gene family  (Fig. 2i). From nine T0 transgenics, we identified an individual (T0-6) that displayed the j2CRej2CR double mutant phenotype with highly branched inflorescences, and one individual (T0-17) with the j2CRej2CRlinCR triple mutant null phenotype of cauliflower-like inflorescences (Fig. 2j-k). Sanger sequencing revealed j2CR and ej2CR mutant alleles in addition to wild-type LIN alleles in the T0–6 individual, while only mutant alleles for all three genes were detected in the T0–17 plant (Additional file 1: Fig. S13g). We then used multiplexed amplicon sequencing and genotyped a population of 184 T1 individuals for edits within protospacer sequences (Fig. 2l-o, see Methods section “Validation of CRISPR-Cas9 editing”). We identified multiple allelic combinations of j2CR, ej2CR and linCR mutations and found that higher read modification rates were associated with increases in inflorescence branching, supporting our previous findings that lower MADS-box gene dosage causes quantitative increases in inflorescence complexity  (Fig. 2p). Together, these results illustrate the effectiveness of S100 as an experimental platform for CRISPR-Cas9 genome editing of individual genes and multigene families for revealing complex genotype-to-phenotype relationships.
In this report, we introduced a new toolset for automated genome assembly scaffolding to generate new reference genomes and elevate a rapid-cycling tomato variety to an effective experimental system for functional genomics and genome-scale editing experiments. Our study outlines a roadmap to rapidly establish multiple personalized reference systems as cornerstones for functional interrogation of the vast genetic variation within and between species.
Plant material, growth conditions, and phenotyping
Seeds of S. lycopersicum cv. M82 (LA3475), Sweet-100 (S100), and Micro-tom (MT) were from our own stocks. Seeds were directly sown and germinated in soil in 96-cell plastic flats. Plants were grown under long-day conditions (16-h light, 8-h dark) in a greenhouse under natural light supplemented with artificial light from high-pressure sodium bulbs (~250 μmol m−2 s−1) at 25°C and 50–60% relative humidity. Seedlings were transplanted to soil to 3.5 L (S100 and MT) or 10 L (M82) pots 3–4 weeks after sowing. Analyses of fruit ripening, flower number, seed number, fruit weight, fruit sugar content (Brix), and inflorescence branching were conducted on mature plants grown in pots. Sugar content (Brix) of fruit juice was quantified using a digital refractometer (Hanna Instruments HI96811). Fruit ripening was quantified by labeling individual flowers at anthesis and counting the days to breaker fruit stage and red fruit stage. The number of replicates is indicated in figures or legends. The source data is included in an Additional file 5.
RagTag supersedes RaGOO as a homology-based genome assembly correction (RagTag “correct”) and scaffolding (RagTag “scaffold”) tool . RagTag implements several general improvements and conveniences for these features but follows the same algorithmic approach as previously reported. RagTag also provides two new tools called “patch” and “merge” for genome assembly improvement. RagTag “patch” uses one genome assembly to “patch” (join contigs and/or fill gaps) sequences in another assembly. RagTag “merge” reconciles two or more distinct scaffolding solutions for the same assembly. Finally, RagTag offers a variety of command-line utilities for calculating assembly statistics, validating AGP files, and manipulating genome assembly file formats. RagTag is open source (distributed under the MIT license) and is available on GitHub: https://github.com/malonge/RagTag.
RagTag whole-genome alignment filtering and merging
Most RagTag tools rely on pairwise (a “query” vs. a “reference/target”) whole-genome alignments. RagTag supports the use of Minimap2, Unimap, or Nucmer for whole-genome alignment, though any alignments in PAF or MUMmer’s delta format can be used [29, 30]. RagTag filters and merges whole-genome alignments to extract useful scaffolding information. To remove repetitive alignments, RagTag uses an integrated version of unique anchor filtering introduced by Assemblytics . RagTag can also remove alignments based on mapping quality score, when available. Filtered alignments are then merged to identify macro-synteny blocks. For each query sequence, alignments are sorted by reference position. Consecutive alignments within 100 kbp (configured using the “-d” parameter) of each other and on the same strand are merged together, taking the minimum coordinate as the new start position and the maximum coordinate as the new end position. Consequently, unmerged alignments are either far apart on the same reference sequence, on different reference sequences, or on different strands. Finally, merged alignments contained within other merged alignments (with respect to the query position) are removed.
Following the approach we developed for RaGOO, RagTag “correct” uses pairwise whole-genome sequence homology to identify and correct putative misassemblies. First, RagTag generates filtered and merged whole-genome alignments between a “query” and a “reference” assembly. The “query” assembly will be corrected and the “reference” assembly will be used to inform correction. Any query sequence with more than one merged alignment is considered for correction. RagTag breaks these query sequences at merged alignment boundaries provided that the boundaries are not within 5 kbp (-b) from either sequence terminus. Users may optionally choose to only break between alignments to the same or different reference sequences (--intra and --inter). If a GFF file is provided to annotate features in the query assembly, the query assembly will never be broken within a defined feature.
When the query and reference assemblies do not represent the same genotypes, unmerged alignments within a contig can indicate genuine structural variation. To help distinguish between structural variation and misassemblies, users can optionally provide Whole Genome Shotgun (WGS) sequencing reads from the same query genotype, such as short accurate reads or long error-corrected reads, to validate putative query breakpoints. RagTag aligns these reads to the query assembly with Minimap2 and computes the read coverage for each position in the query assembly. For each proposed query breakpoint, RagTag will flag exceptionally low (below --min-cov) or high (above --max-cov) coverage within 10 kbp (-v) of the proposed breakpoint. If exceptionally low or high coverage is not observed, the merged alignment boundaries are considered to be caused by true variation, and the query assembly is not broken at this position.
RagTag “scaffold” uses pairwise whole-genome sequence homology to scaffold a genome assembly. First, RagTag generates filtered and merged whole-genome alignments between a “query” and a “reference” assembly. The “query” assembly will be scaffolded and the “reference” assembly will be used to inform scaffolding. The merged alignments are used to compute a clustering, location, and orientation “confidence” score, just as is done in RaGOO, and sequences with confidence scores below certain thresholds are excluded (as set with parameters “-i”, “-a”, and “-s”). For each query sequence, the longest merged alignment is designated as the “primary” alignment. Primary alignments contained within other primary alignments (with respect to the reference coordinates) are removed. Primary alignments are then used to order and orient query sequences. To order query sequences, sequences are assigned to the reference chromosome to which they primarily align. Then, for each reference sequence, primary alignments are sorted by reference coordinate, establishing an order of query sequences. To orient query sequences, the sequence is assigned the same orientation as its primary alignment. Query sequences with no filtered alignments to the reference assembly (“unplaced” sequences) are output without modification or are optionally concatenated together.
By default, 100 bp gaps are placed between adjacent scaffolded query sequences, indicating an “unknown” gap size according to the AGP specification (https://www.ncbi.nlm.nih.gov/assembly/agp/AGP_Specification/). Optionally, RagTag can infer the gap size based on the whole-genome pairwise alignments. Let seq1 (upstream) and seq2 (downstream) be adjacent query sequences, and let aln1 and aln2 be their respective primary alignments. Let rs, re, qs, and qe denote the alignment reference start position, reference end position, query start position, and query end position, respectively. The following function computes the inferred gap length between seq1 and seq2:
where len(seq1) is the length of seq1. All inferred gap sizes must be at least 1 bp, and if the inferred gap size is too small (-g or less than 1) or too large (-m), it is replaced with an “unknown” gap size of 100 bp.
The new RagTag “patch” tool uses pairwise whole-genome sequence homology to make joins between contigs, without introducing gaps, and fill gaps in a “target” genome assembly using sequences from a “query” genome assembly. First, RagTag breaks all target sequences at gaps and generates filtered and merged whole-genome alignments between the query and target assemblies. Merged alignments that are not close (-i) to a target sequence terminus or are shorter than 50,000 bp (-s) are removed. If an alignment is not close to both query sequence termini yet it is not close to either target sequence terminus, meaning the target sequence should be contained within the query sequence, yet large portions of the target sequence do not align to the query sequence, the alignment is discarded.
To ultimately patch the target assembly, RagTag employs a directed version of a “scaffold graph” [18, 32]. Nodes in the graph are target sequence termini (two per target sequence), and edges connect termini of distinct target sequences observed to be adjacent in the input candidate scaffolds. The graph is initialized with the known target sequence adjacencies originally separated by gaps in the target assembly. Next, merged and filtered alignments are processed to identify new target sequence adjacencies. For each query sequence that aligns to more than one target sequence, alignments are sorted by query position. For each pair of adjacent target sequences, an edge is created in the scaffold graph. The edge stores metadata such as query sequence coordinates in order to continuously join the adjacent target sequences. If an edge already exists due to an existing gap, the gap metadata is replaced with the query sequence metadata so that the gap can be replaced with sequence. If an adjacency is supported by more than one alignment, the corresponding edge is discarded. To find a solution to this graph and output a patched assembly, a maximal weight matching is computed with networkx and if there are any cycles, they are broken . RagTag then iterates through each connected component and iteratively builds a sequence from adjacent target sequences. When target sequences are not overlapping, they are connected with sequence from the supporting query sequence. Unpatched target sequences are output without modification.
RagTag “merge” is a new implementation and extension of CAMSA, a tool to reconcile two or more distinct scaffolding solutions for a genome assembly . Input scaffolding solutions must be in valid AGP format, and they must order and orient the same set of genome assembly AGP “components.” RagTag iteratively builds a scaffold graph to store adjacency evidence provided by each AGP file. First, each AGP file is assigned a weight (1 by default). Then, for each AGP file and for each pair of adjacent components, an edge is added to the scaffold graph, and the edge weight is incremented by the weight of the AGP file, just as is done in CAMSA. After the scaffold graph is created, users can optionally replace native edge weights with Hi-C weights. To do this, Hi-C alignments are used to compute h(), the scaffold graph weights according to the SALSA2 algorithm, which uses the same underlying scaffold graph data structure . To find a solution to this graph and to output a merged AGP file, a maximal weight matching is computed with networkx and if there are any cycles, they are broken. RagTag then iterates through each connected component and iteratively builds AGP objects. Unmerged components are output without modification. While RagTag “merge” accepts any arbitrary number of input scaffolds, we advise that users only use the minimal number of informative scaffolds.
Patching a human genome assembly
The CHM13v1.1 assembly is the first-ever published complete sequence of a human genome . Though the original draft assembly was built exclusively from HiFi reads, it was manually inspected and patched at 25 loci, mostly at HiFi coverage dropouts, with sequence from the previously published, ONT-based CHM13v0.7 assembly. Using these 25 manual patches as a benchmark (Additional file 3), we evaluated the ability of RagTag to automatically patch the CHM13 draft assembly with the CHM13v0.7 assembly. RagTag made all 25 patches (Additional file 4), 19 of which were identical to the manual patches. The remaining six patches had slightly shifted patch coordinates, with a median Euclidean distance of 66.4 bp using the start and end genomic coordinates for the two joined sequences and the sequence used for patching. The slight differences in coordinates are due to locally repetitive sequences that cause aligner-specific coordinates to be selected when transitioning from the query and target sequences. RagTag made one false join connecting chr18 and chr10, though this was caused by a misassembly in CHM13v.07 caused by a long, high-identity repeat shared between these chromosomes. Patching was performed with RagTag v2.1.0 (--aligner minimap2).
Patching and merging multiple A. thaliana assemblies of varying quality
We performed patching and merging on several A. thaliana Columbia-0 draft genome assemblies to assess the impact of input genome assembly quality on RagTag accuracy. We used several assemblies of varying quality including published assemblies (“GCA_927323615” (https://www.ncbi.nlm.nih.gov/assembly/GCA_927323615.1/), “GCA_900243935”  (Additional file 1: Fig. S10), and “GCA_902825305” (https://www.ncbi.nlm.nih.gov/assembly/GCA_902825305.1)) and assemblies generated in-house from public data (“VERKKO”, “HIFIASM_L0_10X”, and “HIFIASM_L0_5X”) . Both Hifiasm assemblies used the “-L0” parameter, and the “10X” assembly and the “5X” assembly were derived from a random 10× and 5× subset of reads, respectively, to explore the outcomes for lower contiguity and lower accuracy input assemblies. All assemblies were screened using the same method described for the tomato assemblies. For patching, we first made a modified reference genome, breaking each Col-CEN v1.2 chromosome sequence into arms, excluding the centromeres, to promote contiguous unique alignments. We then ran RagTag “correct” with this reference to correct any potential misassemblies. Then, for patching, we used this modified reference assembly to patch the input “target” assemblies using the “--remove-small -f 75000 –aligner minimap2” parameters .
For merging, for each input assembly, we performed homology-based scaffolding multiple times using several different reference genomes. We used TAIR10, and the An-1 and C24 “1001 Genomes” assemblies as reference genomes [37, 38]. The accuracy of the patches was assessed by aligning the patch sequences including the neighboring 500 bp sequence window to the Col-CEN v1.2 reference genome. The number of simple repeats, satellite repeats, and transposable elements in patched sequences was quantified using EDTA  (Additional file 2: Table S1). For each query assembly, the individual homology-based scaffolding solutions were merged with RagTag “merge” using default parameters. For merging the low-contiguity HIFIASM_L0_5X assembly, we reduced the parameter for minimum contig-length input from the default of 100 kbp to 10 kbp in steps of 10 kbp to accommodate the smaller contigs present in this assembly, which enables it to reach a scaffold N50 more comparable to the high coverage assemblies (Additional file 2: Table S2 and Additional file 1: Fig. S3-S10). For the 10X coverage assembly, we noted the merging produced chromosome-scale scaffolds, but also propagated a few large mis-assemblies that were present in the initial contigs. To demonstrate how to address these errors with RagTag, we used the “correct” module to scan the input assemblies for mis-assemblies based on the alignment to the An-1 reference genome. As expected, this reduced the contiguity of the input assembly from a contig N50 of 2.8 Mbp to 2.4 Mbp. However, the merged assembly after correction achieves high scaffold N50 (14 Mbp) with noticeably fewer mis-assemblies in the final dotplot (Additional file 1: Fig. S5b). Note in nearly all cases the final merged assemblies have lower contiguity than the scaffolding results when using a single reference genome (e.g., a merged result of 14 Mbp vs 23 Mbp for the 10× corrected assembly when scaffolding to a single reference). This is the expected outcome since the scaffolding is conservative in regions where the input reference genomes disagree, such as a few large inversions present in these three reference genomes (Additional file 1: Fig. S10).
Benchmarking of several genome assembly patching tools
We compared RagTag “patch” to DENTIST (v3.0.0, read-coverage: 1, ploidy: 2, allow-single-reads: true, best-pile-up-margin: 1.5, existing-gap-bonus: 3.0, join-policy: contigs, min-reads-per-pile-up: 1, min-spanning-reads: 1, proper-alignment-allowance: 500), SAMBA (MaSuRCA-4.0.9, parameters -d asm -t 40 -m 5000), and Quickmerge (v0.3). Patching was performed as described in “Sweet-100 genome assembly” and “M82 genome assembly” where the respective M82/S100 HiFi contigs were patched with the M82/S100 ONT contigs. We used QUAST to evaluate the results, comparing the patched assemblies to the M82v1.0 reference, the S100v2.0 reference and the SL4.0 reference.
Extraction of high molecular weight DNA and sequencing
Extraction of high molecular weight genomic DNA, construction of Oxford Nanopore Technology libraries, and sequencing were described previously . Libraries for PacBio HiFi sequencing were constructed and sequenced at the Genome Technology Center at UNIL and Genome Center at CSHL. High molecular-weight DNA was sheared with a Megaruptor (Diagenode) to obtain 15–20 kbp fragments. After shearing, the DNA size distribution was evaluated with a Fragment Analyzer (Agilent) and 5–10 μg of the DNA was used to prepare a SMRTbell library with the PacBio SMRTbell Express Template Prep Kit 2.0 (Pacific Biosciences) according to the manufacturer's instructions. The library was size-selected on a BluePippin system (Sage Science) for molecules larger than 12.5 kbp and sequenced on one SMRT cell 8M with v2.0/v2.0 chemistry on a PacBio Sequel II instrument (Pacific Biosciences) at 30 hours movie length. Hi-C experiments were conducted with 2 g of flash-frozen leaf tissue using the Arima high-coverage Hi-C Service at Arima Genomics (San Diego, CA).
BLAST databases for screening contigs
We built each BLAST database with makeblastdb (v2.5.0+, -dbtype nucl) . We used all RefSeq bacterial genomes (downloaded on February 11th, 2021) for the bacterial genomes database. We used a collection of Solanum chloroplast sequences for the chloroplast database, and their GenBank accession IDs are as follows:
MN218076.1, MN218077.1, MN218078.1, MN218079.1, MN218091.1, MN218088.1, MN218089.1, NC_039611.1, NC_035724.1, KX792501.2, NC_041604.1, MH283721.1, NC_039605.1, NC_039600.1, NC_007898.3, MN218081.1, NC_039606.1, NC_030207.1, MT120858.1, MN635796.1, MN218090.1, MT120855.1, MT120856.1, NC_050206.1, MN218087.1, NC_008096.2
We used a collection of Solanum mitochondrial sequences for the mitochondria database, and their GenBank accession IDs are as follows:
MT122954.1, MT122955.1, MT122966.1, MT122969.1, MT122973.1, MT122974.1, MT122977.1, MT122988.1, NC_050335.1, MT122980.1, MT122981.1, MT122982.1, MT122983.1, MF989960.1, MF989961.1, NC_035963.1, MT122970.1, MT122971.1, NC_050334.1, MW122958.1, MW122959.1, MW122960.1, MT122964.1, MT122965.1, MW122949.1, MW122950.1, MW122951.1, MW122952.1, MW122953.1, MW122954.1, MW122961.1, MW122962.1, MW122963.1, MT122978.1, MT122979.1, MF989953.1, MF989957.1, MN114537.1, MN114538.1, MN114539.1, MT122958.1, MT122959.1
We used a collection of Solanum rDNA sequences for the rDNA database, and their GenBank accession IDs are as follows:
X55697.1, AY366528.1, AY366529.1, KF156909.1, KF156910.1, KF156911.1, KF156912.1, KF156913.1, KF156914.1, KF156915.1, KF156916.1, KF156917.1, KF156918.1, KF156919.1, KF156920.1, KF156921.1, KF156922.1, KF603895.1, KF603896.1, X65489.1, X82780.1, AF464863.1, AF464865.1, AY366530.1, AY366531.1, AY875827.1
Sweet-100 genome assembly
The following describes the methods used to produce SollycSweet-100_v2.0 assembly. We independently assembled all HiFi reads (33,815,891,985 bp) with Hifiasm (v0.13-r308, -l0) and we assembled ONT reads at least 30 kbp long (a total of 28,595,007,408 bp) with Flye (v2.8.2-b1689, --genome-size 1g) [41, 42]. The Hifiasm primary contigs were screened to remove contaminant or organellar contigs using the databases described above. Next, we used WindowMasker to mask repeats in the primary contigs (v1.0.0, -mk_counts -sformat obinary -genome_size 882654037) . We then aligned each contig to the bacterial, chloroplast, mitochondria, and rDNA BLAST databases with blastn (v2.5.0+, -task megablast). We only included the WindowMasker file for alignments to the bacterial database (-window_masker_db). For each contig, we counted the percentage of base pairs covered by alignments to each database. If more than 10% of a contig aligned to the rDNA database, we deemed it to be a putative rDNA contig. We then removed any contigs not identified as rDNA contigs that met any of the following criteria: 1) More than 10% of the contig was covered by alignments to the bacterial database; 2) More than 20% of the contig was covered by alignments to the mitochondria database and the contig was less than 1 Mbp long; or 3) More than 20% of the contig was covered by alignments to the chloroplast database and the contig was less than 0.5 Mbp long. In total, we removed 1015 contigs (35,481,360 bp) with an average length of 34,957.005 bp, most of which contained chloroplast sequence.
Even though Sweet-100 is an inbred line, to ensure that the assembly did not contain haplotypic duplication, we aligned all HiFi reads to the screened Hifiasm contigs with Winnowmap2 (v2.0, k=15, --MD -ax map-pb) . We then used purge_dups to compute and visualize the contig coverage distribution, and we determined that haplotypic duplication was not evident in the screened contigs .
We used RagTag “patch” to patch the screened Hifiasm contigs with sequences from the ONT flye contigs, and we manually excluded three incorrect patches caused by a missassembly in the Flye contigs. We then scaffolded the patched contigs using three separate approaches producing three separate AGP files. For the first two approaches, we used RagTag for homology-based scaffolding, once using the SL4.0 reference genome and once using the LA2093 v1.5 reference genome (v2.0.1, --aligner=nucmer --nucmer-params="--maxmatch -l 100 -c 500") [3, 4]. In both cases, only contigs at least 100 kbp long were considered for scaffolding, and the reference chromosome 0 sequences were not used for scaffolding. For the third scaffolding approach, we used Juicebox Assembly Tools to manually scaffold contigs with Hi-C data (using “arima” as the restriction enzyme), and we used a custom script to convert the “.assembly” file to an AGP file. We also separately generated Hi-C alignments by aligning the Hi-C reads to the screened contigs with bwa mem (v0.7.17-r1198-dirty) and processing the alignments with the Arima mapping pipeline (https://github.com/ArimaGenomics/mapping_pipeline) which employs Picard Tools (https://broadinstitute.github.io/picard/) . We merged the three AGP files with RagTag “merge” (v2.0.1, -r 'GATC,GA[ATCG]TC,CT[ATCG]AG,TTAA'), using Hi-C alignments to weight the Scaffold Graph (-b). Finally, using the merged scaffolds as a template, we made four manual scaffolding corrections in Juicebox Assembly tools. The final assembly contained 12 scaffolds corresponding to 12 chromosomes totaling 805,184,690 bp of sequence and 918 unplaced nuclear sequences totaling 40,749,555 bp.
VecScreen did not identify any “strong” or “moderate” hits to the adaptor contamination database (ftp://ftp.ncbi.nlm.nih.gov/pub/kitts/adaptors_for_screening_euks.fa) (https://www.ncbi.nlm.nih.gov/tools/vecscreen/). We packaged the assembly according to the pan-sol v0 specification (https://github.com/pan-sol/pan-sol-spec), and chromosomes were renamed and oriented to match the SL4.0 reference genome. The tomato chloroplast (GenBank accession NC_007898.3) and mitochondria (GenBank accession NC_035963.1) reference genomes were added to the final assembly.
To identify potential misassemblies and heterozygous Structural Variants (SVs), we aligned all HiFi reads (v2.0, k=15, --MD -ax map-pb) and ONT reads longer than 30 kbp (v2.0, k=15, --MD -ax map-ont) to the final assembly with Winnowmap2 and we called structural variants with Sniffles (v1.0.12, -d 50 -n -1 -s 5) . We removed any SVs with less than 30% of reads supporting the ALT allele and we merged the filtered SV calls (317 in total) with Jasmine (v1.0.10, max_dist=500 spec_reads=5 --output_genotypes) .
Sweet-100 gene and repeat annotation
We used Liftoff to annotate the Sweet 100 v2.0 assembly using ITAG4.0 gene models and tomato pan-genome genes as evidence (v1.5.1, -copies) [1, 3, 49]. Chloroplast and mitochondria annotations were replaced with their original GenBank annotation. Transcript, coding sequence, and protein sequences were extracted using gffread (v0.12.3, -y -w -x) . We annotated transposable elements with EDTA (v1.9.6, --cds --overwrite 1 --sensitive 1 --anno 1 --evaluate 1) .
M82 genome assembly
The M82 genome was assembled following the approach used for the Sweet-100 assembly, with the following distinctions. First, Hifiasm v0.15-r327 was used for assembling HiFi reads. Also, the M82 ONT assembly was polished before patching. M82 Illumina short-reads  were aligned to the draft Flye ONT assembly with BWA-MEM (v0.7.17-r1198-dirty) and alignments were sorted and compressed with samtools (v1.10) [46, 51]. Small variants were called with freebayes (v1.3.2-dirty, --skip-coverage 480), and polishing edits were incorporated into the assembly with bcftools “consensus” (v1.10.2, -i'QUAL>1 && (GT="AA" || GT="Aa")' -Hla) . In total, two iterative rounds of polishing were used. RagTag “merge” was also used for scaffolding, though the input scaffolding solutions used different methods than the Sweet-100 assembly. First, homology-based scaffolds were generated with RagTag “scaffold,” using the SL4.0 reference genome (v2.0.1, --aligner=nucmer --nucmer-params="--maxmatch -l 100 -c 500"). Contigs smaller than 300 kbp were not scaffolded (-j), and the reference chromosome 0 was not used to inform scaffolding (-e). Next, SALSA2 was used to derive Hi-C-based scaffolds. Hi-C reads were aligned to the assembly with the pipeline described for Sweet-100. We then produced scaffolds with SALSA2 (-c 300000 -p yes -e GATC -m no) and manually corrected false scaffolding joins in Juicebox Assembly Tools. We reconciled the homology-based and Hi-C-based scaffolds with RagTag “merge” using Hi-C alignments to re-weight the scaffold graph (-b). Finally, we made four manual corrections in Juicebox Assembly Tools. Cooler and HiGlass were used to visualize Hi-C heatmaps [53, 54]. Merqury was used to calculate QV and k-mer completeness metrics using 21-mers from the HiFi data .
Design of CRISPR-Cas9 gRNAs and cloning of constructs
CRISPR-Cas9 mutagenesis was performed as described previously . Briefly, guide RNAs (gRNAs) were designed based on the Sweet 100 v2.0 assembly and the CRISPRdirect tool (https://crispr.dbcls.jp/). Binary vectors for plant transformation were assembled using the Golden Gate cloning system as previously described .
Final vectors were transformed into the tomato cultivar S100 by Agrobacterium tumefaciens-mediated transformation according to Gupta and Van Eck (2016) with minor modifications . Briefly, seeds were sterilized for 15 min in 1.3% bleach followed by 10 min in 70% ethanol and rinsed four times with sterile water before sowing on MS media (4.4 g/L MS salts, 1.5 % sucrose, 0.8 % agar, pH 5.9) in Magenta boxes. Cotyledons were excised 7–8 days after sowing and incubated on 2Z- media  at 25°C in the dark for 24 h before transformation. A. tumefaciens were grown in LB media and washed in MS-0.2% media (4.4 g/L MS salts, 2% sucrose, 100 mg/L myo-inositol, 0.4 mg/L thiamine, 2 mg/L acetosyringone, pH5.8). Explants were co-cultivated with A. tumefaciens on 2Z- media supplemented with 100 μg/L IAA for 48 h at 25°C in the dark and transferred to 2Z selection media (supplemented with 150 mg/L kanamycin). Explants were transferred every two weeks to fresh 2Z selection media until shoot regeneration. Shoots were excised and transferred to selective rooting media  (supplemented with 150 mg/L kanamycin) in Magenta boxes. Well-rooted shoots were transplanted to soil and acclimated in a Percival growth chamber (~50 μmol m−2 s−1, 25°C, 50% humidity) before transfer to the greenhouse.
Validation of CRISPR-Cas9 editing
Genomic DNA was extracted from T0 plants using a quick genomic DNA extraction protocol. Briefly, small pieces of leaf tissue were flash frozen in liquid nitrogen and ground in a bead mill (Qiagen). Tissue powder was incubated in extraction buffer (100 mM Tris-HCl pH9.5, 250 mM KCl, 10 mM EDTA) for 10 min at 95°C followed by 5 min on ice. Extracts were combined with one volume of 3% BSA, vigorously vortexed, and spun at 13,000 rpm for 1 min. One microliter supernatant was used as template for PCR using primers flanking the gRNA target sites. PCR products were separated on agarose gels and purified for Sanger Sequencing (Microsynth) using ExoSAP-IT reagent (Thermo Fisher Scientific). Chimeric PCR products were subcloned before sequencing using StrataClone PCR cloning kits (Agilent).
High-throughput genotyping of T1 individuals was conducted by barcoded amplicon sequencing according to Liu et al. (2021) with minor modifications . Briefly, gene-specific amplicons were diluted ten-fold before barcoding and pools of barcoded amplicons were gel-purified (NEB Monarch DNA gel extraction) before Illumina library preparation and sequencing (Amplicon-EZ service at Genewiz). Editing efficiencies were quantified from a total of 194,590 aligned reads using the CRISPResso2 software (--min_frequency_alleles_around_cut_to_plot 0.1 --quantification_window_size 50) . All oligos used in this study are listed in Additional file 2: Tables S6-S8.
Availability of data and materials
Genome assemblies and annotations are available at https://github.com/pan-sol/pan-sol-data (https://doi.org/10.5281/zenodo.6814693) , and the Solanaceae Genomics Network (https://solgenomics.net/ftp/genomes/M82/;https://solgenomics.net/ftp/genomes/Sweet-100/). Raw sequence data is available on SRA under the BioProject PRJNA779684 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA779684) . Seeds are available on request from S. Soyk. The RagTag software is freely available on GitHub under the MIT license: https://github.com/malonge/RagTag (https://doi.org/10.5281/zenodo.5634263) [61, 62].
Gao L, Gonda I, Sun H, Ma Q, Bao K, Tieman DM, et al. The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor. Nat Genet. 2019;51:1044–51.
Alonge M, Wang X, Benoit M, Soyk S, Pereira L, Zhang L, et al. Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato. Cell. 2020;182:145–61.e23 Elsevier.
Hosmani PS, Flores-Gonzalez M, van de Geest H, Maumus F, Bakker LV, Schijlen E, et al. An improved de novo assembly and annotation of the tomato reference genome using single-molecule sequencing, Hi-C proximity ligation and optical maps. bioRxiv. 2019:767764 biorxiv.org. https://doi.org/10.1101/767764.
Wang X, Gao L, Jiao C, Stravoravdis S, Hosmani PS, Saha S, et al. Genome of Solanum pimpinellifolium provides insights into structural variants during tomato breeding. Nat Commun. 2020;11:5817 nature.com.
Rengs WMJ, Schmidt MHW, Effgen S, Le DB, Wang Y, Zaidan MWAM, et al. A chromosome scale tomato genome built from complementary PacBio and Nanopore sequences alone reveals extensive linkage drag during breeding. Plant J. 2022;110:572–88. Available from: https://onlinelibrary.wiley.com/doi/10.1111/tpj.15690.
Meissner R, Jacobson Y, Melamed S, Levyatuv S, Shalev G, Ashri A, et al. A new model system for tomato genetics. Plant J. 1997;12:1465–72 Wiley.
Eshed Y, Zamir D. An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL. Genetics. 1995;141:1147–62 academic.oup.com.
Menda N, Semel Y, Peled D, Eshed Y, Zamir D. In silico screening of a saturated mutation library of tomato. Plant J. 2004;38:861–72 Wiley Online Library.
Logsdon GA, Vollger MR, Eichler EE. Long-read human genome sequencing and its applications. Nat Rev Genet. 2020;21:597–614.
Wenger AM, Peluso P, Rowell WJ, Chang P-C, Hall RJ, Concepcion GT, et al. Accurate circular consensus long-read sequencing improves variant detection and assembly of a human genome. Nat Biotechnol. 2019;37:1155–62.
Nurk S, Koren S, Rhie A, Rautiainen M, Bzikadze AV, Mikheenko A, et al. The complete sequence of a human genome. Science. 2022;376:44–53 American Association for the Advancement of Science (AAAS).
Zimin AV, Smith DR, Sutton G, Yorke JA. Assembly reconciliation. Bioinformatics. 2007;24:42–5 Oxford Academic.
Alhakami H, Mirebrahim H, Lonardi S. A comparative evaluation of genome assembly reconciliation tools. Genome Biol. 2017;18:1–14 BioMed Central.
Rhie A, McCarthy SA, Fedrigo O, Damas J, Formenti G, Koren S, et al. Towards complete and error-free genome assemblies of all vertebrate species. Nature. 2021;592:737–46.
Howe K, Chow W, Collins J, Pelan S, Pointon D-L, Sims Y, et al. Significantly improving the quality of genome assemblies through curation. Gigascience. 2021:10. https://doi.org/10.1093/gigascience/giaa153academic.oup.com.
Soyk S, Müller NA, Park SJ, Schmalenbach I, Jiang K, Hayama R, et al. Variation in the flowering gene SELF PRUNING 5G promotes day-neutrality and early yield in tomato. Nat Genet. 2017;49:162–8 nature.com.
Alonge M, Soyk S, Ramakrishnan S, Wang X, Goodwin S, Sedlazeck FJ, et al. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol. 2019;20:224.
Aganezov SS, Alekseyev MA. CAMSA: a tool for comparative analysis and merging of scaffold assemblies. BMC Bioinformatics. 2017;18:496.
Ghurye J, Rhie A, Walenz BP, Schmitt A, Selvaraj S, Pop M, et al. Integrating Hi-C links with assembly graphs for chromosome-scale assembly. PLoS Comput Biol. 2019;15:e1007273.
Ludwig A, Pippel M, Myers G, Hiller M. DENTIST—using long reads for closing assembly gaps at high accuracy. Gigascience. 2022;11:giab100 Oxford Academic.
Zimin AV, Salzberg SL. The SAMBA tool uses long reads to improve the contiguity of genome assemblies. PLoS Comput Biol. 2022;18:e1009860 Public Library of Science.
Chakraborty M, Baldwin-Brown JG, Long AD, Emerson JJ. Contiguous and accurate de novo assembly of metazoan genomes with modest long read coverage. Nucleic Acids Res. 2016;44:e147 Oxford Academic.
Dudchenko O, Shamim MS, Batra S, Durand NC. The Juicebox Assembly Tools module facilitates de novo assembly of mammalian genomes with chromosome-length scaffolds for under $1000. bioRxiv. 2018. https://doi.org/10.1101/254797.
Gupta S, Van Eck J. Modification of plant regeneration medium decreases the time for recovery of Solanum lycopersicum cultivar M82 stable transgenic lines. Plant Cell Tiss Org Cult. 2016;127:417–23 Springer Nature.
Jack T, Brockman LL, Meyerowitz EM. The homeotic gene APETALA3 of Arabidopsis thaliana encodes a MADS box and is expressed in petals and stamens. Cell. 1992;68:683–97 Elsevier.
El-Brolosy MA, Kontarakis Z, Rossi A, Kuenne C, Günther S, Fukuda N, et al. Genetic compensation triggered by mutant mRNA degradation. Nature. 2019;568:193–7 nature.com.
Lippman ZB, Cohen O, Alvarez JP, Abu-Abied M, Pekker I, Paran I, et al. The making of a compound inflorescence in tomato and related nightshades. PLoS Biol. 2008;6:e288 journals.plos.org.
Soyk S, Lemmon ZH, Oved M, Fisher J, Liberatore KL, Park SJ, et al. Bypassing Negative Epistasis on Yield in Tomato Imposed by a Domestication Gene. Cell. 2017;169:1142–55.e12.
Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34:3094–100.
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. Versatile and open software for comparing large genomes. Genome Biol. 2004;5:R12.
Nattestad M, Schatz MC. Assemblytics: a web analytics tool for the detection of variants from an assembly. Bioinformatics. 2016;32:3021–3.
Ghurye J, Pop M, Koren S, Bickhart D, Chin C-S. Scaffolding of long read assemblies using long range contact information. BMC Genomics. 2017;18:527.
Galil Z. Efficient algorithms for finding maximum matching in graphs. ACM Comput Surv, vol. 18. New York: Association for Computing Machinery; 1986. p. 23–38.
Jupe F, Rivkin AC, Michael TP, Zander M, Timothy Motley S, Sandoval JP, et al. The complex architecture and epigenomic impact of plant T-DNA insertions. PLoS Genet. 2019;15:e1007819 Public Library of Science.
High-quality Arabidopsis thaliana Genome Assembly with Nanopore and HiFi Long Reads. Genomics Proteomics Bioinformatics. Elsevier; 2021 [cited 2022 Jul 8]; https://doi.org/10.1016/j.gpb.2021.08.003
Naish M, Alonge M, Wlodzimierz P, Tock AJ, Abramson BW, Schmücker A, et al. The genetic and epigenetic landscape of the centromeres. Science. 2021;374:eabi7489.
1001 Genomes Consortium. Electronic address: email@example.com, 1001 Genomes Consortium. 1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana. Cell. 2016;166:481–91.
Berardini TZ, Reiser L, Li D, Mezheritsky Y, Muller R, Strait E, et al. The Arabidopsis information resource: Making and mining the “gold standard” annotated reference plant genome. Genesis. 2015;53:474–85.
Ou S, Su W, Liao Y, Chougule K, Agda JRA, Hellinga AJ, et al. Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline. Genome Biol. 2019;20:275.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Cheng H, Concepcion GT, Feng X, Zhang H, Li H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021;18:170–5 nature.com.
Kolmogorov M, Yuan J, Lin Y, Pevzner PA. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 2019;37:540–6.
Morgulis A, Gertz EM, Schäffer AA, Agarwala R. WindowMasker: window-based masker for sequenced genomes. Bioinformatics. 2006;22:134–41.
Jain C, Rhie A, Hansen NF, Koren S, Phillippy AM. Long-read mapping to repetitive reference sequences using Winnowmap2. Nat Methods. 2022. https://doi.org/10.1038/s41592-022-01457-8.
Guan D, McCarthy SA, Wood J, Howe K, Wang Y, Durbin R. Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics. 2020;36:2896–8.
Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv [q-bio.GN]. 2013. Available from: http://arxiv.org/abs/1303.3997
Sedlazeck FJ, Rescheneder P, Smolka M, Fang H, Nattestad M, von Haeseler A, et al. Accurate detection of complex structural variations using single-molecule sequencing. Nat Methods. 2018;15:461–8 nature.com.
Kirsche M, Prabhu G, Sherman R, Ni B, Aganezov S, Schatz MC. Jasmine: Population-scale structural variant comparison and analysis. bioRxiv. 2021:2021.05.27.445886 [cited 2021 Sep 28]. Available from: https://www.biorxiv.org/content/10.1101/2021.05.27.445886v1.abstract.
Shumate A, Salzberg SL. Liftoff: accurate mapping of gene annotations. Bioinformatics. 2020. https://doi.org/10.1093/bioinformatics/btaa1016.
Pertea G, Pertea M. GFF Utilities: GffRead and GffCompare. F1000Res. 2020;9:304.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, et al. Twelve years of SAMtools and BCFtools. Gigascience. 2021;10. https://doi.org/10.1093/gigascience/giab008.
Kerpedjiev P, Abdennur N, Lekschas F, McCallum C, Dinkla K, Strobelt H, et al. HiGlass: web-based visual exploration and analysis of genome interaction maps. Genome Biol. 2018;19:125 genomebiology.biomedcentral.com.
Abdennur N, Mirny LA. Cooler: scalable storage for Hi-C data and other genomically labeled arrays. Bioinformatics. 2020;36:311–6 academic.oup.com.
Rhie A, Walenz BP, Koren S, Phillippy AM. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol. 2020;21:245.
Brooks C, Nekrasov V, Lippman ZB, Van Eck J. Efficient gene editing in tomato in the first generation using the clustered regularly interspaced short palindromic repeats/CRISPR-associated9 system. Plant Physiol. 2014;166:1292–7 academic.oup.com.
Liu L, Chen R, Fugina CJ, Siegel B, Jackson D. High-throughput and low-cost genotyping method for plant genome editing. Curr Protoc. 2021;1:e100 Wiley.
Clement K, Rees H, Canver MC, Gehrke JM, Farouni R, Hsu JY, et al. CRISPResso2 provides accurate and rapid genome editing sequence analysis. Nat Biotechnol. 2019;37:224–6 nature.com.
Alonge M, Lebeigle L, Kirsche M, Jenike K, Ou S, Aganezov S, et al. Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing: Zenodo; 2022. https://doi.org/10.5281/zenodo.6814693.
Alonge M, Lebeigle L, Kirsche M, Jenike K, Ou S, Aganezov S, et al. Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing: Datasets. NCBI SRA; 2022. https://www.ncbi.nlm.nih.gov/bioproject/PRJNA779684
Alonge M, Lebeigle L, Kirsche M, Jenike K, Ou S, Aganezov S, et al. Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing: Github; 2022. https://github.com/malonge/RagTag
Alonge M, Lebeigle L, Kirsche M, Jenike K, Ou S, Aganezov S, et al. Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing: Zenodo; 2022. https://doi.org/10.5281/zenodo.5634263.
We thank all members of the Lippman, Schatz, and Soyk labs for helpful discussions. We thank B. Tissot-Dit-Sanfin, A. Robadey, and V. Vashanthakumar from UNIL for assistance with plant care. We thank J. Marquis, J. Weber, and M. Dupasquier from the UNIL Genome Technology Facility (GTF) for sequencing support and René Dreos for bioinformatic support. We thank Sara Goodwin, Melissa Kramer, and the Cold Spring Harbor Laboratory Woodbury Genome Center for sequencing support. We thank Lei Liu and Penelope Lindsay for helpful discussions about amplicon sequencing. We thank Steven Salzberg, Aleksey Zimin, Luca Venturini, and Alexander Leonard for their help with software features and improvements.
Peer review information
Kevin Pang was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
The review history is available as Additional file 6.
This work was performed with assistance from the US National Institutes of Health Grant S10OD028632-01. M.A. and M.C.S. were supported by the National Science Foundation (NSF) DBI-1350041, and NSF PGRP grant IOS-1732253 to Z.B.L. and M.C.S. K.J. and M.C.S. were supported in part by NSF PGRP grant IOS-1543901 to M.C.S. S.O. and M.C.S. were supported in part by NFS PGRP grant IOS-1758800 and Human Frontiers Scientific Program (HFSP) award RGP0025/2021 to M.C.S. This work was supported by funding from the Howard Hughes Medical Institutes to Z.B.L., and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (ERC Starting Grant “EPICROP” Grant No. 802008) and a Swiss National Science Foundation (SNSF) Eccellenza Professorial Fellowship (Grant No. PCEFP3_181238) to S.S.
Ethics approval and consent to participate
The authors declared that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Supplementary Figures for Automated Assembly Scaffolding using RagTag elevates a new tomato system for high-throughput genome editing.
Additional file 2. Supplementary Tables for Automated Assembly Scaffolding using RagTag elevates a new tomato system for high-throughput genome editing.
About this article
Cite this article
Alonge, M., Lebeigle, L., Kirsche, M. et al. Automated assembly scaffolding using RagTag elevates a new tomato system for high-throughput genome editing. Genome Biol 23, 258 (2022). https://doi.org/10.1186/s13059-022-02823-7
- Genome sequencing
- Assembly scaffolding
- Genome editing