PacBio Research Articles

Construction and evaluation of a new rat reference genome assembly, GRCr8, from long reads and long-range scaffolding

We report the construction and analysis of a new reference genome assembly forRattus norvegicus, the laboratory rat, a widely used experimental animal model organism. The assembly has been adopted as the rat reference assembly by the Genome Reference Consortium and is named GRCr8. The assembly has employed 40× Pacific Biosciences (PacBio) HiFi sequencing coverage and scaffolding using optical mapping and Hi-C. We used genomic DNA from a male BN/NHsdMcwi (BN) rat of the same strain and from the same colony as the prior reference assembly, mRatBN7.2. The assembly is at chromosome level with 98.7% of the sequence assigned to chromosomes. All chromosomes have increased in size compared with the prior assembly andk-mer analysis indicates that the subject animal is fully inbred and that the genome is represented as a single haploid assembly. Notable increases are observed in Chromosomes 3, 11, and 12 in the prospective rDNA regions. In addition, Chr Y has increased threefold in size and is more consistent with the rat karyotype than previous assemblies. Several other chromosomes have grown by the incorporation of sizable discrete new blocks. These contain highly repetitive sequences and encode numerous previously unannotated genes. In addition, centromeric sequences are incorporated in most chromosomes. Genome annotation has been performed by NCBI RefSeq, which confirms improvement in assembly quality and adds more than 1100 new protein coding genes. PacBio Iso-Seq data have been acquired from multiple tissues of the subject animal and are released concurrently with the new assembly to aid further analyses.

Gapless assembly of complete human and plant chromosomes using only nanopore sequencing.

The combination of ultra-long (UL) Oxford Nanopore Technologies (ONT) sequencing reads with long, accurate Pacific Bioscience (PacBio) High Fidelity (HiFi) reads has enabled the completion of a human genome and spurred similar efforts to complete the genomes of many other species. However, this approach for complete, "telomere-to-telomere" genome assembly relies on multiple sequencing platforms, limiting its accessibility. ONT "Duplex" sequencing reads, where both strands of the DNA are read to improve quality, promise high per-base accuracy. To evaluate this new data type, we generated ONT Duplex data for three widely studied genomes: human HG002, Solanum lycopersicum Heinz 1706 (tomato), and Zea mays B73 (maize). For the diploid, heterozygous HG002 genome, we also used "Pore-C" chromatin contact mapping to completely phase the haplotypes. We found the accuracy of Duplex data to be similar to HiFi sequencing, but with read lengths tens of kilobases longer, and the Pore-C data to be compatible with existing diploid assembly algorithms. This combination of read length and accuracy enables the construction of a high-quality initial assembly, which can then be further resolved using the UL reads, and finally phased into chromosome-scale haplotypes with Pore-C. The resulting assemblies have a base accuracy exceeding 99.999% (Q50) and near-perfect continuity, with most chromosomes assembled as single contigs. We conclude that ONT sequencing is a viable alternative to HiFi sequencing for de novo genome assembly, and provides a multirun single-instrument solution for the reconstruction of complete genomes.

A genome assembly for the California endemic liverwort Calasterella californica.

Calasterella californica belongs to a monotypic genus of liverworts endemic to the west coast of North America, primarily distributed in California. This dioicous species occurs in a variety of ecosystems from deserts to redwood forest; little is known about how this species is adapted to live in those seemingly contrasting environments. In this paper, we report the assembly of the nuclear genome of Calasterella californica. As part of the California Conservation Genomics Project (CCGP), we used Pacific Biosciences HiFi long-read sequencing data to produce a de novo assembly that consists of 772 contigs, with a total length of 517 Mbp and a BUSCO complete score of 95%. C. californica is only the sixth species of liverworts - a group with more than 7200 described species - to have a nuclear reference genome. The availability of this reference genome will facilitate the study of the unique features of C. californica and other liverworts, pave the road towards a comparative understanding of liverwort genomes, and add an important starting point for studies of the geographic variation of this species within the CCGP project.

A Chromosome-Level Genome Assembly of Chiton Acanthochiton rubrolineatus (Chitonida, Polyplacophora, Mollusca)

(1) Background: Chitons (Mollusca, Polyplacophora) are relatively primitive species in Mollusca that allow the study of biomineralization. Although mitochondrial genomes have been isolated from Polyplacophora, there is no genomic information at the chromosomal level; (2) Methods: Here we report a chromosome-level genome assembly for Acanthochiton rubrolineatus using PacBio (Pacific Biosciences, United States) reads and high-throughput chromosome conformation capture (Hi-C) data; (3) Results: The assembly spans 1.08 Gb with a contig N50 of 3.63 Mb and 99.97% of the genome assigned to eight chromosomes. Among the 32,291 predicted genes, 76.32% had functional predictions. The divergence time of Brachiopoda and Mollusca was ~550.8 Mya (million years ago), and that of A. rubrolineatus and other mollusks was ~548.5 Mya; (4) Conclusions: This study not only offers high-quality reference sequences for the Acanthochiton rubrolineatus genome, but also establishes groundwork for investigating the mechanisms of Polyplacophora biomineralization and its evolutionary history. This research will aid in uncovering the genetic foundations of molluscan adaptations across diverse environments.

Decoding the Genetic Code: Scientific Exploration of the Telomere-totelomere (T2T) Genome

Abstract: With the development of third-generation sequencing technology, genome assembly has entered a new era. The combination of multiple sequencing methods can combine the strengths of each, resulting in higher-quality assembly results. Telomere-to-Telomere (T2T) genome refers to a zero gap genome that is assembled at the T2T level of one or more chromosomes through the combination of multiple sequencing technologies, such as Pacific Biosciences High-Fidelity (PacBio HiFi), Oxford Nanopore Technologies (ONT) Ultra-long, High-throughput Chromosome Conformation Capture (Hi-C), and others. High-quality reference genomes are the basis for genomics research, and the T2T genome has enabled the exploration of unknown areas of the genome, such as telomeres and centromeres, which has provided a more in-depth direction for research. This review provides a comprehensive overview of the T2T genome, reviewing the development of sequencing technologies and then outlining sequencing strategies, assembly methods, quality assessment processes, and analysis software for the T2T genome with practical examples. Representative T2T genomic species of plants, animals, and microorganisms (e.g., human, Arabidopsis, brewer's yeast, and so on) are presented separately. A summary of the potential and challenges of current T2T genomic applications is provided, along with an outlook for future developments.

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements.

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth-resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into the Atabeyarchaeum deiterrae-1 chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

Characterizing tandem repeat complexities across long-read sequencing platforms with TREAT and otter.

Tandem repeats (TRs) play important roles in genomic variation and disease risk in humans. Long-read sequencing allows for the accurate characterization of TRs; however, the underlying bioinformatics perspectives remain challenging. We present otter and TREAT: otter is a fast targeted local assembler, cross-compatible across different sequencing platforms. It is integrated in TREAT, an end-to-end workflow for TR characterization, visualization, and analysis across multiple genomes. In a comparison with existing tools based on long-read sequencing data from both Oxford Nanopore Technology (ONT, Simplex and Duplex) and Pacific Bioscience (PacBio, Sequel II and Revio), otter and TREAT achieve state-of-the-art genotyping and motif characterization accuracy. Applied to clinically relevant TRs, TREAT/otter significantly identify individuals with pathogenic TR expansions. When applied to a case-control setting, we replicate previously reported associations of TRs with Alzheimer's disease, including those near or within APOC1 (P = 2.63 × 10-9), SPI1 (P = 6.5 × 10-3), and ABCA7 (P = 0.04) genes. Finally, we use TREAT/otter to systematically evaluate potential biases when genotyping TRs using diverse ONT and PacBio long-read sequencing data sets. We show that, in rare cases (0.06%), long-read sequencing from coverage drops in TRs, including the disease-associated TRs in ABCA7 and RFC1 genes. Such coverage drops can lead to TR misgenotyping, hampering the accurate characterization of TR alleles. Taken together, our tools can accurately genotype TRs across different sequencing technologies and with minimal requirements, allowing end-to-end analysis and comparisons of TRs in human genomes, with broad applications in research and clinical fields.

Global Invasion History and Genomic Signatures of Adaptation of the Highly Invasive Sycamore Lace Bug.

Invasive species cause massive economic and ecological damage. Climate change has resulted in an unprecedented increase in the number and impact of invasive species; however, the mechanisms underlying these invasions are unclear. The sycamore lace bug, Corythucha ciliata, is a highly invasive species originating from North America and has expanded across the Northern Hemisphere since the 1960s. In this study, we assembled the C. ciliata genome using high-coverage Pacific Biosciences (PacBio), Illumina, and high-throughput chromosome conformation capture (Hi-C) sequencing. A total of 15,278 protein-coding genes were identified, and expansions of gene families with oxidoreductase and metabolic activities were observed. In-depth resequencing of 402 samples from native and nine invaded countries across three continents revealed 2.74 million single nucleotide polymorphisms. Two major invasion routes of C. ciliata were identified from North America to Europe and Japan, with a contact zone forming in East Asia. Genomic signatures of selection associated with invasion and long-term balancing selection in native ranges were identified. These genomic signatures overlapped with expanded genes, suggesting improvements in the oxidative stress and thermal tolerance of C. ciliata. These findings offer valuable insights into the genomic architecture and adaptive evolution underlying the invasive capabilities of species during rapid environmental changes.

Barcode 100K Specimens: In a Single Nanopore Run.

It is a global priority to better manage the biosphere, but action must be informed by comprehensive data on the abundance and distribution of species. The acquisition of such information is currently constrained by high costs. DNA barcoding can speed the registration of unknown animal species, the most diverse kingdom of eukaryotes, as the BIN system automates their recognition. However, inexpensive sequencing protocols are critical as the census of all animal species is likely to require the analysis of a billion or more specimens. Barcoding involves DNA extraction followed by PCR and sequencing with the last step dominating costs until 2017. By enabling the sequencing of highly multiplexed samples, the Sequel platforms from Pacific BioSciences slashed costs by 90%, but these instruments are only deployed in core facilities because of their expense. Sequencers from Oxford Nanopore Technologies provide an escape from high capital and service costs, but their low sequence fidelity has, until recently, constrained adoption. However, the improved performance of its latest flow cells (R10.4.1) erases this barrier. This study demonstrates that a MinION flow cell can characterise an amplicon pool derived from 100,000 specimens while a Flongle flow cell can process one derived from several thousand. At $0.01 per specimen, DNA sequencing is now the least expensive step in the barcode workflow.

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.

Fine mapping of RNA isoform diversity using an innovative targeted long-read RNA sequencing protocol with novel dedicated bioinformatics pipeline.

Solving the structure of mRNA transcripts is a major challenge for both research and molecular diagnostic purposes. Current approaches based on short-read RNA sequencing and RT-PCR techniques cannot fully explore the complexity of transcript structure. The emergence of third-generation long-read sequencing addresses this problem by solving this sequence directly. However, genes with low expression levels are difficult to study with the whole transcriptome sequencing approach. To fix this technical limitation, we propose a novel method to capture transcripts of a gene panel using a targeted enrichment approach suitable for Pacific Biosciences and Oxford Nanopore Technologies platforms. We designed a set of probes to capture transcripts of a panel of genes involved in hereditary breast and ovarian cancer syndrome. We present SOSTAR (iSofOrmS annoTAtoR), a versatile pipeline to assemble, quantify and annotate isoforms from long read sequencing using a new tool specially designed for this application. The significant enrichment of transcripts by our capture protocol, together with the SOSTAR annotation, allowed the identification of 1,231 unique transcripts within the gene panel from the eight patients sequenced. The structure of these transcripts was annotated with a resolution of one base relative to a reference transcript. All major alternative splicing events of the BRCA1 and BRCA2 genes described in the literature were found. Complex splicing events such as pseudoexons were correctly annotated. SOSTAR enabled the identification of abnormal transcripts in the positive controls. In addition, a case of unexplained inheritance in a family with a history of breast and ovarian cancer was solved by identifying an SVA retrotransposon in intron 13 of the BRCA1 gene. We have validated a new protocol for the enrichment of transcripts of interest using probes adapted to the ONT and PacBio platforms. This protocol allows a complete description of the alternative structures of transcripts, the estimation of their expression and the identification of aberrant transcripts in a single experiment. This proof-of-concept opens new possibilities for RNA structure exploration in both research and molecular diagnostics.

Long-read genome sequencing and variant reanalysis increase diagnostic yield in neurodevelopmental disorders.

Variant detection from long-read genome sequencing (lrGS) has proven to be more accurate and comprehensive than variant detection from short-read genome sequencing (srGS). However, the rate at which lrGS can increase molecular diagnostic yield for rare disease is not yet precisely characterized. We performed lrGS using Pacific Biosciences HiFi technology on 96 short-read-negative probands with rare diseases that were suspected to be genetic. We generated hg38-aligned variants and de novo phased genome assemblies, and subsequently annotated, filtered, and curated variants using clinical standards. New disease-relevant or potentially relevant genetic findings were identified in 16/96 (16.7%) probands, nine of which (8/96, ~9.4%) harbored pathogenic or likely pathogenic variants. Nine probands (~9.4%) had variants that were accurately called in both srGS and lrGS and represent changes to clinical interpretation, mostly from recently published gene-disease associations. Seven cases included variants that were only correctly interpreted in lrGS, including copy-number variants, an inversion, a mobile element insertion, two low-complexity repeat expansions, and a 1 bp deletion. While evidence for each of these variants is, in retrospect, visible in srGS, they were either not called within srGS data, were represented by calls with incorrect sizes or structures, or failed quality-control and filtration. Thus, while reanalysis of older srGS data clearly increases diagnostic yield, we find that lrGS allows for substantial additional yield (7/96, 7.3%) beyond srGS. We anticipate that as lrGS analysis improves, and as lrGS datasets grow allowing for better variant frequency annotation, the additional lrGS-only rare disease yield will grow over time.

Assembly and annotation of a chromosome-level reference genome for the endangered Colorado pikeminnow (Ptychocheilus lucius).

Advancements in genome sequencing technology have brought unprecedented accessibility of high-throughput sequencing to species of conservation interest. The potential knowledge gained from application of these techniques is maximized by availability of high-quality, annotated reference genomes for endangered species. However, these vital resources are often lacking for endangered minnows of North America (Cypriniformes: Leuciscidae). One such endangered species, Colorado pikeminnow (Ptychocheilus lucius), is the largest North American minnow and the top-level native aquatic predator in the Colorado River Basin of the southwestern United States and northwestern Mexico. Over the past century, Colorado pikeminnow has suffered habitat loss and population declines due to anthropogenic habitat modifications and invasive species introductions. The lack of genetic resources for Colorado pikeminnow has hindered conservation genomic study of this unique organism. This study seeks to remedy this issue by presenting a high-quality reference genome for Colorado pikeminnow developed from Pacific Biosciences HiFi sequencing and Hi-C scaffolding. The final assembly was a 1.1 Gb genome comprised of 305 contigs including 25 chromosome-sized scaffolds. Measures of quality, contiguity, and completeness met or exceeded those observed for Danio rerio (Danionidae) and 2 other Colorado River Basin leuciscids (Meda fulgida and Tiaroga cobitis). Comparative genomic analyses identified enrichment of gene families for growth, development, immune activity, and gene transcription; all of which are important for a large-bodied piscivorous fish living in a dynamic environment. This reference genome will provide a basis for important conservation genomic study of Colorado pikeminnow and help efforts to better understand the evolution of desert fishes.

De novo Whole-Genome Assembly of the 10-Gigabase Fokienia Hodginsii Genome to Reveal Differential Epigenetic Events Between Callus and Xylem.

Fokienia hodginsii (F. hodginsii), belonging to the genus Fokienia of the Cupressaceae. F. hodginsii has significant application value due to its wood properties and great research value in evolutionary studies as a gymnosperm. However, the genome of F. hodginsii remains unknown due to the large size of gymnosperms genome. Pacific Bioscience sequencing, Hi-C mapping, whole-genome Bisulfite Sequencing (BS-Seq), long-read isoform sequencing (Iso-Seq), direct RNA sequencing (DRS), quantitative proteomics, and metabonomics analysis are employed to facilitate genome assembly, gene annotation, and investigation into epigenetic mechanisms. In this study, the 10G F. hodginsii genome is assembled into 11 chromosomes. Furthermore, 50521 protein-coding genes are annotated and determined that 65% of F. hodginsii genome comprises repetitive sequences. It is discovered that transposable element (TE)-including introns is associated with higher expression. The DNA methylome of F. hodginsii reveals that xylem has a higher DNA methylation level compared to callus. Moreover, DRS reveals the significant alterations in RNA full-length ratio, which potentially associated with poly(A) length (PAL) and alternative polyadenylation (APA). Finally, the morphology measurement and metabonomics analysis revealed the difference of 14 cultivars. In summary, the genomes and epigenetics datasets provide a molecular basis for callus formation in the gymnosperm family.

Annotation and Characterization of the Zacco platypus Genome

The pale chub Zacco platypus (Cypriniformes; Xenocyprididae; Jordan & Evermann, 1902) is widely distributed across freshwater ecosystems in East Asia and has been recognized as a potential model fish species for ecotoxicology and environmental monitoring. Here, a high-quality de novo genome assembly of Z. platypus was constructed through the integration of a combination of long-read Pacific Bioscience (PacBio) sequencing, short-read Illumina sequencing, and Hi-C sequencing technologies. Z. platypus has the smallest genome size compared to other species belonging to the order Cypriniformes. The assembled genome encompasses 41.45% repeat sequences. As shown in other fish, a positive correlation was observed between genome size and the composition of transposable elements (TE) in the genome. Among TEs, a relatively higher rate of DNA transposon was observed, which is a common pattern in the members of the order Cypriniformes. Functional annotation was processed using four representative databases, identifying a core set of 12,907 genes shared among them. Orthologous gene family analysis revealed that Z. platypus has experienced more gene family contraction rather than expansion compared to other Cypriniformes species. Among the uniquely expanded gene families in Z. platypus, detoxification and stress-related gene families were identified, suggesting that this species could represent a promising model for ecotoxicology and environmental monitoring. Taken together, the Z. platypus genome assembly will provide valuable data for omics-based health assessments in aquatic ecosystems, offering further insights into the environmental and ecological facets within this species.

Discrete measurements of RNA polymerase and reverse transcriptase fidelity reveal evolutionary tuning.

Direct methods for determining the fidelity of DNA polymerases are robust, with relatively little sample manipulation before sequencing. In contrast, methods for measuring RNA polymerase and reverse transcriptase fidelities are complicated by additional preparation steps that introduce ambiguity and error. Here, we describe a sequencing method, termed Roll-Seq, for simultaneously determining the individual fidelities of RNA polymerases and reverse transcriptases (RT) using Pacific Biosciences single molecule real-time sequencing. By using reverse transcriptases with high rolling-circle activity, Roll-Seq generates long concatemeric cDNA from a circular RNA template. To discern the origin of a mutation, errors are recorded and determined to occur within a single concatemer (reverse transcriptase error) or all concatemers (RNA polymerase error) over the cDNA strand. We used Roll-Seq to measure the fidelities of T7 RNA polymerases, a Group II intron-encoded RT (Induro), and two LINE RTs (Fasciolopsis buski R2-RT and human LINE-1). Substitution rates for Induro and R2-RT are the same for cDNA and second-strand synthesis while LINE-1 has 2.5-fold lower fidelity when performing second-strand synthesis. Deletion and insertion rates increase for all RTs during second-strand synthesis. In addition, we find that a structured RNA template impacts fidelity for both RNA polymerase and RT. The accuracy and precision of Roll-Seq enable this method to be applied as a complementary analysis to structural and mechanistic characterization of RNA polymerases and reverse transcriptases or as a screening method for RNAP and RT fidelity.

Chromosome-scale reference genome of an ancient landrace: unveiling the genetic basis of seed weight in the food legume crop pigeonpea (Cajanus cajan).

Pigeonpea (Cajanus cajan) is a nutrient-rich and versatile food legume crop of tropical and subtropical regions. In this study, we describe the de novo assembly of a high-quality genome for the ancient pigeonpea landrace 'D30', achieved through a combination of Pacific Biosciences high-fidelity (PacBio HiFi) and high-throughput chromatin conformation capture (Hi-C) sequencing technologies. The assembled 'D30' genome has a size of 813.54Mb, with a contig N50 of 10.74Mb, a scaffold N50 of 73.07Mb, and a GC content of 35.67%. Genomic evaluation revealed that the 'D30' genome contains 99.2% of Benchmarking Universal Single-Copy Orthologs (BUSCO) and achieves a 29.06 long terminal repeat (LTR) assembly index (LAI). Genome annotation indicated that 'D30' encompasses 431.37Mb of repeat elements (53.02% of the genome) and 37 977 protein-coding genes. Identification of single-nucleotide polymorphisms (SNPs), insertions/deletions (indels), and structural variations between 'D30' and the published genome of pigeonpea cultivar 'Asha' suggests that genes affected by these variations may play important roles in biotic and abiotic stress responses. Further investigation of genomic regions under selection highlights genes enriched in starch and sucrose metabolism, with 42.11% of these genes highly expressed in seeds. Finally, we conducted genome-wide association studies (GWAS) to facilitate the identification of 28 marker-trait associations for six agronomic traits of pigeonpea. Notably, we discovered a calmodulin-like protein (CcCML) that harbors a dominant haplotype associated with the 100-seed weight of pigeonpea. Our study provides a foundational resource for developing genomics-assisted breeding programs in pigeonpea.

Local read haplotagging enables accurate long-read small variant calling

Long-read sequencing technology has enabled variant detection in difficult-to-map regions of the genome and enabled rapid genetic diagnosis in clinical settings. Rapidly evolving third-generation sequencing platforms like Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) are introducing newer platforms and data types. It has been demonstrated that variant calling methods based on deep neural networks can use local haplotyping information with long-reads to improve the genotyping accuracy. However, using local haplotype information creates an overhead as variant calling needs to be performed multiple times which ultimately makes it difficult to extend to new data types and platforms as they get introduced. In this work, we have developed a local haplotype approximate method that enables state-of-the-art variant calling performance with multiple sequencing platforms including PacBio Revio system, ONT R10.4 simplex and duplex data. This addition of local haplotype approximation simplifies long-read variant calling with DeepVariant.

Comparing Gold-Standard Sanger Sequencing with Two Next-Generation Sequencing Platforms of HIV-1 gp160 Single Genome Amplicons.

Our goal was to assess the accuracy of next generation sequencing (NGS) compared with Sanger. We performed single genome amplification (SGA) of HIV-1 gp160 on extracted tissue DNA from two HIV+ individuals. Amplicons (n = 30) were sequenced with Sanger or reamplified with barcoded primers and pooled before sequencing using Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PB). For each amplicon, a consensus sequence for NGS reads was obtained by (1) mapping reads to the Sanger sequence when available ("reference-based") or (2) mapping reads to a "pseudo-reference" sequence, i.e., a consensus sequence of a subset of NGS reads ("reference-free"). PB reads were clustered based on genetic similarity. A Sanger consensus sequence was obtained for 23/30 amplicons, for which all NGS consensus sequences were identical (n = 9) or nearly identical (n = 14) compared with Sanger. For the nine mismatches between Sanger/NGS, the nucleotide in the NGS sequence matched all other sequences from that patient. Of the 7/30 amplicons without a Sanger sequence, NGS sequences had ≥35 ambiguous calls in five amplicons and 0 ambiguitiesin two amplicons. Analysis of the electropherograms showed failure of a single sequencing primer for the latter two amplicons (consistent with a single template) and overlapping peaks for the other five (consistent with multiple templates). Clustering results closely followed the Sanger/NGS consensus results, where amplicons derived from a single template also had a single cluster and vice versa (with one exception, which could be the result of barcode misidentification). Representative sequences from the clusters contained 2-13 differences compared with Sanger/NGS. In summary, we show that both ONT and PB can produce amplicon consensus sequences with similar or higher accuracy compared with Sanger and, importantly, without the need for a known reference sequence. Clustering could be useful in some circumstances to predict or confirm the presence of multiple starting templates.

Telomere-to-telomere genome assembly of the goose Anser cygnoides

Our study presents the assembly of a high-quality Taihu goose genome at the Telomere-to-Telomere (T2T) level. By employing advanced sequencing technologies, including Pacific Biosciences HiFi reads, Oxford Nanopore long reads, Illumina short reads, and chromatin conformation capture (Hi-C), we achieved an exceptional assembly. The T2T assembly encompasses a total length of 1,197,991,206 bp, with contigs N50 reaching 33,928,929 bp and scaffold N50 attaining 81,007,908 bp. It consists of 73 scaffolds, including 38 autosomes and one pair of Z/W sex chromosomes. Importantly, 33 autosomes were assembled without any gap, resulting in a contiguous representation. Furthermore, gene annotation efforts identified 34,898 genes, including 436,162 RNA transcripts, encompassing 806,158 exons, 743,910 introns, 651,148 coding sequences (CDS), and 135,622 untranslated regions (UTR). The T2T-level chromosome-scale goose genome assembly provides a vital foundation for future genetic improvement and understanding the genetic mechanisms underlying important traits in geese.