A CRISPR-enhanced metagenomic NGS test to improve pandemic preparedness.

The lack of preparedness for detecting and responding to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogen (i.e., COVID-19) has caused enormous harm to public health and the economy. Testing strategies deployed on a population scale at day zero, i.e., the time of the first reported case, would be of significant value. Next-generation sequencing (NGS) has such capabilities; however, it has limited detection sensitivity for low-copy-number pathogens. Here, we leverage the CRISPR-Cas9 system to effectively remove abundant sequences not contributing to pathogen detection and show that NGS detection sensitivity of SARS-CoV-2 approaches that of RT-qPCR. The resulting sequence data can also be used for variant strain typing, co-infection detection, and individual human host response assessment, all in a single molecular and analysis workflow. This NGS work flow is pathogen agnostic and, therefore, has the potential to transform how large-scale pandemic response and focused clinical infectious disease testing are pursued in the future.

T cell subtype profiling measures exhaustion and predicts anti-PD-1 response.

Anti-PD-1 therapy can provide long, durable benefit to a fraction of patients. The on-label PD-L1 test, however, does not accurately predict response. To build a better biomarker, we created a method called T Cell Subtype Profiling (TCSP) that characterizes the abundance of T cell subtypes (TCSs) in FFPE specimens using five RNA models. These TCS RNA models are created using functional methods, and robustly discriminate between naïve, activated, exhausted, effector memory, and central memory TCSs, without the reliance on non-specific, classical markers. TCSP is analytically valid and corroborates associations between TCSs and clinical outcomes. Multianalyte biomarkers based on TCS estimates predicted response to anti-PD-1 therapy in three different cancers and outperformed the indicated PD-L1 test, as well as Tumor Mutational Burden. Given the utility of TCSP, we investigated the abundance of TCSs in TCGA cancers and created a portal to enable researchers to discover other TCSP-based biomarkers.

Analytical Performance of an Immunoprofiling Assay Based on RNA Models.

As immuno-oncology drugs grow more popular in the treatment of cancer, better methods are needed to quantify the tumor immune cell component to determine which patients are most likely to benefit from treatment. Methods such as flow cytometry can accurately assess the composition of infiltrating immune cells; however, they show limited use in formalin-fixed, paraffin-embedded (FFPE) specimens. This article describes a novel hybrid-capture RNA sequencing assay, ImmunoPrism, that estimates the relative percentage abundance of eight immune cell types in FFPE solid tumors. Immune health expression models were generated using machine learning methods and used to uniquely identify each immune cell type using the most discriminatively expressed genes. The analytical performance of the assay was assessed using 101 libraries from 40 FFPE and 32 fresh-frozen samples. With defined samples, ImmunoPrism had a precision of ±2.72%, a total error of 2.75%, and a strong correlation (r2 = 0.81; P < 0.001) to flow cytometry. ImmunoPrism had similar performance in dissociated tumor cell samples (total error of 8.12%) and correlated strongly with immunohistochemistry (CD8: r2 = 0.83; P < 0.001) in FFPE samples. Other performance metrics were determined, including limit of detection, reportable range, and reproducibility. The approach used for analytical validation is shared here so that it may serve as a helpful framework for other laboratories when validating future complex RNA-based assays.

Genomic heterogeneity of ALK fusion breakpoints in non-small-cell lung cancer.

In lung adenocarcinoma, canonical EML4-ALK inversion results in a fusion protein with a constitutively active ALK kinase domain. Evidence of ALK rearrangement occurs in a minority (2-7%) of lung adenocarcinoma, and only ~60% of these patients will respond to targeted ALK inhibition by drugs such as crizotinib and ceritinib. Clinically, targeted anti-ALK therapy is often initiated based on evidence of an ALK genomic rearrangement detected by fluorescence in situ hybridization (FISH) of interphase cells in formalin-fixed, paraffin-embedded tissue sections. At the genomic level, however, ALK rearrangements are heterogeneous, with multiple potential breakpoints in EML4, and alternate fusion partners. Using next-generation sequencing of DNA and RNA together with ALK immunohistochemistry, we comprehensively characterized genomic breakpoints in 33 FISH-positive lung adenocarcinomas. Of these 33 cases, 29 (88%) had detectable DNA level ALK rearrangements involving EML4, KIF5B, or non-canonical partners including ASXL2, ATP6V1B1, PRKAR1A, and SPDYA. A subset of 12 cases had material available for RNA-Seq. Of these, eight of eight (100%) cases with DNA rearrangements showed ALK fusion transcripts from RNA-Seq; three of four cases (75%) without detectable DNA rearrangements were similarly negative by RNA-Seq, and one case was positive by RNA-Seq but negative by DNA next-generation sequencing. By immunohistochemistry, 17 of 19 (89%) tested cases were clearly positive for ALK protein expression; the remaining cases had no detectable DNA level rearrangement or had a non-canonical rearrangement not predicted to form a fusion protein. Survival analysis of patients treated with targeted ALK inhibitors demonstrates a significant difference in mean survival between patients with next-generation sequencing confirmed EML4-ALK rearrangements, and those without (20.6 months vs 5.4 months, P<0.01). Together, these data demonstrate abundant genomic heterogeneity among ALK-rearranged lung adenocarcinoma, which may account for differences in treatment response with targeted ALK inhibitors.

Hybrid capture and next-generation sequencing identify viral integration sites from formalin-fixed, paraffin-embedded tissue.

Although next-generation sequencing (NGS) has been the domain of large genome centers, it is quickly becoming more accessible to general pathology laboratories. In addition to finding single-base changes, NGS allows for the detection of larger structural variants, including insertions/deletions, translocations, and viral insertions. We describe the use of targeted NGS on DNA extracted from formalin-fixed, paraffin-embedded (FFPE) tissue, and show that the short read lengths of NGS are ideally suited to fragmented DNA obtained from FFPE tissue. Further, we describe a novel method for performing hybrid-capture target enrichment using PCR-generated capture probes. As a model, we captured the 5.3-kb Merkel cell polyomavirus (MCPyV) genome in FFPE cases of Merkel cell carcinoma using inexpensive, PCR-derived capture probes, and achieved up to 37,000-fold coverage of the MCPyV genome without prior virus-specific PCR amplification. This depth of coverage made it possible to reproducibly detect viral genome deletions and insertion sites anywhere within the human genome. Out of four cases sequenced, we identified the 5' insertion sites in four of four cases and the 3' sites in three of four cases. These findings demonstrate the potential for an inexpensive gene targeting and NGS method that can be easily adapted for use with FFPE tissue to identify large structural rearrangements, opening up the possibility for further discovery from archival tissue.

SLOPE: a quick and accurate method for locating non-SNP structural variation from targeted next-generation sequence data.

MOTIVATION: Targeted 'deep' sequencing of specific genes or regions is of great interest in clinical cancer diagnostics where some sequence variants, particularly translocations and indels, have known prognostic or diagnostic significance. In this setting, it is unnecessary to sequence an entire genome, and target capture methods can be applied to limit sequencing to important regions, thereby reducing costs and the time required to complete testing. Existing 'next-gen' sequencing analysis packages are optimized for efficiency in whole-genome studies and are unable to benefit from the particular structure of targeted sequence data. RESULTS: We developed SLOPE to detect structural variants from targeted short-DNA reads. We use both real and simulated data to demonstrate SLOPE's ability to rapidly detect insertion/deletion events of various sizes as well as translocations and viral integration sites with high sensitivity and low false discovery rate. AVAILABILITY: Binary code available at http://www-genepi.med.utah.edu/suppl/SLOPE/index.html

contribution of genetically restricted, methicillin-susceptible strains to the ongoing epidemic of community-acquired Staphylococcus aureus infections.

BACKGROUND: Within the current worldwide epidemic of community-acquired Staphylococcus aureus infections, attention has focused on the role of methicillin-resistant strains. We characterize methicillin-susceptible strains that also contribute to this epidemic. METHODS: We tracked cultures from abscess specimens submitted to the microbiology laboratory at St. Louis Children's Hospital and examined Panton-Valentine leukocidin (PVL) genes in methicillin-susceptible S. aureus (MSSA) isolates. We further characterized some isolates by multilocus sequence typing, pulsed-field gel electrophoresis, antibiotic susceptibility, accessory gene regulator (agr) allele, and presence of the arcA gene of the arginine catabolic mobile element. RESULTS: From 1999 to 2007, we detected a 250-fold increase in cultures of abscesses yielding methicillin-resistant S. aureus (MRSA) and a 5-fold increase in abscess cultures yielding MSSA. MSSA isolates from abscesses and wounds were more likely to encode PVL than isolates from other sources. In contrast to PVL-negative isolates of MSSA, which were genetically diverse, PVL-positive isolates were predominantly multilocus sequence typing type 8 and agr type 1. More than half of PVL-positive MSSA isolates were resistant to erythromycin and susceptible to clindamycin with the absence of inducible resistance, a pattern uncommon in PVL-negative MSSA but frequent in the USA300 clone of MRSA. In addition, pulsed-field gel electrophoresis of PVL-positive MSSA strains revealed the USA300 pattern. CONCLUSIONS: In addition to methicillin-resistant strains, the current epidemic of S. aureus infections includes infections caused by methicillin-susceptible strains that are closely related genetically and share phenotypic characteristics other than susceptibility to methicillin. These findings suggest that factors other than methicillin resistance are driving the epidemic.

Panton-Valentine Leukocidin-positive methicillin-resistant Staphylococcus aureus lung infection in patients with cystic fibrosis.

BACKGROUND: Panton-Valentine Leukocidin-expressing (PVL+) methicillin-resistant Staphylococcus aureus (MRSA) is an emerging pathogen worldwide causing fatal necrotizing pneumonias in otherwise healthy individuals but has not been described in patients with cystic fibrosis (CF). Following two cases of patients with CF admitted with lung abscesses in association with PVL+ MRSA, we examined the incidence and the clinical characteristics of MRSA acquisition in our CF patient population. METHODS: Newly acquired MRSA isolates from patients with CF followed up at St. Louis Children's Hospital were analyzed for the presence of Panton-Valentine leukocidin coding region, clindamycin susceptibility, staphylococcal cassette chromosome (SCC) mec type, and multilocus sequence type. Medical records and pulmonary function studies at the time of MRSA isolation were reviewed. RESULTS: MRSA isolates from 40 CF patients were available for analysis. Six children (15%) had PVL+ MRSA infection. All PVL+ organisms were clindamycin susceptible. Patients who acquired a PVL+ organism were more likely to have a focal pulmonary infiltrate on chest radiograph, including cavitary lung lesions in two patients (p = 0.04), a markedly greater decline in FEV1 at the time of MRSA detection (p = 0.01), and a significantly higher WBC count (p = 0.04) and absolute neutrophil count (p = 0.04). These patients were more likely to be admitted for IV antibiotic therapy for respiratory illnesses (p < 0.01). CONCLUSIONS: We describe the emergence of PVL+ MRSA in our CF population in association with development of invasive lung infections including lung abscesses. Early identification and treatment of CF patients with newly acquired PVL+ MRSA may be crucial.

Glutamate regulates neurite outgrowth of cultured descending brain neurons from larval lamprey.

In larval lamprey, descending brain neurons, which regenerate their axons following spinal cord injury, were isolated and examined in cell culture to identify some of the factors that regulate neurite outgrowth. Focal application of 5 mM or 25 mM L-glutamate to single growth cones inhibited outgrowth of the treated neurite, but other neurites from the same neuron were not inhibited, an effect that has not been well studied for neurons in other systems. Glutamate-induced inhibition of neurite outgrowth was abolished by 10 mM kynurenic acid. Application of high potassium media to growth cones inhibited neurite outgrowth, an effect that was blocked by 2 mM cobalt or 100 microM cadmium, suggesting that calcium influx via voltage-gated channels contributes to glutamate-induced regulation of neurite outgrowth. Application of glutamate to growth cones in the presence of 2 microM omega-conotoxin MVIIC (CTX) still inhibited neurite outgrowth, while CTX blocked high potassium-induced inhibition of neurite outgrowth. Thus, CTX blocked virtually all of the calcium influx resulting from depolarization. To our knowledge, this is the first direct demonstration that calcium influx via ligand-gated ion channels can contribute to regulation of neurite outgrowth. Finally, focal application of glutamate to the cell bodies of descending brain neurons inhibited outgrowth of multiple neurites from the same neuron, and this is the first demonstration that multiple neurites can be regulated in this fashion. Signaling mechanisms involving intracellular calcium, similar to those shown here, may be important for regulating axonal regeneration following spinal cord injury in the lamprey.

Expression of a bcr-1 isoform of RARalpha-PML does not affect the penetrance of acute promyelocytic leukemia or the acquisition of an interstitial deletion on mouse chromosome 2.

Expression of a bcr-3 isoform of retinoic acid receptor alpha-promyelocytic leukemia (RARalpha-PML) in mice expressing a bcr-1 isoform of PML-RARalpha is associated with increased penetrance of murine acute promyelocytic leukemia (APL) and the frequent acquisition of an interstitial deletion of one copy of mouse chromosome 2 (del(2)). To determine whether the isoform of RARalpha-PML is important for these effects, we created mice that expressed a bcr-1 isoform of RARalpha-PML. Coexpression with the bcr-1 isoform of PML-RARalpha did not increase the penetrance of APL (7 of 45 animals developed APL with PML-RARalpha alone vs 12 of 44 with both transgenes; P=.19). Furthermore, the frequency of del(2) in APL cells from doubly transgenic mice was not different from that of mice expressing PML-RARalpha alone (3 of 6 vs 6 of 12, respectively-P=1.38-compared with 11 of 11 for mice coexpressing PML-RARalpha and bcr-3 RARalpha-PML). The bcr-1 and bcr-3 isoforms of RARalpha-PML, therefore, have different biological activities that may be relevant for the pathogenesis of murine APL.

Blinded comparison of repetitive-sequence PCR and multilocus sequence typing for genotyping methicillin-resistant Staphylococcus aureus isolates from a children's hospital in St. Louis, Missouri.

We performed a blinded study to compare repetitive-sequence PCR and multilocus sequence typing for genotyping hospital- and community-acquired methicillin-resistant Staphylococcus aureus (MRSA). The MRSA strains that were sequence type 8 (ST8), staphylococcal cassette chromosome mec (SCCmec) type IV, and Panton-Valentine leukocidin-positive clustered separately from those that were ST5 and SCCmec type II.

Identification of genes subject to positive selection in uropathogenic strains of Escherichia coli: a comparative genomics approach.

Escherichia coli is a model laboratory bacterium, a species that is widely distributed in the environment, as well as a mutualist and pathogen in its human hosts. As such, E. coli represents an attractive organism to study how environment impacts microbial genome structure and function. Uropathogenic E. coli (UPEC) must adapt to life in several microbial communities in the human body, and has a complex life cycle in the bladder when it causes acute or recurrent urinary tract infection (UTI). Several studies designed to identify virulence factors have focused on genes that are uniquely represented in UPEC strains, whereas the role of genes that are common to all E. coli has received much less attention. Here we describe the complete 5,065,741-bp genome sequence of a UPEC strain recovered from a patient with an acute bladder infection and compare it with six other finished E. coli genome sequences. We searched 3,470 ortholog sets for genes that are under positive selection only in UPEC strains. Our maximum likelihood-based analysis yielded 29 genes involved in various aspects of cell surface structure, DNA metabolism, nutrient acquisition, and UTI. These results were validated by resequencing a subset of the 29 genes in a panel of 50 urinary, periurethral, and rectal E. coli isolates from patients with UTI. These studies outline a computational approach that may be broadly applicable for studying strain-specific adaptation and pathogenesis in other bacteria.

A Drosophila model of multiple endocrine neoplasia type 2.

Dominant mutations in the Ret receptor tyrosine kinase lead to the familial cancer syndrome multiple endocrine neoplasia type 2 (MEN2). Mammalian tissue culture studies suggest that RetMEN2 mutations significantly alter Ret-signaling properties, but the precise mechanisms by which RetMEN2 promotes tumorigenesis remain poorly understood. To determine the signal transduction pathways required for RetMEN2 activity, we analyzed analogous mutations in the Drosophila Ret ortholog dRet. Overexpressed dRetMEN2 isoforms targeted to the developing retina led to aberrant cell proliferation, inappropriate cell fate specification, and excessive Ras pathway activation. Genetic analysis indicated that dRetMEN2 acts through the Ras-ERK, Src, and Jun kinase pathways. A genetic screen for mutations that dominantly suppress or enhance dRetMEN2 phenotypes identified new genes that are required for the phenotypic outcomes of dRetMEN2 activity. Finally, we identified human orthologs for many of these genes and examined their status in human tumors. Two of these loci showed loss of heterozygosity (LOH) within both sporadic and MEN2-associated pheochromocytomas, suggesting that they may contribute to Ret-dependent oncogenesis.