Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “PacBio”

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 109 records · Page 6Linked to original sources
PubMed2026

A chromosome-level genome of the Nicobar pigeon, Caloenas nicobarica.

The Nicobar pigeon (Caloenas nicobarica), the closest living relative of the extinct Dodo (Raphus cucullatus), is endemic to Southeast Asia with a fragmented distribution across numerous small islands. It suffers from habitat loss, hunting, and predation from invasive species, resulting in its classification as Near Threatened by the International Union for the Conservation of Nature. We have generated a haplotype-resolved and chromosome-level genome assembly of the Nicobar pigeon using a combination of PacBio HiFi long-read sequencing and Arima Hi-C chromatin interaction mapping. This assembly includes two haplotypes, each spanning approximately 1.2 Gb. Haplotype 1 has a contig N50 of 25.2 Mb and a scaffold N50 of 79.7 Mb, whereas haplotype 2 has a contig N50 of 24.7 Mb and a scaffold N50 of 107.9 Mb. As the first high-quality genome assembly of any bird in the Columbidae Indo-Pacific clade, this resource provides valuable insights for phylogenetic studies. Furthermore, the phylogenetic proximity of the Nicobar pigeon to the Dodo (R. cucullatus) and the Rodrigues Solitaire (Pezophaps solitaria) offers a unique opportunity to study these extinct species, making this assembly a critical resource for evolutionary studies. It also offers a unique model for studying genetic diversity, adaptation, and speciation in island environments. This genomic resource will not only enhance our understanding of the evolutionary history of the Nicobar pigeon but also serve as a valuable tool for future conservation efforts aimed at preserving this unique species and its fragile island ecosystem.

Animals↗
PubMed2026

Chromosome-Level Reference Genome of the Desert Night Lizard Xantusia vigilis.

We present a reference-quality genome assembly for the desert night lizard (Xantusia vigilis). The night lizards (Xantusiidae) are a family of small-bodied lizards found in North America (Xantusia), Central America (Lepidophyma), and Cuba (Cricosaura). The night lizard family has an independent evolutionary history of at least 80 million years from its sister taxa within Scincoidea. The Xantusiids have several unique ecological, behavioral and evolutionary characteristics. For instance, the family contains the only squamate species that form diploid, unisexual, parthenogenic lineages. In addition, most night lizards are viviparous and form stable kin groups that are maintained over multiple years, an unusual life history strategy among lizards. Combining PacBio long-read sequencing, Hi-C, and RNAseq data we developed a reference-quality genome for the desert night lizard, X. vigilis. We assembled a complete mitochondrion and ~ 2.2 Gb nuclear genome, with 20 scaffolds that correlate in size to the X. vigilis karyotype. In addition, we found that X. vigilis chromosome 1 aligns with gene content of both of macrochromosome 1 and microchromosome 9 from a genome assembly of a species in the sister family Cordylidae (Hemicordylus capensis).

Xantusia↗
PubMed2026

Improved Genomic Resources for the swordtail cricket, Laupala kohalensis Otte 1994.

Advances in genetic tools such as next and third generation sequencing, paired with a focus on representative clades, provide insight into how processes including adaptation, admixture, and genome structure shape the evolution and maintenance of species. However, our understanding of the genomics of speciation is dominated by systems where ecological adaptations are thought to cause initial barriers to gene exchange. In contrast to other model systems, the 38 species of the genus Laupala constitute a very rapid radiation, where evolution of reproductive barriers and speciation is thought to be driven by sexual selection. Here, with novel PacBio HiFi reads and RNA- and Iso-Seq data, we provide a highly contiguous, chromosome-level genome and markedly improved annotation of the endemic Hawaiian cricket, Laupala kohalensis Otte, 1994. Our new resources advance previous efforts, placing 99% of 47 scaffolds on 7 autosomes and 1 sex chromosome in the 1.67 Gb assembly, with a 98.8% BUSCO score (insecta_db10), N50 of ~268 Mb, and L50 of 3. Using a custom repeat library, we estimate the genome to have 46.09% repeat content, and the new annotation includes an increased estimate of 17,670 genes, which coincides with that known from other Orthopterans. Notably, we find a large nuclear DNA segment of mitochondrial origin on chromosome 7. This new resource provides a powerful tool to identify and compare genomic causes of phenotypic diversification in a system characterized by strong signatures of sexual differentiation, representing an underappreciated but potentially widespread cause of speciation.

Hawaii↗
PubMed2024

DeepSomatic: Accurate somatic small variant discovery for multiple sequencing technologies.

Somatic variant detection is an integral part of cancer genomics analysis. While most methods have focused on short-read sequencing, long-read technologies now offer potential advantages in terms of repeat mapping and variant phasing. We present DeepSomatic, a deep learning method for detecting somatic SNVs and insertions and deletions (indels) from both short-read and long-read data, with modes for whole-genome and exome sequencing, and able to run on tumor-normal, tumor-only, and with FFPE-prepared samples. To help address the dearth of publicly available training and benchmarking data for somatic variant detection, we generated and make openly available a dataset of five matched tumor-normal cell line pairs sequenced with Illumina, PacBio HiFi, and Oxford Nanopore Technologies, along with benchmark variant sets. Across samples and technologies (short-read and long-read), DeepSomatic consistently outperforms existing callers, particularly for indels.

Journal Article↗
PubMed2025

Perplexity as a Metric for Isoform Diversity in the Human Transcriptome.

Long-read sequencing (LRS) has revealed a far greater diversity of RNA isoforms than earlier technologies, increasing the critical need to determine which, and how many, isoforms per gene are biologically meaningful. To define the space of relevant isoforms from LRS, many existing analysis pipelines rely on arbitrary expression cutoffs, but a single threshold cannot accommodate the broad variability in isoform complexity across genes, cell-types, and disease states captured by LRS. To address this, we propose using perplexity-an interpretable measure derived from entropy-that quantifies the effective number of isoforms per gene based on the full, unfiltered isoform ratio distribution. Calculating perplexity for 124 ENCODE4 PacBio LRS datasets spanning 55 human cell types, we show that it provides intuitive assessments of isoform diversity and captures uncertainty across genes with varying complexity. Perplexity can be calculated at multiple gene regulatory levels-from transcript to protein-to compare how isoform diversity is reduced across stages of gene expression. On average, genes have an ORF-level perplexity of 2.1, indicating production of two distinct protein isoforms. We extended this analysis to evaluate expression variation across tissues and identified 4,593 ORFs across 3,102 genes with moderate to extreme tissue-specificity. We propose perplexity as a consistent, quantitative metric for interpreting isoform diversity across genes, cell types, and disease states. All results are compiled into a community resource to enable cross-study comparisons of novel isoforms.

Journal Article↗
PubMed2025

Genome Report: De novo genome assembly of the greater Bermuda land snail, Poecilozonites bermudensis (Mollusca: Gastropoda), confirms ancestral genome duplication.

Poecilozonites bermudensis, the greater Bermuda land snail, is a critically endangered species and one of only two extant members in its genus. These snails are one of Bermuda's few endemic animal clades and their rich fossil record was the basis for the punctuated equilibria model of speciation. Once thought extinct, recent conservation efforts have focused on the recovery of the species, yet no genomic information or other molecular sequences have been available to inform these initiatives. We present a high-quality, annotated genome for P. bermudensis generated using PacBio long read and Omni-C short read sequencing. The resulting assembly is approximately 1.36 Gb with a scaffold N50 of 44.t Mb and 31 chromosome-length scaffolds. Nearly 43 percent of the genome was identified as repeat content. This assembly will serve as a resource for the conservation and study of P. bermudensis, and its only close extant and also critically endangered relative, P. circumfirmatus. Additionally, this genome adds to the growing body of data needed for a more complete understanding of gastropod evolution and for evolutionary processes in general.

Annotation↗
PubMed2025

Germline stem cell isolation, lineage tracing, and aging in a protochordate.

Germline stem cells (GSCs), the source of gametes, are the only stem cells capable of passing genes to future generations and are therefore considered units of natural selection. Yet, the factors that influence GSC fitness, and thus govern GSC competition, which exist in both protochordates and mammals, remain poorly understood. We studied how aging affects GSC fitness in the protochordate Botryllus schlosseri, an evolutionary crosspoint between invertebrates and vertebrates. GSCs were isolated and distinguished from developing and mature gametes using flow cytometry and scRNA-Seq, facilitated by a new PacBio genome assembly. Moreover, their function was validated through a novel lineage tracing approach that combines membrane-labeled GSC transplantation with scRNA-Seq. Leveraging our method to isolate them, single-cell transcriptomics showed significant age-related changes between young and old GSCs. Spermatids and sperm, however, showed minimal changes, suggesting that reproductive aging is governed by GSCs rather than by gametes. Reduced expressions of markers like DDX4 and PIWIL1 in aged GSCs mirrored trends in mammalian datasets, pointing to a conserved GSC-driven aging mechanism across chordate evolution. This study provides new techniques that lay the foundation to investigate further drivers of GSC fitness and highlights fertility-related genes as promising targets for therapies to preserve reproductive health.

Journal Article↗
PubMed2025

High-resolution metagenome assembly for modern long reads with myloasm.

Long-read metagenome assembly promises complete genomic recovery from microbiomes. However, the complexity of metagenomes poses challenges. We present myloasm, a metagenome assembler for PacBio HiFi and Oxford Nanopore Technologies (ONT) R10.4 long reads. Myloasm uses polymorphic k-mers to construct a high-resolution string graph and then leverages differential abundance for graph simplification. On real-world ONT metagenomes, myloasm assembled three times more complete circular contigs than the next-best assembler. Myloasm can make ONT and HiFi comparable for assembly: for a jointly sequenced gut metagenome, myloasm with ONT assembled more complete circular genomes than any assembler with HiFi. Myloasm recovers previously inaccessible within-species diversity; we recovered six complete Prevotella copri single-contig genomes from a gut metagenome and eight complete TM7 (Saccharibacteria) contigs with > 93% similarity from an oral metagenome. With this improved resolution, we resolved two 98% similar ermF antibiotic resistance genes spreading through distinct strain-specific mobile genetic elements in a human gut.

Journal Article↗
PubMed2025

Comprehensive benchmarking of somatic structural variant detection at ultra-low allele fractions.

Postzygotic mosaicism gives rise to somatic structural variants (SVs) at ultra-low variant allele fractions (VAFs), which pose challenges for detection due to the high-coverage sequencing required and noise introduced by sequencing artifacts. Although somatic SV detection has been extensively studied in cancer, these studies are not directly applicable to the study of tissue mosaicism, as they rely on matched normals, target higher VAF ranges, and are enriched for different types of SVs. We present comprehensive benchmark data and best practices for non-cancer somatic SV detection. We created a synthetic mosaic sample by combining six HapMap individuals at varying proportions, generating allele fractions as low as 0.25%. This sample was sequenced to ~2,300x total coverage using Illumina, PacBio, and Nanopore technologies across multiple sequencing centers. A high-confidence benchmark SV set containing over 21,000 pseudo-somatic insertions and deletions ≥50bp was derived from haplotype-resolved assemblies. We evaluated 12 SV discovery pipelines and identified caller-specific strengths and sequencing platform-specific shortcomings. We find that short read-based approaches show reduced recall for insertions and repeat-associated SVs, whereas long-read sequencing achieves high accuracy throughout the genome, increasing linearly with coverage. The best algorithm's sensitivity exceeded 80% for VAFs ≥4% and 15% for VAFs of 0.5-1% with 60x coverage. The publicly available benchmarking data and comparative analysis of current methods provide a foundation for robust discovery of SV mosaicism in non-cancer tissues..

Journal Article↗
PubMed2025

The UTRs of Leishmania donovani vary in length and are enriched in potential regulatory structures.

Leishmania spp. regulate gene expression largely post-transcriptionally, yet untranslated regions (UTRs) remain poorly delineated. We generated high-quality genome and transcriptome datasets for Leishmania donovani strain 1S2D (Ld1S) by combining PacBio HiFi de novo assembly with Oxford Nanopore direct RNA sequencing of promastigotes and axenic amastigotes. The genome assembly consists of 65 scaffolds totaling ~33.3 Mb. Structural comparisons to LdBPK282A1 revealed numerous rearrangements, including some reshuffling genes among polycistronic transcription units and validated by polycistronic reads from RNA sequencing. Promastigote and amastigote RNA sequencing produced 469,010 and 46,729 monocistronic reads containing a spliced-leader and a polyA tail sequences, defining 8,479 transcripts and supporting 7,415 of the 7,969 annotated protein coding genes, as well as 604 putative long non-coding RNAs. We annotated UTRs for 4,921 genes and observed that putative RNA G-quadruplexes were markedly enriched in UTRs. We also noted that 31.9% and 11.5% were expressed into multiple isoforms in promastigotes and amastigotes, respectively. Collectively, these data provide a comprehensive annotation of L. donovani genes and their UTRs and reveal widespread and stage-specific UTR length polymorphisms, and, overall, points to an important role of 3' UTR in post-transcriptional regulation in L. donovani.

Journal Article↗
PubMed2026

Assessing the readiness of Oxford Nanopore sequencing for clinical genomics applications.

Long-read sequencing (LRS) technologies, namely, Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), have emerged as promising solutions to overcome the limitations of short-read sequencing (SRS). Nevertheless, the still higher sequencing error rates compared with SRS, need for customized pipelines, rapidly updating software, and incipient scalability continue to present challenges for adopting ONT in standard clinical practice. Here we assess the performance of ONT (R9 and R10 chemistries) in comparison to Illumina and MGI across 17 well-characterized reference samples with 11 clinical variants representing nine different genetic diseases. To enable this, we have implemented a production-ready pipeline including SNV, indel, STR, SV, and CNV detection, alongside reporting key summary metrics to ensure high-quality data at the production sequencing level. Our results show high accuracy of ONT across SNVs (F-score 0.978-0.983) and SVs (F-score = 0.75) but still weaknesses across indels (F-score 0.659-0.758). However, we highlight that ONT accurately detected all four pathogenic indels as well as the performance improvement in exons and with the newer R10 chemistry. We further demonstrated the importance of long reads to detect clinically impactful variants such as a FMR1 pathogenic expansion, often misclassified by SRS as being in the premutation range. Our multiplatform analysis and Sanger validation uncovered a 1 bp error in the Coriell annotation for a cystic fibrosis-causing indel in GM07829. This work underscores the growing readiness of ONT for clinical applications, highlighting both its advancements and its potential for broader adoption in clinical genomics and large-scale operations.

Humans↗
PubMed2026

Post-Hoc Long-Read Sequencing Links Leukemic Mutation Status to Single-Cell Transcriptomes.

Single-cell RNA-sequencing-based characterization of cells that belong to the neoplastic clone is a major challenge in hematologic neoplasms, where malignant and normal cells coexist. Confident molecular profiling requires simultaneous analysis of gene expression and genetic mutations in individual cells, an ability that is not supported by the standard 10X Genomics workflow. Here, we systematically evaluated the potential and limitations of repurposing amplified cDNA generated during the 10X Genomics 3' workflow for post hoc genotyping of individual cells. We first established a mixed leukemic cell line system comprising one cell line with KIT point mutations and another with the BCR::ABL1 fusion gene. Targeted long-read PacBio sequencing enabled post hoc assignment of mutation data to transcriptionally profiled cells, but recovery differed between targets. Consistent with ambient RNA in microfluidics-based single-cell workflows, mutation-associated transcripts were detected in cells not expected to carry the corresponding mutations, illustrating how transcript recovery complicates cell-level genotype assignment. Target-specific thresholds mitigated this source of misclassification. In primary chronic myeloid leukemia samples, the post hoc approach detected BCR::ABL1-positive cells at diagnosis, but not during imatinib treatment. Together, we present a framework for adding mutation status to cells already profiled using the 10X Genomics workflow and highlight broader considerations for transcript-based single-cell genotyping.

BCR::ABL1↗
PubMed2026

Integrative Long-Read Multi-Omics of a Patient With GPI Deficiency: A Molecular Case Study of a Candidate Dual-Effect GPI Variant.

The molecular determinants of phenotypic severity in red cell enzymopathies are often obscured by the disconnect between coding sequence variants and their regulatory landscapes. Here we present a single-patient molecular case study that uses an integrative multi-omic approach-combining short-read WGS, PacBio HiFi long-read sequencing, native CpG methylation profiling, and Iso-Seq full-length transcriptomics-to characterize a severe, transfusion-dependent hemolytic anaemia. We identified a compound heterozygous state in the glucose-6-phosphate isomerase (GPI) gene, with no wild-type allele present. One allele (Haplotype 1) carried a missense variant (p.His191Arg); the other (Haplotype 2) carried a distinct missense variant, c.1414C>T (p.Arg472Cys), previously reported as biochemically unstable. Long-read phasing placed the two variants in trans. Allele-resolved transcript counts showed a directionally consistent but statistically non-significant trend toward higher expression of Haplotype 2 across two Iso-Seq replicates. Notably, the c.1414C>T transition abolishes a local CpG dinucleotide; in a small number of haplotype-2 reads spanning this position, the corresponding cytosine on the wild-type/Haplotype-1 background was methylated. We did not measure GPI protein abundance, enzymatic activity, or stability in this patient, and we do not establish that methylation at this site regulates GPI transcription. On the basis of these correlative observations in a single patient, we propose-as a hypothesis for future testing-that a coding variant might simultaneously perturb protein stability and disrupt a local epigenetic mark, and we outline the experiments required to test whether such a dual effect contributes to disease. This case illustrates the value of integrative long-read multi-omics for generating mechanistic hypotheses about variants of uncertain significance, while underscoring that causal claims require dedicated functional validation.

Humans↗
PubMed2025

Genetic Adaptation to Brackish Water and Spawning Season in European Cisco.

How species adapt to diverse environmental conditions is essential for understanding evolution and the maintenance of biodiversity. The European cisco (Coregonus albula) is a salmonid that occurs in both fresh and brackish water, and this together with the presence of sympatric spring- and autumn-spawning lacustrine populations provides an opportunity for studying the genetics of adaptation in relation to salinity and timing of reproduction. Here, we present a high-quality reference genome of the European cisco based on PacBio HiFi long read sequencing and HiC-directed scaffolding. We generated low-coverage whole-genome sequencing data from 336 individuals across 12 population samples to explore population structure and genetics of ecological adaptation. We found a major subdivision between two groups of populations most likely reflecting colonisation from different glacial refugia. Within the two major groups, we detected further genetic differentiation between spring- and autumn-spawning populations and between populations from freshwater lakes, rivers and brackish water (Bothnian Bay). A genome-wide screen for genetic differentiation among populations identified a set of outlier SNPs strongly correlated with spawning timing and salinity. Several of the genes associated with spawning time, including BHLHE40, TIMELESS and CPT1A, have previously been shown to have a role in circadian rhythm biology. As many as 17 loci were associated with genetic differentiation between populations reproducing in fresh and brackish water. This study provides insights into the genomic basis of ecological adaptation in European cisco with implications for sustainable fishery management.

Animals↗
PubMed2026

Draft genome sequence of the almond red leaf blotch pathogen Polystigma amygdalinum assembled from infected almond leaves collected in California, USA.

We report a draft genome assembly of Polystigma amygdalinum, the causal agent of almond red leaf blotch. DNA extracted from infected leaves was sequenced using PacBio HiFi, and host-derived reads were removed bioinformatically. The 238.7-Mb assembly (90.8% BUSCO completeness) is highly repetitive (82.3%) and unusually large for an ascomycete.

Polystigma amygdalinum↗
PubMed2026

Complete genome sequence of an MDR and hypermucoviscous Klebsiella pneumoniae strain KP99 isolated from a 62-year-old patient in Dalian, China.

We employed a PacBio-Illumina hybrid sequencing approach to sequence the strain KP99 isolated from the sputum of a clinical patient. Classified as ST592, KP99 carries multiple resistance genes (including blaSHV-26, fosA5, and oqxA/B) and harbors a hypervirulent factor system, including the rmpA gene linked to CPS and siderophore-associated genes.

Klebsiella pneumoniae↗