An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

BACKGROUND: Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. RESULTS: Here, we generate single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We develop an unsupervised machine learning approach ("automatic consensus nonnegative matrix factorization" (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirm a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly, however, this weak-mesenchymal-like program is maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 h, suggesting an uncharacterized therapy-escape mechanism. CONCLUSIONS: Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Isolation and Characterization of a Novel Alpha-Hemolytic Streptococcus spp. from the Oral Cavity and Blood of Septicemic Periparturient Immunodeficient Mice.

MISTRG is an immunodeficient mouse strain that expresses multiple human cytokines that support hematopoietic stem cell maintenance and myelopoiesis. While establishing a breeding colony of MISTRG mice in a dedicated barrier room, 6 cases of death or disease occurred in pregnant or postpartum mice. Clinically, this manifested as hunched posture, dyspnea, and 1 case of emaciation with ataxia. Pathologic analysis of 7 mice revealed multisystemic necrosuppurative inflammation variably affecting the uterus and placenta, joints, meninges, inner and middle ears, kidneys, and small intestine. Bacteria cultured from the blood of septic mice were identified with 89% probability by the Vitek 2 identification system as Streptococcus sanguinus with atypical biochemical parameters; the API 20E/NE system fully differentiated the isolates as a novel Streptococcus species. MALDI Biotyper-based mass spectrometry also indicated that the phenotype represented a novel Streptococcus spp. Sequencing revealed that the full-length 16S rRNA gene identity was below 97% with known Streptococcus species, including the 2 closest species Streptococcus acidominimus and Streptococcus azizii. We propose the name Streptococcus murisepticum spp. nov to our novel isolates. All male mice in this colony remained healthy despite their association with diseased female mice. Overall, 19% of the colony carried the novel Streptococcus in their oral cavity, but it could not be detected in feces. The organism was sensitive to amoxicillin, which was administered via drinking water throughout pregnancy and weaning to establish a colony of pathogen-negative future breeders. The colony remained disease-free and culture-negative for Streptococcus murisepticum spp. nov after treatment with amoxicillin. We suspect that oral colonization of MISTRG mice with the novel Streptococcus species and its associated unique pathology in periparturient mice is potentially the principal cause of loss of this strain at several institutions. Therefore, screening the oral cavity for α-hemolytic streptococci followed by targeted antibiotic treatment may be necessary when establishing MISTRG and allied immunodeficient mouse strains.

Measuring Influenza Virus Infection Using Bioluminescent Reporter Viruses for In Vivo Imaging and In Vitro Replication Assays.

To streamline standard virological assays, we developed bioluminescent replication-competent influenza reporter viruses that mimic their parental counterparts. These reporter viruses provide a rapid and quantitative readout of viral infection and replication. Moreover, they permit real-time in vivo measures of viral load, tissue distribution, and transmission in the same cohort of animals over the entire course of infection-measurements that were not previously possible. Here we provide detailed protocols using bioluminescent reporter viruses for in vivo imaging in mice and ferrets. We also describe cell culture-based techniques using reporter viruses for quantification of viral titers and performing microneutralization assays. The ease, speed, and adaptability of these approaches have the potential to accelerate multiple areas of influenza virus research.

Visualizing real-time influenza virus infection, transmission and protection in ferrets.

Influenza transmission efficiency in ferrets is vital for risk-assessment studies. However, the inability to monitor viral infection and transmission dynamics in real time only provides a glimpse into transmissibility. Here we exploit a replication-competent influenza reporter virus to investigate dynamics of infection/transmission in ferrets. Bioluminescent imaging of ferrets infected with A/California/04/2009 H1N1 virus (CA/09) encoding NanoLuc (NLuc) luciferase provides the first real-time snapshot of influenza infection/transmission. Luminescence in the respiratory tract and in less well-characterized extra-pulmonary sites is observed, and imaging identifies infections in animals that would have otherwise been missed by traditional methods. Finally, the reporter virus significantly increases the speed and sensitivity of virological and serological assays. Thus, bioluminescent imaging of influenza infections rapidly determines intra-host dissemination, inter-host transmission and viral load, revealing infection dynamics and pandemic potential of the virus. These results have important implications for antiviral drug susceptibility, vaccine efficacy, transmissibility and pathogenicity studies.

Protective roles for caspase-8 and cFLIP in adult homeostasis.

Caspase-8 or cellular FLICE-like inhibitor protein (cFLIP) deficiency leads to embryonic lethality in mice due to defects in endothelial tissues. Caspase-8(-/-) and receptor-interacting protein kinase-3 (RIPK3)(-/-), but not cFLIP(-/-) and RIPK3(-/-), double-knockout animals develop normally, indicating that caspase-8 antagonizes the lethal effects of RIPK3 during development. Here, we show that the acute deletion of caspase-8 in the gut of adult mice induces enterocyte death, disruption of tissue homeostasis, and inflammation, resulting in sepsis and mortality. Likewise, acute deletion of caspase-8 in a focal region of the skin induces local keratinocyte death, tissue disruption, and inflammation. Strikingly, RIPK3 ablation rescues both phenotypes. However, acute loss of cFLIP in the skin produces a similar phenotype that is not rescued by RIPK3 ablation. TNF neutralization protects from either acute loss of caspase-8 or cFLIP. These results demonstrate that caspase-8-mediated suppression of RIPK3-induced death is required not only during development but also for adult homeostasis. Furthermore, RIPK3-dependent inflammation is dispensable for the skin phenotype.

Sub-Doppler infrared spectroscopy of CH2D radical in a slit supersonic jet: isotopic symmetry breaking in the CH stretching manifold.

First high-resolution infrared absorption spectra in the fundamental symmetric/asymmetric CH stretching region of isotopically substituted methyl radical, CH(2)D, are reported and analyzed. These studies become feasible in the difference frequency spectrometer due to (i) high density radical generation via dissociative electron attachment to CH(2)DI in a discharge, (ii) low rotational temperatures (23 K) from supersonic cooling in a slit expansion, (iii) long absorption path length (64 cm) along the slit axes, and (iv) near shot noise limited absorption sensitivity (5 × 10(-7)/√(Hz)). The spectra are fully rovibrationally resolved and fit to an asymmetric top rotational Hamiltonian to yield rotational/centrifugal constants and vibrational band origins. In addition, the slit expansion collisionally quenches the transverse velocity distribution along the laser probe direction, yielding sub-Doppler resolution of spin-rotation structure and even partial resolution of nuclear hyperfine structure for each rovibrational line. Global least-squares fits to the line shapes provide additional information on spin-rotation and nuclear hyperfine constants, which complement and clarify previous FTIR studies [K. Kawaguchi, Can. J. Phys. 79, 449 (2001)] of CH(2)D in the out-of-plane bending region. Finally, analysis of the spectral data from the full isotopomeric CH(m)D(3-m) series based on harmonically coupled Morse oscillators establishes a predictive framework for describing the manifold of planar stretching vibrations in this fundamental combustion radical.

Quantum deconstruction of the infrared spectrum of CH5+.

We present two quantum calculations of the infrared spectrum of protonated methane (CH5+) using full-dimensional, ab initio-based potential energy and dipole moment surfaces. The calculated spectra compare well with a low-resolution experimental spectrum except below 1000 cm(-1), where the experimental spectrum shows no absorption. The present calculations find substantial absorption features below 1000 cm(-1), in qualitative agreement with earlier classical calculations of the spectrum. The major spectral bands are analyzed in terms of the molecular motions. Of particular interest is an intense feature at 200 cm(-1), which is due to an isomerization mode that connects two equivalent minima. Very recent high-resolution jet-cooled spectra in the CH stretch region (2825 to 3050 cm(-1)) are also reported, and assignments of the band origins are made, based on the present quantum calculations.