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 55 records · Page 3Linked to original sources

Utilization of long-read sequencing for the detection of structural rearrangements with AgileStructure.

MOTIVATION: Changes in genome organisation contribute to genetic disease when they disrupt gene function or regulation. Structural rearrangements may interrupt coding sequence or alter expression through promoter loss or gain, chromatin changes, copy-number variation, or disruption of short-range regulatory elements. Although short-read sequencing excels at detecting small variants, it performs poorly at resolving breakpoints of large rearrangements, especially in repetitive or low-complexity regions. Long-read sequencing overcomes these limitations, but analytical tools have not kept pace, making accurate identification and annotation of large structural variants challenging. RESULTS: We developed AgileStructure, a desktop application for locating and annotating large‑scale genomic rearrangements using aligned long‑read data. The software enables user‑guided exploration of breakpoint‑spanning reads, supporting accurate interpretation of complex events and filling a key gap in current structural variant analysis workflows. AVAILABILITY AND IMPLEMENTATION: Source code, binaries, user guide, and example aligned read data, are available on GitHub: https://github.com/msjimc/AgileStructure. An archived version is also available on Zenodo at https://doi.org/10.5281/zenodo.18610110.

Software↗

Long-read DNA sequencing resolves a rare case of alloimmune hemolysis mimicking autoimmune hemolysis.

BACKGROUND: Immune hemolytic anemia poses a significant challenge in transfusion medicine, as identification of underlying alloantibodies can be masked by warm and/or cold autoantibodies. This increases the risk of transfusing incompatible blood, which can precipitate or exacerbate hemolysis. Identifying alloantibodies in the presence of autoantibodies remains difficult with standard serologic and genotypic methods, often delaying accurate diagnosis and appropriate transfusion strategies. CASE REPORT: We describe a 63-year-old woman with autoimmune hemolytic anemia who suffered near-fatal hemolysis following transfusion. Despite extensive serologic and genotypic testing, the cause of her hemolytic transfusion reactions remained elusive. Given her clinical course and transfusion history, we hypothesized that her acute hemolytic transfusion reactions could be due to immune sensitization to a high-incidence RBC antigen. Research whole-genome long-read sequencing (LRS) revealed homozygosity for a rare KEL*02N.16 allele, consistent with a rare Ko phenotype, which was validated by Sanger sequencing. Retrospective serologic testing with Ko RBCs further confirmed alloimmunization within the Kell system. CONCLUSION: This case highlights the limitations of conventional serologic and genotypic methods in detecting rare blood group phenotypes, and emphasizes the diagnostic power of long-read sequencing in transfusion medicine. Early molecular testing in complex hemolytic cases can facilitate targeted transfusion strategies, reduce the risk of severe hemolysis, and improve patient outcomes. As sequencing technologies become more accessible, they have the potential to revolutionize blood group typing and alloimmunization risk assessment in clinical practice.

Humans↗

LD-PCR coupled to long-read direct sequencing: an approach for mutation detection in genes with compact genomic structures.

A number of techniques have been developed as primary screens to scan for DNA sequence variants, including denaturing gradient gel electrophoresis, denaturing high-performance liquid chromatography, single-strand conformation polymorphism and heteroduplex analysis. Variant alleles detected by these assays are subsequently characterised by DNA sequencing. Sequencing itself is rarely used as a primary screen because of labour intensity, cost, and, upon automation, occasional inaccuracy in identifying heterozygous sites. We have previously developed an approach based on coupling long-distance PCR (LD-PCR) to long-read direct sequencing to allow the detection of mutations in the approximately 1.1 kb exon 3 of MECP2. Our use of dye-labelled primers generated high-quality bi-directional sequence runs > 650 bp and allowed easy discrimination of heterozygous bases. We now describe the application of this approach to the detection of mutations in a considerably larger 6.35 kb LD-PCR fragment spanning 10 exons (exons 32-41) of the structurally complex, but genomically compact, TSC2 gene. In a blind analysis, 15/15 previously characterised mutations were successfully identified using seven overlapping bi-directional sequencing reactions. Our approach of long-read sequencing of long-distance PCR products may allow rapid sequencing of multiple exons of compact genes and may be appropriate as a highly sensitive primary screen for mutations.

Base Sequence↗

Long-read, high-coverage reference genome of the nymphalid butterfly Catonephele acontius (Nymphalidae: Biblidinae).

Catonephele acontius (Nymphalidae:Biblidinae:Epicalinii) is a butterfly species with a wide distribution across the Neotropics including the Amazon. Here, we present a long-read high-coverage reference genome for this species to serve as a genomic resource for future studies on Biblidinae butterflies, a group that is the subject of ongoing studies of seasonal adaptation under climate change. We used PacBio HiFi and IsoSeq reads to generate a highly contiguous and well-annotated reference genome. Five libraries were constructed, 4 using RNA from different tissues and 1 using high molecular weight (HMW) DNA from a wild-caught female. The DNA was sequenced using PacBio HiFi technology, and the RNA was sequenced using long read PacBio IsoSeq technology. About 20 Gb of raw HiFi data were generated and assembled to an initial size of 520.7 Mb (39 × homozygous coverage) in 90 contigs. The assembly was then polished and decontaminated into 40 contigs with an N50 of 19.927 Mb (BUSCO completeness: 99.0%; duplication: 0.5%; fragmentation: 0.7%; and missing: 0.3%). Final assembly size was 519.2 Mb. Repeats were annotated, showing that the genome consisted of 40.4% transposable elements. IsoSeq transcriptome data from antennae, leg, ovary, and digestive tissue was then used to structurally and functionally annotate gene models for the softmasked genome, uncovering ∼18,500 genes, with 70% of them given functional annotation. This reference assembly joins many published genomes in the Nymphalidae family but represents one of the first high-quality genomes from the Biblidinae subfamily. It provides a valuable resource to study the evolution of plastic and seasonal traits and will help investigate the genetic processes that may influence these species' responses to rapid climate change.

Animals↗

FuFiHLA: a tool for full-field HLA typing from long-read data.

MOTIVATION: Allele typing for Human Leukocyte Antigen (HLA) genes has many important clinical applications. Popular short-read typing can only accurately distinguish alleles at the coding sequence level, which potentially limit our understanding of the effect of variants in non-coding region. Long read data has been proved to be useful in typing HLA alleles in full resolution, but only a few tools are publicly available and with significant limitations in practical application. RESULTS: We developed FuFiHLA, a lightweight open-source software, to type HLA alleles. Currently it supports typing alleles of six HLA genes (HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, and HLA-DQB1) from long reads. Evaluation using 233 PacBio HiFi WGS samples from HPRC shows that FuFiHLA achieves 99.6% accuracy in the full field allele typing and QV as 51.8 for consensus allele sequence construction. Additional testing on four Nanopore R10 reads demonstrates slightly reduced accuracy in the fourth field. AVAILABILITY: FuFiHLA is available at https://github.com/jingqing-hu/FuFiHLA under MIT License.

Humans↗

Novel bacterial hosts and mobile genetic structure of tet(X) variants in tetracycline-contaminated aquatic environment uncovered by culture and long-read metagenomics.

Clinically important tigecycline (3rd-generation tetracycline) resistance tet(X) variants were inferred to have evolutionarily originated from environmental bacteria, and have been recognized among environment, human and animals. However, genetic basis for environmental proliferation and dissemination of tet(X) variants remains ambiguous. This study profiled tet(X) variants at gene, contig, isolate, and community levels in environmental community subjected to long-term stepwise increasing oxytetracycline (1st-generation tetracycline) or tigecycline pressure using long-term microcosm experiments, quantitative PCR, bacterial isolation, whole-genome sequencing, and Nanopore-based long-read metagenomics. We confirmed that both oxytetracycline and tigecycline enriched the abundance of tetracycline resistance genes especially oxytetracycline-enriched tet(X3). Unexpectedly diverse bacterial hosts and genetic structure of tet(X)-positive mobile elements in the environment microbiome were identified using bacterial isolation and long-read Nanopore metagenomics. Pseudomonas defluvii was first reported to carry tet(X3) in the chromosome, forming IS26-tet(X3)-res-ISCR2 circular intermediate to transfer between different DNA molecules. Database mining revealed similar mobile segments have prevailed among animal-derived Acinetobacter species. Unlike the widely reported ISCR2-mediated transfer of tet(X6), we identified a novel mobile multidrug transposon TnAs3 where tet(X6) and class 1 integron co-transferred as its passenger region. Mobile tet(X2)-ere(D)-aadS-erm(F)-blaOXA-347 segment was annotated in Runella, and co-occurrences of tet(X2) and ere(D), aadS, blaOXA-347 were also found in Flavobacterium, Arsenicibacter, Chryseobacterium and Pedobacter. Overall, tetracycline-contaminated aquatic microbiome harboured diverse mobile tet(X)-positive segments which have not yet been acquired by clinical pathogens, and thus served as the genetic pool of tet(X) variants together with indigenous bacterial hosts, especially the newly reported Pseudomonas defluvii. Reducing pollution of older-generation tetracyclines would be a proactive way to mitigate environmental evolution and possible clinical effects of tet(X) variants.

Metagenomics↗

Long-read sequencing to interrogate strain-level variation among adherent-invasive Escherichia coli isolated from human intestinal tissue.

Adherent-invasive Escherichia coli (AIEC) is a pathovar linked to inflammatory bowel diseases (IBD), especially Crohn's disease, and colorectal cancer. AIEC are genetically diverse, and in the absence of a universal molecular signature, are defined by in vitro functional attributes. The relative ability of difference AIEC strains to colonize, persist, and induce inflammation in an IBD-susceptible host is unresolved. To evaluate strain-level variation among tissue-associated E. coli in the intestines, we develop a long-read sequencing approach to identify AIEC by strain that excludes host DNA. We use this approach to distinguish genetically similar strains and assess their fitness in colonizing the intestine. Here we have assembled complete genomes using long-read nanopore sequencing for a model AIEC strain, NC101, and seven strains isolated from the intestinal mucosa of Crohn's disease and non-Crohn's tissues. We show these strains can colonize the intestine of IBD susceptible mice and induce inflammatory cytokines from cultured macrophages. We demonstrate that these strains can be quantified and distinguished in the presence of 99.5% mammalian DNA and from within a fecal population. Analysis of global genomic structure and specific sequence variation within the ribosomal RNA operon provides a framework for efficiently tracking strain-level variation of closely-related E. coli and likely other commensal/pathogenic bacteria impacting intestinal inflammation in experimental settings and IBD patients.

Animals↗

Novel insights into tomato leaf curl New Delhi virus introduction and evolution in Southeastern France using an advanced long-read sequencing workflow.

The Mediterranean population of tomato leaf curl New Delhi virus (ToLCNDV-ES) is characterized by a high genetic uniformity, distinguishing it from its Asian counterparts. ToLCNDV-ES is thought to have a monophyletic origin, likely resulting from a single recombination event, prior to its spread throughout the Mediterranean region. Following its first detection in southeastern France in 2020, ToLCNDV-ES re-emerged in France in 2022. Our analysis based on advanced long-read sequencing, circular DNA profiling, and phylogeny indicates both local persistence of French ToLCNDV-ES and multiple independent introduction events. Signatures of positive selection were identified in French ToLCNDV-ES populations, whereas no clear evidence of recombination was found. Bayesian time-structured phylogenetic analyses suggest that introductions in France occurred between 2018 and 2021 from the major ToLCNDV-ES clade, while several Italian ToLCNDV-ES isolates diverged prior to the virus introduction in the Mediterranean basin. Overall, this study demonstrates the value of an optimized long-read sequencing approach for resolving circular DNA virus diversity, and sheds light on the complex evolutionary history of ToLCNDV-ES in the Mediterranean Basin, particularly in southeastern France.

France↗

Evaluating detection of Histophilus somni immunoglobulin-binding protein A DR2 Fic: A species-specific gene target for recombinase polymerase amplification relative to long-read sequencing of respiratory samples from feedlot calves.

Histophilosis is an important cause of morbidity and mortality as well as antimicrobial use in feedlot cattle across North America. Detection of Histophilus somni by culture is challenging, and there is no standardized tool for distinguishing isolates that carry virulence factors most likely to contribute to disease. The DR2 repeat of H. somni-associated virulence factor 'immunoglobulin-binding protein A' (ibpA DR2) harbors a Fic domain that mediates host cell cytotoxicity and is essential for histophilosis. For rapid detection of ibpA DR2 in extracted DNA, we developed a real-time recombinase polymerase amplification (RPA) assay with a runtime of 24&#xa0;min at 39&#xa0;&#xb0;C. DNA from H. somni-RPA-positive respiratory swabs (n&#xa0;=&#xa0;73) was screened for ibpA DR2 using the novel RPA assay and long-read metagenomic sequencing, as well as nanopore whole-genome sequencing (WGS) of H. somni isolated from the same samples. IbpA DR2 was identified in 71% and 70% of tested samples using RPA and WGS, respectively, and in &#x2264;41% of samples using metagenomic sequencing. The likelihood of detection by RPA did not differ (OR 1.1, 95% CI (0.42, 2.9), P&#xa0;>&#xa0;0.99) from WGS; however, agreement between these assays was only fair (&#x3ba;&#xa0;=&#xa0;0.31). Conversely, RPA (OR 3.4, 95% CI (1.6, 8.2)) and WGS (OR 8.0, 95% CI (2.4, 42)) were more likely (P&#xa0;<&#xa0;0.001) to detect ibpA DR2 than metagenomic sequencing, likely reflecting limited coverage of H. somni by metagenomics. This study demonstrated that RPA and long-read WGS detected ibpA DR2 with similar frequencies in extracted DNA and H. somni isolates, respectively. Further testing of non-target isolates confirmed the analytical specificity of ibpA DR2 to H. somni. Further investigation of the diagnostic validity for RPA-based ibpA DR2 detection is required in a larger cohort of field samples, as a rapid screening tool for H. somni most likely to contribute to disease.

Animals↗

Optimized Amplicon Strategy for Long-Read Sequencing of the Chikungunya Virus Genome.

Chikungunya virus (CHIKV) is a positive-sense RNA alphavirus transmitted to humans primarily by Aedes aegypti and Aedes albopictus mosquitoes. Its global circulation and significant public health impact underscore the need to better understand the molecular mechanisms driving CHIKV pathogenesis and transmission. Although robust molecular biology methods exist for CHIKV genome sequencing, a major limitation for surveillance and research is the inability to determine whether two nucleotide variations co-occur within the same viral genome when they are separated beyond the span of typical short-read designs. Here, we describe an optimized approach for processing CHIKV RNA samples that generates large amplicons suitable for long-read nanopore sequencing. This protocol enables amplification of the complete CHIKV genome in only two or three amplicons and facilitates detection of co-occurring nucleotide variations across 4-7.5&#x2009;kb within the same molecule, thereby simplifying sequencing workflows and improving resolution in studies of viral evolution.

Chikungunya virus↗

Sex differences in MAGEL2 gene promoter methylation in high functioning autism - trends from a pilot study using nanopore Cas9 targeted long read sequencing.

BACKGROUND: MAGEL2 is an autism susceptibility gene whose deficiency has been associated with autism-related behaviors in animal models and in syndromic human autism spectrum disorders (ASDs) such as Schaaf-Yang syndrome, but has not been studied in the broader autism spectrum. Given the capabilities of long-read sequencing technologies, this pilot study used a targeted nanopore sequencing approach to simultaneously examine MAGEL2 DNA sequence and methylation in adults with high-functioning autism (HFA) compared to neurotypical controls (NC). METHODS: Using DNA extracted from peripheral blood, Cas9-targeted nanopore DNA sequencing was used to analyze MAGEL2, including its entire regulatory construct (chr15:23639316-23651466), for sequence variation and 5-methyl-cytosine (5mC) modification in a cohort of adults with HFA compared to sex- and age-matched NC. Given the known sex differences in ASD and MAGEL2 KO animal models, results were further analyzed by sex. RESULTS: 20 adults with HFA (10 males, 10 females) and 20 NC were included. While there were no overall differences in MAGEL2 DNA sequence and 5mC modification between HFA and NC, we found a significant difference in MAGEL2 gene promoter methylation between males and females with HFA and NC of both sexes, with HFA males tending to show hypomethylation in a 300&#xa0;bp long differentially methylated region (chr15:23647640-23647939) around the MAGEL2 transcription start site. CONCLUSIONS: In this pilot study utilizing nanopore Cas9 targeted DNA sequencing, significant sex-specific differences in MAGEL2 gene promoter methylation were identified in male adults with HFA in comparison to control groups, suggesting the potential for sex-specific epigenetic differences. However, further replication in larger cohorts is required to validate these findings.

Humans↗

Prospective clinical validation of targeted long-read sequencing for preimplantation genetic testing of &#x3b1;-thalassaemia.

BACKGROUND: Preimplantation genetic testing for monogenic disorders (PGT-M) can prevent transmission of severe &#x3b1;-thalassaemia, but conventional workflows remain limited by family-specific assay design for direct variant detection, dependence on additional family samples for haplotype construction, and labour-intensive multi-step procedures across several platforms. Targeted long-read sequencing-based PGT-M for &#x3b1;-thalassaemia (tlrPGT-&#x3b1;-thal) integrates direct variant detection and haplotype linkage analysis within a single assay, but prospective clinical validation is lacking. METHODS: This prospective clinical study enrolled 103 families at high risk of transmitting &#x3b1;-thalassaemia at a reproductive medicine centre between August 2024 and March 2025. All families underwent blinded parallel analysis using both conventional NGS-based PGT-M (comparator) and tlrPGT-&#x3b1;-thal. RESULTS: In the primary concordance analysis, tlrPGT-&#x3b1;-thal was fully concordant with conventional NGS-based PGT-M (507/507, 100.0%; exact 95% CI, 99.3-100.0). Direct variant detection was successful in 501/507 embryos (98.82%; 95% CI, 97.4-99.6), haplotype linkage was established in 505/507 embryos (99.61%; 95% CI, 98.6-100.0), and one meiotic recombination event was identified. Among 93 families proceeding to embryo transfer, 57 pregnancies underwent invasive prenatal diagnosis, and all were concordant with the corresponding tlrPGT-&#x3b1;-thal results. Of the 26 comparator-inconclusive embryos, tlrPGT-&#x3b1;-thal resolved 6 complex cases, including cases with incomplete pedigrees or insufficient informative SNPs. Among the remaining 20 embryos with HBA-region aneuploidies, genotype and parental origin could be determined in 12. CONCLUSIONS: The findings show that tlrPGT-&#x3b1;-thal enables direct detection of diverse &#x3b1;-thalassaemia-causing variants together with efficient haplotype linkage analysis within a single workflow, without requiring family-specific assay design or additional family samples. The method demonstrated high diagnostic accuracy while providing added value in complex scenarios. Taken together, tlrPGT-&#x3b1;-thal represents a simplified and broadly applicable strategy for &#x3b1;-thalassaemia PGT-M.

Humans↗

An isoform-resolution transcriptomic atlas of colorectal cancer from long-read single-cell sequencing.

Colorectal cancer (CRC) ranks as the second leading cause of cancer deaths globally. In recent years, short-read single-cell RNA sequencing (scRNA-seq) has been instrumental in deciphering tumor heterogeneities. However, these studies only enable gene-level quantification but neglect alterations in transcript structures arising from alternative end processing or splicing. In this study, we integrated short- and long-read scRNA-seq of CRC samples to build an isoform-resolution CRC transcriptomic atlas. We identified 394 dysregulated transcript structures in tumor epithelial cells, including 299 resulting from various combinations of splicing events. Second, we characterized genes and isoforms associated with epithelial lineages and subpopulations exhibiting distinct prognoses. Among 31,935 isoforms with novel junctions, 330 were supported by The Cancer Genome Atlas RNA-seq and mass spectrometry data. Finally, we built an algorithm that integrated novel peptides derived from open reading frames of recurrent tumor-specific transcripts with mass spectrometry data and identified recurring neoepitopes that may aid the development of cancer vaccines.

Colorectal Neoplasms↗

Identifying inversions with breakpoints in the Dystrophin gene through long-read sequencing: report of two cases.

BACKGROUND: Duchenne Muscular Dystrophy (DMD) is an X-linked disorder caused by mutations in the DMD gene, with large deletions being the most common type of mutation. Inversions involving the DMD gene are a less frequent cause of the disorder, largely because they often evade detection by standard diagnostic methods such as multiplex ligation probe amplification (MLPA) and whole exome sequencing (WES). CASE PRESENTATION: Our research identified two intrachromosomal inversions involving the dystrophin gene in two unrelated families through Long-read sequencing (LRS). These variants were subsequently confirmed via Sanger sequencing. The first case involved a pericentric inversion extending from DMD intron 47 to Xq27.3. The second case featured a paracentric inversion between DMD intron 42 and Xp21.1, inherited from the mother. In both cases, simple repeat sequences (SRS) were present at the breakpoints of these inversions. CONCLUSIONS: Our findings demonstrate that LRS is an effective tool for detecting atypical mutations. The identification of SRS at the breakpoints in DMD patients enhances our understanding of the mechanisms underlying structural variations, thereby facilitating the exploration of potential treatments.

Humans↗

CLN3 transcript complexity revealed by long-read RNA sequencing analysis.

BACKGROUND: Batten disease is a group of rare inherited neurodegenerative diseases. Juvenile CLN3 disease is the most prevalent type, and the most common pathogenic variant shared by most patients is the "1-kb" deletion which removes two internal coding exons (7 and 8) in CLN3. Previously, we identified two transcripts in patient fibroblasts homozygous for the 1-kb deletion: the 'major' and 'minor' transcripts. To understand the full variety of disease transcripts and their role in disease pathogenesis, it is necessary to first investigate CLN3 transcription in "healthy" samples without juvenile CLN3 disease. METHODS: We leveraged PacBio long-read RNA sequencing datasets from ENCODE to investigate the full range of CLN3 transcripts across various tissues and cell types in human control samples. Then we sought to validate their existence using data from different sources. RESULTS: We found that a readthrough gene affects the quantification and annotation of CLN3. After taking this into account, we detected over 100 novel CLN3 transcripts, with no dominantly expressed CLN3 transcript. The most abundant transcript has median usage of 42.9%. Surprisingly, the known disease-associated 'major' transcripts are detected. Together, they have median usage of 1.5% across 22 samples. Furthermore, we identified 48 CLN3 ORFs, of which 26 are novel. The predominant ORF that encodes the canonical CLN3 protein isoform has median usage of 66.7%, meaning around one-third of CLN3 transcripts encode protein isoforms with different stretches of amino acids. The same ORFs could be found with alternative UTRs. Moreover, we were able to validate the translational potential of certain transcripts using public mass spectrometry data. CONCLUSION: Overall, these findings provide valuable insights into the complexity of CLN3 transcription, highlighting the importance of studying both canonical and non-canonical CLN3 protein isoforms as well as the regulatory role of UTRs to fully comprehend the regulation and function(s) of CLN3. This knowledge is essential for investigating the impact of the 1-kb deletion and rare pathogenic variants on CLN3 transcription and disease pathogenesis.

Humans↗

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↗

Comparison of a long-read amplicon sequencing approach to short-read amplicons for microbiome analysis.

Most microbiome studies to date rely on sequencing short amplicons of the 16S rRNA gene on Illumina's platforms. Because of the short read length, sequences often can be identified reliably only to the family or genus levels. Long read sequencing with whole-length 16S rRNA sequencing can improve taxonomic resolution, but often only to the species level. StrainID is an alternative approach that amplifies a large segment of the ribosomal operon, including the entire 16S rRNA gene, internal transcribed spacer, and a portion of the 23S rRNA gene. This longer amplicon is designed to allow ribotype-level classification. Although studies have demonstrated the utility of StrainID for several sample types, it has not yet been validated for saliva. Here, we compared the performance of StrainID to short read amplicons with saliva samples as well as a synthetic mock DNA community and human and mouse fecal samples. Short reads were amplified with primer pairs appropriate for the corresponding sample type, and were classified with two different taxonomic databases. For both saliva and fecal samples, we found that StrainID performed similarly to short reads overall and demonstrated a key benefit with phylogenetic-based beta diversity tests and taxonomic classification. Our results further build on establishing StrainID as a valid method and specifically validate its use with saliva samples.

Journal Article↗