Urine Single-Cell RNA Sequencing in Focal Segmental Glomerulosclerosis Reveals Inflammatory Signatures.

INTRODUCTION: Individuals with focal segmental glomerular sclerosis (FSGS) typically undergo kidney biopsy only once, which limits the ability to characterize kidney cell gene expression over time. METHODS: We used single-cell RNA sequencing (scRNA-seq) to explore disease-related molecular signatures in urine cells from subjects with FSGS. We collected 17 urine samples from 12 FSGS subjects and captured these as 23 urine cell samples. The inflammatory signatures from renal epithelial and immune cells were evaluated in bulk gene expression data sets of FSGS and minimal change disease (MCD) (The Nephrotic Syndrome Study Network [NEPTUNE] study) and an immune single-cell data set from lupus nephritis (Accelerating Medicines Partnership). RESULTS: We identified immune cells, predominantly monocytes, and renal epithelial cells in the urine. Further analysis revealed 2 monocyte subtypes consistent with M1 and M2 monocytes. Shed podocytes in the urine had high expression of marker genes for epithelial-to-mesenchymal transition (EMT). We selected the 16 most highly expressed genes from urine immune cells and 10 most highly expressed EMT genes from urine podocytes as immune signatures and EMT signatures, respectively. Using kidney biopsy transcriptomic data from NEPTUNE, we found that urine cell immune signature and EMT signature genes were more highly expressed in FSGS biopsies compared with MCD biopsies. CONCLUSION: The identification of monocyte subsets and podocyte expression signatures in the urine samples of subjects with FSGS suggests that urine cell profiling might serve as a diagnostic and prognostic tool in nephrotic syndrome. Furthermore, this approach may aid in the development of novel biomarkers and identifying personalized therapies targeting particular molecular pathways in immune cells and podocytes.

ATAC Sequencing Protocol For Cryopreserved Mammalian Cells.

ATAC-seq (assay for transposase-accessible chromatin with high-throughput sequencing) is a powerful method to evaluate chromatin accessibility and nucleosome positioning at a genome-wide scale. This assay uses a hyperactive Tn5 transposase, to simultaneously cut open chromatin and insert adapter sequences. After sequencing, the reads generated through this technique are generally indicative of transcriptional regulatory elements that are located in accessible chromatin. This method was originally developed by Buenrostro et al. (2013), and since then it has been improved by the same authors several times, until their last update called OMNI ATAC-seq ( Corces et al., 2017 ). Here, we describe an ATAC-seq protocol based on the OMNI-ATAC method, with a special focus on the initial steps of thawing cryopreserved cells, and the final steps of library purification using magnetic beads. This protocol can be of interest for laboratories working in a fast-paced environment. Graphic abstract: Flowchart of the protocol.

Establishing community reference samples, data and call sets for benchmarking cancer mutation detection using whole-genome sequencing.

The lack of samples for generating standardized DNA datasets for setting up a sequencing pipeline or benchmarking the performance of different algorithms limits the implementation and uptake of cancer genomics. Here, we describe reference call sets obtained from paired tumor-normal genomic DNA (gDNA) samples derived from a breast cancer cell line-which is highly heterogeneous, with an aneuploid genome, and enriched in somatic alterations-and a matched lymphoblastoid cell line. We partially validated both somatic mutations and germline variants in these call sets via whole-exome sequencing (WES) with different sequencing platforms and targeted sequencing with >2,000-fold coverage, spanning 82% of genomic regions with high confidence. Although the gDNA reference samples are not representative of primary cancer cells from a clinical sample, when setting up a sequencing pipeline, they not only minimize potential biases from technologies, assays and informatics but also provide a unique resource for benchmarking 'tumor-only' or 'matched tumor-normal' analyses.

T Cell-Mediated Antitumor Immunity Cooperatively Induced By TGFßR1 Antagonism and Gemcitabine Counteracts Reformation of the Stromal Barrier in Pancreatic Cancer.

The desmoplastic stroma of pancreatic cancers forms a physical barrier that impedes intratumoral drug delivery. Attempts to modulate the desmoplastic stroma to increase delivery of administered chemotherapy have not shown positive clinical results thus far, and preclinical reports in which chemotherapeutic drugs were coadministered with antistromal therapies did not universally demonstrate increased genotoxicity despite increased intratumoral drug levels. In this study, we tested whether TGFß antagonism can break the stromal barrier, enhance perfusion and tumoral drug delivery, and interrogated cellular and molecular mechanisms by which the tumor prevents synergism with coadministered gemcitabine. TGFß inhibition in genetically engineered murine models (GEMM) of pancreas cancer enhanced tumoral perfusion and increased intratumoral gemcitabine levels. However, tumors rapidly adapted to TGFß-dependent stromal modulation, and intratumoral perfusion returned to pre-treatment levels upon extended TGFß inhibition. Perfusion was governed by the phenotypic identity and distribution of cancer-associated fibroblasts (CAF) with the myelofibroblastic phenotype (myCAFs), and myCAFs which harbored unique genomic signatures rapidly escaped the restricting effects of TGFß inhibition. Despite the reformation of the stromal barrier and reversal of initially increased intratumoral exposure levels, TGFß inhibition in cooperation with gemcitabine effectively suppressed tumor growth via cooperative reprogramming of T regulatory cells and stimulation of CD8 T cell-mediated antitumor activity. The antitumor activity was further improved by the addition of anti-PD-L1 immune checkpoint blockade to offset adaptive PD-L1 upregulation induced by TGFß inhibition. These findings support the development of combined antistroma anticancer therapies capable of impacting the tumor beyond the disruption of the desmoplastic stroma as a physical barrier to improve drug delivery.

A multicenter study benchmarking single-cell RNA sequencing technologies using reference samples.

Comparing diverse single-cell RNA sequencing (scRNA-seq) datasets generated by different technologies and in different laboratories remains a major challenge. Here we address the need for guidance in choosing algorithms leading to accurate biological interpretations of varied data types acquired with different platforms. Using two well-characterized cellular reference samples (breast cancer cells and B cells), captured either separately or in mixtures, we compared different scRNA-seq platforms and several preprocessing, normalization and batch-effect correction methods at multiple centers. Although preprocessing and normalization contributed to variability in gene detection and cell classification, batch-effect correction was by far the most important factor in correctly classifying the cells. Moreover, scRNA-seq dataset characteristics (for example, sample and cellular heterogeneity and platform used) were critical in determining the optimal bioinformatic method. However, reproducibility across centers and platforms was high when appropriate bioinformatic methods were applied. Our findings offer practical guidance for optimizing platform and software selection when designing an scRNA-seq study.

Optimization for Sequencing and Analysis of Degraded FFPE-RNA Samples.

Gene expression analysis by RNA sequencing (RNA-seq) enables unique insights into clinical samples that can potentially lead to mechanistic understanding of the basis of various diseases as well as resistance and/or susceptibility mechanisms. However, FFPE tissues, which represent the most common method for preserving tissue morphology in clinical specimens, are not the best sources for gene expression profiling analysis. The RNA obtained from such samples is often degraded, fragmented, and chemically modified, which leads to suboptimal sequencing libraries. In turn, these generate poor quality sequence data that may not be reliable for gene expression analysis and mutation discovery. In order to make the most of FFPE samples and obtain the best possible data from low quality samples, it is important to take certain precautions while planning experimental design, preparing sequencing libraries, and during data analysis. This includes the use of appropriate metrics for precise sample quality control (QC), identifying the best methods for various steps during the sequencing library generation, and careful library QC. In addition, applying correct software tools and parameters for sequence data analysis is critical in order to identify artifacts in RNA-seq data, filter out contamination and low quality reads, assess uniformity of gene coverage, and measure the reproducibility of gene expression profiles among biological replicates. These steps can ensure high accuracy and reproducibility for profiling of very heterogeneous RNA samples. Here we describe the various steps for sample QC, library preparation and QC, sequencing, and data analysis that can help to increase the amount of useful data obtained from low quality RNA, such as that obtained from FFPE-RNA tissues.

Mannose receptor (CD206) activation in tumor-associated macrophages enhances adaptive and innate antitumor immune responses.

Solid tumors elicit a detectable immune response including the infiltration of tumor-associated macrophages (TAMs). Unfortunately, this immune response is co-opted into contributing toward tumor growth instead of preventing its progression. We seek to reestablish an antitumor immune response by selectively targeting surface receptors and endogenous signaling processes of the macrophage subtypes driving cancer progression. RP-182 is a synthetic 10-mer amphipathic analog of host defense peptides that selectively induces a conformational switch of the mannose receptor CD206 expressed on TAMs displaying an M2-like phenotype. RP-182-mediated activation of this receptor in human and murine M2-like macrophages elicits a program of endocytosis, phagosome-lysosome formation, and autophagy and reprograms M2-like TAMs to an antitumor M1-like phenotype. In syngeneic and autochthonous murine cancer models, RP-182 suppressed tumor growth, extended survival, and was an effective combination partner with chemo- or immune checkpoint therapy. Antitumor activity of RP-182 was also observed in CD206high patient-derived xenotransplantation models. Mechanistically, via selective reduction of immunosuppressive M2-like TAMs, RP-182 improved adaptive and innate antitumor immune responses, including increased cancer cell phagocytosis by reprogrammed TAMs.

T cell receptor and cytokine signal integration in CD8+ T cells is mediated by the protein Themis.

T cell homeostasis and functional responsiveness require signals from self-peptide-major histocompatibility complex (self-pMHC) and cytokines, but the mechanisms controlling this signal integration are unknown. Using a conditional deletion of the T cell lineage-specific protein Themis, we show that Themis is required for the maintenance of peripheral CD8+ T cells and for proliferative CD8+ T cell responses to low-affinity pMHC aided by cytokines. Themis-deficient peripheral T cells show a phenotype indicative of reduced tonic signaling from self-pMHC, strongly suggesting that Themis is a positive regulator of T cell receptor signal strength in response to low-affinity self-pMHC in peripheral T cells. Signals from low-affinity pMHC and cytokines synergistically induce phosphorylation of the kinase Akt, metabolic changes and c-Myc transcription factor induction in CD8+ T cells only in the presence of Themis. This function of Themis is mediated through Shp1 phosphatase, as peripheral Themis and Shp1 double deletion rescues the peripheral CD8+ T cell maintenance.

Single-Cell Profiling Defines Transcriptomic Signatures Specific to Tumor-Reactive versus Virus-Responsive CD4+ T Cells.

Most current tumor immunotherapy strategies leverage cytotoxic CD8+ T cells. Despite evidence for clinical potential of CD4+ tumor-infiltrating lymphocytes (TILs), their functional diversity limits our ability to harness their activity. Here, we use single-cell mRNA sequencing to analyze the response of tumor-specific CD4+ TILs and draining lymph node (dLN) T cells. Computational approaches to characterize subpopulations identify TIL transcriptomic patterns strikingly distinct from acute and chronic anti-viral responses and dominated by diversity among T-bet-expressing T helper type 1 (Th1)-like cells. In contrast, the dLN response includes T follicular helper (Tfh) cells but lacks Th1 cells. We identify a type I interferon-driven signature in Th1-like TILs and show that it is found in human cancers, in which it is negatively associated with response to checkpoint therapy. Our study provides a proof-of-concept methodology to characterize tumor-specific CD4+ T cell effector programs. Targeting these programs should help improve immunotherapy strategies.

Tumor Cell Biodiversity Drives Microenvironmental Reprogramming in Liver Cancer.

Cellular diversity in tumors is a key factor for therapeutic failures and lethal outcomes of solid malignancies. Here, we determined the single-cell transcriptomic landscape of liver cancer biospecimens from 19 patients. We found varying degrees of heterogeneity in malignant cells within and between tumors and diverse landscapes of tumor microenvironment (TME). Strikingly, tumors with higher transcriptomic diversity were associated with patient's worse overall survival. We found a link between hypoxia-dependent vascular endothelial growth factor expression in tumor diversity and TME polarization. Moreover, T cells from higher heterogeneous tumors showed lower cytolytic activities. Consistent results were found using bulk genomic and transcriptomic profiles of 765 liver tumors. Our results offer insight into the diverse ecosystem of liver cancer and its impact on patient prognosis.

Robustness of RNA sequencing on older formalin-fixed paraffin-embedded tissue from high-grade ovarian serous adenocarcinomas.

Formalin-fixed paraffin-embedded (FFPE) tissues are among the most widely available clinical specimens. Their potential utility as a source of RNA for transcriptome studies would greatly enhance population-based cancer studies. Although preliminary studies suggest FFPE tissue may be used for RNA sequencing, the effect of storage time on these specimens needs to be determined. We conducted this study to determine whether RNA in archived FFPE high-grade ovarian serous adenocarcinomas from Surveillance, Epidemiology and End Results (SEER) registries was present in sufficient quantity and quality for RNA-Seq analysis. FFPE tissues, stored from 7 to 32 years, were obtained from three SEER sites. RNA was extracted, quantified, quality assessed, and subjected to RNA-Seq (a whole transcriptome sequencing technology). FFPE specimens stored for longer periods of time had poorer RNA sample quality as indicated by negative correlations between specimen storage time and fragment distribution values (DV). In addition, sample contamination was a common issue among the RNA, with 41 of 67 samples having 5% to 48% bacterial contamination. However, regardless of specimen storage time and bacterial contamination, 60% of the samples yielded data that enabled gene expression quantification, identifying more than 10,000 genes, with the correlations among most biological replicates above 0.7. This study demonstrates that FFPE high-grade ovarian serous adenocarcinomas specimens stored in repositories for up to 32 years and under varying storage conditions are a promising source of RNA for RNA-Seq. We also describe certain caveats to be considered when designing RNA-Seq studies using archived FFPE tissues.

The Emergence and Functional Fitness of Memory CD4+ T Cells Require the Transcription Factor Thpok.

Memory CD4+ T cells mediate long-term immunity, and their generation is a key objective of vaccination strategies. However, the transcriptional circuitry controlling the emergence of memory cells from early CD4+ antigen-responders remains poorly understood. Here, using single-cell RNA-seq to study the transcriptome of virus-specific CD4+ T cells, we identified a gene signature that distinguishes potential memory precursors from effector cells. We found that both that signature and the emergence of memory CD4+ T cells required the transcription factor Thpok. We further demonstrated that Thpok cell-intrinsically protected memory cells from a dysfunctional, effector-like transcriptional program, similar to but distinct from the exhaustion pattern of cells responding to chronic infection. Mechanistically, Thpok- bound genes encoding the transcription factors Blimp1 and Runx3 and acted by antagonizing their expression. Thus, a Thpok-dependent circuitry promotes both memory CD4+ T cells' differentiation and functional fitness, two previously unconnected critical attributes of adaptive immunity.

Genome Assembly and Annotation of the Trichoplusia ni Tni-FNL Insect Cell Line Enabled by Long-Read Technologies.

BACKGROUND: Trichoplusiani derived cell lines are commonly used to enable recombinant protein expression via baculovirus infection to generate materials approved for clinical use and in clinical trials. In order to develop systems biology and genome engineering tools to improve protein expression in this host, we performed de novo genome assembly of the Trichoplusiani-derived cell line Tni-FNL. METHODS: By integration of PacBio single-molecule sequencing, Bionano optical mapping, and 10X Genomics linked-reads data, we have produced a draft genome assembly of Tni-FNL. RESULTS: Our assembly contains 280 scaffolds, with a N50 scaffold size of 2.3 Mb and a total length of 359 Mb. Annotation of the Tni-FNL genome resulted in 14,101 predicted genes and 93.2% of the predicted proteome contained recognizable protein domains. Ortholog searches within the superorder Holometabola provided further evidence of high accuracy and completeness of the Tni-FNL genome assembly. CONCLUSIONS: This first draft Tni-FNL genome assembly was enabled by complementary long-read technologies and represents a high-quality, well-annotated genome that provides novel insight into the complexity of this insect cell line and can serve as a reference for future large-scale genome engineering work in this and other similar recombinant protein production hosts.

Themis-associated phosphatase activity controls signaling in T cell development.

Thymocyte-expressed molecule involved in selection (Themis) has been shown to be important for T cell selection by setting the threshold for positive versus negative selection. Themis interacts with the protein tyrosine phosphatase (PTP) Src-homology domain containing phosphatase-1 (Shp1), a negative regulator of the T cell receptor (TCR) signaling cascade. However, how Themis regulates Shp1 is still not clear. Here, using a very sensitive phosphatase assay on ex vivo thymocytes, we have found that Themis enhances Shp1 phosphatase activity by increasing its phosphorylation. This positive regulation of Shp1 activity by Themis is found in thymocytes, but not in peripheral T cells. Shp1 activity is modulated by different affinity peptide MHC ligand binding in thymocytes. Themis is also associated with phosphatase activity, due to its constitutive interaction with Shp1. In the absence of Shp1 in thymocytes, Themis interacts with Shp2, which leads to almost normal thymic development in Shp1 conditional knockout (cKO) mice. Double deletion of both Themis and Shp1 leads to a thymic phenotype similar to that of Themis KO. These findings demonstrate unequivocally that Themis positively regulates Shp1 phosphatase activity in TCR-mediated signaling in developing thymocytes.

Single-cell analysis reveals cancer stem cell heterogeneity in hepatocellular carcinoma.

Intratumor molecular heterogeneity of hepatocellular carcinoma is partly attributed to the presence of hepatic cancer stem cells (CSCs). Different CSC populations defined by various cell surface markers may contain different oncogenic drivers, posing a challenge in defining molecularly targeted therapeutics. We combined transcriptomic and functional analyses of hepatocellular carcinoma cells at the single-cell level to assess the degree of CSC heterogeneity. We provide evidence that hepatic CSCs at the single-cell level are phenotypically, functionally, and transcriptomically heterogeneous. We found that different CSC subpopulations contain distinct molecular signatures. Interestingly, distinct genes within different CSC subpopulations are independently associated with hepatocellular carcinoma prognosis, suggesting that a diverse hepatic CSC transcriptome affects intratumor heterogeneity and tumor progression. CONCLUSION: Our work provides unique perspectives into the biodiversity of CSC subpopulations, whose molecular heterogeneity further highlights their role in tumor heterogeneity, prognosis, and hepatic CSC therapy. (Hepatology 2018;68:127-140).

SHP1-ing thymic selection.

Thymocyte development and maintenance of peripheral T-cell numbers and functions are critically dependent on T-cell receptor (TCR) signal strength. SHP1 (Src homology region 2 domain-containing phosphatase-1), a tyrosine phosphatase, acts as a negative regulator of TCR signal strength. Moreover, germline SHP1 knockout mice have shown impaired thymic development. However, this has been recently questioned by an analysis of SHP1 conditional knockout mice, which reported normal thymic development of SHP1 deficient thymocytes. Using this SHP1 conditional knockout mice, in this issue of the European Journal of Immunology, Martinez et al. [Eur. J. Immunol. 2016. 46: 2103-2110] show that SHP1 indeed does have a role in the negative regulation of TCR signal strength in positively selected thymocytes, and in the final maturation of single positive thymocytes. They report that thymocyte development in such mice shows loss of mature, post-selection cells. This is due to increased TCR signal transduction in thymocytes immediately post positive-selection, and increased cell death in response to weak TCR ligands. Thus, SHP1-deficiency shows strong similarities to deficiency in the T-cell specific SHP1-associated protein Themis.

Drosophila TAP/p32 is a core histone chaperone that cooperates with NAP-1, NLP, and nucleophosmin in sperm chromatin remodeling during fertilization.

Nuclear DNA in the male gamete of sexually reproducing animals is organized as sperm chromatin compacted primarily by sperm-specific protamines. Fertilization leads to sperm chromatin remodeling, during which protamines are expelled and replaced by histones. Despite our increased understanding of the factors that mediate nucleosome assembly in the nascent male pronucleus, the machinery for protamine removal remains largely unknown. Here we identify four Drosophila protamine chaperones that mediate the dissociation of protamine-DNA complexes: NAP-1, NLP, and nucleophosmin are previously characterized histone chaperones, and TAP/p32 has no known function in chromatin metabolism. We show that TAP/p32 is required for the removal of Drosophila protamine B in vitro, whereas NAP-1, NLP, and Nph share roles in the removal of protamine A. Embryos from P32-null females show defective formation of the male pronucleus in vivo. TAP/p32, similar to NAP-1, NLP, and Nph, facilitates nucleosome assembly in vitro and is therefore a histone chaperone. Furthermore, mutants of P32, Nlp, and Nph exhibit synthetic-lethal genetic interactions. In summary, we identified factors mediating protamine removal from DNA and reconstituted in a defined system the process of sperm chromatin remodeling that exchanges protamines for histones to form the nucleosome-based chromatin characteristic of somatic cells.

Mouse BAZ1A (ACF1) is dispensable for double-strand break repair but is essential for averting improper gene expression during spermatogenesis.

ATP-dependent chromatin remodelers control DNA access for transcription, recombination, and other processes. Acf1 (also known as BAZ1A in mammals) is a defining subunit of the conserved ISWI-family chromatin remodelers ACF and CHRAC, first purified over 15 years ago from Drosophila melanogaster embryos. Much is known about biochemical properties of ACF and CHRAC, which move nucleosomes in vitro and in vivo to establish ordered chromatin arrays. Genetic studies in yeast, flies and cultured human cells clearly implicate these complexes in transcriptional repression via control of chromatin structures. RNAi experiments in transformed mammalian cells in culture also implicate ACF and CHRAC in DNA damage checkpoints and double-strand break repair. However, their essential in vivo roles in mammals are unknown. Here, we show that Baz1a-knockout mice are viable and able to repair developmentally programmed DNA double-strand breaks in the immune system and germ line, I-SceI endonuclease-induced breaks in primary fibroblasts via homologous recombination, and DNA damage from mitomycin C exposure in vivo. However, Baz1a deficiency causes male-specific sterility in accord with its high expression in male germ cells, where it displays dynamic, stage-specific patterns of chromosomal localization. Sterility is caused by pronounced defects in sperm development, most likely a consequence of massively perturbed gene expression in spermatocytes and round spermatids in the absence of BAZ1A: the normal spermiogenic transcription program is largely intact but more than 900 other genes are mis-regulated, primarily reflecting inappropriate up-regulation. We propose that large-scale changes in chromatin composition that occur during spermatogenesis create a window of vulnerability to promiscuous transcription changes, with an essential function of ACF and/or CHRAC chromatin remodeling activities being to safeguard against these alterations.

Rewiring of genetic networks in response to DNA damage.

Although cellular behaviors are dynamic, the networks that govern these behaviors have been mapped primarily as static snapshots. Using an approach called differential epistasis mapping, we have discovered widespread changes in genetic interaction among yeast kinases, phosphatases, and transcription factors as the cell responds to DNA damage. Differential interactions uncover many gene functions that go undetected in static conditions. They are very effective at identifying DNA repair pathways, highlighting new damage-dependent roles for the Slt2 kinase, Pph3 phosphatase, and histone variant Htz1. The data also reveal that protein complexes are generally stable in response to perturbation, but the functional relations between these complexes are substantially reorganized. Differential networks chart a new type of genetic landscape that is invaluable for mapping cellular responses to stimuli.

Individual lysine acetylations on the N terminus of Saccharomyces cerevisiae H2A.Z are highly but not differentially regulated.

The multi-functional histone variant Htz1 (Saccharomyces cerevisiae H2A.Z) is acetylated on up to four N-terminal lysines at positions 3, 8, 10, and 14. It has thus been posited that specific acetylated forms of the histone could regulate distinct roles. Antibodies against Htz1-K8(Ac), -K10(Ac), and -K14(Ac) show that all three modifications are added by Esa1 acetyltransferase and removed by Hda1 deacetylase. Completely unacetylatable htz1 alleles exhibit widespread interactions in genome scale genetic screening. However, singly mutated (e.g. htz1-K8R) or singly acetylable (e.g. the triple mutant htz1-K3R/K10R/K14R) alleles show no significant defects in these analyses. This suggests that the N-terminal acetylations on Htz1 are internally redundant. Further supporting this proposal, each acetylation decays with similar kinetics when Htz1 transcription is repressed, and proteomic screening did not find a single condition in which one Htz1(Ac) was differentially regulated. However, whereas the individual acetylations on Htz1 may be redundant, they are not dispensable. Completely unacetylatable htz1 alleles display genetic interactions and phenotypes in common with and distinct from htz1Δ. In addition, each Htz1 N-terminal lysine is deacetylated by Hda1 in response to benomyl and reacetylated when this agent is removed. Such active regulation suggests that acetylation plays a significant role in Htz1 function.