Utility Analyses of AVITI Sequencing Chemistry.

BACKGROUND: DNA sequencing is a critical tool in modern biology. Over the last two decades, it has been revolutionized by the advent of massively parallel sequencing, leading to significant advances in the genome and transcriptome sequencing of various organisms. Nevertheless, challenges with accuracy, lack of competitive options and prohibitive costs associated with high throughput parallel short-read sequencing persist. RESULTS: Here, we conduct a comparative analysis using matched DNA and RNA short-reads assays between Element Biosciences AVITI chemistry and Illumina NextSeq 550. Similar comparisons were evaluated for synthetic long-read sequencing for RNA and targeted single-cell transcripts between the AVITI and Illumina NovaSeq 6000. For both DNA and RNA short-read applications, the study found that the AVITI produced significantly higher per sequence quality scores. For PCR-free DNA libraries, we observed up to a 10-fold lower experimentally determined error rate for using the AVITI chemistry compared to the NextSeq 550. For short-read RNA quantification, both AVITI and the NextSeq 550 demonstrated comparable accuracy. With regards to synthetic long-read mRNA and targeted synthetic long read single cell mRNA sequencing, both platforms respective chemistries performed comparably in quantification of genes and isoforms. The AVITI displayed a marginally lower error rate for long reads, with fewer chemistry-specific errors and a higher mutation detection rate. CONCLUSION: These results point to the potential of the AVITI platform as a competitive candidate in high-throughput short read sequencing analyses when juxtaposed with the Illumina NextSeq 550.

Long-read single-cell sequencing reveals expressions of hypermutation clusters of isoforms in human liver cancer cells.

The protein diversity of mammalian cells is determined by arrays of isoforms from genes. Genetic mutation is essential in species evolution and cancer development. Accurate long-read transcriptome sequencing at single-cell level is required to decipher the spectrum of protein expressions in mammalian organisms. In this report, we developed a synthetic long-read single-cell sequencing technology based on LOOPSeq technique. We applied this technology to analyze 447 transcriptomes of hepatocellular carcinoma (HCC) and benign liver from an individual. Through Uniform Manifold Approximation and Projection analysis, we identified a panel of mutation mRNA isoforms highly specific to HCC cells. The evolution pathways that led to the hyper-mutation clusters in single human leukocyte antigen molecules were identified. Novel fusion transcripts were detected. The combination of gene expressions, fusion gene transcripts, and mutation gene expressions significantly improved the classification of liver cancer cells versus benign hepatocytes. In conclusion, LOOPSeq single-cell technology may hold promise to provide a new level of precision analysis on the mammalian transcriptome.

Long-read single-cell sequencing reveals expressions of hypermutation clusters of isoforms in human liver cancer cells.

The protein diversity of mammalian cells is determined by arrays of isoforms from genes. Genetic mutation is essential in species evolution and cancer development. Accurate Long-read transcriptome sequencing at single-cell level is required to decipher the spectrum of protein expressions in mammalian organisms. In this report, we developed a synthetic long-read single-cell sequencing technology based on LOOPseq technique. We applied this technology to analyze 447 transcriptomes of hepatocellular carcinoma (HCC) and benign liver from an individual. Through Uniform Manifold Approximation and Projection (UMAP) analysis, we identified a panel of mutation mRNA isoforms highly specific to HCC cells. The evolution pathways that led to the hyper-mutation clusters in single human leukocyte antigen (HLA) molecules were identified. Novel fusion transcripts were detected. The combination of gene expressions, fusion gene transcripts, and mutation gene expressions significantly improved the classification of liver cancer cells versus benign hepatocytes. In conclusion, LOOPseq single-cell technology may hold promise to provide a new level of precision analysis on the mammalian transcriptome.

Genomic analyses of pneumococci from children with sickle cell disease expose host-specific bacterial adaptations and deficits in current interventions.

Sickle cell disease (SCD) patients are at high risk of contracting pneumococcal infection. To address this risk, they receive pneumococcal vaccines, and antibiotic prophylaxis and treatment. To assess the impact of SCD and these interventions on pneumococcal genetic architecture, we examined the genomes of more than 300 pneumococcal isolates from SCD patients over 20 years. Modern SCD strains retained invasive capacity but shifted away from the serotypes used in vaccines. These strains had specific genetic changes related to antibiotic resistance, capsule biosynthesis, metabolism, and metal transport. A murine SCD model coupled with Tn-seq mutagenesis identified 60 noncapsular pneumococcal genes under differential selective pressure in SCD, which correlated with aspects of SCD pathophysiology. Further, virulence determinants in the SCD context were distinct from the general population, and protective capacity of potential antigens was lost over time in SCD. This highlights the importance of understanding bacterial pathogenesis in the context of high-risk individuals.

Control of virulence by small RNAs in Streptococcus pneumoniae.

Small noncoding RNAs (sRNAs) play important roles in gene regulation in both prokaryotes and eukaryotes. Thus far, no sRNA has been assigned a definitive role in virulence in the major human pathogen Streptococcus pneumoniae. Based on the potential coding capacity of intergenic regions, we hypothesized that the pneumococcus produces many sRNAs and that they would play an important role in pathogenesis. We describe the application of whole-genome transcriptional sequencing to systematically identify the sRNAs of Streptococcus pneumoniae. Using this approach, we have identified 89 putative sRNAs, 56 of which are newly identified. Furthermore, using targeted genetic approaches and Tn-seq transposon screening, we demonstrate that many of the identified sRNAs have important global and niche-specific roles in virulence. These data constitute the most comprehensive analysis of pneumococcal sRNAs and provide the first evidence of the extensive roles of sRNAs in pneumococcal pathogenesis.

Chlamydia psittaci genetic variants differ in virulence by modulation of host immunity.

BACKGROUND: Psittacosis is a zoonosis caused by Chlamydia psittaci and is characterized by severe pneumonia and systemic infection. We sought to determine the basis of the 1000-fold difference in lethal dose of 2 C. psittaci 6BC strains in mice. METHODS: Genomes of the strains were sequenced. Mice were infected intraperitoneally and the growth kinetics, immune responses, and pathology were compared. RESULTS: The 2 strains differed by the presence of a 7.5-kb plasmid in the attenuated strain and 7 nonsynonomous single-nucleotide polymorphisms between the chromosomes, including a serine/threonine protein kinase gene pkn5. The plasmid was cured from the attenuated strain, but it remained nonlethal. Strains did not differ in growth kinetics in vitro or in vivo. Infection with the attenuated strain led to influx of activated macrophages with relatively minor organ damage. In contrast, the virulent strain caused an influx of nonactivated macrophages, neutrophils, and significant end organ damage. Mice infected with the virulent strain survived challenge when coinfected with either the plasmid-positive or plasmid-negative attenuated strain, indicating that an active process elicited by the attenuated strain reduces inflammation and disease. CONCLUSIONS: C. psittaci modulates virulence by alteration of host immunity, which is conferred by small differences in the chromosome.

Identification of errors introduced during high throughput sequencing of the T cell receptor repertoire.

BACKGROUND: Recent advances in massively parallel sequencing have increased the depth at which T cell receptor (TCR) repertoires can be probed by >3log10, allowing for saturation sequencing of immune repertoires. The resolution of this sequencing is dependent on its accuracy, and direct assessments of the errors formed during high throughput repertoire analyses are limited. RESULTS: We analyzed 3 monoclonal TCR from TCR transgenic, Rag-/- mice using Illumina® sequencing. A total of 27 sequencing reactions were performed for each TCR using a trifurcating design in which samples were divided into 3 at significant processing junctures. More than 20 million complementarity determining region (CDR) 3 sequences were analyzed. Filtering for lower quality sequences diminished but did not eliminate sequence errors, which occurred within 1-6% of sequences. Erroneous sequences were pre-dominantly of correct length and contained single nucleotide substitutions. Rates of specific substitutions varied dramatically in a position-dependent manner. Four substitutions, all purine-pyrimidine transversions, predominated. Solid phase amplification and sequencing rather than liquid sample amplification and preparation appeared to be the primary sources of error. Analysis of polyclonal repertoires demonstrated the impact of error accumulation on data parameters. CONCLUSIONS: Caution is needed in interpreting repertoire data due to potential contamination with mis-sequence reads. However, a high association of errors with phred score, high relatedness of erroneous sequences with the parental sequence, dominance of specific nt substitutions, and skewed ratio of forward to reverse reads among erroneous sequences indicate approaches to filter erroneous sequences from repertoire data sets.

Examination of type IV pilus expression and pilus-associated phenotypes in Kingella kingae clinical isolates.

Kingella kingae is a gram-negative bacterium that is being recognized increasingly as a cause of septic arthritis and osteomyelitis in young children. Previous work established that K. kingae expresses type IV pili that mediate adherence to respiratory epithelial and synovial cells. PilA1 is the major pilus subunit in K. kingae type IV pili and is essential for pilus assembly. To develop a better understanding of the role of K. kingae type IV pili during colonization and invasive disease, we examined a collection of clinical isolates for pilus expression and in vitro adherence. In addition, in a subset of isolates we performed nucleotide sequencing to assess the level of conservation of PilA1. The majority of respiratory and nonendocarditis blood isolates were piliated, while the majority of joint fluid, bone, and endocarditis blood isolates were nonpiliated. The piliated isolates formed either spreading/corroding or nonspreading/noncorroding colonies and were uniformly adherent, while the nonpiliated isolates formed domed colonies and were nonadherent. PilA1 sequence varied significantly from strain to strain, resulting in substantial variability in antibody reactivity. These results suggest that type IV pili may confer a selective advantage on K. kingae early in infection and a selective disadvantage on K. kingae at later stages in the pathogenic process. We speculate that PilA1 is immunogenic during natural infection and undergoes antigenic variation to escape the immune response.

Emergence of H5N1 avian influenza viruses with reduced sensitivity to neuraminidase inhibitors and novel reassortants in Lao People's Democratic Republic.

Pandemic influenza viruses can emerge through continuous evolution and the acquisition of specific mutations or through reassortment. This study assessed the pandemic potential of H5N1 viruses isolated from poultry outbreaks occurring from July 2006 to September 2008 in the Lao People's Democratic Republic (PDR). We analyzed 29 viruses isolated from chickens and ducks and two from fatal human cases in 2007. Prior to 2008, all H5N1 isolates in Lao PDR were from clade 2.3.4; however, clade 2.3.2 was introduced in September 2008. Of greatest concern was the circulation of three isolates that showed reduced sensitivity to the neuraminidase (NA) inhibitor oseltamivir in an enzyme inhibition assay, each with different NA mutations - V116A, I222L and K150N, and a previously unreported S246N mutation. In addition, six isolates had an S31N mutation in the M2 protein, which conferred resistance to amantadine not previously reported in clade 2.3.4 viruses. Two H5N1 reassortants were isolated whose polymerase genes, PB1 and PB2, were homologous to those of Eurasian viruses giving rise to a novel H5N1 genotype, genotype P. All H5N1 viruses retained avian-like receptor specificity, but four had altered affinities for alpha2,3-linked sialic acid. This study shows that, in a genetically similar population of H5N1 viruses in Lao PDR, mutants emerged with natural resistance to antivirals and altered affinities for alpha2,3-linked sialic acids, together with reassortants with polymerase genes homologous to Eurasian viruses. These changes may contribute to the emergence of a pandemic influenza strain and are critical in devising surveillance strategies.

Genomic analysis reveals few genetic alterations in pediatric acute myeloid leukemia.

Pediatric de novo acute myeloid leukemia (AML) is an aggressive malignancy with current therapy resulting in cure rates of only 60%. To better understand the cause of the marked heterogeneity in therapeutic response and to identify new prognostic markers and therapeutic targets a comprehensive list of the genetic mutations that underlie the pathogenesis of AML is needed. To approach this goal, we examined diagnostic leukemic samples from a cohort of 111 children with de novo AML using single-nucleotide-polymorphism microarrays and candidate gene resequencing. Our data demonstrate that, in contrast to pediatric acute lymphoblastic leukemia (ALL), de novo AML is characterized by a very low burden of genomic alterations, with a mean of only 2.38 somatic copy-number alterations per leukemia, and less than 1 nonsynonymous point mutation per leukemia in the 25 genes analyzed. Even more surprising was the observation that 34% of the leukemias lacked any identifiable copy-number alterations, and 28% of the leukemias with recurrent translocations lacked any identifiable sequence or numerical abnormalities. The only exception to the presence of few mutations was acute megakaryocytic leukemias, with the majority of these leukemias being characterized by a high number of copy-number alterations but rare point mutations. Despite the low overall number of lesions across the patient cohort, novel recurring regions of genetic alteration were identified that harbor known, and potential new cancer genes. These data reflect a remarkably low burden of genomic alterations within pediatric de novo AML, which is in stark contrast to most other human malignancies.

Zoonotic potential of highly pathogenic avian H7N3 influenza viruses from Pakistan.

H5 and H7 avian influenza viruses can become highly pathogenic in chickens after interspecies transmission. These viruses have transmitted directly to humans from birds in Eurasia and Africa (H5N1), the Netherlands (H7N7), and Canada (H7N3). Here we report antigenic, sequence, and phylogenetic analyses of H7N3 viruses isolated from chickens in Pakistan from 1995 to 2002. We compared the pathogenic and zoonotic potential of the Pakistani viruses in avian and mammalian hosts. In chickens, all of the isolates showed high pathogenicity with poor transmissibility to contact birds. Viral shedding from the trachea and cloaca was equivalent, but cloacal shedding occurred longer; dissemination of virus into the tissues was widespread. In contrast, the viruses replicated poorly in 6-week-old mallard ducks. In mammalian hosts, of the two Pakistani H7N3/02 viruses that caused weight loss, one also caused 40% mortality in mice without prior adaptation, and preliminary experiments in ferrets showed significant virus multiplication in the lungs, intestine, and conjunctiva. We conclude that the H7N3/02 isolates from Pakistan show limited antigenic drift and have evolved slowly during their 8-year circulation in chickens; however, these viruses have the potential to infect mammals.

Influenza in migratory birds and evidence of limited intercontinental virus exchange.

Migratory waterfowl of the world are the natural reservoirs of influenza viruses of all known subtypes. However, it is unknown whether these waterfowl perpetuate highly pathogenic (HP) H5 and H7 avian influenza viruses. Here we report influenza virus surveillance from 2001 to 2006 in wild ducks in Alberta, Canada, and in shorebirds and gulls at Delaware Bay (New Jersey), United States, and examine the frequency of exchange of influenza viruses between the Eurasian and American virus clades, or superfamilies. Influenza viruses belonging to each of the subtypes H1 through H13 and N1 through N9 were detected in these waterfowl, but H14 and H15 were not found. Viruses of the HP Asian H5N1 subtypes were not detected, and serologic studies in adult mallard ducks provided no evidence of their circulation. The recently described H16 subtype of influenza viruses was detected in American shorebirds and gulls but not in ducks. We also found an unusual cluster of H7N3 influenza viruses in shorebirds and gulls that was able to replicate well in chickens and kill chicken embryos. Genetic analysis of 6,767 avian influenza gene segments and 248 complete avian influenza viruses supported the notion that the exchange of entire influenza viruses between the Eurasian and American clades does not occur frequently. Overall, the available evidence does not support the perpetuation of HP H5N1 influenza in migratory birds and suggests that the introduction of HP Asian H5N1 to the Americas by migratory birds is likely to be a rare event.

Assessment of molecular typing methods to determine invasiveness and to differentiate clones of Streptococcus pneumoniae.

In the United States, Streptococcus pneumoniae is the leading cause of community-acquired pneumonia and invasive bacterial disease. As antimicrobial resistance increases, it will become critical to determine if strains circulating in a population are likely to cause invasive pneumococcal disease (IPD). This is possible by comparison of an isolate's genotype to strains known to be invasive. In this work, we compared pulse-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST), comparative genomic hybridization (CGH) and multi-invasive-locus sequence typing (MILST) for their ability to distinguish between known IPD causing and carrier strains using phylogenetic analyses. In addition, we assess the ability of these techniques to determine true clones from highly related strains. The resulting trees suggest that despite similar overall topologies, the clearest picture of invasiveness and genetic relatedness can be viewed when typing methods are used collectively.

Quinupristin-dalfopristin nonsusceptibility in pneumococci from sickle cell disease patients.

Sickle cell disease (SCD) is a risk factor for fatal pneumococcal infection. Nonsusceptibilty to quinupristin-dalfopristin (Q-D) was absent from 105 non-SCD-associated pneumococcal isolates but was present in 33/148 (22%) SCD-associated isolates. One-third of the isolates harbored a known resistance mechanism. Q-D is not optimal for use for the treatment of pneumococcal infection in SCD patients.

Identification of a Candidate Streptococcus pneumoniae core genome and regions of diversity correlated with invasive pneumococcal disease.

Streptococcus pneumoniae is a leading cause of community-acquired pneumonia and gram-positive sepsis. While multiple virulence determinants have been identified, the combination of features that determines the propensity of an isolate to cause invasive pneumococcal disease (IPD) remains unknown. In this study, we determined the genetic composition of 42 invasive and 30 noninvasive clinical isolates of serotypes 6A, 6B, and 14 by comparative genomic hybridization. Comparison of the present/absent gene matrix (i.e., comparative genomic analysis [CGA]) identified a candidate core genome consisting of 1,553 genes (73% of the TIGR4 genome), 154 genes whose presence correlated with the ability to cause IPD, and 176 genes whose presence correlated with the noninvasive phenotype. Genes identified by CGA were cross-referenced with the published signature-tagged mutagenesis studies, which served to identify core and IPD-correlated genes required for in vivo passage. Among these, two pathogenicity islands, region of diversity 8a (RD8a), which encodes a neuraminidase and V-type sodium synthase, and RD10, which encodes PsrP, a protein homologous to the platelet adhesin GspB in Streptococcus gordonii, were identified. Mice infected with a PsrP mutant were delayed in the development of bacteremia and demonstrated reduced mortality versus wild-type-infected controls. Finally, the presence of seven RDs was determined to correlate with the noninvasive phenotype, a finding that suggests some RDs may contribute to asymptomatic colonization. In conclusion, RDs are unequally distributed between invasive and noninvasive isolates, RD8a and RD10 are correlated with the propensity of an isolate to cause IPD, and PsrP is required for full virulence in mice.

Vancomycin stress response in a sensitive and a tolerant strain of Streptococcus pneumoniae.

The vancomycin stress response was studied in Streptococcus pneumoniae strains T4 (TIGR4) and Tupelo. Vancomycin affected the expression of 175 genes, including genes encoding transport functions and enzymes involved in aminosugar metabolism. The two-component systems TCS03, TCS11, and CiaRH also responded to antibiotic treatment. We hypothesize that the three regulons are an important part of the bacterium's response to vancomycin stress.

Cephalosporin-resistant pneumococci and sickle cell disease.

Sickle cell anemia patients have 600 times the risk for invasive pneumococcal disease than their healthy peers. High-level cephalosporin resistance was described in the 1990s in healthy children from Tennessee, but its prevalence in sickle cell disease patients is unknown. Pneumococcal isolates from sickle cell disease patients from Tennessee were subjected to multilocus sequence typing to characterize antimicrobial drug-resistant strains. Twenty-one percent of strains were resistant to cefotaxime and penicillin. Of the 14 cephalosporin-resistant strains, 9 were sequence types previously described as highly cephalosporin resistant, while resistance was found for the first time in 3 clones: Maryland6B, ST660, and a novel clone, ST1753. High-level cephalosporin resistance exists in more settings than initially recognized, and its high prevalence in sickle cell disease patients may decrease the efficacy of third-generation cephalosporins in invasive pneumococcal disease.