SARS-CoV-2 variants of concern surveillance including Omicron using RT-PCR-based genotyping offers comparable performance to whole genome sequencing.

Known SARS-CoV-2 variants of concern (VOCs) can be detected and differentiated using an RT-PCR-based genotyping approach, which offers quicker time to result, lower cost, higher flexibility, and use of the same laboratory instrumentation for detection of SARS-CoV-2 when compared with whole genome sequencing (WGS). In the current study, we demonstrate how we applied a genotyping approach for identification of all VOCs and that such technique can offer comparable performance to WGS for identification of known SARS-CoV-2 VOCs, including more recent strains, Omicron BA.1 and BA.2.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

Developmental validation of the Quantifiler(®) HP and Trio Kits for human DNA quantification in forensic samples.

The quantification of human genomic DNA is a necessary first step in the DNA casework sample analysis workflow. DNA quantification determines optimal sample input amounts for subsequent STR (short tandem repeat) genotyping procedures, as well as being a useful screening tool to identify samples most likely to provide probative genotypic evidence. To better mesh with the capabilities of newest-generation STR analysis assays, the Quantifiler(®) HP and Quantifiler(®) Trio DNA Quantification Kits were designed for greater detection sensitivity and more robust performance with samples that contain PCR inhibitors or degraded DNA. The new DNA quantification kits use multiplex TaqMan(®) assay-based fluorescent probe technology to simultaneously quantify up to three human genomic targets, allowing samples to be assessed for total human DNA, male contributor (i.e., Y-chromosome) DNA, as well as a determination of DNA degradation state. The Quantifiler HP and Trio Kits use multiple-copy loci to allow for significantly improved sensitivity compared to earlier-generation kits that employ single-copy target loci. The kits' improved performance provides better predictive ability for results with downstream, newest-generation STR assays, and their shortened time-to-result allows more efficient integration into the forensic casework analysis workflow.

Using Whole Genome Analysis to Examine Recombination across Diverse Sequence Types of Staphylococcus aureus.

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

Dominance of multidrug resistant CC271 clones in macrolide-resistant streptococcus pneumoniae in Arizona.

BACKGROUND: Rates of resistance to macrolide antibiotics in Streptococcus pneumoniae are rising around the world due to the spread of mobile genetic elements harboring mef(E) and erm(B) genes and post-vaccine clonal expansion of strains that carry them. RESULTS: Characterization of 592 clinical isolates collected in Arizona over a 10 year period shows 23.6% are macrolide resistant. The largest portion of the macrolide-resistant population, 52%, is dual mef(E)/erm(B)-positive. All dual-positive isolates are multidrug-resistant clonal lineages of Taiwan19F-14, mostly multilocus sequence type 320, carrying the recently described transposon Tn2010. The remainder of the macrolide resistant S. pneumoniae collection includes 31% mef(E)-positive, and 9% erm(B)-positive strains. CONCLUSIONS: The dual-positive, multidrug-resistant S. pneumoniae clones have likely expanded by switching to non-vaccine serotypes after the heptavalent pneumococcal conjugate vaccine release, and their success limits therapy options. This upsurge could have a considerable clinical impact in Arizona.

Prospective genomic characterization of the German enterohemorrhagic Escherichia coli O104:H4 outbreak by rapid next generation sequencing technology.

An ongoing outbreak of exceptionally virulent Shiga toxin (Stx)-producing Escherichia coli O104:H4 centered in Germany, has caused over 830 cases of hemolytic uremic syndrome (HUS) and 46 deaths since May 2011. Serotype O104:H4, which has not been detected in animals, has rarely been associated with HUS in the past. To prospectively elucidate the unique characteristics of this strain in the early stages of this outbreak, we applied whole genome sequencing on the Life Technologies Ion Torrent PGM™ sequencer and Optical Mapping to characterize one outbreak isolate (LB226692) and a historic O104:H4 HUS isolate from 2001 (01-09591). Reference guided draft assemblies of both strains were completed with the newly introduced PGM™ within 62 hours. The HUS-associated strains both carried genes typically found in two types of pathogenic E. coli, enteroaggregative E. coli (EAEC) and enterohemorrhagic E. coli (EHEC). Phylogenetic analyses of 1,144 core E. coli genes indicate that the HUS-causing O104:H4 strains and the previously published sequence of the EAEC strain 55989 show a close relationship but are only distantly related to common EHEC serotypes. Though closely related, the outbreak strain differs from the 2001 strain in plasmid content and fimbrial genes. We propose a model in which EAEC 55989 and EHEC O104:H4 strains evolved from a common EHEC O104:H4 progenitor, and suggest that by stepwise gain and loss of chromosomal and plasmid-encoded virulence factors, a highly pathogenic hybrid of EAEC and EHEC emerged as the current outbreak clone. In conclusion, rapid next-generation technologies facilitated prospective whole genome characterization in the early stages of an outbreak.

Design and validation of a novel multiplex real-time PCR assay for Vibrio pathogen detection.

Three species--Vibrio cholerae, Vibrio parahaemolyticus, and Vibrio vulnificus--account for the majority of vibrio infections in humans. Rapid and accurate identification of Vibrio species has been problematic because phenotypic characteristics are variable within species. Additionally, biochemical identification and confirmation require 2 or more days to complete. Rapid and sensitive molecular techniques for the detection of vibrio pathogens would be useful for the surveillance and management of outbreaks. To facilitate the identification of human-pathogenic species, we designed and validated a highly sensitive, specific, and robust multiplex real-time PCR assay to identify V. cholerae, V. parahaemolyticus, and V. vulnificus using a four-dye configuration in a convenient lyophilized format. Multiple Vibrio strains were sequenced to verify candidate target TaqMan sites. Several individual assays within the multiplex contain multiple primers or probes to ensure detection of polymorphic variants. V. cholerae, V. parahaemolyticus, and V. vulnificus were detected either individually or in mixtures at ≤30 genomic copies. V. cholerae was specifically detected in the presence or absence of Vibrio mimicus. The Vibrio multiplex assay showed 100% specificity to all targets analyzed and no detection of nearest neighbor strains. Each assay exhibited 100% ± 10% efficiency. Multiplex real-time PCR can simplify pathogen detection and reduce costs per test since three species can be analyzed in a single reaction tube. Rapid and accurate detection of pathogenic vibrios in shellfish or seawater samples will improve the microbiological safety of seafood for consumers.

Tumor transcriptome sequencing reveals allelic expression imbalances associated with copy number alterations.

Due to growing throughput and shrinking cost, massively parallel sequencing is rapidly becoming an attractive alternative to microarrays for the genome-wide study of gene expression and copy number alterations in primary tumors. The sequencing of transcripts (RNA-Seq) should offer several advantages over microarray-based methods, including the ability to detect somatic mutations and accurately measure allele-specific expression. To investigate these advantages we have applied a novel, strand-specific RNA-Seq method to tumors and matched normal tissue from three patients with oral squamous cell carcinomas. Additionally, to better understand the genomic determinants of the gene expression changes observed, we have sequenced the tumor and normal genomes of one of these patients. We demonstrate here that our RNA-Seq method accurately measures allelic imbalance and that measurement on the genome-wide scale yields novel insights into cancer etiology. As expected, the set of genes differentially expressed in the tumors is enriched for cell adhesion and differentiation functions, but, unexpectedly, the set of allelically imbalanced genes is also enriched for these same cancer-related functions. By comparing the transcriptomic perturbations observed in one patient to his underlying normal and tumor genomes, we find that allelic imbalance in the tumor is associated with copy number mutations and that copy number mutations are, in turn, strongly associated with changes in transcript abundance. These results support a model in which allele-specific deletions and duplications drive allele-specific changes in gene expression in the developing tumor.

An efficient and high-throughput approach for experimental validation of novel human gene predictions.

A highly automated RT-PCR-based approach has been established to validate novel human gene predictions with no prior experimental evidence of mRNA splicing (ab initio predictions). Ab initio gene predictions were selected for high-throughput validation using predicted protein classification, sequence similarity to other genomes, colocalization with an MPSS tag, or microarray expression. Initial microarray prioritization followed by RT-PCR validation was the most efficient combination, resulting in approximately 35% of the ab initio predictions being validated by RT-PCR. Of the 7252 novel genes that were prioritized and processed, 796 constituted real transcripts. In addition, high-throughput RACE successfully extended the 5' and/or 3' ends of >60% of RT-PCR-validated genes. Reevaluation of these transcripts produced 574 novel transcripts using RefSeq as a reference. RT-PCR sequencing in combination with RACE on ab initio gene predictions could be used to define the transcriptome across all species.

In situ-synthesized novel microarray optimized for mouse stem cell and early developmental expression profiling.

Applications of microarray technologies to mouse embryology/genetics have been limited, due to the nonavailability of microarrays containing large numbers of embryonic genes and the gap between microgram quantities of RNA required by typical microarray methods and the miniscule amounts of tissue available to researchers. To overcome these problems, we have developed a microarray platform containing in situ-synthesized 60-mer oligonucleotide probes representing approximately 22,000 unique mouse transcripts, assembled primarily from sequences of stem cell and embryo cDNA libraries. We have optimized RNA labeling protocols and experimental designs to use as little as 2 ng total RNA reliably and reproducibly. At least 98% of the probes contained in the microarray correspond to clones in our publicly available collections, making cDNAs readily available for further experimentation on genes of interest. These characteristics, combined with the ability to profile very small samples, make this system a resource for stem cell and embryogenomics research.