Recovering Gene Interactions from Single-Cell Data Using Data Diffusion.

Single-cell RNA sequencing technologies suffer from many sources of technical noise, including under-sampling of mRNA molecules, often termed "dropout," which can severely obscure important gene-gene relationships. To address this, we developed MAGIC (Markov affinity-based graph imputation of cells), a method that shares information across similar cells, via data diffusion, to denoise the cell count matrix and fill in missing transcripts. We validate MAGIC on several biological systems and find it effective at recovering gene-gene relationships and additional structures. Applied to the epithilial to mesenchymal transition, MAGIC reveals a phenotypic continuum, with the majority of cells residing in intermediate states that display stem-like signatures, and infers known and previously uncharacterized regulatory interactions, demonstrating that our approach can successfully uncover regulatory relations without perturbations.

Single-Cell Map of Diverse Immune Phenotypes in the Breast Tumor Microenvironment.

Knowledge of immune cell phenotypes in the tumor microenvironment is essential for understanding mechanisms of cancer progression and immunotherapy response. We profiled 45,000 immune cells from eight breast carcinomas, as well as matched normal breast tissue, blood, and lymph nodes, using single-cell RNA-seq. We developed a preprocessing pipeline, SEQC, and a Bayesian clustering and normalization method, Biscuit, to address computational challenges inherent to single-cell data. Despite significant similarity between normal and tumor tissue-resident immune cells, we observed continuous phenotypic expansions specific to the tumor microenvironment. Analysis of paired single-cell RNA and T cell receptor (TCR) sequencing data from 27,000 additional T cells revealed the combinatorial impact of TCR utilization on phenotypic diversity. Our results support a model of continuous activation in T cells and do not comport with the macrophage polarization model in cancer. Our results have important implications for characterizing tumor-infiltrating immune cells.

Molecular Analysis of Arthrobacter Myovirus vB_ArtM-ArV1: We Blame It on the Tail.

This is the first report on a myophage that infects Arthrobacter A novel virus, vB_ArtM-ArV1 (ArV1), was isolated from soil using Arthrobacter sp. strain 68b for phage propagation. Transmission electron microscopy showed its resemblance to members of the family Myoviridae: ArV1 has an isometric head (∼74 nm in diameter) and a contractile, nonflexible tail (∼192 nm). Phylogenetic and comparative sequence analyses, however, revealed that ArV1 has more genes in common with phages from the family Siphoviridae than it does with any myovirus characterized to date. The genome of ArV1 is a linear, circularly permuted, double-stranded DNA molecule (71,200 bp) with a GC content of 61.6%. The genome includes 101 open reading frames (ORFs) yet contains no tRNA genes. More than 50% of ArV1 genes encode unique proteins that either have no reliable identity to database entries or have homologues only in Arthrobacter phages, both sipho- and myoviruses. Using bioinformatics approaches, 13 ArV1 structural genes were identified, including those coding for head, tail, tail fiber, and baseplate proteins. A further 6 ArV1 ORFs were annotated as encoding putative structural proteins based on the results of proteomic analysis. Phylogenetic analysis based on the alignment of four conserved virion proteins revealed that Arthrobacter myophages form a discrete clade that seems to occupy a position somewhat intermediate between myo- and siphoviruses. Thus, the data presented here will help to advance our understanding of genetic diversity and evolution of phages that constitute the order CaudoviralesIMPORTANCE Bacteriophages, which likely originated in the early Precambrian Era, represent the most numerous population on the planet. Approximately 95% of known phages are tailed viruses that comprise three families: Podoviridae (with short tails), Siphoviridae (with long noncontractile tails), and Myoviridae (with contractile tails). Based on the current hypothesis, myophages, which may have evolved from siphophages, are thought to have first emerged among Gram-negative bacteria, whereas they emerged only later among Gram-positive bacteria. The results of the molecular characterization of myophage vB_ArtM-ArV1 presented here conform to the aforementioned hypothesis, since, at a glance, bacteriophage vB_ArtM-ArV1 appears to be a siphovirus that possesses a seemingly functional contractile tail. Our work demonstrates that such "chimeric" myophages are of cosmopolitan nature and are likely characteristic of the ecologically important soil bacterial genus Arthrobacter.

Generation in yeast of recombinant virus-like particles of porcine circovirus type 2 capsid protein and their use for a serologic assay and development of monoclonal antibodies.

BACKGROUND: Porcine circovirus type 2 (PCV2) is considered to be an important emerging pathogen associated with a number of different syndromes and diseases in pigs known as PCV2-associated diseases. It has been responsible for significant mortality among pigs and remains a serious economic problem to the swine industry worldwide leading to significant negative impacts on profitability of pork production. RESULTS: In this study we have demonstrated that PCV2 capsid (Cap) protein based virus-like particles (VLPs) were efficiently produced in yeast S. cerevisiae and induced production of monoclonal antibodies (MAbs) reactive with virus-infected cells. Moreover, PCV2 Cap VLPs served as a highly specific recombinant antigen for the development of an indirect IgG PCV2 Cap VLP-based ELISA for the detection of virus-specific IgG antibodies in swine sera. Four hundred-nine serum samples collected from pigs in Lithuania were tested for PCV2-specific IgG to determine the sensitivity and specificity of the newly developed ELISA in parallel using a commercial SERELISA test as a gold standard. From 409 tested serum samples, 297 samples were positive by both assays. Thirty-nine sera from 112 serum samples were determined as negative by SERELISA but were found to be positive both in the newly developed indirect IgG PCV2 Cap VLP-based ELISA and the PCR test. CONCLUSIONS: We have demonstrated that S. cerevisiae expression system is an alternative to insect/baculovirus expression system for production of homogenous in size and shape PCV2 Cap protein-based VLPs similar to native virions. Yeast expression system tolerated native virus genes encoding PCV2 Cap protein variants as well as the codon-optimized gene. Moreover, yeast-derived PCV2 Cap VLPs were capable to induce the generation of PCV2-specific MAbs that did not show any cross-reactivity with PCV1-infected cells. The high sensitivity and specificity of the indirect IgG PCV2 Cap VLP-based ELISA clearly suggested that this assay is potentially useful diagnostic tool for screening PCV2-suspected samples.

Single-cell transcriptional profiling of clear cell renal cell carcinoma reveals a tumor-associated endothelial tip cell phenotype.

Clear cell renal cell carcinoma (ccRCC) is the most prevalent form of renal cancer, accounting for over 75% of cases. The asymptomatic nature of the disease contributes to late-stage diagnoses and poor survival. Highly vascularized and immune infiltrated microenvironment are prominent features of ccRCC, yet the interplay between vasculature and immune cells, disease progression and response to therapy remains poorly understood. Using droplet-based single-cell RNA sequencing we profile 50,236 transcriptomes from paired tumor and healthy adjacent kidney tissues. Our analysis reveals significant heterogeneity and inter-patient variability of the tumor microenvironment. Notably, we discover a previously uncharacterized vasculature subpopulation associated with epithelial-mesenchymal transition. The cell-cell communication analysis reveals multiple modes of immunosuppressive interactions within the tumor microenvironment, including clinically relevant interactions between tumor vasculature and stromal cells with immune cells. The upregulation of the genes involved in these interactions is associated with worse survival in the TCGA KIRC cohort. Our findings demonstrate the role of tumor vasculature and stromal cell populations in shaping the ccRCC microenvironment and uncover a subpopulation of cells within the tumor vasculature that is associated with an angiogenic phenotype.

Multi-omics at single-cell resolution: comparison of experimental and data fusion approaches.

Biological samples are inherently heterogeneous and complex. Tackling this complexity requires innovative technological and analytical solutions. Recent advances in high-throughput single-cell isolation and nucleic acid barcoding methods are rapidly changing the technological landscape of biological sciences and now make it possible to measure the (epi)genomic, transcriptomic, or proteomic state of individual cells. In addition, few experimental approaches enable multi-omics measurements of the same cell. However, merging-omics data collected from different experiments remains a considerable challenge. Although several strategies for merging transcriptomics datasets have recently been introduced, cell-to-cell variability and heterogeneity remains one of the confounding factors limiting data fusion and integration. Here, we focus our discussion on the latest single-cell technological and analytical solutions to achieve high data dimensionality and resolution. Obtaining datasets with a wealth of multi-omics information will undoubtedly provide new avenues for researchers to unravel the complexity of biological samples encountered in modern biological research and molecular diagnostics.

Single-cell barcoding and sequencing using droplet microfluidics.

Single-cell RNA sequencing has recently emerged as a powerful tool for mapping cellular heterogeneity in diseased and healthy tissues, yet high-throughput methods are needed for capturing the unbiased diversity of cells. Droplet microfluidics is among the most promising candidates for capturing and processing thousands of individual cells for whole-transcriptome or genomic analysis in a massively parallel manner with minimal reagent use. We recently established a method called inDrops, which has the capability to index >15,000 cells in an hour. A suspension of cells is first encapsulated into nanoliter droplets with hydrogel beads (HBs) bearing barcoding DNA primers. Cells are then lysed and mRNA is barcoded (indexed) by a reverse transcription (RT) reaction. Here we provide details for (i) establishing an inDrops platform (1 d); (ii) performing hydrogel bead synthesis (4 d); (iii) encapsulating and barcoding cells (1 d); and (iv) RNA-seq library preparation (2 d). inDrops is a robust and scalable platform, and it is unique in its ability to capture and profile >75% of cells in even very small samples, on a scale of thousands or tens of thousands of cells.

Identification of Two Novel Members of the Tentative Genus Wukipolyomavirus in Wild Rodents.

Two novel polyomaviruses (PyVs) were identified in kidney and chest-cavity fluid samples of wild bank voles (Myodes glareolus) and common voles (Microtus arvalis) collected in Germany. All cloned and sequenced genomes had the typical PyV genome organization, including putative open reading frames for early regulatory proteins large T antigen and small T antigen on one strand and for structural late proteins (VP1, VP2 and VP3) on the other strand. Virus-like particles (VLPs) were generated by yeast expression of the VP1 protein of both PyVs. VLP-based ELISA and large T-antigen sequence-targeted polymerase-chain reaction investigations demonstrated signs of infection of these novel PyVs in about 42% of bank voles and 18% of common voles. In most cases only viral DNA, but not VP1-specific antibodies were detected. In additional animals exclusively VP1-specific antibodies, but no viral DNA was detected, indicative for virus clearance. Phylogenetic and clustering analysis including all known PyV genomes placed novel bank vole and common vole PyVs amongst members of the tentative Wukipolymavirus genus. The other known four rodent PyVs, Murine PyV and Hamster PyV, and Murine pneumotropic virus and Mastomys PyV belong to different phylogenetic clades, tentatively named Orthopolyomavirus I and Orthopolyomavirus II, respectively. In conclusion, the finding of novel vole-borne PyVs may suggest an evolutionary origin of ancient wukipolyomaviruses in rodents and may offer the possibility to develop a vole-based animal model for human wukipolyomaviruses.

Purification of recombinant virus-like particles of porcine circovirus type 2 capsid protein using ion-exchange monolith chromatography.

Diseases associated with porcine circovirus type 2 (PCV2) infection are having a severe economic impact on swine-producing countries. The PCV2 capsid (Cap) protein expressed in eukaryotic systems self-assemble into virus-like particles (VLPs) which can serve as antigens for diagnostics or/and as vaccine candidates. In this work, conventional adsorbents as well as a monolithic support with large pore sizes were examined for the chromatographic purification of PCV2 Cap VLPs from clarified yeast lysate. Q Sepharose XL was used for the initial separation of VLPs from residual host nucleic acids and some host cell proteins. For the further purification of PCV2 Cap VLPs, SP Sepharose XL, Heparin Sepharose CL-6B and CIMmultus SO3 monolith were tested. VLPs were not retained on SP Sepharose XL. The purity of VLPs after chromatography on Heparin Sepharose CL-6B was only 4-7% and the recovery of VLPs was 5-7%. Using ion-exchange chromatography on the CIMmultus SO3 monolith, PCV2 Cap VLPs with the purity of about 40% were obtained. The recovery of VLPs after chromatography on the CIMmultus SO3 monolith was 15-18%. The self-assembly of purified PCV2 Cap protein into VLPs was confirmed by electron microscopy. Two-step chromatographic purification procedure of PCV2 Cap VLPs from yeast lysate was developed using Q Sepharose XL and cation-exchange CIMmultus SO3 monolith.

Isolation and characterization of vB_ArS-ArV2 - first Arthrobacter sp. infecting bacteriophage with completely sequenced genome.

This is the first report on a complete genome sequence and biological characterization of the phage that infects Arthrobacter. A novel virus vB_ArS-ArV2 (ArV2) was isolated from soil using Arthrobacter sp. 68b strain for phage propagation. Based on transmission electron microscopy, ArV2 belongs to the family Siphoviridae and has an isometric head (∼63 nm in diameter) with a non-contractile flexible tail (∼194×10 nm) and six short tail fibers. ArV2 possesses a linear, double-stranded DNA genome (37,372 bp) with a G+C content of 62.73%. The genome contains 68 ORFs yet encodes no tRNA genes. A total of 28 ArV2 ORFs have no known functions and lack any reliable database matches. Proteomic analysis led to the experimental identification of 14 virion proteins, including 9 that were predicted by bioinformatics approaches. Comparative phylogenetic analysis, based on the amino acid sequence alignment of conserved proteins, set ArV2 apart from other siphoviruses. The data presented here will help to advance our understanding of Arthrobacter phage population and will extend our knowledge about the interaction between this particular host and its phages.