Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “long read”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

Accelerated long-read variant calling with Clair3 for whole-genome sequencing.

SUMMARY: The rapid growth of genomic data and increasing adoption of long-read sequencing technologies have rendered variant calling one of the most computationally demanding tasks in genomic analysis. Although deep learning-based methods currently outperform conventional approaches in distinguishing true variants from complex sequencing noise, they impose prohibitive computational and time requirements. To address this limitation, we present a computational framework based on Clair3 that integrates parallelized feature generation, enhanced variant phasing, in-memory read haplotagging, and GPU-accelerated neural network inference to accelerate variant calling. By dynamically optimizing the use of both GPU and CPU resources, our method achieves substantial runtime improvements without compromising accuracy. We evaluated our framework across a range of sequencing depths, diverse samples, and multiple hardware configurations. Our results demonstrate that the optimized pipeline completes variant calling for a 30× whole-genome sequence in 12-20 minutes using standard computational resources (32 CPU threads and one NVIDIA GPU), and in 12-15 minutes on an Apple Mac Studio (32 threads), which is ∼10-20-fold speedup compared with its initial release. In addition to exceptional efficiency, our method maintains state-of-the-art accuracy, achieving SNP F1-scores of 99.32% and 99.70% on 30× ONT and PacBio GIAB HG003 datasets, respectively. This work introduces a rapid, accurate, and scalable variant calling framework that effectively supports large-cohort genomic studies and time-sensitive clinical applications. AVAILABILITY AND IMPLEMENTATION: The accelerated implementation of Clair3 is open source and available at: https://github.com/HKU-BAL/Clair3/tree/gpu.

Whole Genome Sequencing↗

needLR: long-read structural variant annotation with population-scale frequency estimation.

SUMMARY: We present needLR, a structural variant (SV) annotation tool that can be used for filtering and prioritization of candidate pathogenic SVs from long-read sequencing data using population allele frequencies, annotations for genomic context, and gene-phenotype associations. When using population data from 500 presumably healthy individuals to evaluate nine test cases with known pathogenic SVs, needLR assigned allele frequencies to over 97.5% of all detected SVs and reduced the average number of novel genic SVs to 121 per case while retaining all known pathogenic variants. AVAILABILITY AND IMPLEMENTATION: needLR is implemented in bash with dependencies including Truvari v4.2.2, BEDTools v2.31.1, and BCFtools v1.19. Source code, documentation, and pre-computed population allele frequency data are freely available at https://github.com/jgust1/needLR under an MIT license and archived on Zenodo at https://zenodo.org/records/19463479.

Software↗

Charting host structural variations in cervical cancer by long-read sequencing pinpoints a functional deletion in PIAS1.

Host structural variations (SVs) are critical in cancer development but their landscape and interaction with HPV integration in cervical carcinogenesis remain unclear. In this study, we performed Nanopore long-read sequencing on five HPV-positive cervical cancer tissues and two cell lines to profile host SVs. We identified thousands of SVs and statistically demonstrated their significant enrichment in genomic windows ±25 to ±50 kb from HPV integration sites. Cross-sample analysis revealed 60 shared SVs, including a recurrent deletion within the PIAS1 gene. Multi-omics integration (Hi-C, H3K27ac ChIP-seq, and TCGA data) showed that this deletion is associated with reduced PIAS1 expression, disruption of local topologically associating domains, advanced pathological tumor stage, and poorer overall survival. Functional assays confirmed that PIAS1 deficiency inhibits cervical cancer cell proliferation and migration. Our findings identify a PIAS1 deletion as a candidate driver event, and underscore the pivotal role of host genomic instability in HPV-associated oncogenesis.

Cervical cancer↗

Long-read proteogenomic atlas of human neuronal differentiation reveals isoform diversity informing neurodevelopmental risk mechanisms.

RNA splicing shapes neuronal identity and disease risk, yet current maps lack the developmental resolution and depth to resolve this complexity. Here, we integrate deep long-read RNA sequencing and proteomics in induced pluripotent stem cell-derived cortical neurons to generate a high-resolution proteogenomic atlas of human neuron development. We identify 182,371 mRNA isoforms (over half previously unknown) and provide direct peptide evidence for the translation of hundreds of novel protein-coding sequences. Population genetics demonstrates that variants affecting novel exons and splice sites are under negative selection, underscoring the potential significance of these isoforms. During neuronal maturation, we observe that autism risk genes undergo dynamic isoform switching, including microexon inclusion and intron retention, that remodel key protein domains and regulatory regions. Furthermore, we uncover widespread, long-range coordination between alternative transcript processing events, including transcription start sites, exon splicing, and polyadenylation. Finally, our atlas enables variant reinterpretation in autism, highlighting the value of an isoform-centric view for interpreting pathogenic variation in neurodevelopment.

Humans↗

Long-read sequencing resolves complex CYP21A2 variants and identifies 2+0 carriers in 21-hydroxylase deficiency.

The complex CYP21A2 variants arising from high homology with its pseudogene CYP21A1P challenge the diagnosis of 21-hydroxylase deficiency (21-OHD). This study systematically evaluated long-read sequencing (LRS) for identifying complex structural variants of the CYP21A2 gene in 21-OHD in comparison with conventional molecular diagnostic methods, including multiplex ligation-dependent probe amplification (MLPA), CNVplex, and SNaPshot. Twenty patients with suspected 21-OHD and defined CYP21A2 structural variants identified via initial MLPA screening were enrolled. Variants were further analyzed using CNVplex and SNaPshot, then all samples underwent LRS for comprehensive variant detection, breakpoint mapping, and haplotype resolution. LRS overcame key limitations of conventional methods. It reliably identified a novel large-fragment deletion and defined its boundaries. Notably, LRS identified "2+0" carriers, where deletions masked by duplications cause false-negatives with standard techniques. Moreover, LRS accurately distinguished CYP21A1P/CYP21A2_CH-4 and CH-9 chimera subtypes which were indistinguishable by the combined conventional assays. Furthermore, LRS enabled the precise identification and characterization of TNXA/TNXB chimeric deletions. These are frequently misclassified as CYP21A1P/CYP21A2 chimeras by conventional methods but are critical for diagnosing associated conditions such as CAH-X syndrome. LRS provides a superior, integrated solution for the molecular diagnosis of 21-OHD, offering precise structural variant characterization, accurate carrier detection, and reliable breakpoint mapping. Its application enhances diagnostic accuracy, supports advanced genetic counseling, and paves the way for genotype-informed clinical management.

Journal Article↗

Long-read sequencing of single cell-derived melanoma subclones reveals divergent and parallel genomic and epigenomic evolutionary trajectories.

Tumor evolution is driven by various mutational processes, ranging from single-nucleotide variants (SNVs) to large structural variants (SVs) to dynamic shifts in DNA methylation. Current short-read sequencing methods struggle to accurately capture the full spectrum of these genomic and epigenomic alterations due to inherent technical limitations. To overcome that, here we introduce an approach for long-read sequencing of single-cell derived subclones, and use it to profile 23 subclones of a mouse melanoma cell line, characterized with distinct growth phenotypes and treatment responses. We develop a computational framework for harmonization and joint analysis of different variant types in the evolutionary context. Uniquely, our framework enables detection of recurrent amplifications of putative driver genes, generated by independent SVs across different lineages, suggesting parallel evolution. In addition, our approach revealed gradual and lineage-specific methylation changes associated with aggressive clonal phenotypes. We also show our set of phylogeny-constrained variant calls along with openly released sequencing data can be a valuable resource for the development of new computational methods.

Journal Article↗

Community-driven updates for comprehensive long-read metagenomics and enhanced binning in nf-core/mag v5.

SUMMARY: nf-core/mag is a reproducible Nextflow pipeline for best-practice metagenomic de novo assembly and binning within the nf-core framework. Here we present a major update that adds support for long-read-only assembly and bin refinement, includes five new binning tools, expands taxonomic classification to viruses and eukaryotes, and improves bin quality evaluation with new tools and latest databases. Through sustained community-driven development spanning seven years and four primary curator teams, nf-core/mag remains actively developed as an open source workflow for metagenomic analysis, benefiting from contributions from across the broader metagenomics, nf-core, and Nextflow ecosystem. AVAILABILITY AND IMPLEMENTATION: The source code of nf-core/mag v5 is available on GitHub (https://github.com/nf-core/mag) under the open source MIT license, with v5.5.0 source code archived on Zenodo (https://zenodo.org/records/21735731). Documentation is viewable on the nf-core website (https://nf-co.re/mag).

Metagenomics↗

Long-read sequencing reveals widespread novel splicing and neojunction-derived neoantigens in nasopharyngeal carcinoma.

The widespread transcriptomic diversity driven by alternative splicing (AS) contributes to all hallmarks of cancer and represents a critical source of neoantigens for personalized immunotherapy. However, unlike other major malignancies, the full repertoire of AS in nasopharyngeal carcinoma (NPC) remains underexplored. Here, we employ long-read sequencing (LR-seq) to generate a high-resolution, isoform-level transcriptomic atlas from a cohort of 14 NPC tumor samples and four immortalized nasopharyngeal epithelial cell lines. We identify a substantial number of full-length novel transcripts (22,687; ∼44.38%), which reveal diverse splicing patterns and previously unannotated splicing events. By integrating short-read RNA-seq data to quantify isoform expression, we discover a subset of novel transcripts that are differentially expressed between tumor samples and immortalized nasopharyngeal epithelial cell lines. Furthermore, LR-seq enables precise identification of chimeric readthrough fusion transcripts, such as CLDN15-FIS1 and FOXRED2-TXN2 Finally, we develop a computational framework, tumor-specific splicing neoantigen detection (TS-SNAD), to predict neoantigens originating from novel exon-exon junctions (neojunctions) in tumor-specific novel transcripts. Using this framework, we identify neojunction-derived neoantigens and experimentally validate the immunogenicity of selected HLA-B*40:01-restricted neoantigens. These neojunction-derived peptides constitute a new class of noncanonical neoantigens with significant potential for developing personalized cancer vaccines for NPC.

Humans↗

Detecting and reconstructing breakage-fusion-bridge cycles from long-read sequencing using BFBArchitect.

MOTIVATION: Focal oncogene amplification is a key driver of tumor progression. Remarkably, the increased pathology depends on the context-whether the amplification is extrachromosomal (ecDNA) or intrachromosomal. EcDNA amplifications promote heterogeneity, therapy resistance, and poor prognosis. Focal intrachromosomal amplifications often arise through breakage-fusion-bridge (BFB) cycles, which produce highly rearranged but stable chromosomes. Distinguishing BFB from ecDNA remains challenging due to overlapping genomic signatures. To address this, we present BFBArchitect, a computational method leveraging long-read Oxford Nanopore data to identify BFB sequences consistent with both copy number and structural variations. RESULTS: We provide a novel combinatorial characterization of BFB, which naturally leads to an integer linear programming (ILP) optimization. The ILP optimization generates a BFB sequence that best explains experimentally observed copy numbers and foldback structural variants. We implement this idea in a tool called BFBArchitect, which achieves near-perfect accuracy in distinguishing BFB from non-BFB structures in extensive simulations as well as on 18 validated tumor samples. Moreover, it generates sequence-level BFB reconstructions that provide mechanistic insights into BFB formation, including repair mechanisms with template switching and other structural variants, and recapture of telomere for stabilization. AVAILABILITY AND IMPLEMENTATION: BFBArchitect is available at https://github.com/AmpliconSuite/BFBArchitect.

Sequence Analysis, DNA↗

Long-read sequence analysis of the MECP2 gene in Rett syndrome patients: correlation of disease severity with mutation type and location.

Mutations in the methyl-CpG-binding protein gene MECP2 at Xq28 cause Rett syndrome (RTT), an X-linked dominant neurodevelopmental disorder characterized by a period of stagnation followed by regression in the development of young girls. Mutations were sought in MECP2 in 48 females with classical sporadic RTT, seven families with possible familial RTT and five sporadic females with features suggestive, but not diagnostic of RTT. Long distance PCR coupled with long-read direct sequencing was employed to sequence the entire MECP2 gene coding region in all cases. Mutations were identified in 44/55 (80%) unrelated classical sporadic and familial RTT patients, but only 1/5 (20%) sporadic cases with suggestive but non-diagnostic features of RTT. Twenty-one different mutations were identified (12 missense, four nonsense and five frame-shift mutations); 14 of these were novel. All missense mutations were located either in the methyl-CpG-binding domain or in the transcription repression domain. Nine recurrent mutations were characterized in a total of 33 unrelated cases (73% of all cases with MECP2 mutations). Significantly milder disease was noted in patients carrying missense mutations as compared with those with truncating mutations ( P = 0. 0023), and milder disease was associated with late as compared with early truncating mutations ( P = 0.0190).

Amino Acid Sequence↗

Dissecting the relationship between haplotypes around ATXN2 CAG repeats and the number of CAA interruptions by long-read sequencing.

BACKGROUND: CAG repeat expansions in ATXN2 are implicated as risk factors for several neurological diseases, including spinocerebellar ataxia type 2 (SCA2) when >=33 CAG repeats are present, and amyotrophic lateral sclerosis (ALS) when 27-33 CAG repeats are present. However, how haplotypes around the repeats and CAA interruptions within the repeats are associated with disease phenotypes remains poorly understood. Previous studies on haplotypes around ATXN2 were limited to SNPs very close to the repeats (<5kb) or were based on statistical inference only. METHODS: Here, we used long-read sequencing on the Oxford Nanopore Technologies (ONT) platform to simultaneously infer haplotypes around ATXN2, the number of CAG repeats, and the number of CAA interruptions, along with NYGC ALS Consortium NGS dataset. We further sequenced 41 individuals (EUR = 39) with neurological diseases with intermediate repeats by ONT. RESULTS: We found that haplotypes around ATXN2 and the number of interruptions show ethnicity-specific and ALS-specific distribution. Three CAA interruptions are present at low prevalence (~1%) in control populations in multiple ancestry groups, but high prevalence (~55%) in ALS individuals with intermediate repeats. Furthermore, we examined 159 individuals with ALS (~90% European ancestry) with intermediate ATXN2 repeats and found a unique haplotype in ALS individuals with three CAA interruptions, which can be tagged by an SNV, rs148019457. We also validated that the rs148019457-G allele is only present in haplotypes with three CAA interruptions. CONCLUSIONS: In summary, our study shows that 3 CAA interruptions are rarely seen in healthy controls but are common in those with expanded ATXN2 CAG repeats who have neurological disorders, and that rs148019457 tags a specific haplotype with 3 CAA interruptions within expanded ATXN2 CAG repeats in individuals of European ancestry. These results have implications for the development of precision genomic medicine for neurological disorders, and the tag SNP may help identify those with interruptions from existing population genotyping data.

ATXN2↗

Rethinking the pathogenicity of intragenic DMD duplications detected by carrier screening: High prevalence of nontandem duplications revealed by long-read sequencing.

PURPOSE: The pathogenicity of intragenic duplications depends on their structural configuration. Tandem duplications often disrupt reading frames and cause gene loss of function, whereas interspersed (nontandem) duplications are largely benign. When the configuration cannot be determined, current guidelines presume a tandem structure, leading to some laboratories automatically classifying such variants as likely pathogenic or pathogenic. This study evaluates the validity of this presumption for DMD, in patients with and without clinical indications of dystrophinopathy. METHODS: We performed high-coverage long-read genome sequencing on 15 patients with intragenic DMD duplications. A total of 4 patients had clinically indicated dystrophinopathy testing, whereas in the remaining 11 patients, the duplications were detected without clear indications of dystrophinopathy (eg, through carrier screening). RESULTS: All 4 patients with clinical indications had tandem duplications. In contrast, 64% (7/11) of the cases without such indications had interspersed duplications, with 4 subsequently reclassified as likely benign, 2 (likely) pathogenic, and 1 uncertain. These duplications were often complex, involving coduplications or codeletions with other regions. CONCLUSION: Our findings challenge the presumption that intragenic DMD duplications are predominantly in tandem. This highlights the need for a cautious variant interpretation approach, particularly in carrier screening and other settings in which variants are identified without indications of dystrophinopathy.

Humans↗

Versatile wastewater monitoring of pathogens and antimicrobial resistance enabled by metatranscriptomics and long-read metagenomics.

Widespread interest in the development of population-wide pathogen and antimicrobial resistance (AMR) monitoring has revealed wastewater's microbial footprint as a marker of public health. Near-source wastewater remains a difficult sample type for microbiome analyses but represents a closer link to human health than the downstream products of its treatment. Few studies integrate methods for non-targeted monitoring applications, and critically, current methods cannot connect AMR genes to species, nor resolve full genomes. We address these challenges by developing a pipeline that enables untargeted metagenomics, metatranscriptomics, and novel long-read metagenomics (LRG). We achieve untargeted pathogen detection, limited by highly abundant resident species, while retaining microbial information with near-source sampling. Furthermore, LRG identifies antibiotic resistance gene-containing microbes and enables assembly of culture-independent genomes with previously unreported AMR genes. We establish an integrated approach to broadly monitor pathogens in wastewater, while demonstrating the importance of LRG to illuminate microbial AMR at the species level.

Journal Article↗

Inferring the demographic history of Chinese and Indian rhesus macaque (Macaca mulatta) populations from PacBio HiFi long-read sequencing data.

The rhesus macaque (Macaca mulatta) is one of the most widely used animal models in biomedical research, both as it resembles humans in key biological aspects and as it is characterized by a broad geographic range. Most of the individuals housed in U.S. research colonies have been sampled from either China or India, though notably the source population of these animals has significantly shifted over time. Given the substantial genetic and immunological differences between these populations, a deeper understanding of the underlying population structure is critically important for biomedical interpretation. Despite this, the demographic histories of these two populations remain poorly resolved. Here, we present an analysis of whole-genome, PacBio HiFi long-read sequencing data from ten unrelated individuals of each population, applying four related model- and non-model based demographic inference approaches, in order to reconstruct their ancestral history. We evaluated the fit of the subsequently estimated models against the empirical data, and incorporated underlying uncertainty in the mutation rates used for scaling. We inferred a well-fitting population history characterized by substantial structure between Chinese and Indian populations, with a split time &#x223c;140,000 generations ago from an ancestral population of &#x223c;65,000 individuals. We additionally inferred the subsequent history of size change within, and gene flow between, these populations, reaching the current estimated sizes of &#x223c;220,000 individuals in the Chinese population and &#x223c;14,000 individuals in the Indian population. The robust baseline demographic model established in this study will serve as a valuable resource for future research on this species, including for improved fine-scale recombination mapping, selection inference, and association studies.

Cercopithecidae↗

Long-read Sequences Mapped to a Complete Reference Genome Uncover Uncaptured Structural Variants across the Beta-globin Cluster in Africans with Sickle Cell Disease.

African genomes are marked by extensive complexity in the number and distribution of variants, yet remain under-represented in genetic databases and the human reference genome. This gap in representation limits the broad application of genomic medicine. Sickle cell disease (SCD) - one of the most common monogenic diseases - has its highest prevalence in Africa, and variation in disease severity has consistently been linked to the beta-globin locus, including levels of fetal hemoglobin (HbF). Modulation of HbF is central to current SCD gene therapies; however, the inherent complexity and variation at the locus in African genomes presents a challenge to translating these advances to Africa. Here, we align long-read single molecule sequences (LRS) targeted to the beta-globin region to the hg38 and T2T-CHM13v2 genome references in 40 individuals with SCD, predominantly recruited from three African countries. We demonstrate that the expanded T2T-CHM13v2 reference sequence at this locus reduces Structural Variant (SV) calls by 70% and uncovers uncaptured single nucleotide variants (SNVs). Across the cluster we report 343 SVs and 196 SNVs that have not been previously reported, including in LRS data from the All of Us project. By including African populations from ethnolinguistic groups that have not been previously surveyed we improve variant resolution and bolster evidence for observed variation. Finally, we identify a common &#x223c;4kb insertion locus overlapping the HBB promoter among individuals with high HbF. These results demonstrate the utility of combining a comprehensive reference genome with LRS in African populations to uncover genomic variation at disease-associated loci.

SNV↗

Long-read direct infrared sequencing of crude PCR products for prediction of resistance to HIV-1 reverse transcriptase and protease inhibitors.

Patients infected with human immunodeficiency virus type 1 (HIV-1) are being treated with a number of different combinations of antiretroviral compounds that target the essential viral enzymes reverse transcriptase and protease. Different sets of HIV-1 mutations that confer drug resistance have been well defined; they allow reasonable prediction of the drug sensitivity pattern from analysis of the HIV-1 genotype in vivo. Since periodical monitoring of genotypic resistance is expected to improve clinical management in a large number of infected patients, practical and cost-effective methods are highly desirable to set at least medium-scale sequencing in clinical diagnostic settings. We present a complete protocol for direct sequencing of HIV-1 reverse transcriptase and protease-coding regions. Features making the system amenable to routine clinical use include: 1. Highly robust presequencing steps (plasma RNA extraction, reverse transcription, and nested PCR); 2. Direct use of the crude unpurified PCR product as the sequencing template; and 3. Use of infrared-labeled sequencing primers consistently allowing long reads, thus obviating the need for sequencing of both DNA strands.

Base Sequence↗

Genome size estimation from long read overlaps.

MOTIVATION: Accurate genome size estimation is an important component of genomic analyses such as assembly and coverage calculation, though existing tools are primarily optimized for short-read data. RESULTS: We present LRGE, a novel tool that uses read-to-read overlap information to estimate genome size in a reference-free manner. LRGE calculates per-read genome size estimates by analysing the expected number of overlaps for each read, considering read lengths and a minimum overlap threshold. The final size is taken as the median of these estimates, ensuring robustness to outliers such as reads with no overlaps. Additionally, LRGE provides an expected confidence range for the estimate. We validate LRGE on a large, diverse bacterial dataset and confirm it generalizes to eukaryotic datasets. On bacterial genomes, LRGE outperforms k-mer-based methods in both accuracy and computational efficiency and produces genome size estimates comparable to those from assembly-based approaches, like Raven, while using significantly less computational resources. AVAILABILITY AND IMPLEMENTATION: Our method, LRGE (Long Read-based Genome size Estimation from overlaps), is implemented in Rust and is available as a precompiled binary for most architectures, a Bioconda package, a prebuilt container image, and a crates.io package as a binary (lrge) or library (liblrge). The source code is available at https://github.com/mbhall88/lrge and an archive at https://doi.org/10.5281/zenodo.17183812 under an MIT license.

Genome Size↗

k-mer-based Upstream Preprocessing of long reads for Isoform Discovery.

Eukaryotic genes can encode multiple protein isoforms based on alternative splicing of their transcribed regions. Most modern novel isoform discovery methods function by identifying and assembling exon splice junctions from an RNA-seq sample. However, splice junctions can only be accurately annotated with time-intensive dynamic programming alignment. This manuscript introduces KuPID, a method for preprocessing long RNA-seq reads with the goal of better identifying novel isoform transcripts. KuPID utilizes k-mer sketching as a prefilter to quickly pseudo-align reads to known reference isoforms. Full alignment need only then be applied to reads that are most relevant to isoform discovery. Not only does KuPID speed up the discovery pipeline, it also increases downstream accuracy by filtering out extraneous reads. KuPID preprocessing simultaneously increases the f1 accuracy of isoform discovery pipelines by up to 11.6 points while decreasing the runtime by a factor of 2-3&#xd7;;. An optional mode permits a KuPID sample to be paired with both isoform discovery and transcript quantification.

Journal Article↗