Single-molecule Real-time Long-read Sequencing Research Articles

HPV genotyping in clinical samples using long-read single-molecule real-time sequencing.

Human papillomavirus (HPV) genotyping is widely used, particularly in combination with high-risk (HR) HPV tests for cervical cancer screening. We developed a genotyping method using sequences of approximately 800 bp in the E6/E7 region obtained by PacBio single molecule real-time sequencing (SMRT) and evaluated its performance against MY09-11 L1 sequencing and after the APTIMA HPV genotyping assay. The levels of concordance of PacBio E6/E7 SMRT sequencing with MY09-11 L1 sequencing and APTIMA HPV genotyping were 100% and 90.8%, respectively. The sensitivity of PacBio E6/EA7 SMRT was slightly greater than that of L1 sequencing and, as expected, lower than that of HR-HPV tests. In the context of cervical cancer screening, PacBio E6/E7 SMRT is then best used after a positive HPV test. PacBio E6/E7 SMRT genotyping is an attractive alternative for HR and LR-HPV genotyping of clinical samples. PacBio SMRT sequencing provides unbiased genotyping and can detect multiple HPV infections and haplotypes within a genotype.

Identification of Hb Lepore, Hb anti-Lepore, and α-globin gene triplications by long-read single-molecule real-time sequencing.

Hemoglobin (Hb) Lepore and Hb anti-Lepore are infrequent fusion gene variants that result from nonhomologous crossovers during meiosis. Conventional molecular testing methods may face challenges in identifying these variants. During Hb analysis using capillary electrophoresis, we encountered 6 cases with unusual Hb variants. Our aim was to identify the alterations in their globin genes. Gap-polymerase chain reaction (PCR), reverse dot-blot assay (RDB), Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA), and long-read single-molecule real-time (SMRT) sequencing were used to confirm the presence of globin gene alterations. The routine thalassemia gene test kit using the gap-PCR and RDB techniques did not detect common gene variations. Direct sequencing failed to identify any known or unknown globin gene alterations. The MLPA analysis, however, revealed the possible presence of α-globin gene triplications as well as 2 types of fusion gene alterations. Further analysis using long-read SMRT sequencing accurately identified 3 rare gene variations: αααanti-3.7, Hb Lepore-Boston-Washington, and Hb anti-Lepore P-India. Conventional methods may overlook rare thalassemias or Hb variants. Long-read SMRT sequencing has the potential to identify breakpoints in fusion genes, demonstrating that it is a promising technique for detecting rare thalassemias.

SARS-CoV-2 Co-Infections and Recombinations Identified by Long-Read Single-Molecule Real-Time Sequencing.

Co-infection with at least 2 strains of virus is the prerequisite for recombination, one of the means of genetic diversification. Little is known about the prevalence of these events in SARS-CoV-2, partly because it is difficult to detect them. We used long-read PacBio single-molecule real-time (SMRT) sequencing technology to sequence whole genomes and targeted regions for haplotyping. We identified 17 co-infections with SARS-CoV-2 strains belonging to different clades in 6829 samples sequenced between January and October, 2022 (prevalence 0.25%). There were 3 Delta/Omicron co-infections and 14 Omicron/Omicron co-infections (4 cases of 21K/21L, 1 case of 21L/22A, 2 cases of 21L/22B, 4 cases of 22A/22B, 2 cases of 22B/22C and 1 case of 22B/22E). Four of these patients (24%) also harbored recombinant minor haplotypes, including one with a recombinant virus that was selected in the viral quasispecies over the course of his chronic infection. While co-infections remain rare among SARS-CoV-2-infected individuals, long-read SMRT sequencing is a useful tool for detecting them as well as recombinant events, providing the basis for assessing their clinical impact, and a precise indicator of epidemic evolution. IMPORTANCE SARS-CoV-2 variants have been responsible for the successive waves of infection over the 3 years of pandemic. While co-infection followed by recombination is one driver of virus evolution, there have been few reports of co-infections, mainly between Delta and Omicron variants or between the first 2 Omicron variants 21K_BA.1 and 21L_BA.2. The 17 co-infections we detected during 2022 included cases with the recent clades of Omicron 22A, 22B, 22C, and 22E; 24% harbored recombinant variants. This study shows that long-read SMRT sequencing is well suited to SARS-CoV-2 genomic surveillance.

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

ABSTRACTObjectiveIn 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.MethodsHematological 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.ResultsLong-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.ConclusionsIdentification 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.

Detection of rare thalassemia mutations using long-read single-molecule real-time sequencing.

Gap- polymerase chain reaction (PCR), reverse dot-blot assay (RDB), real-time PCR based multicolor melting curve analysis (MMCA assay), multiplex ligation-dependent probe amplification (MLPA) and Sanger sequencing are conventional methods to diagnose thalassemia but all of them have limitations. In this study, we applied single-molecule real-time (SMRT) sequencing following multiplex long-range PCR to uncover rare mutations in nine patients and their family members. The patients with different results between Gap-PCR and MMCA assay or with phenotype not matching genotype were included. Using SMRT sequencing, we first identified the carriers with αααanti3.7/HKαα, -α762bpα/αα (chr16:172,648-173,409), ααfusion/αQSα (in a trans configuration), two cases with novel gene rearrangements and another case with a novel 341bp insertion in α-globin gene cluster, respectively. One carrier with --SEA/αααanti4.2, and two carriers with the coexistence of globin variant and an α-globin gene duplication were also found. Most importantly, we could determine two defects in α-globin gene cluster being a cis or trans configuration in a single test. Our results showed that SMRT has great advantages in detection of α-globin gene triplications, rare deletions and determination of a cis or trans configuration. SMRT is a comprehensive and one-step method for thalassemia screening and diagnosis, especially for detection of rare thalassemia mutations.

Mutational spectrum of hepatitis C virus in patients with chronic hepatitis C determined by single molecule real-time sequencing

The emergence of hepatitis C virus (HCV) with resistance-associated substitution (RAS), produced by mutations in the HCV genome, is a major problem in direct acting antivirals (DAA) treatment. This study aimed to clarify the mutational spectrum in HCV-RNA and the substitution pattern for the emergence of RASs in patients with chronic HCV infection. HCV-RNA from two HCV replicon cell lines and the serum HCV-RNA of four non-liver transplant and four post-liver transplant patients with unsuccessful DAA treatment were analyzed using high-accuracy single-molecule real-time long-read sequencing. Transition substitutions, especially A>G and U>C, occurred prominently under DAAs in both non-transplant and post-transplant patients, with a mutational bias identical to that occurring in HCV replicon cell lines during 10-year culturing. These mutational biases were reproduced in natural courses after DAA treatment. RASs emerged via both transition and transversion substitutions. NS3-D168 and NS5A-L31 RASs resulted from transversion mutations, while NS5A-Y93 RASs was caused by transition substitutions. The fidelity of the RNA-dependent RNA polymerase, HCV-NS5B, produces mutational bias in the HCV genome, characterized by dominant transition mutations, notably A>G and U>C substitutions. However, RASs are acquired by both transition and transversion substitutions, and the RASs-positive HCV clones are selected and proliferated under DAA treatment pressure.

CFH and CFHR Copy Number Variations in C3 Glomerulopathy and Immune Complex-Mediated Membranoproliferative Glomerulonephritis.

C3 Glomerulopathy (C3G) and Immune Complex-Mediated Membranoproliferative glomerulonephritis (IC-MPGN) are rare diseases characterized by glomerular deposition of C3 caused by dysregulation of the alternative pathway (AP) of complement. In approximately 20% of affected patients, dysregulation is driven by pathogenic variants in the two components of the AP C3 convertase, complement C3 (C3) and Factor B (CFB), or in complement Factor H (CFH) and Factor I (CFI), two genes that encode complement regulators. Copy number variations (CNVs) involving the CFH-related genes (CFHRs) that give rise to hybrid FHR proteins also have been described in a few C3G patients but not in IC-MPGN patients. In this study, we used multiplex ligation-dependent probe amplification (MLPA) to study the genomic architecture of the CFH-CFHR region and characterize CNVs in a large cohort of patients with C3G (n = 103) and IC-MPGN (n = 96) compared to healthy controls (n = 100). We identified new/rare CNVs resulting in structural variants (SVs) in 5 C3G and 2 IC-MPGN patients. Using long-read single molecule real-time sequencing (SMRT), we detected the breakpoints of three SVs. The identified SVs included: 1) a deletion of the entire CFH in one patient with IC-MPGN; 2) an increased number of CFHR4 copies in one IC-MPGN and three C3G patients; 3) a deletion from CFHR3-intron 3 to CFHR3-3′UTR (CFHR34–6Δ) that results in a FHR3-FHR1 hybrid protein in a C3G patient; and 4) a CFHR31–5-CFHR410 hybrid gene in a C3G patient. This work highlights the contribution of CFH-CFHR CNVs to the pathogenesis of both C3G and IC-MPGN.

Long-read sequencing identified a novel nonsense and a de novo missense of PPA2 in trans in a Chinese patient with autosomal recessive infantile sudden cardiac failure

Background and aimsBiallelic missense variants in PPA2 gene cause infantile sudden cardiac failure (SCFI; OMIM #617222) characterized by sudden cardiac failure, sudden cardiac death in infants. Here, we present an unusual survivor with one inherited plus one de novo variant in PPA2. Since next-generation sequencing (NGS) fails to resolve variant phasing, which require long-read sequencing to clarify the diagnosis. Materials and methodsWhole exome and Sanger sequencing were initially performed to identify the causative variants. PCR-based short tandem repeats (STRs) analysis and long-read single molecule real-time (SMRT) sequencing were further implemented for paternity testing and variant phasing. Pathogenicity evaluation of the biallelic variants in PPA2 was conducted according to the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) guidelines based on VarSome. ResultsWhole exome and Sanger sequencing revealed two variants in PPA2, with one novel nonsense variant (c.524C > G; p.Ser175*) inherited from the mother and one de novo missense variant (c.379C > T; p.Arg127Cys). PCR-based STRs analysis verified the paternity. And long-read SMRT sequencing phased the two variants in trans and identified the paternal origin of the de novo variant. The genetic diagnosis clarified the genetic etiology of the proband and assisted in patient management and counseling. ConclusionWe identified a rare combination of one inherited plus one de novo variant of PPA2 in a patient with autosomal recessive SCFI, which expanded the mutation spectrum of PPA2 and demonstrated the power of target long-read sequencing to make up the diagnostic gap of prevailing NGS.

Phased Haplotype Resolution of the SLC6A4 Promoter Using Long-Read Single Molecule Real-Time (SMRT) Sequencing.

The SLC6A4 gene has been implicated in psychiatric disorder susceptibility and antidepressant response variability. The SLC6A4 promoter is defined by a variable number of homologous 20–24 bp repeats (5-HTTLPR), and long (L) and short (S) alleles are associated with higher and lower expression, respectively. However, this insertion/deletion variant is most informative when considered as a haplotype with the rs25531 and rs25532 variants. Therefore, we developed a long-read single molecule real-time (SMRT) sequencing method to interrogate the SLC6A4 promoter region. A total of 120 samples were subjected to SLC6A4 long-read SMRT sequencing, primarily selected based on available short-read sequencing data. Short-read genome sequencing from the 1000 Genomes (1KG) Project (~5X) and the Genetic Testing Reference Material Coordination Program (~45X), as well as high-depth short-read capture-based sequencing (~330X), could not identify the 5-HTTLPR short (S) allele, nor could short-read sequencing phase any identified variants. In contrast, long-read SMRT sequencing unambiguously identified the 5-HTTLPR short (S) allele (frequency of 0.467) and phased SLC6A4 promoter haplotypes. Additionally, discordant rs25531 genotypes were reviewed and determined to be short-read errors. Taken together, long-read SMRT sequencing is an innovative and robust method for phased resolution of the SLC6A4 promoter, which could enable more accurate pharmacogenetic testing for both research and clinical applications.

Highly accurate mtGenome haplotypes from long-read SMRT sequencing can distinguish between monozygotic twins

Discriminating between monozygotic twins (MZ) remains a challenge in the field of forensics globally. It is very difficult to find sequence variants within MZ twins, despite using ultra-deep next generation sequencing (NGS) for nuclear DNA. However, mitochondrial DNA might be a potential marker owing to its higher mutation rate and easier sequencing via NGS. Here, we aimed to introduce a long-read single molecule real-time sequencing (SMRT) strategy, with better continuity and fewer alignment errors, to obtain more accurate mitochondrial genome (mtGenome) sequence on the Sequel platform. Compared to Ion Torrent Personal Genome Machine (PGM), the long-read SMRT sequencing strategy generated highly accurate and mapped circular consensus sequence (CCS) reads and exhibited robust performance in terms of reliable repeatability, consistent coverage pattern, and balanced strands in mtGenome recovery. Moreover, the long-read SMRT strategy exhibited superior ability to not only identify accurate haplotypes but also discover a total of 785 low-level variants within 16 MZ twin pairs with threshold of 2% and 20 CCS reads with Q30 quality. Taken together, our findings suggested the long-read SMRT technology as an appreciable strategy for obtaining accurate mitotypes and providing a promising solution for distinguishing between MZ twin pairs in forensic genetics.

Discovery of Neuropeptides in the Ganglia of the Pacific Abalone Haliotis discus hannai Using Long-Read Transcriptome Analysis

Advances in the optimization of abalone culture can be achieved through a deeper understanding of the molecules modulating physiological functions including growth, metabolism, and reproduction. To identify neuropeptide molecules that play key roles in the modulation, transcriptome analysis of the ganglia (e.g., cerebral and pleuropedal ganglia) of the abalone Haliotis discus hannai was performed. The long-read single-molecule real-time sequencing using the PacBio RSII platform generated 48,344 full-length transcripts, with an N50 length of 3,020 nucleotides (nt), a GC ratio of 44.76%, and a size range of 309–17,545 nt. Of these transcripts, 32,901 (68.1%), 40,515 (83.8%), 7,591 (15.7%), and 14,945 (30.9%) were annotated using the NCBI nonredundant, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Cluster of Orthologous Groups databases, respectively. Using in silico data mining, 18 putative neuropeptide genes were identified from the transcriptome of H. discus hannai. Comparative analysis revealed that the sequences of identified neuropeptides of H. discus hannai were conserved with other molluscan neuropeptide homologs; however, some features of H. discus hannai neuropeptides were found distinct from those of other molluscan neuropeptides. Interestingly, this is the first study to report a precursor encoding a DYamide peptide in molluscan, which is likely to be a homolog of DLamide peptides, which have been discovered only in annelids to date. This study provides comprehensive neuronal transcriptome data for H. discus hannai, which could help improve aquaculture methods both practically and economically by providing a basis for neuroendocrine-based strategies.

Comparison of different technologies for the decipherment of the whole genome sequence of Campylobacter jejuni BfR-CA-14430

BackgroundCampylobacter jejuni is a zoonotic pathogen that infects the human gut through the food chain mainly by consumption of undercooked chicken meat, raw chicken cross-contaminated ready-to-eat food or by raw milk. In the last decades, C. jejuni has increasingly become the most common bacterial cause for food-born infections in high income countries, costing public health systems billions of euros each year. Currently, different whole genome sequencing techniques such as short-read bridge amplification and long-read single molecule real-time sequencing techniques are applied for in-depth analysis of bacterial species, in particular, Illumina MiSeq, PacBio and MinION.ResultsIn this study, we analyzed a recently isolated C. jejuni strain from chicken meat by short- and long-read data from Illumina, PacBio and MinION sequencing technologies. For comparability, this strain is used in the German PAC-CAMPY research consortium in several studies, including phenotypic analysis of biofilm formation, natural transformation and in vivo colonization models. The complete assembled genome sequence most likely consists of a chromosome of 1,645,980 bp covering 1665 coding sequences as well as a plasmid sequence with 41,772 bp that encodes for 46 genes. Multilocus sequence typing revealed that the strain belongs to the clonal complex CC-21 (ST-44) which is known to be involved in C. jejuni human infections, including outbreaks. Furthermore, we discovered resistance determinants and a point mutation in the DNA gyrase (gyrA) that render the bacterium resistant against ampicillin, tetracycline and (fluoro-)quinolones.ConclusionThe comparison of Illumina MiSeq, PacBio and MinION sequencing and analyses with different assembly tools enabled us to reconstruct a complete chromosome as well as a circular plasmid sequence of the C. jejuni strain BfR-CA-14430. Illumina short-read sequencing in combination with either PacBio or MinION can substantially improve the quality of the complete chromosome and epichromosomal elements on the level of mismatches and insertions/deletions, depending on the assembly program used.

Methicillin-Resistant Staphylococcus aureus Blood Isolates Harboring a Novel Pseudo-staphylococcal Cassette Chromosome mec Element.

The aim of this work was to assess a novel pseudo-staphylococcal cassette chromosome mec (ΨSCCmec) element in methicillin-resistant Staphylococcus aureus (MRSA) blood isolates. Community-associated MRSA E16SA093 and healthcare-associated MRSA F17SA003 isolates were recovered from the blood specimens of patients with S. aureus bacteremia in 2016 and in 2017, respectively. Antimicrobial susceptibility was determined via the disk diffusion method, and SCCmec typing was conducted by multiplex polymerase chain reaction. Whole genome sequencing was carried out by single molecule real-time long-read sequencing. Both isolates belonged to sequence type 72 and agr-type I, and they were negative for Panton-Valentine leukocidin and toxic shock syndrome toxin. The spa-types of E16SA093 and F17SA003 were t324 and t2460, respectively. They had a SCCmec IV-like element devoid of the cassette chromosome recombinase (ccr) gene complex, designated as ΨSCCmecE16SA093. The element was manufactured from SCCmec type IV and the deletion of the ccr gene complex and a 7.0- and 31.9-kb portion of each chromosome. The deficiency of the ccr gene complex in the SCCmec unit is likely resulting in mobility loss, which would be an adaptive evolutionary mechanism. The dissemination of this clone should be monitored closely.

The genome assembly of the fungal pathogen Pyrenochaeta lycopersici from Single-Molecule Real-Time sequencing sheds new light on its biological complexity.

The first draft genome sequencing of the non-model fungal pathogen Pyrenochaeta lycopersici showed an expansion of gene families associated with heterokaryon incompatibility and lacking of mating-type genes, providing insights into the genetic basis of this “imperfect” fungus which lost the ability to produce the sexual stage. However, due to the Illumina short-read technology, the draft genome was too fragmented to allow a comprehensive characterization of the genome, especially of the repetitive sequence fraction. In this work, the sequencing of another P. lycopersici isolate using long-read Single Molecule Real-Time sequencing technology was performed with the aim of obtaining a gapless genome. Indeed, a gapless genome assembly of 62.7 Mb was obtained, with a fraction of repetitive sequences representing 30% of the total bases. The gene content of the two P. lycopersici isolates was very similar, and the large difference in genome size (about 8 Mb) might be attributable to the high fraction of repetitive sequences detected for the new sequenced isolate. The role of repetitive elements, including transposable elements, in modulating virulence effectors is well established in fungal plant pathogens. Moreover, transposable elements are of fundamental importance in creating and re-modelling genes, especially in imperfect fungi. Their abundance in P. lycopersici, together with the large expansion of heterokaryon incompatibility genes in both sequenced isolates, suggest the presence of possible mechanisms alternative to gene re-assorting mediated by sexual recombination. A quite large fraction (~9%) of repetitive elements in P. lycopersici, has no homology with known classes, strengthening this hypothesis. The availability of a gapless genome of P. lycopersici allowed the in-depth analysis of its genome content, by annotating functional genes and TEs. This goal will be an important resource for shedding light on the evolution of the reproductive and pathogenic behaviour of this soilborne pathogen and the onset of a possible speciation within this species.

Cytogenomic identification and long-read single molecule real-time (SMRT) sequencing of a Bardet\u2013Biedl Syndrome 9 (BBS9) deletion

Bardet–Biedl syndrome (BBS) is a recessive disorder characterized by heterogeneous clinical manifestations, including truncal obesity, rod-cone dystrophy, renal anomalies, postaxial polydactyly, and variable developmental delays. At least 20 genes have been implicated in BBS, and all are involved in primary cilia function. We report a 1-year-old male child from Guyana with obesity, postaxial polydactyly on his right foot, hypotonia, ophthalmologic abnormalities, and developmental delay, which together indicated a clinical diagnosis of BBS. Clinical chromosomal microarray (CMA) testing and high-throughput BBS gene panel sequencing detected a homozygous 7p14.3 deletion of exons 1–4 of BBS9 that was encompassed by a 17.5 Mb region of homozygosity at chromosome 7p14.2–p21.1. The precise breakpoints of the deletion were delineated to a 72.8 kb region in the proband and carrier parents by third-generation long-read single molecule real-time (SMRT) sequencing (Pacific Biosciences), which suggested non-homologous end joining as a likely mechanism of formation. Long-read SMRT sequencing of the deletion breakpoints also determined that the aberration included the neighboring RP9 gene implicated in retinitis pigmentosa; however, the clinical significance of this was considered uncertain given the paucity of reported cases with unambiguous RP9 mutations. Taken together, our study characterized a BBS9 deletion, and the identification of this shared haplotype in the parents suggests that this pathogenic aberration may be a BBS founder mutation in the Guyanese population. Importantly, this informative case also highlights the utility of long-read SMRT sequencing to map nucleotide breakpoints of clinically relevant structural variants.

HIV-1 Infection of Primary CD4+ T Cells Regulates the Expression of Specific Human Endogenous Retrovirus HERV-K (HML-2) Elements.

Endogenous retroviruses (ERVs) occupy extensive regions of the human genome. Although many of these retroviral elements have lost their ability to replicate, those whose insertion took place more recently, such as the HML-2 group of HERV-K elements, still retain intact open reading frames and the capacity to produce certain viral RNA and/or proteins. Transcription of these ERVs is, however, tightly regulated by dedicated epigenetic control mechanisms. Nonetheless, it has been reported that some pathological states, such as viral infections and certain cancers, coincide with ERV expression, suggesting that transcriptional reawakening is possible. HML-2 elements are reportedly induced during HIV-1 infection, but the conserved nature of these elements has, until recently, rendered their expression profiling problematic. Here, we provide comprehensive HERV-K HML-2 expression profiles specific for productively HIV-1-infected primary human CD4+ T cells. We combined enrichment of HIV-1 infected cells using a reporter virus expressing a surface reporter for gentle and efficient purification with long-read single-molecule real-time sequencing. We show that three HML-2 proviruses-6q25.1, 8q24.3, and 19q13.42-are upregulated on average between 3- and 5-fold in HIV-1-infected CD4+ T cells. One provirus, HML-2 12q24.33, in contrast, was repressed in the presence of active HIV replication. In conclusion, this report identifies the HERV-K HML-2 loci whose expression profiles differ upon HIV-1 infection in primary human CD4+ T cells. These data will help pave the way for further studies on the influence of endogenous retroviruses on HIV-1 replication.IMPORTANCE Endogenous retroviruses inhabit big portions of our genome. Moreover, although they are mainly inert, some of the evolutionarily younger members maintain the ability to express both RNA and proteins. We have developed an approach using long-read single-molecule real-time (SMRT) sequencing that produces long reads that allow us to obtain detailed and accurate HERV-K HML-2 expression profiles. We applied this approach to study HERV-K expression in the presence or absence of productive HIV-1 infection of primary human CD4+ T cells. In addition to using SMRT sequencing, our strategy also includes the magnetic selection of the infected cells so that levels of background expression due to uninfected cells are kept at a minimum. The results presented here provide a blueprint for in-depth studies of the interactions of the authentic upregulated HERV-K HML-2 elements and HIV-1.

Characterization of the Dynamic Transcriptome of a Herpesvirus with Long-read Single Molecule Real-Time Sequencing

Herpesvirus gene expression is co-ordinately regulated and sequentially ordered during productive infection. The viral genes can be classified into three distinct kinetic groups: immediate-early, early, and late classes. In this study, a massively parallel sequencing technique that is based on PacBio Single Molecule Real-time sequencing platform, was used for quantifying the poly(A) fraction of the lytic transcriptome of pseudorabies virus (PRV) throughout a 12-hour interval of productive infection on PK-15 cells. Other approaches, including microarray, real-time RT-PCR and Illumina sequencing are capable of detecting only the aggregate transcriptional activity of particular genomic regions, but not individual herpesvirus transcripts. However, SMRT sequencing allows for a distinction between transcript isoforms, including length- and splice variants, as well as between overlapping polycistronic RNA molecules. The non-amplified Isoform Sequencing (Iso-Seq) method was used to analyse the kinetic properties of the lytic PRV transcripts and to then classify them accordingly. Additionally, the present study demonstrates the general utility of long-read sequencing for the time-course analysis of global gene expression in practically any organism.

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.

Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data

We present a hierarchical genome-assembly process (HGAP) for high-quality de novo microbial genome assemblies using only a single, long-insert shotgun DNA library in conjunction with Single Molecule, Real-Time (SMRT) DNA sequencing. Our method uses the longest reads as seeds to recruit all other reads for construction of highly accurate preassembled reads through a directed acyclic graph-based consensus procedure, which we follow with assembly using off-the-shelf long-read assemblers. In contrast to hybrid approaches, HGAP does not require highly accurate raw reads for error correction. We demonstrate efficient genome assembly for several microorganisms using as few as three SMRT Cell zero-mode waveguide arrays of sequencing and for BACs using just one SMRT Cell. Long repeat regions can be successfully resolved with this workflow. We also describe a consensus algorithm that incorporates SMRT sequencing primary quality values to produce de novo genome sequence exceeding 99.999% accuracy.