A high-quality telomere-to-telomere LSDV genome assembly

This study presents a high-quality, telomere-to-telomere assembly of the Lumpy skin disease virus (LSDV) Oman 2009 isolate using a hybrid short- and long-read sequencing approach, which resolves complex telomeric structures, corrects previous misassemblies, and establishes a robust reference for future genomic surveillance and evolutionary analysis.

The Gist

Imagine a virus as a complex instruction manual for building a tiny, dangerous machine. For years, scientists have been trying to read the Lumpy Skin Disease Virus (LSDV) manual, but they were only reading the middle chapters. The beginning and the end of the book—the parts that tell the virus how to start, stop, and trick the animal's immune system—were a jumbled mess of repeated sentences that standard reading tools couldn't decipher.

This paper is like a team of expert editors who finally managed to reconstruct the entire manual from cover to cover, including the messy, repetitive pages at the very start and finish.

Here is the story of how they did it, explained simply:

1. The Problem: The "Blurry" Book

LSDV is a virus that hurts cows, causing painful lumps on their skin and costing farmers billions of dollars. To fight it, scientists need to understand its genetic code (its DNA).

Previously, scientists tried to read this code using short-read sequencing. Think of this like trying to assemble a 1,000-piece puzzle where every piece is a tiny, identical blue square. You can see the blue color, but you can't tell which piece goes where, especially at the edges. This led to "blurry" maps of the virus, with missing pieces and errors at the ends.

2. The Solution: The "Long-Read" Camera

The researchers used a new, high-tech approach called hybrid assembly. They combined two types of data:

• Short Reads (Illumina): These are like high-resolution photos. They are very accurate but short.
• Long Reads (Nanopore): These are like a wide-angle drone video. They might be slightly grainier, but they capture long stretches of the landscape in one go.

By using the "drone video" (long reads) to lay out the big picture and the "high-res photos" (short reads) to fill in the tiny details, they created a Telomere-to-Telomere (T2T) assembly. This means they didn't just get the middle; they got the very first page to the very last page, perfectly connected.

3. The Big Discovery: The "Bookends"

The most exciting part of this new map is the ITRs (Inverted Terminal Repeats). Imagine the virus genome as a book. The ITRs are the front and back covers, which are identical but flipped. These covers are crucial because they contain the "on/off switches" for the virus.

• The Old Map: Was missing about 22% of these covers. It was like trying to understand a book but having the first 50 pages torn off and the last 50 pages crumpled.
• The New Map: Shows the covers perfectly. The researchers found that the virus has a specific "hairpin" structure at the ends (like a folded paperclip) that acts as a key to unlock the virus's ability to replicate. They also found a specific sequence called the CRS (Concatemer Resolution Sequence), which is the "scissors" the virus uses to cut its long chains of DNA into individual virus packets.

4. What Was Hidden in the Shadows?

Because the old maps were so blurry, scientists missed some important details:

• Truncated Genes: They found that two specific genes (LSDV019 and LSDV026) are "broken" or shortened in the wild version of the virus. It's like finding a car engine with two missing pistons. This explains why the wild virus behaves differently than the vaccine version (which has the full, working engines).
• A New Character: They discovered a tiny, previously invisible gene (LSDV042.5). It's so small it was hiding in the shadows, but it turns out to be a crucial tool the virus uses to sneak into cow cells.

5. Why Does This Matter?

Think of this new genome as the definitive "Master Blueprint" for Lumpy Skin Disease.

• Better Vaccines: Now that we know exactly where the "broken" parts are, scientists can design vaccines that target the real weaknesses of the virus.
• Faster Detection: When a new outbreak happens, doctors can use this perfect map to quickly scan samples and spot mutations (typos in the manual) that might make the virus stronger or resistant to drugs.
• Future Proofing: This paper proves that for complex viruses with repetitive ends, you need the "drone video" (long reads) to get the full picture. It sets a new standard for how we study not just this virus, but many others in the future.

In a nutshell: The scientists stopped guessing at the edges of the virus's instruction manual. They used a mix of high-tech tools to read the whole thing, revealing hidden secrets about how the virus starts, stops, and tricks its host. This gives us a much better chance of stopping it in its tracks.

Technical Summary

1. Problem Statement

Lumpy Skin Disease Virus (LSDV), a capripoxvirus (CaPV), causes significant economic losses in livestock across Africa, Asia, and Europe. While genomic variation influences virulence and transmission, existing understanding of the LSDV genome structure is limited by previous assembly methods.

• Limitations of Short-Read Sequencing: Most existing LSDV genomes were assembled using short-read Illumina technology. This approach fails to resolve highly repetitive regions, specifically the Inverted Terminal Repeats (ITRs) at the 5' and 3' ends of the linear genome.
• Consequences: Short-read assemblies result in fragmented contigs, structural ambiguities, and misassemblies in the terminal regions. These errors obscure true biological variation, hinder accurate annotation of virulence genes located at the termini (e.g., LSDV001 and LSDV156), and confound phylogenetic and evolutionary analyses.
• Knowledge Gap: The precise architecture of LSDV ITRs, including hairpin structures and concatemer resolution sequences (CRS), remains poorly defined compared to other poxviruses like Vaccinia virus.

2. Methodology

The authors generated a Telomere-to-Telomere (T2T) reference genome for the LSDV Oman 2009 isolate using a hybrid sequencing strategy that integrates long-read and short-read data.

• Sample Source: The virus was isolated from a cattle skin nodule and passaged in MDBK cells.
• Sequencing Strategy:
  • Long-Reads (ONT): Oxford Nanopore Technology (ONT) R10.4.1 flow cells were used to generate long reads (median 533 bp, max ~81 kb) without target enrichment. This provided the primary scaffold for the assembly.
  • Short-Reads (Illumina): Illumina NextSeq 2000 (paired-end 150 bp) was used with a custom myBaits probe panel to enrich specifically for LSDV DNA, achieving 99% viral read content.
• Assembly Pipeline:
  1. De Novo Assembly: The Flye assembler (v2.9.1) was used on the ONT reads to generate the initial contig.
  2. Polishing: A multi-step polishing process was applied to correct errors:
     • Medaka: Corrected large indels and homopolymer errors using ONT reads.
     • Pilon: Fixed SNPs and small indels using Illumina reads.
     • Polypolish: Corrected errors in repetitive regions using Illumina reads.
     • Pypolca: Final correction step using Illumina reads.
  3. Quality Control: Read depth was visualized using Minimap2 and SAMtools; assembly quality was assessed via QUAST and CheckV.
  4. Annotation: Genes were annotated using Prokka (guided by reference KX894508) and GATU (mapped from AF325528.1), followed by manual inspection in IGV.

3. Key Contributions

• First T2T LSDV Assembly: This is the first complete, telomere-to-telomere assembly of an LSDV genome, fully resolving the complex, repeat-rich ITRs.
• Hybrid Assembly Validation: The study demonstrates that a hybrid approach (ONT scaffolding + Illumina polishing) effectively resolves structural ambiguities that short-read-only or even PacBio-only assemblies miss.
• Correction of Reference Errors: The new assembly corrects misassemblies found in previous genomes, particularly regarding the length and sequence of the terminal repeats.
• Discovery of Novel Features: The high-resolution assembly allowed for the identification of previously unannotated or truncated genes and confirmed clade-specific structural variations.

4. Key Results

• Genome Statistics:
  • Length: 151,091 bp (single contig, N50 = 151,091 bp).
  • GC Content: 25.9%.
  • ORFs: 157 annotated open reading frames.
• Resolution of ITRs:
  • The assembly revealed 5' and 3' ITRs spanning approximately 2.5 kb each (positions 1–2,547 and 148,557–151,091).
  • Short-read Bias: Short-read mapping showed overrepresentation (22% overestimation) in the ITRs due to repetitive structures, whereas long-reads provided uniform coverage, confirming the true length.
  • CRS Identification: The study identified specific Concatemer Resolution Sequences (CRS) at the hairpin termini (123 bp hairpin), matching known GTPV sequences.
• Comparative Analysis:
  • vs. PacBio Assembly (PX492334): The Oman 2009 assembly showed only 12 mismatches across the whole genome but revealed that the PacBio assembly was missing the first 67 bp at both termini. The hybrid approach resolved an additional 56 bp at each end.
• Gene Annotation Findings:
  • Truncated Genes: Confirmed truncation of LSDV019 and LSDV026 in wild-type clade 1.1 strains (Oman 2009), contrasting with full-length versions in clade 1.2 vaccine strains.
  • New ORF: Identified a novel, small ORF named LSDV042.5 (29 amino acids), an orthologue of Vaccinia virus O3L, which is essential for host cell penetration. It was previously missed due to size cutoffs in older annotations.

5. Significance

• Improved Genomic Surveillance: The T2T reference provides a robust template for accurate read mapping, mutation detection, and evolutionary inference, which are critical for tracking the spread of LSDV.
• Vaccine and Diagnostic Development: Accurate characterization of terminal genes (involved in immune modulation and virulence) and clade-specific truncations aids in the development of better diagnostic assays and safer, more effective vaccines.
• Methodological Benchmark: The study validates that long-read sequencing (specifically ONT R10) is sufficient for high-quality CaPV assembly, with short-reads serving best for polishing and resequencing. This sets a new standard for resolving complex viral genome structures.
• Data Availability: All raw reads (SRA: PRJNA1366067) and the assembly (Accession: PV877838) are publicly available, facilitating immediate community use.