Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

High accuracy of single-molecule real-time sequencing in detecting a rare α-globin fusion gene in carrier screening population.

INTRODUCTION: The α-globin fusion gene between the HBA2 and HBAP1 genes becomes clinically important in thalassemia screening because this fusion gene can cause severe hemoglobin (Hb) H disease when combining with α0 -thalassemia (α0 -thal). Due to its uncommon rearrangement in the α gene cluster without dosage changes, this fusion gene is undetectable by common molecular testing approaches used for α-thal diagnosis. METHODS: In this study, we used the single-molecule real-time (SMRT) sequencing technique to detect this fusion gene in 23 carriers identified by next-generation sequencing (NGS) among 16,504 screened individuals. Five primers for α and ß thalassemia were utilized. RESULTS: According to the NGS results, the 23 carriers include 14 pure heterozygotes, eight compound heterozygotes with common α-thal alleles, and one homozygote. By using SMRT, the fusion mutant was successfully detected in all 23 carriers. Furthermore, SMRT corrected the diagnosis in two "pure" heterozygotes: one was compound heterozygote with anti-3.7 triplication, and the other was homozygote. CONCLUSION: Our results indicate that SMRT is a superior method compared to NGS in detecting the α fusion gene, attributing to its efficient, accurate, and one-step properties.

Context-dependent DNA polymerization effects can masquerade as DNA modification signals.

BACKGROUND: Single molecule measurements of DNA polymerization kinetics provide a sensitive means to detect both secondary structures in DNA and deviations from primary chemical structure as a result of modified bases. In one approach to such analysis, deviations can be inferred by monitoring the behavior of DNA polymerase using single-molecule, real-time sequencing with zero-mode waveguide. This approach uses a Single Molecule Real Time (SMRT)-sequencing measurement of time between fluorescence pulse signals from consecutive nucleosides incorporated during DNA replication, called the interpulse duration (IPD). RESULTS: In this paper we present an analysis of loci with high IPDs in two genomes, a bacterial genome (E. coli) and a eukaryotic genome (C. elegans). To distinguish the potential effects of DNA modification on DNA polymerization speed, we paired an analysis of native genomic DNA with whole-genome amplified (WGA) material in which DNA modifications were effectively removed. Adenine modification sites for E. coli are known and we observed the expected IPD shifts at these sites in the native but not WGA samples. For C. elegans, such differences were not observed. Instead, we found a number of novel sequence contexts where IPDs were raised relative to the average IPDs for each of the four nucleotides, but for which the raised IPD was present in both native and WGA samples. CONCLUSION: The latter results argue strongly against DNA modification as the underlying driver for high IPD segments for C. elegans, and provide a framework for separating effects of DNA modification from context-dependent DNA polymerase kinetic patterns inherent in underlying DNA sequence for a complex eukaryotic genome.

Genome Resource of a Hypervirulent Strain C9-3 of Xanthomonas oryzae pv. oryzae Causing Bacterial Blight of Rice.

Xanthomonas oryzae pv. oryzae is the causal agent of bacterial blight, one of the most devastating diseases of rice. Here, a hypervirulent strain, C9-3, defeating Xa1, Xa10, xa13, and Xa23 resistance genes, was used to extract genomic DNA for single molecule real-time (SMRT) sequencing. After assembly, the genome consists of a single-circular chromosome with the size of 4,924,298 bp with G+C content of 63.7% and contains 4,715 genes. Annotation and analysis of the TALE genes using a suite of applications named AnnoTALE suggested that 17 transcription activator-like effectors, including 15 typical TALEs and 2 iTALEs/truncTALEs, were encoded in the genome. The approach and genome resource will contribute to the discovery of new virulence effectors and understanding on rice-X. oryzae pv. oryzae interactions.

Pathogenic 12-kb copy-neutral inversion in syndromic intellectual disability identified by high-fidelity long-read sequencing.

We report monozygotic twin girls with syndromic intellectual disability who underwent exome sequencing but with negative pathogenic variants. To search for variants that are unrecognized by exome sequencing, high-fidelity long-read genome sequencing (HiFi LR-GS) was applied. A 12-kb copy-neutral inversion was precisely identified by HiFi LR-GS after trio-based variant filtering. This inversion directly disrupted two genes, CPNE9 and BRPF1, the latter of which attracted our attention because pathogenic BRPF1 variants have been identified in autosomal dominant intellectual developmental disorder with dysmorphic facies and ptosis (IDDDFP), which later turned out to be clinically found in the twins. Trio-based HiFi LR-GS together with haplotype phasing revealed that the 12-kb inversion occurred de novo on the maternally transmitted chromosome. This study clearly indicates that submicroscopic copy-neutral inversions are important but often uncharacterized culprits in monogenic disorders and that long-read sequencing is highly advantageous for detecting such inversions involved in genetic diseases.

A Novel 5 kb Deletion in the ß-Globin Gene Cluster Identified in a Chinese Patient.

ß-Thalassemia (ß-thal), a highly prevalent disease in tropical and subtropical regions of Southern China, is caused mainly by point mutations in the ß-globin gene cluster. However, large deletions have also been found to contribute to some types of ß-thal. We identified a novel 5 kb deletion in the ß-globin cluster in a Chinese patient using multiplex ligation-dependent probe amplification (MLPA), and characterized it with single molecule real-time (SMRT) sequencing, gap-polymerase chain reaction (gap-PCR) and Sanger sequencing. The deletion was located between positions 5226189 and 5231091 on chromosome 11 (GRCh38), extending from 4 kb upstream of the 5' untranslated region (5'UTR) to the second intron of the ß-globin gene. The patient with this deletion presented with microcytosis and hypochromic red cells, as well as relatively high Hb F and Hb A2 levels. Our research indicated that SMRT sequencing is a useful tool for accurate detection of large deletions. Our study broadens the spectrum of deletional ß-thalassemias and provides a perspective for further study of the function of the ß-globin cluster.

De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth.

Casuarina equisetifolia (C. equisetifolia), a conifer-like angiosperm with resistance to typhoon and stress tolerance, is mainly cultivated in the coastal areas of Australasia. C. equisetifolia, making it a valuable model to study secondary growth associated genes and stress-tolerance traits. However, the genome sequence is unavailable and therefore wood-associated growth rate and stress resistance at the molecular level is largely unexplored. We therefore constructed a high-quality draft genome sequence of C. equisetifolia by a combination of Illumina second-generation sequencing reads and Pacific Biosciences single-molecule real-time (SMRT) long reads to advance the investigation of this species. Here, we report the genome assembly, which contains approximately 300 megabases (Mb) and scaffold size of N50 is 1.06 Mb. Additionally, gene annotation, assisted by a combination of prediction and RNA-seq data, generated 29 827 annotated protein-coding genes and 1983 non-coding genes, respectively. Furthermore, we found that the total number of repetitive sequences account for one-third of the genome assembly. Here we also construct the genome-wide map of DNA modification, such as two novel forms N6 -adenine (6mA) and N4-methylcytosine (4mC) at the level of single-nucleotide resolution using single-molecule real-time (SMRT) sequencing. Interestingly, we found that 17% of 6mA modification genes and 15% of 4mC modification genes also included alternative splicing events. Finally, we investigated cellulose, hemicellulose, and lignin-related genes, which were associated with secondary growth and contained different DNA modifications. The high-quality genome sequence and annotation of C. equisetifolia in this study provide a valuable resource to strengthen our understanding of the diverse traits of trees.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds.

BACKGROUND: The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map. RESULTS: Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features. CONCLUSIONS: The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Investigating the bacterial microbiota of traditional fermented dairy products using propidium monoazide with single-molecule real-time sequencing.

Traditional fermented dairy foods have been the major components of the Mongolian diet for millennia. In this study, we used propidium monoazide (PMA; binds to DNA of nonviable cells so that only viable cells are enumerated) and single-molecule real-time sequencing (SMRT) technology to investigate the total and viable bacterial compositions of 19 traditional fermented dairy foods, including koumiss from Inner Mongolia (KIM), koumiss from Mongolia (KM), and fermented cow milk from Mongolia (CM); sample groups treated with PMA were designated PKIM, PKM, and PCM. Full-length 16S rRNA sequencing identified 195 bacterial species in 121 genera and 13 phyla in PMA-treated and untreated samples. The PMA-treated and untreated samples differed significantly in their bacterial community composition and α-diversity values. The predominant species in KM, KIM, and CM were Lactobacillus helveticus, Streptococcus parauberis, and Lactobacillus delbrueckii, whereas the predominant species in PKM, PKIM, and PCM were Enterobacter xiangfangensis, Lactobacillus helveticus, and E. xiangfangensis, respectively. Weighted and unweighted principal coordinate analyses showed a clear clustering pattern with good separation and only minor overlapping. In addition, a pure culture method was performed to obtain lactic acid bacteria resources in dairy samples according to the results of SMRT sequencing. A total of 102 LAB strains were identified and Lb. helveticus (68.63%) was the most abundant, in agreement with SMRT sequencing results. Our results revealed that the bacterial communities of traditional dairy foods are complex and vary by type of fermented dairy product. The PMA treatment induced significant changes in bacterial community structure.

Comprehensive profiling of rhizome-associated alternative splicing and alternative polyadenylation in moso bamboo (Phyllostachys edulis).

Moso bamboo (Phyllostachys edulis) represents one of the fastest-spreading plants in the world, due in part to its well-developed rhizome system. However, the post-transcriptional mechanism for the development of the rhizome system in bamboo has not been comprehensively studied. We therefore used a combination of single-molecule long-read sequencing technology and polyadenylation site sequencing (PAS-seq) to re-annotate the bamboo genome, and identify genome-wide alternative splicing (AS) and alternative polyadenylation (APA) in the rhizome system. In total, 145 522 mapped full-length non-chimeric (FLNC) reads were analyzed, resulting in the correction of 2241 mis-annotated genes and the identification of 8091 previously unannotated loci. Notably, more than 42 280 distinct splicing isoforms were derived from 128 667 intron-containing full-length FLNC reads, including a large number of AS events associated with rhizome systems. In addition, we characterized 25 069 polyadenylation sites from 11 450 genes, 6311 of which have APA sites. Further analysis of intronic polyadenylation revealed that LTR/Gypsy and LTR/Copia were two major transposable elements within the intronic polyadenylation region. Furthermore, this study provided a quantitative atlas of poly(A) usage. Several hundred differential poly(A) sites in the rhizome-root system were identified. Taken together, these results suggest that post-transcriptional regulation may potentially have a vital role in the underground rhizome-root system.

Adaptive Strategies of Bacillus thuringiensis Isolated from Acid Mine Drainage Site in Sabah, Malaysia.

The adaptive process in bacteria is driven by specific genetic elements which regulate phenotypic characteristics such as tolerance to high metal ion concentrations and the secretion of protective biofilms. Extreme environments such as those associated with heavy metal pollution and extremes of acidity offer opportunities to study the adaptive mechanisms of microorganisms. This study focused on the genome analysis of Bacillus thuringiensis (Bt MCMY1), a gram positive rod shaped bacterium isolated from an acid mine drainage site in Sabah, Malaysia by using a combination of Single Molecule Real Time DNA Sequencing, Scanning Electron Microscopy (SEM) and Fourier Transform Infrared Spectroscopy (FTIR). The genome size of Bt MCMY1 was determined to be 5,458,152 bases which was encoded on a single chromosome. Analysis of the genome revealed genes associated with resistance to Copper, Mercury, Arsenic, Cobalt, Zinc, Cadmium and Aluminum. Evidence from SEM and FTIR indicated that the bacterial colonies form distinct films which bear the signature of polyhydroxyalkanoates (PHA) and this finding was supported by the genome data indicating the presence of a genetic pathway associated with the biosynthesis of PHAs. This is the first report of a Bacillus sp. isolated from an acid mine drainage site in Sabah, Malaysia and the genome sequence will provide insights into the manner in which B. thuringiensis adapts to acid mine drainage.

Long-Read Single Molecule Real-Time Full Gene Sequencing of Cytochrome P450-2D6.

The cytochrome P450-2D6 (CYP2D6) enzyme metabolizes ∼25% of common medications, yet homologous pseudogenes and copy number variants (CNVs) make interrogating the polymorphic CYP2D6 gene with short-read sequencing challenging. Therefore, we developed a novel long-read, full gene CYP2D6 single molecule real-time (SMRT) sequencing method using the Pacific Biosciences platform. Long-range PCR and CYP2D6 SMRT sequencing of 10 previously genotyped controls identified expected star (*) alleles, but also enabled suballele resolution, diplotype refinement, and discovery of novel alleles. Coupled with an optimized variant-calling pipeline, CYP2D6 SMRT sequencing was highly reproducible as triplicate intra- and inter-run nonreference genotype results were completely concordant. Importantly, targeted SMRT sequencing of upstream and downstream CYP2D6 gene copies characterized the duplicated allele in 15 control samples with CYP2D6 CNVs. The utility of CYP2D6 SMRT sequencing was further underscored by identifying the diplotypes of 14 samples with discordant or unclear CYP2D6 configurations from previous targeted genotyping, which again included suballele resolution, duplicated allele characterization, and discovery of a novel allele and tandem arrangement. Taken together, long-read CYP2D6 SMRT sequencing is an innovative, reproducible, and validated method for full-gene characterization, duplication allele-specific analysis, and novel allele discovery, which will likely improve CYP2D6 metabolizer phenotype prediction for both research and clinical testing applications.

High-resolution community profiling of arbuscular mycorrhizal fungi.

Community analyses of arbuscular mycorrhizal fungi (AMF) using ribosomal small subunit (SSU) or internal transcribed spacer (ITS) DNA sequences often suffer from low resolution or coverage. We developed a novel sequencing based approach for a highly resolving and specific profiling of AMF communities. We took advantage of previously established AMF-specific PCR primers that amplify a c. 1.5-kb long fragment covering parts of SSU, ITS and parts of the large ribosomal subunit (LSU), and we sequenced the resulting amplicons with single molecule real-time (SMRT) sequencing. The method was applicable to soil and root samples, detected all major AMF families and successfully discriminated closely related AMF species, which would not be discernible using SSU sequences. In inoculation tests we could trace the introduced AMF inoculum at the molecular level. One of the introduced strains almost replaced the local strain(s), revealing that AMF inoculation can have a profound impact on the native community. The methodology presented offers researchers a powerful new tool for AMF community analysis because it unifies improved specificity and enhanced resolution, whereas the drawback of medium sequencing throughput appears of lesser importance for low-diversity groups such as AMF.

High-accuracy de novo assembly and SNP detection of chloroplast genomes using a SMRT circular consensus sequencing strategy.

A circular consensus sequencing (CCS) strategy involving single molecule, real-time (SMRT) DNA sequencing technology was applied to de novo assembly and single nucleotide polymorphism (SNP) detection of chloroplast genomes. Chloroplast DNA was purified from enriched chloroplasts of pooled individuals to construct a shotgun library for each species. The sequencing reactions were performed on a PacBio RS platform. CCS sub-reads were generated from polymerase reads that passed the native dumbbell-shaped DNA templates multiple times. The complete chloroplast genome sequence was generated by mapping all reads to the draft sequence constructed in a step-by-step manner. The full-chain, PCR-free approach eliminates the possible context-specific biases in library construction and sequencing reaction. The chloroplast genome was easily and completely assembled using the data generated from one SMRT Cell without requiring a reference genome. Comparisons of the three assembled Fritillaria genomes to 34.1 kb of validation Sanger sequences revealed 100% concordance, and the detected intraspecies SNPs at a minimum variant frequency of 15% were all confirmed. This simple approach with potential for parallel sequencing yields high-quality chloroplast genomes for sensitive SNP detection and comparative analyses. We recommend this approach for its powerful applicability for evolutionary genetics and genomics studies in plants based on the sequences of chloroplast genomes.

Unveiling genomic features linked to traits of plant growth-promoting bacterial communities from sugarcane.

Microorganisms are ubiquitous, and those inhabiting plants have been the subject of several studies. Plant-associated bacteria exhibit various biological mechanisms that enable them to colonize host plants and, in some cases, enhance their fitness. In this study, we describe the genomic features predicted to be associated with plant growth-promoting traits in six bacterial communities isolated from sugarcane. The use of highly accurate single-molecule real-time sequencing technology for metagenomic samples from these bacterial communities allowed us to recover 17 genomes. The taxonomic assignments for the binned genomes were performed, revealing taxa distributed across three main phyla: Bacillota, Bacteroidota, and Pseudomonadota, with the latter being the most representative. Subsequently, we functionally annotated the metagenome-assembled genomes (MAGs) to characterize their metabolic pathways related to plant growth-promoting traits. Our study successfully identified the enrichment of important functions related to phosphate and potassium acquisition, modulation of phytohormones, and mechanisms for coping with abiotic stress. These findings could be linked to the robust colonization of these sugarcane endophytes.

Optimized SMRT-UMI protocol produces highly accurate sequence datasets from diverse populations-Application to HIV-1 quasispecies.

Pathogen diversity resulting in quasispecies can enable persistence and adaptation to host defenses and therapies. However, accurate quasispecies characterization can be impeded by errors introduced during sample handling and sequencing, which can require extensive optimizations to overcome. We present complete laboratory and bioinformatics workflows to overcome many of these hurdles. The Pacific Biosciences single molecule real-time platform was used to sequence polymerase-chain reaction (PCR) amplicons derived from cDNA templates tagged with unique molecular identifiers (SMRT-UMI). Optimized laboratory protocols were developed through extensive testing of different sample preparation conditions to minimize between-template recombination during PCR. The use of UMI allowed accurate template quantitation as well as removal of point mutations introduced during PCR and sequencing to produce a highly accurate consensus sequence from each template. Production of highly accurate sequences from the large datasets produced from SMRT-UMI sequencing is facilitated by a novel bioinformatic pipeline, Probabilistic Offspring Resolver for Primer IDs (PORPIDpipeline). PORPIDpipeline automatically filters and parses circular consensus reads by sample, identifies and discards reads with UMIs likely created from PCR and sequencing errors, generates consensus sequences, checks for contamination within the dataset, and removes any sequence with evidence of PCR recombination, heteroduplex formation, or early cycle PCR errors. The optimized SMRT-UMI sequencing and PORPIDpipeline methods presented here represent a highly adaptable and established starting point for accurate sequencing of diverse pathogens. These methods are illustrated through characterization of human immunodeficiency virus quasispecies in a virus transmitter-recipient pair of individuals.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features.

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Hb Q-Thailand heterozygosity unlinked with the (-α4.2/) α+-thalassemia deletion allele identified by long-read SMRT sequencing: hematological and molecular analyses.

OBJECTIVE: In the present study, two unrelated cases of Hb Q-Thailand heterozygosity unlinked with the (-α4.2/) α+-thalassemia deletion allele were identified by long-read single molecule real-time (SMRT) sequencing in southern China. The aim of this study was to report the hematological and molecular features as well as diagnostic aspects of the rare manifestation. METHODS: Hematological parameters and hemoglobin analysis results were recorded. A suspension array system for routine thalassemia genetic analysis and long-read SMRT sequencing were applied in parallel for thalassemia genotyping. Traditional methods, including Sanger sequencing, multiplex gap-polymerase chain reaction (gap-PCR) and multiplex ligation-dependent probe amplification (MLPA), were used together to confirm the thalassemia variants. RESULTS: Long-read SMRT sequencing was used to diagnose two Hb Q-Thailand heterozygous patients for whom the hemoglobin variant was unlinked to the (-α4.2/) allele for the first time. The hitherto undescribed genotypes were verified by traditional methods. Hematological parameters were compared with those of Hb Q-Thailand heterozygosity linked with the (-α4.2/) deletion allele in our study. For the positive control samples, long-read SMRT sequencing revealed a linkage relationship between the Hb Q-Thailand allele and the (-α4.2/) deletion allele. CONCLUSIONS: Identification of the two patients confirms that the linkage relationship between the Hb Q-Thailand allele and the (-α4.2/) deletion allele is a common possibility but not a certainty. Remarkably, as it is superior to traditional methods, SMRT technology may eventually serve as a more comprehensive and precise method that holds promising prospects in clinical practice, especially for rare variants.

The Long and Short of Genome Sequencing: Using a Hybrid Sequencing Strategy to Sequence Oral Microbial Genomes.

Since our chapter on genome sequencing using the GS-FLX pyrosequencer in the First Edition of this book, significant advances have been made in next-generation DNA sequencing (NGS) technology. Not only has the GS-FLX become extinct, but the more recent introduction and establishment of the so-called third-generation DNA sequencers by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) has revolutionized genomics yet again by generating ultra-long (>100,000 basepair) sequence reads concomitant with an incredible reduction in cost per sequenced basepair. Unfortunately, the ultra-high sequence yields of third-generation sequencers are compromised by their inherent sequencing error rates, prompting an alternative sequencing strategy, i.e., a hybrid sequencing strategy, which combines PacBio/ONT primary datasets with complementary datasets generated by mainstream short-read NGS platforms, e.g., Illumina or Ion Torrent. Although the concept of a hybrid sequencing strategy is not new, existing yields and accuracy of ultra-long and short-read sequencing technologies makes such a strategy achievable, resulting in complete genome sequences in one hit. In this chapter, we describe our updated laboratory and bioinformatic protocols that will allow the average research group to obtain complete oral microbial genome sequences assembled from a combination of DNA sequence data generated by NGS and third-generation platforms.

A novel 15.8 kb deletion α-thalassemia confirmed by long-read single-molecule real-time sequencing: Hematological phenotypes and molecular characterization.

BACKGROUND: Thalassemia is the most frequent recessive Mendelian inherited monogenic disease worldwide, and is characterized by the impaired synthesis of globin chains due to disease-causing variants in α- or ß-globin genes. There are many conventional methods to diagnose thalassemia but all of them have limitations. CASE REPORT: We present the case of a 37-year-old female with abnormal values of routine hematological indices who was admitted for genetic screening of thalassemia. Genomic DNA was extracted and used for genetic assays covering the known and potential novel genotypes in HBA and HBB genes using a suspension-array system, gap-polymerase chain reaction (Gap-PCR), PCR-reverse dot blot (PCR-RDB) and multiplex ligation-dependent probe amplification (MLPA). Finally, using long-read single-molecule real-time (SMRT) sequencing, we first confirmed the case with a novel 15.8 kb deletion located in the HBA gene (Chr16:163886-179768, GRch38/hg38). CONCLUSIONS: Our results showed that long-read SMRT sequencing has great advantages in the detection of rare α-globin gene variants. This study may provide a reference protocol for the use of long-read SMRT sequencing for the detection of known and potential novel genotypes of thalassemia in the population and improve the accuracy of genetic counseling and prenatal diagnosis.