Long-read Sequencing Technologies Research Articles

Chromosome-level genome assembly of the smallscale yellowfin (Plagiognathops microlepis).

The small-scale yellowfin (Plagiognathops microlepis) is a highly valued species in East Asian aquaculture due to its adaptability and high yield. However, the lack of genomic data has impeded genetic research and breeding efforts. In this study, we utilize PacBio Hifi long-read sequencing and Hi-C technologies to construct a highly detailed genome of P. microlepis at the chromosomal level. The assembly encompasses 976.41 Mb, with an exceptional 99.84% distribution across 24 chromosomes. Notably, the contig N50 was 34.41 Mb and scaffold N50 was 38.38 Mb. The completeness of the P. microlepis genome assembly is underscored by a BUSCO score of 98.08%. A total of 25,389 protein-coding genes were identified, with a BUSCO score of 96.98%, and 99.85% of these genes were functionally annotated. Synteny relationships at the chromosome level with Danio rerio and Chanodichthys erythropterus genomes uncover small-scale chromosomal rearrangements. This high-fidelity genome assembly serves as a pivotal resource for forthcoming endeavors such as the genome structure, functional elements, comparative genomics, and evolutionary characteristics of P. microlepis and its relative species.

Abstract A023: Advancing RNA liquid biopsy technology via long-read sequencing

Abstract Recent studies have begun to highlight the advantages of RNA-based liquid biopsy technology for cancer early detection, and we sought to further advance the development of RNA liquid biopsies by incorporating single-molecule, long-read sequencing technology. Cell-free RNAs are dynamically secreted by the cells that comprise our tissues and organs but are thought to be largely degraded in the blood. To determine the true length of cell-free RNAs, we performed long-read nanopore sequencing on cell-free RNAs isolated from pancreatic and esophageal cancer patient plasma and found that thousands of genes are detectable as full-length cell-free RNAs in vivo. By creating a custom annotation of the transcriptome for cell-free RNA quantification, we found that cancer patients exhibit high levels of repetitive cell-free RNAs in plasma, with full-length, SINE-derived cell-free RNAs being among the most abundant repetitive RNAs. Moreover, we found that utilizing repetitive cell-free RNA features substantially improved the performance of our machine learning models for cancer classification. Notably, nanopore sequencing enabled us to discover over 250,000 novel cell-free RNAs that have not been previously annotated in the human genome. Our findings highlight the value and utility of long-read nanopore sequencing of cell-free RNAs for discovering and leveraging novel RNA biomarkers for cancer early detection via RNA liquid biopsy technology. Citation Format: Daniel H. Kim. Advancing RNA liquid biopsy technology via long-read sequencing [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr A023.

Real-time plasmid transmission detection pipeline.

Plasmid-mediated spread of antimicrobial resistance is a major challenge for clinical microbiology, and monitoring of potential plasmid transmissions is essential to combat further dissemination. Whole-genome sequencing is often used to surveil nosocomial transmissions but usually limited to the detection of clonal transmissions (based on chromosomal markers). Recent advances in long-read sequencing technologies enable full reconstruction of plasmids and the detection of very similar plasmids, but so far, easy-to-use bioinformatic tools for this purpose have been missing. Here, we present an evaluation of an innovative real-time plasmid transmission detection pipeline. It is integrated into the GUI-based SeqSphere+ software, which already offers core-genome multi-locus sequence typing-based pathogen outbreak detection. It requires very limited bioinformatics knowledge, and its database, automated analyses, and alert system make it well suited for prospective clinical application.

Complete genome sequences of five bacteria isolated from rice plants in a paddy field in Sekinchan, Selangor, Malaysia.

This study examines the genome sequences of five endophytic bacterial isolates from the Oryza sativa microbiome to assess their potential as plant bio-inoculants. The five complete bacterial genomes from the genera Pseudomonas, Burkholderia, Sphingobacterium, Stenotrophomonas, and Pantoea were sequenced using Nanopore long-read sequencing technology.

Long-read RNA sequencing reveals allele-specific N6-methyladenosine modifications.

Long-read sequencing technology enables highly accurate detection of allele-specific RNA expression, providing insights into the effects of genetic variation on splicing and RNA abundance. Furthermore, the ability to directly sequence RNA using the Oxford Nanopore technology promises the detection of RNA modifications in tandem with ascertaining the allelic origin of each molecule. Here, we leverage these advantages to determine allele-biased patterns of N6-methyladenosine (m6A) modifications in native mRNA. We utilized human and mouse cells with known genetic variants to assign allelic origin of each mRNA molecule combined with a supervised machine learning model to detect read-level m6A modification ratios. Our analyses revealed the importance of sequences adjacent to the DRACH-motif in determining m6A deposition, in addition to allelic differences that directly alter the motif. Moreover, we discovered allele-specific m6A modification (ASM) events with no genetic variants in close proximity to the differentially modified nucleotide, demonstrating the unique advantage of using long reads and surpassing the capabilities of antibody-based short-read approaches. This technological advancement promises to advance our understanding of the role of genetics in determining mRNA modifications.

An optimized protocol for quality control of gene therapy vectors using nanopore direct RNA sequencing.

Despite recent advances made toward improving the efficacy of lentiviral gene therapies, a sizeable proportion of produced vector contains an incomplete and thus potentially nonfunctional RNA genome. This can undermine gene delivery by the lentivirus as well as increase manufacturing costs and must be improved to facilitate the widespread clinical implementation of lentiviral gene therapies. Here, we compare three long-read sequencing technologies for their ability to detect issues in vector design and determine nanopore direct RNA sequencing to be the most powerful. We show how this approach identifies and quantifies incomplete RNA caused by cryptic splicing and polyadenylation sites, including a potential cryptic polyadenylation site in the widely used Woodchuck Hepatitis Virus Posttranscriptional Regulatory Element (WPRE). Using artificial polyadenylation of the lentiviral RNA, we also identify multiple hairpin-associated truncations in the analyzed lentiviral vectors (LVs), which account for most of the detected RNA fragments. Finally, we show that these insights can be used for the optimization of LV design. In summary, nanopore direct RNA sequencing is a powerful tool for the quality control and optimization of LVs, which may help to improve lentivirus manufacturing and thus the development of higher quality lentiviral gene therapies.

Generation and analysis of a mouse multitissue genome annotation atlas.

Generating an accurate and complete genome annotation for an organism is complex because the cells within each tissue can express a unique set of transcript isoforms from a unique set of genes. A comprehensive genome annotation should contain information on what tissues express what transcript isoforms at what level. This tissue-level isoform information can then inform a wide range of research questions as well as experiment designs. Long-read sequencing technology combined with advanced full-length cDNA library preparation methods has now achieved throughput and accuracy where generating these types of annotations is achievable. Here, we show this by generating a genome annotation of the mouse (Mus musculus). We used the nanopore-based R2C2 long-read sequencing method to generate 64 million highly accurate full-length cDNA consensus reads-averaging 5.4 million reads per tissue for a dozen tissues. Using the Mandalorion tool, we processed these reads to generate the Tissue-level Atlas of Mouse Isoforms which is available as a trackhub for the UCSC Genome Browser and contains at least one full-length isoform for the vast majority of expressed genes in each tissue.

GCI: a continuity inspector for complete genome assembly.

Recent advances in long-read sequencing technologies have significantly facilitated the production of high-quality genome assembly. The telomere-to-telomere (T2T) gapless assembly has become the new golden standard of genome assembly efforts. Several recent efforts have claimed to produce T2T-level reference genomes. However, a universal standard is still missing to qualify a genome assembly to be at T2T standard. Traditional genome assembly assessment metrics (N50 and its derivatives) have no capacity in differentiating between nearly T2T assembly and the truly T2T assembly in continuity either globally or locally. Additionally, these metrics are independent of raw reads, making them inflated easily by artificial operations. Therefore, a gaplessness evaluation tool at single-nucleotide resolution to reflect true completeness is urgently needed in the era of complete genomes. Here, we present a tool called Genome Continuity Inspector (GCI), designed to assess genome assembly continuity at single-base resolution, and evaluate how close an assembly is to the T2T level. GCI utilizes multiple aligners to map long reads from various sequencing platforms back to the assembly. By incorporating curated mapping coverage of high-confidence read alignments, GCI identifies potential assembly issues. Meanwhile, it provides GCI scores that quantify overall assembly continuity on the whole genome or chromosome scales. The open-source GCI code is freely available on Github (https://github.com/yeeus/GCI) under the MIT license.

The Ribosomal Operon Database: A Full-Length rDNA Operon Database Derived From Genome Assemblies.

Current rDNA reference sequence databases are tailored towards shorter DNA markers, such as parts of the 16/18S marker or the internally transcribed spacer (ITS) region. However, due to advances in long-read DNA sequencing technologies, longer stretches of the rDNA operon are increasingly used in environmental sequencing studies to increase the phylogenetic resolution. There is, therefore, a growing need for longer rDNA reference sequences. Here, we present the ribosomal operon database (ROD), which includes eukaryotic full-length rDNA operons fished from publicly available genome assemblies. Full-length operons were detected in 34.1% of the 34,701 examined eukaryotic genome assemblies from NCBI. In most cases (53.1%), more than one operon variant was detected, which can be due to intragenomic operon copy variability, allelic variation in non-haploid genomes, or technical errors from the sequencing and assembly process. The highest copy number found was 5947 in Zea mays. In total, 453,697 unique operons were detected, with 69,480 operon variant clusters remaining after intragenomic clustering at 99% sequence identity. The operon length varied extensively across eukaryotes, ranging from 4136 to 16,463 bp, which will lead to considerable polymerase chain reaction (PCR) bias during amplification of the entire operon. Clustering the full-length operons revealed that the different parts (i.e., 18S, 28S, and the hypervariable regions V4 and V9 of 18S) provide divergent taxonomic resolution, with 18S, the V4 and V9 regions being the most conserved. The ROD will be updated regularly to provide an increasing number of full-length rDNA operons to the scientific community.

A chromosome-scale genome assembly of the peanut beetle, Ulomoides dermestoides (Coleptera: Tenebrionidae)

Ulomoides dermestoides, a member of the Tenebrionidae family, is known by various common names such as Chinese beetle, Chinese weevil, peanut beetle, cancer beetle, or asthma beetle. Despite its origin in Asia, it has achieved global distribution owing to its perceived medicinal properties. In this study, we employed a combination of Illumina short-read, PacBio long-read, and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to produce a high-quality, chromosome-scale genome assembly of U. dermestoides. The assembled genome spans 253.38 Mb across 11 chromosomes, including the two sex chromosomes, X and Y. Key metrics include contig and scaffold N50 values of 9.83 Mb and 28.62 Mb, respectively, with an estimated genome completeness of 99.62%. Annotation efforts identified 15,553 protein-coding genes and revealed that 48.11% of the genome comprises repetitive sequences. This meticulously assembled genome serves as a pivotal resource for elucidating the genetic determinants underlying the species’ ecological habits and potential medicinal properties.

NANOPORE-BASED INTRA-OPERATIVE CLASSIfiCATION OF CENTRAL NERVOUS SYSTEM TUMOURS IN THE NHS

Abstract AIMS The 2021 WHO classification of CNS tumours necessitates the use of molecular testing in order to reach a confident diagnosis of many tumour types. At present, this is achieved using array-based technology delivered by a centralised model. This approach has inherent delays that mean it is often several weeks before the treating clinical team is informed of diagnosis, delaying the start of adjuvant therapy and causing significant distress for patients. As such, there is an urgent need to rethink this pathway and improve time to diagnosis. Nanopore longread sequencing technology has the potential to enable a paradigm shift in diagnostic and therapeutic pathways, offering rapid, comprehensive molecular diagnosis that can be performed in the intra-operative setting. We aimed to test the hypothesis that intra-operative methylation-based classification can be performed using nanopore sequencing and thereby expediently inform surgical and oncological decision-making. METHOD Oxford Nanopore Technology (ONT) is a native-strand, long-read sequencing platform that can perform copy- number profiling along with parallel single-gene methylation, structural variant, mutational and methylation analyses. Capital and consumable costs are relatively low, and testing can be performed in the local setting on single samples. We have optimised DNA extraction and library preparation for rapid tissue processing. We have developed a robust bioinformatic pipeline which is able to produce intraoperative classification within minutes of sequencing and full molecular profiling after 24 hours. RESULTS We demonstrate the feasibility of using ONT to deliver intra-operative methylation-based classification of CNS tumours in ‘realtime’ within an NHS hospital setting and will discuss the challenges of implementing this testing and the opportunities to inform neuro-oncological practice. CONCLUSION By giving a summary of our novel protocol of rapid DNA extraction, library preparation and bioinformatic analysis, we will show that this approach can revolutionise the diagnosis of brain tumours and has the potential for rapid adoption by non-expert laboratories.

7030 Molecular Analysis Of Congenital Adrenal Hyperplasia-X Syndrome And Development Of A Pilot Panel For Analyzing Structural Variations With Long-Read Sequencing

Abstract Disclosure: J. Kim: None. S. Park: None. J. Yoon: None. D. Kim: None. S. Hwang: None. G. Kim: None. J. Kim: None. J. Choi: None. H. Yoo: None. Purpose: Genomic recombination events between the CYP21A2 and CYP21A1P genes are frequent due to their high sequence homology. Deletions in both the CYP21A2 and TNXB genes result in a chimeric TNXA/TNXB gene, leading to congenital adrenal hyperplasia (CAH)-X syndrome, characterized by a hypermobile form of Ehlers-Danlos syndrome. The aim of study was to determine the frequency of CAH-X syndrome and to develop a pilot panel for structural variation analysis using long-read sequencing technology. Methods: Among the 165 unrelated families with confirmed genetic mutations in the CYP21A2 gene, 44 probands had large deletions in one or both alleles. Long-range PCR was performed in 12 patients out of 44 probands using a forward primer for the 5’-UTR of the CYP21A2 and CYP21A1P and a reverse primer specific for exon 32 of TNXB. The amplified PCR products were then subjected to Sanger sequencing. We constructed a custom panel of 140 kb containing the RCCX modules arranged as STK19-C4A-CYP21A1P-TNXA-STK19B-C4B-CYP21A2-TNXB. In 4 patients with deletions in both the CYP21A2 and TNXB genes, long-read sequencing was performed using the custom panel on the PacBio Sequel II platform. Results: Out of 12 probands, eight harbored TNXA/TNXB chimeras. One patient was homozygous for CAH-X CH1 with a 120 bp deletion in exon 35 and had a Beighton score of 6 for joint hypermobility and multiple diverticula in the ascending colon. Two patients were heterozygous for CAH-X CH1, while three were heterozygous for the CAH-X CH2 with an intact exon 35 and the c.P4058W in exon 40. Two patients carried heterozygous mutations p.S4175N in exon 43, which is an unclassified type. Using long-read sequencing with a custom panel, missense mutations and TNXA/TNXB chimeras were detected in four patients previously identified by conventional direct sequencing and MLPA. Conclusions: This study identified 8 patients with CAH-X syndrome carrying four different TNXA/TNXB chimeras by long-range PCR. Long-read sequencing emerged as a comprehensive and efficient method for the molecular analysis of CAH, especially for the identification of structural variations. Biallelic deletion of TNXB was associated with severe clinical manifestations. Therefore, screening for CAH-X is essential for clinical management and prevention of the long-term consequences of CAH-X syndrome. Presentation: 6/1/2024

9266 Molecular Analysis Of Congenital Adrenal Hyperplasia-X Syndrome And Development Of A Pilot Panel For Analyzing Structural Variations With Long-Read Sequencing

Abstract Disclosure: J. Kim: None. S. Park: None. J. Yoon: None. D. Kim: None. S. Hwang: None. G. Kim: None. J. Kim: None. J. Choi: None. H. Yoo: None. Purpose: Genomic recombination events between the CYP21A2 and CYP21A1P genes are frequent due to their high sequence homology. Deletions in both the CYP21A2 and TNXB genes result in a chimeric TNXA/TNXB gene, leading to congenital adrenal hyperplasia (CAH)-X syndrome, characterized by a hypermobile form of Ehlers-Danlos syndrome. The aim of study was to determine the frequency of CAH-X syndrome and to develop a pilot panel for structural variation analysis using long-read sequencing technology. Methods: Among the 165 unrelated families with confirmed genetic mutations in the CYP21A2 gene, 44 probands had large deletions in one or both alleles. Long-range PCR was performed in 12 patients out of 44 probands using a forward primer for the 5’-UTR of the CYP21A2 and CYP21A1P and a reverse primer specific for exon 32 of TNXB. The amplified PCR products were then subjected to Sanger sequencing. We constructed a custom panel of 140 kb containing the RCCX modules arranged as STK19-C4A-CYP21A1P-TNXA-STK19B-C4B-CYP21A2-TNXB. In 4 patients with deletions in both the CYP21A2 and TNXB genes, long-read sequencing was performed using the custom panel on the PacBio Sequel II platform. Results: Out of 12 probands, eight harbored TNXA/TNXB chimeras. One patient was homozygous for CAH-X CH1 with a 120 bp deletion in exon 35 and had a Beighton score of 6 for joint hypermobility and multiple diverticula in the ascending colon. Two patients were heterozygous for CAH-X CH1, while three were heterozygous for the CAH-X CH2 with an intact exon 35 and the c.P4058W in exon 40. Two patients carried heterozygous mutations p.S4175N in exon 43, which is an unclassified type. Using long-read sequencing with a custom panel, missense mutations and TNXA/TNXB chimeras were detected in four patients previously identified by conventional direct sequencing and MLPA. Conclusions: This study identified 8 patients with CAH-X syndrome carrying four different TNXA/TNXB chimeras by long-range PCR. Long-read sequencing emerged as a comprehensive and efficient method for the molecular analysis of CAH, especially for the identification of structural variations. Biallelic deletion of TNXB was associated with severe clinical manifestations. Therefore, screening for CAH-X is essential for clinical management and prevention of the long-term consequences of CAH-X syndrome. Presentation: 6/1/2024

AAVolve: Concatenated long-read deep sequencing enables whole capsid tracking during shuffled AAV library selection

Gene therapies using recombinant adeno-associated virus (AAV) vectors have demonstrated considerable clinical success in the treatment of genetic disorders. Improved vectors with favourable tropism profiles, decreased immunogenicity and enhanced manufacturability are poised to further improve the state of gene therapies. Such vectors can be identified through directed evolution, a process of subjecting a diverse capsid library to a selection pressure to identify individual variants with a desired trait. Currently, libraries that involve changes distributed throughout the AAV capsid coding region, such as DNA family shuffled libraries, are largely characterised using low-throughput Sanger sequencing of individual clones. However, improvements in long-read sequencing technologies have increased their applicability to capsid libraries and evaluation of the selection process. Here we explore the application of Oxford Nanopore Technologies refined by a concatemeric consensus method for initial library characterization and monitoring selection of a shuffled AAV capsid library. Furthermore, we present AAVolve, a bioinformatic pipeline for processing long-read data from AAV directed evolution experiments. Our approach allows high throughput characterisation of AAV capsids in a streamlined manner, facilitating deeper insights into library composition through multiple rounds of selection, and generalization through training of machine learning models.

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research.

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation.

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level.

Chromosome-level genome assembly of Huai pig (Sus scrofa)

Although advances in long-read sequencing technology and genome assembly techniques have facilitated the study of genomes, little is known about the genomes of unique Chinese indigenous breeds, including the Huai pig. Huai pig is an ancient domestic pig breed and is well-documented for its redder meat color and high forage tolerance compared to European domestic pigs. In the present study, we sequenced and assembled the Huai pig genome using PacBio, Hi-C, and Illumina sequencing technologies. The final highly contiguous chromosome-level Huai pig genome spans 2.53 Gb with a scaffold N50 of 138.92 Mb. The Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score for the assembled genome was 95.33%. Remarkably, 23,389 protein-coding genes were annotated in the Huai-pig genome, along with 45.87% repetitive sequences. Overall, this study provided new foundational resources for future genetic research on Chinese domestic pigs.

Characterisation of the Novel HLA-DRB1*16:79 Allele Using Short- and Long-Read Sequencing Technologies.

The novel HLA-DRB1*16:79 allele, first described in an individual from Brazil.

Comparison of short-read and long-read next-generation sequencing technologies for determining HIV-1 drug resistance.

Accurate HIV-1 genome sequencing is necessary to identify drug resistance mutations (DRMs) in people with HIV-1 (PWH). Next-generation-sequencing (NGS) allows the detection of minor variants and is now available in many laboratories. Our study aimed to compare two NGS approaches, a "short read" sequencing protocol using DeepChek® Whole Genome HIV-1 Assay on Illumina, and a "long read" sequencing protocol of HIV-1 pol and env single-molecule real-time sequencing (SMRT) on Pacific Biosciences (PacBio). We analyzed 16 plasma samples and 13 cellular samples from PWH. HIV-1 whole genome was amplified into five amplicons using DeepChek® Whole Genome HIV-1 Assay and sequenced on an iSeq. 100. In parallel, HIV-1 pol and env genes were separately amplified and sequenced using PacBio SMRT system with the circular consensus sequencing mode on a Sequel IIe. Concordance rates for determining DRMs with both approaches varied depending on the HIV-1 region, with higher concordance in the integrase region compared to the reverse transcriptase and protease regions. DeepChek® Whole Genome HIV-1 Assay exhibited better sensitivity in HIV-1 RNA sequencing of plasmas with lower viral loads. In cell HIV-1 DNA sequencing, the DeepChek® Whole Genome HIV-1 Assay performed better in pol and env sequencing but detected more APOBEC-induced DRMs, which can represent defective proviruses. Our findings indicate that both DeepChek® Whole Genome HIV-1 Assay and PacBio SMRT sequencing exhibit good performance for subtype determination, detection, and quantification of DRMs of the HIV-1 genome. However, some discrepancies were found in cellular samples, highlighting the challenges of interpreting HIV-1 DNA DRMs.

Performance of somatic structural variant calling in lung cancer using Oxford Nanopore sequencing technology

BackgroundLung cancer is a heterogeneous disease and the primary cause of cancer-related mortality worldwide. Somatic mutations, including large structural variants, are important biomarkers in lung cancer for selecting targeted therapy. Genomic studies in lung cancer have been conducted using short-read sequencing. Emerging long-read sequencing technologies are a promising alternative to study somatic structural variants, however there is no current consensus on how to process data and call somatic events. In this study, we preformed whole genome sequencing of lung cancer and matched non-tumour samples using long and short read sequencing to comprehensively benchmark three sequence aligners and seven structural variant callers comprised of generic callers (SVIM, Sniffles2, DELLY in generic mode and cuteSV) and somatic callers (Severus, SAVANA, nanomonsv and DELLY in somatic modes).ResultsDifferent combinations of aligners and variant callers influenced somatic structural variant detection. The choice of caller had a significant influence on somatic structural variant detection in terms of variant type, size, sensitivity, and accuracy. The performance of each variant caller was assessed by comparing to somatic structural variants identified by short-read sequencing. When compared to somatic structural variants detected with short-read sequencing, more events were detected with long-read sequencing. The mean recall of somatic variant events identified by long-read sequencing was higher for the somatic callers (72%) than generic callers (53%). Among the somatic callers when using the minimap2 aligner, SAVANA and Severus achieved the highest recall at 79.5% and 79.25% respectively, followed by nanomonsv with a recall of 72.5%.ConclusionLong-read sequencing can identify somatic structural variants in clincal samples. The longer reads have the potential to improve our understanding of cancer development and inform personalized cancer treatment.