Daniel Schulz

Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast

To obtain rates of mRNA synthesis and decay in yeast, we established dynamic transcriptome analysis (DTA). DTA combines non‐perturbing metabolic RNA labeling with dynamic kinetic modeling. DTA reveals that most mRNA synthesis rates are around several transcripts per cell and cell cycle, and most mRNA half‐lives range around a median of 11 min. DTA can monitor the cellular response to osmotic stress with higher sensitivity and temporal resolution than standard transcriptomics. In contrast to monotonically increasing total mRNA levels, DTA reveals three phases of the stress response. During the initial shock phase, mRNA synthesis and decay rates decrease globally, resulting in mRNA storage. During the subsequent induction phase, both rates increase for a subset of genes, resulting in production and rapid removal of stress‐responsive mRNAs. During the recovery phase, decay rates are largely restored, whereas synthesis rates remain altered, apparently enabling growth at high salt concentration. Stress‐induced changes in mRNA synthesis rates are predicted from gene occupancy with RNA polymerase II. DTA‐derived mRNA synthesis rates identified 16 stress‐specific pairs/triples of cooperative transcription factors, of which seven were known. Thus, DTA realistically monitors the dynamics in mRNA metabolism that underlie gene regulatory systems.

Comparative dynamic transcriptome analysis (cDTA) reveals mutual feedback between mRNA synthesis and degradation

To monitor eukaryotic mRNA metabolism, we developed comparative dynamic transcriptome analysis (cDTA). cDTA provides absolute rates of mRNA synthesis and decay in Saccharomyces cerevisiae (Sc) cells with the use of Schizosaccharomyces pombe (Sp) as an internal standard. cDTA uses nonperturbing metabolic labeling that supersedes conventional methods for mRNA turnover analysis. cDTA reveals that Sc and Sp transcripts that encode orthologous proteins have similar synthesis rates, whereas decay rates are fivefold lower in Sp, resulting in similar mRNA concentrations despite the larger Sp cell volume. cDTA of Sc mutants reveals that a eukaryote can buffer mRNA levels. Impairing transcription with a point mutation in RNA polymerase (Pol) II causes decreased mRNA synthesis rates as expected, but also decreased decay rates. Impairing mRNA degradation by deleting deadenylase subunits of the Ccr4-Not complex causes decreased decay rates as expected, but also decreased synthesis rates. Extended kinetic modeling reveals mutual feedback between mRNA synthesis and degradation that may be achieved by a factor that inhibits synthesis and enhances degradation.

An Immune Atlas of Clear Cell Renal Cell Carcinoma

Summary Immune cells in the tumor microenvironment modulate cancer progression and are attractive therapeutic targets. Macrophages and T cells are key components of the microenvironment, yet their phenotypes and relationships in this ecosystem and to clinical outcomes are ill defined. We used mass cytometry with extensive antibody panels to perform in-depth immune profiling of samples from 73 clear cell renal cell carcinoma (ccRCC) patients and five healthy controls. In 3.5 million measured cells, we identified 17 tumor-associated macrophage phenotypes, 22 T cell phenotypes, and a distinct immune composition correlated with progression-free survival, thereby presenting an in-depth human atlas of the immune tumor microenvironment in this disease. This study revealed potential biomarkers and targets for immunotherapy development and validated tools that can be used for immune profiling of other tumor types.

Rpb4 Subunit Functions Mainly in mRNA Synthesis by RNA Polymerase II

Background: The RNA polymerase II (Pol II) subunit Rpb4 was reported to function in mRNA degradation. Results: Rpb4 deletion causes defects in mRNA synthesis and compensatory changes in mRNA degradation. Covalent attachment of Rpb4 to the Pol II core largely restores mRNA metabolism. Conclusion: Rpb4 is an integral part of Pol II and functions mainly in Pol II transcription. Significance: Rpb4 functions mainly in transcription. RNA polymerase II (Pol II) is the central enzyme that carries out eukaryotic mRNA transcription and consists of a 10-subunit catalytic core and a subcomplex of subunits Rpb4 and Rpb7 (Rpb4/7). Rpb4/7 has been proposed to dissociate from Pol II, enter the cytoplasm, and function there in mRNA translation and degradation. Here we provide evidence that Rpb4 mainly functions in nuclear mRNA synthesis by Pol II, as well as evidence arguing against an important cytoplasmic role in mRNA degradation. We used metabolic RNA labeling and comparative Dynamic Transcriptome Analysis to show that Rpb4 deletion in Saccharomyces cerevisiae causes a drastic defect in mRNA synthesis that is compensated by down-regulation of mRNA degradation, resulting in mRNA level buffering. Deletion of Rpb4 can be rescued by covalent fusion of Rpb4 to the Pol II core subunit Rpb2, which largely restores mRNA synthesis and degradation defects caused by Rpb4 deletion. Thus, Rpb4 is a bona fide Pol II core subunit that functions mainly in mRNA synthesis.

Measurement of genome-wide RNA synthesis and decay rates with Dynamic Transcriptome Analysis (DTA)

Standard transcriptomics measures total cellular RNA levels. Our understanding of gene regulation would be greatly improved if we could measure RNA synthesis and decay rates on a genome-wide level. To that end, the Dynamic Transcriptome Analysis (DTA) method has been developed. DTA combines metabolic RNA labeling with standard transcriptomics to measure RNA synthesis and decay rates in a precise and non-perturbing manner. Here, we present the open source R/Bioconductor software package DTA. It implements all required bioinformatics steps that allow the accurate absolute quantification and comparison of RNA turnover.

histoCAT: analysis of cell phenotypes and interactions in multiplex image cytometry data

Single-cell, spatially resolved omics analysis of tissues is poised to transform biomedical research and clinical practice. We have developed an open-source, computational histology topography cytometry analysis toolbox (histoCAT) to enable interactive, quantitative, and comprehensive exploration of individual cell phenotypes, cell–cell interactions, microenvironments, and morphological structures within intact tissues. We highlight the unique abilities of histoCAT through analysis of highly multiplexed mass cytometry images of human breast cancer tissues.

Simultaneous Multiplexed Imaging of mRNA and Proteins with Subcellular Resolution in Breast Cancer Tissue Samples by Mass Cytometry

Summary To build comprehensive models of cellular states and interactions in normal and diseased tissue, genetic and proteomic information must be extracted with single-cell and spatial resolution. Here, we extended imaging mass cytometry to enable multiplexed detection of mRNA and proteins in tissues. Three mRNA target species were detected by RNAscope-based metal in situ hybridization with simultaneous antibody detection of 16 proteins. Analysis of 70 breast cancer samples showed that HER2 and CK19 mRNA and protein levels are moderately correlated on the single-cell level, but that only HER2, and not CK19, has strong mRNA-to-protein correlation on the cell population level. The chemoattractant CXCL10 was expressed in stromal cell clusters, and the frequency of CXCL10-expressing cells correlated with T cell presence. Our flexible and expandable method will allow an increase in the information content retrieved from patient samples for biomedical purposes, enable detailed studies of tumor biology, and serve as a tool to bridge comprehensive genomic and proteomic tissue analysis.

GenoGAM: genome-wide generalized additive models for ChIP-Seq analysis

Motivation: Chromatin immunoprecipitation followed by deep sequencing (ChIP‐Seq) is a widely used approach to study protein‐DNA interactions. Often, the quantities of interest are the differential occupancies relative to controls, between genetic backgrounds, treatments, or combinations thereof. Current methods for differential occupancy of ChIP‐Seq data rely however on binning or sliding window techniques, for which the choice of the window and bin sizes are subjective. Results: Here, we present GenoGAM (Genome‐wide Generalized Additive Model), which brings the well‐established and flexible generalized additive models framework to genomic applications using a data parallelism strategy. We model ChIP‐Seq read count frequencies as products of smooth functions along chromosomes. Smoothing parameters are objectively estimated from the data by cross‐validation, eliminating ad hoc binning and windowing needed by current approaches. GenoGAM provides base‐level and region‐level significance testing for full factorial designs. Application to a ChIP‐Seq dataset in yeast showed increased sensitivity over existing differential occupancy methods while controlling for type I error rate. By analyzing a set of DNA methylation data and illustrating an extension to a peak caller, we further demonstrate the potential of GenoGAM as a generic statistical modeling tool for genome‐wide assays. Availability and Implementation: Software is available from Bioconductor: https://www.bioconductor.org/packages/release/bioc/html/GenoGAM.html. Contact: gagneur@in.tum.de Supplementary information: Supplementary information is available at Bioinformatics online.

Bringing together what belongs together: Optimizing murine infection models by using mouse-adapted Staphylococcus aureus strains

Staphylococcus (S.) aureus is a leading cause of bacterial infection world-wide, and currently no vaccine is available for humans. Vaccine development relies heavily on clinically relevant infection models. However, the suitability of mice for S. aureus infection models has often been questioned, because experimental infection of mice with human-adapted S. aureus requires very high infection doses. Moreover, mice were not considered to be natural hosts of S. aureus. The latter has been disproven by our recent findings, showing that both laboratory mice, as well as wild small mammals including mice, voles, and shrews, are naturally colonized with S. aureus. Here, we investigated whether mouse-and vole-derived S. aureus strains show an enhanced virulence in mice as compared to the human-adapted strain Newman. Using a step-wise approach based on the bacterial genotype and in vitro assays for host adaptation, we selected the most promising candidates for murine infection models out of a total of 254 S. aureus isolates from laboratory mice as well as wild rodents and shrews. Four strains representing the clonal complexes (CC) 8, 49, and 88 (n = 2) were selected and compared to the human-adapted S. aureus strain Newman (CC8) in murine pneumonia and bacteremia models. Notably, a bank vole-derived CC49 strain, named DIP, was highly virulent in BALB/c mice in pneumonia and bacteremia models, whereas the other murine and vole strains showed virulence similar to or lower than that of Newman. At one tenth of the standard infection dose DIP induced disease severity, bacterial load and host cytokine and chemokine responses in the murine bacteremia model similar to that of Newman. In the pneumonia model, DIP was also more virulent than Newman but the effect was less pronounced. Whole genome sequencing data analysis identified a pore-forming toxin gene, lukF-PV(P83)/lukM, in DIP but not in the other tested S. aureus isolates. To conclude, the mouse-adapted S. aureus strain DIP allows a significant reduction of the inoculation dose in mice and is hence a promising tool to develop clinically more relevant infection models.

Natural Colonization of Laboratory Mice with Staphylococcus aureus Primes a Systemic Immune Response

Background Whether mice are an appropriate model for S. aureus infection and vaccination studies is a matter of debate, because they are not considered as natural hosts of S. aureus. Sparked by an outbreak of S. aureus infections in laboratory mice, we investigated whether laboratory mice are commonly colonized with S. aureus and whether this might impact on infection experiments. Methods We characterized 99 S. aureus isolates from laboratory mice (spa typing, virulence gene PCR), and quantified murine antibodies using FlexMap technology. Results Specific-pathogen-free mice from various vendors were frequently colonized with S. aureus (0-21%). S. aureus was readily transmitted from murine parents to offspring, which became persistently colonized. Most murine isolates belonged to the lineage CC88 (54%). Murine strains showed features of host adaptation, such as absence of hlb-converting phages and superantigen genes, as well as enhanced coagulation of murine plasma. Importantly, S. aureus colonization induced a systemic IgG response specific for numerous S. aureus proteins, including several vaccine candidates. Conclusion Laboratory mice are natural hosts of S. aureus and, therefore, provide better infection models than previously assumed. Pre-exposure to S. aureus is a possible confounder in S. aureus infection and vaccination studies.