Long-read Sequencing Techniques Research Articles

Haplotype determination of the Bombyx mori nucleopolyhedrovirus by Nanopore sequencing and linkage of single nucleotide variants.

Naturally occurring isolates of baculoviruses, such as the Bombyx mori nucleopolyhedrovirus (BmNPV), usually consist of numerous genetically different haplotypes. Deciphering the different haplotypes of such isolates is hampered by the large size of the dsDNA genome, as well as the short read length of next generation sequencing (NGS) techniques that are widely applied for baculovirus isolate characterization. In this study, we addressed this challenge by combining the accuracy of NGS to determine single nucleotide variants (SNVs) as genetic markers with the long read length of Nanopore sequencing technique. This hybrid approach allowed the comprehensive analysis of genetically homogeneous and heterogeneous isolates of BmNPV. Specifically, this allowed the identification of two putative major haplotypes in the heterogeneous isolate BmNPV-Ja by SNV position linkage. SNV positions, which were determined based on NGS data, were linked by the long Nanopore reads in a Position Weight Matrix. Using a modified Expectation-Maximization algorithm, the Nanopore reads were assigned according to the occurrence of variable SNV positions by machine learning. The cohorts of reads were de novo assembled, which led to the identification of BmNPV haplotypes. The method demonstrated the strength of the combined approach of short- and long-read sequencing techniques to decipher the genetic diversity of baculovirus isolates.

Complete Genomic Assembly of Mauritian Cynomolgus Macaque Killer Ig-like Receptor and Natural Killer Group 2 Haplotypes.

Mauritian-origin cynomolgus macaques (MCMs) serve as a powerful nonhuman primate model in biomedical research due to their unique genetic homogeneity, which simplifies experimental designs. Despite their extensive use, a comprehensive understanding of crucial immune-regulating gene families, particularly killer Ig-like receptors (KIR) and NK group 2 (NKG2), has been hindered by the lack of detailed genomic reference assemblies. In this study, we employ advanced long-read sequencing techniques to completely assemble eight KIR and seven NKG2 genomic haplotypes, providing an extensive insight into the structural and allelic diversity of these immunoregulatory gene clusters. Leveraging these genomic resources, we prototype a strategy for genotyping KIR and NKG2 using short-read, whole-exome capture data, illustrating the potential for cost-effective multilocus genotyping at colony scale. These results mark a significant enhancement for biomedical research in MCMs and underscore the feasibility of broad-scale genetic investigations.

Abstract P160: Prevalence of Genetic Dilated Cardiomyopathy: A Systematic Literature Review

Background: Dilated cardiomyopathy (DCM) is a leading cause of heart failure and cardiac transplantation worldwide. Many cases of DCM are due to genetic or idiopathic etiology. Two prominent contributors to genetic DCM (gDCM) are the TTN and MYH7 genes, which encode proteins associated with DCM etiology. Purpose: This study aimed to synthesize the published literature on the global prevalence of adult and pediatric gDCM within the clinical DCM population and the general population. Methods: A systematic review was conducted by searching MEDLINE® and Embase from database inception to 09/19/2022. We included observational studies reporting the prevalence of gDCM within any population; prevalence of TTN and MYH7 mutations were also extracted. Excluded were studies using family history as a proxy for gDCM. Median with interquartile range (IQR) or mean with standard deviation (SD) weighted by sample size (in cases of calculations with <10 studies) were used to summarize the data. Results: Of 2,720 identified abstracts, 57 studies were included in the review. Among the adult or mixed (adults and children) DCM population in any country, median prevalence was 21.0% for overall DCM-related genetic mutations, 11.4% for TTN mutations, and 3.2% for MYH7 mutations ( Table ). The prevalence of pediatric gDCM within the DCM population was similar (weighted mean: 21.3%). Data on prevalence of gDCM within the general population were scarce. Conclusion: This study highlights the variable prevalence of gDCM reported across different populations and settings. The current evidence may underestimate the genetic causes due to technical limitations in genetic testing and limited screening. Accurately determining the prevalence of gDCM is critical for identifying high-risk populations, developing targeted prevention strategies for disease complications, and improving outcomes; epidemiological studies using long-read sequencing techniques and large cohort datasets with genotype-phenotype correlation analyses are needed.

A chromosome-level genome of mango exclusively from long-read sequence data.

Improvements in long-read sequencing techniques have greatly accelerated plant genome sequencing. Current de novo assemblies are routinely achieved by assembling long-read sequence data into contigs that are assembled to chromosome level by chromatin conformation capture. We report here a chromosome-level mango genome using only PacBio high-fidelity (HiFi) long reads. HiFi reads at high coverage (204x) resulted in the assembly of 17 chromosomes, each as a single contig with telomeres at both ends. The remaining three chromosomes were represented each by two contigs, with telomeres at one end and ribosomal repeats at the other end. Analyzing contig ends allowed them to be paired and linked to generate the remaining three complete chromosomes, telomere-to-telomere but with ribosomal repeats of uncertain length. The assembled genome was 365 Mb with 100% completeness as assessed by Benchmarking Universal Single-Copy Orthologs analysis. The haplotypes assembled demonstrated extensive structural differences. This approach using very high genome coverage may be useful for assembling high-quality genomes for many other plants.

Characterisation of a LINE-1 Insertion in the RP1 Gene by Targeted Adaptive Nanopore Sequencing in a Family with Retinitis Pigmentosa

Retinitis pigmentosa (RP) is a group of inherited degenerative retinal disorders affecting more than 1.5 million people worldwide. For 30-50% of individuals with RP, the genetic cause remains unresolved by current clinical diagnostic gene panels. It is likely explained by variants in novel RP-associated genes or noncoding regulatory regions, or by complex genetic alterations such as large structural variants. Recent developments in long-read sequencing techniques have opened an opportunity for efficient analysis of complex genetic variants. We analysed a Finnish family with dominantly inherited RP affecting six individuals in three generations. Two affected individuals underwent a comprehensive clinical examination in combination with a clinical diagnostic gene panel, followed by whole exome sequencing in our laboratory. They exhibited typical signs of RP, yet initial sequence analysis found no causative variants. Reanalysis of the sequencing data detected a LINE-1 (L1) retrotransposon insertion of unknown size in exon 4 of the RP1 axonemal microtubule-associated (RP1) gene. The large chimeric L1 insertion that segregated with the disease was further characterised using targeted adaptive nanopore sequencing of RP1, allowing us to identify a 5.6 kb L1 transposable element insertion in RP1 as the cause of RP in this family with dominantly inherited RP.

Improved sensitivity for detection of pathogenic variants in familial NF2-related schwannomatosis

PurposeTo determine the impact of additional genetic screening techniques on the rate of detection of pathogenic variants leading to familial NF2-related schwannomatosis.MethodsWe conducted genetic screening of a cohort of 168...

Gross Chromosomal Rearrangement at Centromeres.

Centromeres play essential roles in the faithful segregation of chromosomes. CENP-A, the centromere-specific histone H3 variant, and heterochromatin characterized by di- or tri-methylation of histone H3 9th lysine (H3K9) are the hallmarks of centromere chromatin. Contrary to the epigenetic marks, DNA sequences underlying the centromere region of chromosomes are not well conserved through evolution. However, centromeres consist of repetitive sequences in many eukaryotes, including animals, plants, and a subset of fungi, including fission yeast. Advances in long-read sequencing techniques have uncovered the complete sequence of human centromeres containing more than thousands of alpha satellite repeats and other types of repetitive sequences. Not only tandem but also inverted repeats are present at a centromere. DNA recombination between centromere repeats can result in gross chromosomal rearrangement (GCR), such as translocation and isochromosome formation. CENP-A chromatin and heterochromatin suppress the centromeric GCR. The key player of homologous recombination, Rad51, safeguards centromere integrity through conservative noncrossover recombination between centromere repeats. In contrast to Rad51-dependent recombination, Rad52-mediated single-strand annealing (SSA) and microhomology-mediated end-joining (MMEJ) lead to centromeric GCR. This review summarizes recent findings on the role of centromere and recombination proteins in maintaining centromere integrity and discusses how GCR occurs at centromeres.

The highly improved genome of Ixodes scapularis with X and Y pseudochromosomes.

Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the ∼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.

FLIBase: a comprehensive repository of full-length isoforms across human cancers and tissues.

Regulatory processes at the RNA transcript level play a crucial role in generating transcriptome diversity and proteome composition in human cells, impacting both physiological and pathological states. This study introduces FLIBase (www.FLIBase.org), a specialized database that focuses on annotating full-length isoforms using long-read sequencing techniques. We collected and integrated long-read (351 samples) and short-read (12469 samples) RNA sequencing data from diverse normal and cancerous human tissues and cells. The current version of FLIBase comprises a total of 983789 full-length spliced isoforms, identified through long-read sequences and verified using short-read exon-exon splice junctions. Of these, 188248 isoforms have been annotated, while 795541 isoforms remain unannotated. By overcoming the limitations of short-read RNA sequencing methods, FLIBase provides an accurate and comprehensive representation of full-length transcripts. These comprehensive annotations empower researchers to undertake various downstream analyses and investigations. Importantly, FLIBase exhibits a significant advantage in identifying a substantial number of previously unannotated isoforms and tumor-specific RNA transcripts. These tumor-specific RNA transcripts have the potential to serve as a source of immunogenic recurrent neoantigens. This remarkable discovery holds tremendous promise for advancing the development of tailored RNA-based diagnostic and therapeutic strategies for various types of human cancer.

Population genomics-guided engineering of phenazine biosynthesis in Pseudomonas chlororaphis

The emergence of next-generation sequencing (NGS) technologies has made it possible to not only sequence entire genomes, but also identify metabolic engineering targets across the pangenome of a microbial population. This study leverages NGS data as well as existing molecular biology and bioinformatics tools to identify and validate genomic signatures for improving phenazine biosynthesis in Pseudomonas chlororaphis. We sequenced a diverse collection of 34 Pseudomonas isolates using short- and long-read sequencing techniques and assembled whole genomes using the NGS reads. In addition, we assayed three industrially relevant phenotypes (phenazine production, biofilm formation, and growth temperature) for these isolates in two different media conditions. We then provided the whole genomes and phenazine production data to a unitig-based microbial genome-wide association study (mGWAS) tool to identify novel genomic signatures responsible for phenazine production in P. chlororaphis. Post-processing of the mGWAS analysis results yielded 330 significant hits influencing the biosynthesis of one or more phenazine compounds. Based on a quantitative metric (called the phenotype score), we elucidated the most influential hits for phenazine production and experimentally validated them in vivo in the most optimal phenazine producing strain. Two genes significantly increased phenazine-1-carboxamide (PCN) production: a histidine transporter (ProY_1), and a putative carboxypeptidase (PS__04251). A putative MarR-family transcriptional regulator decreased PCN titer when overexpressed in a high PCN producing isolate. Overall, this work seeks to demonstrate the utility of a population genomics approach as an effective strategy in enabling the identification of targets for metabolic engineering of bioproduction hosts.

De Novo Structural Variations of Escherichia coli Detected by Nanopore Long-Read Sequencing.

Spontaneous mutations power evolution, whereas large-scale structural variations (SVs) remain poorly studied, primarily because of the lack of long-read sequencing techniques and powerful analytical tools. Here, we explore the SVs of Escherichia coli by running 67 wild-type (WT) and 37 mismatch repair (MMR)-deficient (ΔmutS) mutation accumulation lines, each experiencing more than 4,000 cell divisions, by applying Nanopore long-read sequencing and Illumina PE150 sequencing and verifying the results by Sanger sequencing. In addition to precisely repeating previous mutation rates of base-pair substitutions and insertion and deletion (indel) mutation rates, we do find significant improvement in insertion and deletion detection using long-read sequencing. The long-read sequencing and corresponding software can particularly detect bacterial SVs in both simulated and real data sets with high accuracy. These lead to SV rates of 2.77 × 10-4 (WT) and 5.26 × 10-4 (MMR-deficient) per cell division per genome, which is comparable with previous reports. This study provides the SV rates of E. coli by applying long-read sequencing and SV detection programs, revealing a broader and more accurate picture of spontaneous mutations in bacteria.

Decomposing mosaic tandem repeats accurately from long reads.

Over the past 30 years, extended tandem repeats (TRs) have been correlated with approximately 60 diseases with high odds ratios, and most known TRs consist of single repeat units. However, in the last few years, mosaic TRs composed of different units have been found to be associated with several brain disorders by long-read sequencing techniques. Mosaic TRs are difficult-to-characterize sequence configurations that are usually confirmed by manual inspection. Widely used tools are not designed to solve the mosaic TR problem and often fail to properly decompose mosaic TRs. We propose an efficient algorithm that can decompose mosaic tandem repeats in the input string with high sensitivity. Using synthetic benchmark data, we demonstrate that our program named uTR outperforms TRF and RepeatMasker in terms of prediction accuracy, this is especially true when mosaic TRs are more complex, and uTR is faster than TRF and RepeatMasker in most cases. The software program uTR that implements the proposed algorithm is available at https://github.com/morisUtokyo/uTR.

Symbiont-screener: A reference-free tool to separate host sequences from symbionts for error-prone long reads

Metagenomic sequencing facilitates large-scale constitutional analysis and functional characterization of complex microbial communities without cultivation. Recent advances in long-read sequencing techniques utilize long-range information to simplify repeat-aware metagenomic assembly puzzles and complex genome binning tasks. However, it remains methodologically challenging to remove host-derived DNA sequences from the microbial community at the read resolution due to high sequencing error rates and the absence of reference genomes. We here present Symbiont-Screener (https://github.com/BGI-Qingdao/Symbiont-Screener), a reference-free approach to identifying high-confidence host’s long reads from symbionts and contaminants and overcoming the low sequencing accuracy according to a trio-based screening model. The remaining host’s sequences are then automatically grouped by unsupervised clustering. When applied to both simulated and real long-read datasets, it maintains higher precision and recall rates of identifying the host’s raw reads compared to other tools and hence promises the high-quality reconstruction of the host genome and associated metagenomes. Furthermore, we leveraged both PacBio HiFi and nanopore long reads to separate the host’s sequences on a real host-microbe system, an algal-bacterial sample, and retrieved an obvious improvement of host assembly in terms of assembly contiguity, completeness, and purity. More importantly, the residual symbiotic microbiomes illustrate improved genomic profiling and assemblies after the screening, which elucidates a solid basis of data for downstream bioinformatic analyses, thus providing a novel perspective on symbiotic research.

De novel heterozygous copy number deletion on 7q31.31-7q31.32 involving TSPAN12 gene with familial exudative vitreoretinopathy in a Chinese family.

To investigate the genetic and clinical characteristics of patients with a large heterozygous copy number deletion on 7q31.31-7q31.32. A family with familial exudative vitreoretinopathy (FEVR) phenotype was included in the study. Whole-exome sequencing (WES) was initially used to locate copy number variations (CNVs) on 7q31.31-31.32, but failed to detect the precise breakpoint. The long-read sequencing, Oxford Nanopore sequencing Technology (ONT) was used to get the accurate breakpoint which is verified by quantitative real-time polymerase chain reaction (QPCR) and Sanger Sequencing. The proband, along with her father and younger brother, were found to have a heterozygous 4.5 Mb CNV deletion located on 7q31.31-31.32, which included the FEVR-related gene TSPAN12. The specific deletion was confirmed as del(7)(q31.31q31.32)chr7:g.119451239_123956818del. The proband exhibited a phase 2A FEVR phenotype, characterized by a falciform retinal fold, macular dragging, and peripheral neovascularization with leaking of fluorescence. These symptoms led to a significant decrease in visual acuity in both eyes. On the other hand, the affected father and younger brother showed a milder phenotype. The heterozygous CNV deletion located on 7q31.31-7q31.32 is associated with the FEVR phenotype. The use of long-read sequencing techniques is essential for accurate molecular diagnosis of genetic disorders.

Genome characterization and CRISPR-Cas9 editing of a human neocentromere

The maintenance of genome integrity is ensured by proper chromosome inheritance during mitotic and meiotic cell divisions. The chromosomal counterpart responsible for chromosome segregation to daughter cells is the centromere, at which the spindle apparatus attaches through the kinetochore. Although all mammalian centromeres are primarily composed of megabase-long repetitive sequences, satellite-free human neocentromeres have been described. Neocentromeres and evolutionary new centromeres have revolutionized traditional knowledge about centromeres. Over the past 20 years, insights have been gained into their organization, but in spite of these advancements, the mechanisms underlying their formation and evolution are still unclear. Today, through modern and increasingly accessible genome editing and long-read sequencing techniques, research in this area is undergoing a sudden acceleration. In this article, we describe the primary sequence of a previously described human chromosome 3 neocentromere and observe its possible evolution and repair results after a chromosome breakage induced through CRISPR-Cas9 technologies. Our data represent an exciting advancement in the field of centromere/neocentromere evolution and chromosome stability.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing.

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

Evaluation of CircRNA Sequence Assembly Methods Using Long Reads.

The functional study on circRNAs has been increasing in the past decade due to its important roles in micro RNA sponge, protein coding, the initiation, and progression of diseases. The study of circRNA functions depends on the full-length sequences of circRNA, and current sequence assembly methods based on short reads face challenges due to the existence of linear transcript. Long reads produced by long-read sequencing techniques such as Nanopore technology can cover full-length sequences of circRNA and therefore can be used to evaluate the correctness and completeness of circRNA full sequences assembled from short reads of the same sample. Using long reads of the same samples, one from human and the other from mouse, we have comprehensively evaluated the performance of several well-known circRNA sequence assembly algorithms based on short reads, including circseq_cup, CIRI_full, and CircAST. Based on the F1 score, the performance of CIRI-full was better in human datasets, whereas in mouse datasets CircAST was better. In general, each algorithm was developed to handle special situations or circumstances. Our results indicated that no single assembly algorithm generated better performance in all cases. Therefore, these assembly algorithms should be used together for reliable full-length circRNA sequence reconstruction. After analyzing the results, we have introduced a screening protocol that selects out exonic circRNAs with full-length sequences consisting of all exons between back splice sites as the final result. After screening, CIRI-full showed better performance for both human and mouse datasets. The average F1 score of CIRI-full over four circRNA identification algorithms increased from 0.4788 to 0.5069 in human datasets, and it increased from 0.2995 to 0.4223 in mouse datasets.

Using structural variants to understand the ecological and evolutionary dynamics of fungal plant pathogens.

Deletions, duplications, insertions, inversions and translocations are commonly referred to as structural variants (SVs). Fungal plant pathogens have compact genomes, facilitating the generation of accurate maps of SVs for these species in recent studies. Structural variants have been found to constitute a significant proportion of the standing genetic variation in fungal plant pathogen populations, potentially leading to the generation of accessory genes, regions or chromosomes enriched in pathogenicity factors. Structural variants are involved in the rapid adaptation and ecological traits of pathogens, including host specialization and mating. Long-read sequencing techniques coupled with theoretical and experimental approaches have considerable potential for elucidating the phenotypic effects of SVs and deciphering the evolutionary and genomic mechanisms underlying the formation of SVs in fungal plant pathogens.

Comparative analysis of alignment tools for application on Nanopore sequencing data

Abstract INTRODUCTION: Long-read sequencing techniques such as Oxford Nanopore sequencing, are representing a promising novel approach in molecular-biological methodology, enabling potential facilitation in mapping and de novo assembly. In comparison to conventional sequencing methods, novel alignment tools are mandated to compensate differing data structures (especially high error rate) to achieve acceptably accurate analysis results. METHODS: In this study, benchmarking for long read aligners BLASR, GraphMap, LAST, minimap2, NGMLR and the short-read aligner BWA MEM on three experimental datasets was conducted. Obtained alignment results were compared for various quality and performance criteria, such as match rate, mismatch rate, error rate, working memory usage and computational time. RESULTS: The comparison yielded differences in alignment quality and performance of tools under test. Tool LAST showed the largest differences among all tools. Minimap2 achieved constant quality with good performance. BLASR, GraphMap, BWA MEM and NGMLR showed slight differences only. CONCLUSION: Differences among the tools could be reasoned with dataset characteristics and algorithm approaches of individual tools. All tools except BLASR seem applicable for Nanopore sequencing data. Therefore, selection of the tool should be done under consideration of the experimental design and the further downstream analysis

Meta-analytic approach for transcriptome profiling of herpes simplex virus type 1

In this meta-analysis, we re-analysed and compared herpes simplex virus type 1 transcriptomic data generated by eight studies using various short- and long-read sequencing techniques and different library preparation methods. We identified a large number of novel mRNAs, non-coding RNAs and transcript isoforms, and validated many previously published transcripts. Here, we present the most complete HSV-1 transcriptome to date. Furthermore, we also demonstrate that various sequencing techniques, including both cDNA and direct RNA sequencing approaches, are error-prone, which can be circumvented by using integrated approaches. This work draws attention to the need for using multiple sequencing approaches and meta-analyses in transcriptome profiling studies to obtain reliable results.