The Poisson distribution model fits UMI-based single-cell RNA-sequencing data.

BACKGROUND: Modeling of single cell RNA-sequencing (scRNA-seq) data remains challenging due to a high percentage of zeros and data heterogeneity, so improved modeling has strong potential to benefit many downstream data analyses. The existing zero-inflated or over-dispersed models are based on aggregations at either the gene or the cell level. However, they typically lose accuracy due to a too crude aggregation at those two levels. RESULTS: We avoid the crude approximations entailed by such aggregation through proposing an independent Poisson distribution (IPD) particularly at each individual entry in the scRNA-seq data matrix. This approach naturally and intuitively models the large number of zeros as matrix entries with a very small Poisson parameter. The critical challenge of cell clustering is approached via a novel data representation as Departures from a simple homogeneous IPD (DIPD) to capture the per-gene-per-cell intrinsic heterogeneity generated by cell clusters. Our experiments using real data and crafted experiments show that using DIPD as a data representation for scRNA-seq data can uncover novel cell subtypes that are missed or can only be found by careful parameter tuning using conventional methods. CONCLUSIONS: This new method has multiple advantages, including (1) no need for prior feature selection or manual optimization of hyperparameters; (2) flexibility to combine with and improve upon other methods, such as Seurat. Another novel contribution is the use of crafted experiments as part of the validation of our newly developed DIPD-based clustering pipeline. This new clustering pipeline is implemented in the R (CRAN) package scpoisson.

Interleukin-1 Receptor-Associated Kinase (IRAK) Signaling in Kaposi Sarcoma-Associated Herpesvirus-Induced Primary Effusion Lymphoma.

Kaposi sarcoma-associated herpesvirus (KSHV) is necessary but not sufficient for primary effusion lymphoma (PEL) development. Alterations in cellular signaling pathways are also a characteristic of PEL. Other B cell lymphomas have acquired an oncogenic mutation in the myeloid differentiation primary response 88 (MYD88) gene. The MYD88 L265P mutant results in the activation of interleukin-1 receptor associated kinase (IRAK). To probe IRAK/MYD88 signaling in PEL, we employed CRISPR/Cas9 technology to generate stable deletion clones in BCBL-1Cas9 and BC-1Cas9 cells. To look for off-target effects, we determined the complete exome of the BCBL-1Cas9 and BC-1Cas9 cells. Deletion of either MYD88, IRAK4, or IRAK1 abolished interleukin-1 beta (IL-1ß) signaling; however, we were able to grow stable subclones from each population. Transcriptome sequencing (RNA-seq) analysis of IRAK4 knockout cell lines (IRAK4 KOs) showed that the IRAK pathway induced cellular signals constitutively, independent of IL-1ß stimulation, which was abrogated by deletion of IRAK4. Transient complementation with IRAK1 increased NF-κB activity in MYD88 KO, IRAK1 KO, and IRAK4 KO cells even in the absence of IL-1ß. IL-10, a hallmark of PEL, was dependent on the IRAK pathway, as IRAK4 KOs showed reduced IL-10 levels. We surmise that, unlike B cell receptor (BCR) signaling, MYD88/IRAK signaling is constitutively active in PEL, but that under cell culture conditions, PEL rapidly became independent of this pathway.IMPORTANCE One hundred percent of primary effusion lymphoma (PEL) cases are associated with Kaposi sarcoma-associated herpesvirus (KSHV). PEL cell lines, such as BCBL-1, are the workhorse for understanding this human oncovirus and the host pathways that KSHV dysregulates. Understanding their function is important for developing new therapies as well as identifying high-risk patient groups. The myeloid differentiation primary response 88 (MYD88)/interleukin-1 receptor associated kinase (IRAK) pathway, which has progrowth functions in other B cell lymphomas, has not been fully explored in PEL. By performing CRISPR/Cas9 knockout (KO) studies targeting the IRAK pathway in PEL, we were able to determine that established PEL cell lines can circumvent the loss of IRAK1, IRAK4, and MYD88; however, the deletion clones are deficient in interleukin-10 (IL-10) production. Since IL-10 suppresses T cell function, this suggests that the IRAK pathway may serve a function in vivo and during early-stage development of PEL.

Extracellular vesicles from Kaposi Sarcoma-associated herpesvirus lymphoma induce long-term endothelial cell reprogramming.

Extracellular signaling is a mechanism that higher eukaryotes have evolved to facilitate organismal homeostasis. Recent years have seen an emerging interest in the role of secreted microvesicles, termed extracellular vesicles (EV) or exosomes in this signaling network. EV contents can be modified by the cell in response to stimuli, allowing them to relay information to neighboring cells, influencing their physiology. Here we show that the tumor virus Kaposi's Sarcoma-associated herpesvirus (KSHV) hijacks this signaling pathway to induce cell proliferation, migration, and transcriptome reprogramming in cells not infected with the virus. KSHV-EV activates the canonical MEK/ERK pathway, while not alerting innate immune regulators, allowing the virus to exert these changes without cellular pathogen recognition. Collectively, we propose that KSHV establishes a niche favorable for viral spread and cell transformation through cell-derived vesicles, all while avoiding detection.

Epstein-Barr virus enhances genome maintenance of Kaposi sarcoma-associated herpesvirus.

Primary effusion lymphoma (PEL) is a B cell lymphoma that is always associated with Kaposi's sarcoma-associated herpesvirus (KSHV) and in many cases also with Epstein-Barr virus (EBV); however, the requirement for EBV coinfection is not clear. Here, we demonstrate that adding exogenous EBV to KSHV+ single-positive PEL leads to increased KSHV genome maintenance and KSHV latency-associated nuclear antigen (LANA) expression. To show that EBV was necessary for naturally coinfected PEL, we nucleofected KSHV+/EBV+ PEL cell lines with an EBV-specific CRISPR/Cas9 plasmid to delete EBV and observed a dramatic decrease in cell viability, KSHV genome copy number, and LANA expression. This phenotype was reversed by expressing Epstein-Barr nuclear antigen 1 (EBNA-1) in trans, even though EBNA-1 and LANA do not colocalize in infected cells. This work reveals that EBV EBNA-1 plays an essential role in the pathogenesis of PEL by increasing KSHV viral load and LANA expression.

CD81 fusion alters SARS-CoV-2 Spike trafficking.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic caused the biggest public health crises in recent history. Many expect future coronavirus introductions into the human population. Hence, it is essential to understand the basic biology of these viruses. In natural infection, the SARS-CoV-2 Spike (S) glycoprotein is co-expressed with all other viral proteins, which modify cellular compartments to maximize virion assembly. By comparison, most of S is degraded when the protein is expressed in isolation, as in current molecular vaccines. To probe the maturation pathway of S, we redirected its maturation by fusing S to the tetraspanin protein CD81. CD81 is a defining constituent of extracellular vesicles (EVs) or exosomes. EVs are generated in large numbers by all cells, extruded into blood and lymph, and transfer cargo between cells and systemically (estimated 1012 EVs per mL plasma). EVs, like platelets, can be transfused between unrelated donors. When fusing the proline-stabilized form of strain Delta S into the flexible, large extracellular loop of CD81 rather than being degraded in the lysosome, S was extruded into EVs. CD81-S fusion containing EVs were produced in large numbers and could be isolated to high purity. Purified CD81::S EVs bound ACE2, and S displayed on individual EV was observed by cryogenic electron microscopy (EM). The CD81::S-fusion EVs were non-toxic and elicited an anti-S trimer and anti-RBD antibody response in mice. This report shows a design path to maximize viral glycoprotein assembly and release without relying on the co-expression of potentially pathogenic nonstructural viral proteins. IMPORTANCE: The severe acute respiratory syndrome coronavirus 2 pandemic caused the biggest public health crises in recent history. To understand the maturation pathway of S, we fused S to the tetraspanin protein CD81. The resulting molecule is secreted in extracellular vesicles and induces antibodies in mice. This may be a general design path for viral glycoprotein vaccines.

The Poisson distribution model fits UMI-based single-cell RNA-sequencing data.

Background: Modeling of single cell RNA-sequencing (scRNA-seq) data remains challenging due to a high percentage of zeros and data heterogeneity, so improved modeling has strong potential to benefit many downstream data analyses. The existing zero-inflated or over-dispersed models are based on aggregations at either the gene or the cell level. However, they typically lose accuracy due to a too crude aggregation at those two levels. Results: We avoid the crude approximations entailed by such aggregation through proposing an Independent Poisson Distribution (IPD) particularly at each individual entry in the scRNA-seq data matrix. This approach naturally and intuitively models the large number of zeros as matrix entries with a very small Poisson parameter. The critical challenge of cell clustering is approached via a novel data representation as Departures from a simple homogeneous IPD (DIPD) to capture the per-gene-per-cell intrinsic heterogeneity generated by cell clusters. Our experiments using real data and crafted experiments show that using DIPD as a data representation for scRNA-seq data can uncover novel cell subtypes that are missed or can only be found by careful parameter tuning using conventional methods. Conclusions: This new method has multiple advantages, including (1) no needfor prior feature selection or manual optimization of hyperparameters; (2) flexibility to combine with and improve upon other methods, such as Seurat. Another novel contribution is the use of crafted experiments as part of the validation of our newly developed DIPD-based clustering pipeline. This new clustering pipeline is implemented in the R (CRAN) package scpoisson .

Large-scale heparin-based bind-and-elute chromatography identifies two biologically distinct populations of extracellular vesicles.

Purifying extracellular vesicles (EVs) has been challenging because EVs are heterogeneous in cargo yet share similar sizes and densities. Most surface marker-based affinity separation methods are limited to research or diagnostic scales. We report that heparin chromatography can separate purified EVs into two distinct subpopulations as ascertained by MS/MS: a non-heparin-binding (NHB) fraction that contains classical EV markers such as tetraspanins and a heparin-binding (HB) fraction enriched in fibronectins and histones. Both fractions were similarly fusogenic but induced different transcriptional responses in endothelial cells. While EVs that were purified by conventional, non-affinity methods alone induced ERK1/2 phosphorylation and Ki67, the NHB fraction did not. This result suggests heparin chromatography as an additional novel fractionation step that is inherently scalable, does not lead to loss of material, and separates inflammatory and pyrogenic EVs from unreactive EVs, which will improve clinical applications.

Intra-Host Evolution Provides for the Continuous Emergence of SARS-CoV-2 Variants.

Variants of concern (VOC) in SARS-CoV-2 refer to viruses whose viral genomes differ from the ancestor virus by ≥3 single-nucleotide variants (SNVs) and that show the potential for higher transmissibility and/or worse clinical progression. VOC have the potential to disrupt ongoing public health measures and vaccine efforts. Still, too little is known regarding how frequently new viral variants emerge and under what circumstances. We report a study to determine the degree of SARS-CoV-2 sequence evolution in 94 patients and to estimate the frequency at which highly diverse variants emerge. Two cases accumulated ≥9 SNVs over a 2-week period and one case accumulated 23 SNVs over 3 weeks, including three nonsynonymous mutations in the spike protein (D138H, E554D, D614G). The remainder of the infected patients did not show signs of intra-host evolution. We estimate that in as much as 2% of hospitalized COVID-19 cases, variants with multiple mutations in the spike glycoprotein emerge in as little as 1 month of persistent intra-host virus replication. This suggests the continued local emergence of variants with multiple nonsynonymous SNVs, even in patients without overt immune deficiency. Surveillance by sequencing for (i) viremic COVID-19 patients, (ii) patients suspected of reinfection, and (iii) patients with diminished immune function may offer broad public health benefits. IMPORTANCE New SARS-CoV-2 variants can potentially disrupt ongoing public health measures and vaccine efforts. Still, little is known regarding how frequently new viral variants emerge and under what circumstances. Based on this study, we estimate that in hospitalized COVID-19 cases, variants with multiple mutations may emerge locally in as little as 1 month, even in patients without overt immune deficiency. Surveillance by sequencing for continuously shedding patients, patients suspected of reinfection, and patients with diminished immune function may offer broad public health benefits.

Evidence for Multiple Subpopulations of Herpesvirus-Latently Infected Cells.

Kaposi's sarcoma-associated herpesvirus (KSHV)-associated primary effusion lymphomas (PEL) are traditionally viewed as homogenous regarding viral transcription and lineage of origin, but so far this contention has not been explored at the single-cell level. Single-cell RNA sequencing of latently infected PEL supports the existence of multiple subpopulations even within a single cell line. At most 1% of the cells showed evidence of near-complete lytic transcription. The majority of cells only expressed the canonical viral latent transcripts: those originating from the latency locus, the viral interferon regulatory factor locus, and the viral lncRNA nut-1/Pan/T1.1; however, a significant fraction of cells showed various degrees of more permissive transcription, and some showed no evidence of KSHV transcripts whatsoever. Levels of viral interleukin-6 (IL-6)/K2 mRNA emerged as the most distinguishing feature to subset KSHV-infected PEL. One newly uncovered phenotype is the existence of BCBL-1 cells that readily adhered to fibronectin and that displayed mesenchymal lineage-like characteristics. IMPORTANCE Latency is the defining characteristic of the Herpesviridae and central to the tumorigenesis phenotype of Kaposi's sarcoma-associated herpesvirus (KSHV). KSHV-driven primary effusion lymphomas (PEL) rapidly develop resistance to therapy, suggesting tumor instability and plasticity. At any given time, a fraction of PEL cells spontaneously reactivate KSHV, suggesting transcriptional heterogeneity even within a clonal cell line under optimal growth conditions. This study employed single-cell mRNA sequencing to explore the within-population variability of KSHV transcription and how it relates to host cell transcription. Individual clonal PEL cells exhibited differing patterns of viral transcription. Most cells showed the canonical pattern of KSHV latency (LANA, vCyc, vFLIP, Kaposin, and vIRFs), but a significant fraction evidenced extended viral gene transcription, including of the viral IL-6 homolog, open reading frame K2. This study suggests new targets of intervention for PEL. It establishes a conceptual framework to design KSHV cure studies analogous to those for HIV.

Whole-genome sequencing of Kaposi sarcoma-associated herpesvirus (KSHV/HHV8) reveals evidence for two African lineages.

Kaposi sarcoma (KS)-associated herpesvirus (KSHV/HHV-8) was first sequenced from the body cavity (BC) lymphoma cell line, BC-1, in 1996. Few other KSHV genomes have been reported. Our knowledge of sequence variation for this virus remains spotty. This study reports additional genomes from historical US patient samples and from African KS biopsies. It describes an assay that spans regions of the virus that cannot be covered by short read sequencing. These include the terminal repeats, the LANA repeats, and the origins of replication. A phylogenetic analysis, based on 107 genomes, identified three distinct clades; one containing isolates from USA/Europe/Japan collected in the 1990s and two of Sub-Saharan Africa isolates collected since 2010. This analysis indicates that the KSHV strains circulating today differ from the isolates collected at the height of the AIDS epidemic. This analysis helps experimental designs and potential vaccine studies.

Imaging of surface microdomains on individual extracellular vesicles in 3-D.

Extracellular vesicles (EVs) are secreted from all cell types and are intimately involved in tissue homeostasis. They are being explored as vaccine and gene therapy platforms, as well as potential biomarkers. As their size is below the diffraction limit of light microscopy, direct visualizations have been daunting and single-particle studies under physiological conditions have been hampered. Here, direct stochastic optical reconstruction microscopy (dSTORM) was employed to visualize EVs in three-dimensions and to localize molecule clusters such as the tetraspanins CD81 and CD9 on the surface of individual EVs. These studies demonstrate the existence of membrane microdomains on EVs. These were confirmed by Cryo-EM. Individual particle visualization provided insights into the heterogeneity, structure, and complexity of EVs not previously appreciated.

High-Density Amplicon Sequencing Identifies Community Spread and Ongoing Evolution of SARS-CoV-2 in the Southern United States.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is constantly evolving. Prior studies focused on high-case-density locations, such as the northern and western metropolitan areas of the United States. This study demonstrates continued SARS-CoV-2 evolution in a suburban southern region of the United States by high-density amplicon sequencing of symptomatic cases. 57% of strains carry the spike D614G variant, which is associated with higher genome copy numbers, and its prevalence expands with time. Four strains carry a deletion in a predicted stem loop of the 3' UTR. The data are consistent with community spread within local populations and the larger continental United States. The data instill confidence in current testing sensitivity and validate "testing by sequencing" as an option to uncover cases, particularly nonstandard coronavirus disease 2019 (COVID-19) clinical presentations. This study contributes to the understanding of COVID-19 through an extensive set of genomes from a non-urban setting and informs vaccine design by defining D614G as a dominant and emergent SARS-CoV-2 isolate in the United States.

Large-scale, cross-flow based isolation of highly pure and endocytosis-competent extracellular vesicles.

Isolation of extracellular vesicles (EVs) from cell culture supernatant or plasma can be accomplished in a variety of ways. Common measures to quantify relative success are: concentration of the EVs, purity from non-EVs associated protein, size homogeneity and functionality of the final product. Here, we present an industrial-scale workflow for isolating highly pure and functional EVs using cross-flow based filtration coupled with high-molecular weight Capto Core size exclusion. Through this combination, EVs loss is kept to a minimum. It outperforms other isolation procedures based on a number of biochemical and biophysical assays. Moreover, EVs isolated through this method can be further concentrated down or directly immunopurified to obtain discreet populations of EVs. From our results, we propose that cross-flow/Capto Core isolation is a robust method of purifying highly concentrated, homogenous, and functionally active EVs from industrial-scale input volumes with few contaminants relative to other methods.