Implementierung und Anwendung von bioinformatischen Methoden zur Analyse und Visualisierung von Einzelzell-Sequenzierungsdaten

RNA sequencing (RNA-seq) has become a transformative method to profile genome-wide gene expression and whole transcriptome analysis over the last decade. In recent years, with the development of new technologies, it has become possible to study gene expression at single-cell level. This new advances in single-cell RNA-sequencing has revolutionized the way scientists study biological processes. Single-cell RNA-sequencing has been used in different areas to better understand the underlying mechanisms of biological processes. In particular, single-RNA-sequencing is a suitable method to study infectious diseases. Infection is composed of heterogeneous mechanisms on either the host or pathogen side and the best way to understand the heterogeneity of these mechanisms and how they interact with each other is to study infectious diseases at the single-cell level. Studying infection processes at the single-cell level can reveal not only the heterogeneity but also the dynamics of infection and the interplay between the host and pathogen at the molecular level. In this thesis, we implemented and applied different single-cell RNA-seq technologies to better understand infectious diseases. In the present work, we conducted four independent but related research works to shed light on different aspects of infection biology: ● We took advantage of this novel technology to study the consequences of RSV infection on primary human epithelial cells. The primary human epithelial cells were collected from six donors and cultured in air liquid interface (ALI) cell culture inoculated with respiratory syncytial virus (RSV). In this project, we discovered ciliated cells as the susceptible cell types in RSV infection. We applied viral load as an indicator of infection progression and used it to reconstruct the dynamics of host response to RSV infection. Reconstruction of the dynamics of infection revealed many host genes and pathways that were suppressed or induced as a result of RSV infection. Pathways related to innate immune response and interferon response were suppressed during the progression of infection and on the other hand pathways like protein targeting to endoplasmic reticulum and apoptosis were induced. ● We developed a new method which is capable of sequencing the transcriptome of a bacterium at the single-cell level and potentially can help us to characterize the bacterial heterogeneity during the course of infection. In this research project, bacteria were cultured in three different culture conditions namely Late stationary phase, Anaerobic shock and NaCl shock and we used a poly(A)-independent single-cell RNA-sequencing protocol to sequence bacteria at the single-cell level. In this work, we report the faithful capture of growth-dependent gene expression patterns in individual Salmonella and Pseudomonas bacteria. The results of our analysis showed that not only we could capture transcripts across different RNA classes but also our method is capable of discerning the transcriptome of bacteria across different culture conditions. ● We used single-cell RNA-sequencing technology to characterize the immune cells landscape over the course of atherosclerosis. Atherosclerosis is considered a cardiac disease which is highly related to infections and previous infections with bacteria or viruses is considered as a risk factor for atherosclerosis. We performed single-cell RNA sequencing of aortic CD45+ cells extracted from healthy and atherosclerotic aorta of mice. We managed to find certain cell populations which were specifically present in atherosclerotic mice. One of the atheroschelorotic populations was previously undescribed TREM2high macrophages showing enrichment in Trem2 gene expression. This population of macrophages seemed to be involved in functions like lipid metabolism and catabolism and lesion calcification. This work revealed the phenotypic heterogeneity and immune cells landscape of different immune cell populations at different stages of atherosclerosis. Our work paves the way to better describe the relation between different infectious diseases and cardiovascular diseases. ● We developed a web-based platform called Infection Atlas to browse and visualize single-cell RNA-sequencing data. Infection Atlas platform provides a user-friendly interface to study different aspects of infectious diseases at the single-cell level and can potentially promote targeted approaches to intervene in infectious diseases. This platform which is available at infection-atlas.org in the short term provides a user-friendly interface to browse and visualize different aspects of infectious diseases and in the long-term is expected to be a comprehensive atlas of infection in human and mouse across different tissues and different pathogens. Overall, in this thesis we provide a framework to study infectious diseases at the single cell level with providing novel data analysis methods and this thesis paves the way for future studies to study host-pathogen encounters at the single-cell level.RNA-Sequenzierung (RNA-Seq) ist in den letzten zehn Jahren zu einer revolutionären Technik für genomweite Genexpressionsanalysen, sowie für Gesamt-Transkriptom-Analysen geworden. In den letzten Jahren ist es mit der Entwicklung neuer Technologien möglich geworden die Genexpression auf Einzelzell-Niveau zu untersuchen. Diese Fortschritte in der Einzelzell-RNA- Sequenzierung haben die Art wie Wissenschaftler biologische Prozesse betrachten von Grund auf verändert. Einzelzell-Sequenzierung wird in unterschiedlichen Bereichen angewendet, um die grundlegenden Mechanismen biologischer Prozesse besser zu verstehen. Besonders Einzelzell-Sequenzierung ist eine geeignete Methode, um Infektionskrankheiten zu untersuchen. Infektionen sind durch heterogene Mechanismen auf Wirts- und Erreger Seite gekennzeichnet. Der beste Weg die Heterogenität dieser Mechanismen zu verstehen und wie sie interagieren ist die Analyse von Infektionskrankheiten auf Einzelzell-Niveau. Untersuchungen von Infektionsprozessen auf Einzelzell-Ebene können nicht nur die Heterogenität, sondern auch die Dynamik einer Infektion und das Wechselspiel zwischen Wirt und Pathogen auf molekularer Stufe aufzeigen. In dieser Dissertation wurden unterschiedliche Einzelzell-RNA-Sequenzierung Technologien implementiert und angewandt um ein besseres Verständnis von Infektionskrankheiten zu erlangen. In der vorliegenden Arbeit haben wir vier unabhängige, aber verwandte Forschungsarbeiten durchgeführt, um unterschiedliche Aspekte von Infektionsbiologie näher zu betrachten. ● Wir nutzten die Vorteile dieser neuen Technologie, um die Konsequenzen einer RSV Infektion bei primären humanen Epithelzellen zu untersuchen. Die primären humanen Epithelzellen stammten von sechs Spendern und wurden in Luft-Flüssigkeits-Grenzflächen (ALI) Zellkultur mit dem Respiratorischen Syncytial-Virus (kurz RS-Virus) infiziert. In diesem Projekt konnten wir ciliierte Zellen als anfällige Zelltypen einer RSV Infektion zeigen. Wir haben die Viruslast als Indikator für den Fortschritt der Infektion herangezogen, als auch für die Rekonstruktion der Wirtsantwort Dynamik gegenüber einer RSV Infektion. Die Rekonstruktion der Infektionsdynamik zeigte viele Wirtsgene und Signalwege, die durch die RSV Infektion unterdrückt oder induziert wurden. Signalwege, die mit der angeborenen Immunantwort und der Interferonantwort assoziiert waren, wurden durch die fortschreitende Infektion unterdrückt und andererseits waren Signalwege, wie die Zielsteuerung von Proteinen zum endoplasmatischen Retikulum und Apoptose induziert. ● Wir haben eine neue Methode entwickelt, die es ermöglicht das Transkriptom eines Bakteriums auf Einzelzell-Niveau zu sequenzieren und potenziell helfen könnte die bakterielle Heterogenität während des Verlaufs einer Infektion zu charakterisieren. In diesem Forschungsprojekt wurden Bakterien unter folgenden drei unterschiedlichen Konditionen angezogen: Späte stationäre Phase, anaerober Schock und Natriumchlorid Schock. Anschließend wendeten wir ein poly(A) unabhängiges Einzelzell-RNA Sequenzier-Protokoll an, um Bakterien auf Einzelzell-Niveau zu sequenzieren. In dieser Arbeit berichten wir die von wachstumsabhängigen Genexpressionsmustern in einzelnen Salmonellen und Pseudomonaden. Das Ergebnis unserer Analyse zeigte, dass wir nicht nur Transkripte unterschiedlicher RNA-Klassen, sondern auch das Transkriptom von Bakterien in unterschiedlichen Wachstumsbedingungen erfassen können. ● Wir haben Einzelzell-RNA Sequenzierungs-Technologien verwendet, um die Immunzellen Zusammensetzung während des Verlaufs der Athereosklerose zu betrachten. Die Atherosklerose wird als Herzkrankheit betrachtet, die eng mit Infektionen in Zusammenhang gebracht wird. Vorherige Infektionen mit Bakterien oder Viren werden als Risikofaktor für Atherosklerose angenommen. Wir haben für aortische CD45 Zellen von der gesunden und atherosklerotischen Aorta von Mäusen Einzelzell-RNA-Sequenzierungen durchgeführt. Hierbei konnten wir bestimmte Zellpopulationen identifizieren, die spezifisch in atherosklerotischen Mäusen vorkommen. Eine der athereosklerotischen Populationen war eine zuvor unbeschriebene TREM2high Makrophagen Population, die eine erhöhte Trem2 Genexpression zeigte. Diese Population von Makrophagen schien in Funktionen wie Lipid Metabolismus, Katabolismus, sowie Kalzifizierung von Verletzungen involviert zu sein. Diese Arbeit hat die phänotypische Heterogenität und das Feld unterschiedlicher Immunzellpopulationen in unterschiedlichen Stadien der Atherosklerose aufgezeigt. Unsere Arbeit bereitet den Weg, um die Beziehung zwischen unterschiedlichen Infektionskrankheiten und kardiovaskulären Krankheiten besser zu beschreiben. ● Wir haben eine webbasierte Plattform namens „Infektionsatlas“ entwickelt, um Einzelzell-RNA-Sequenzierungsdaten zu visualisieren und zu durchsuchen. Die „Infektionsatlas“ Plattform stellt eine nutzerfreundliche Oberfläche zur Untersuchung von unterschiedlichen Aspekten von Infektionskrankheiten auf Einzelzell-Niveau bereit und kann möglicherweise zielgerichtete Ansätze voranzutreiben, um Infektionskrankheiten zu verhindern. Diese Plattform, die unter „infection-atlas.org“ verfügbar ist, bietet im Moment eine nutzerfreundliche Oberfläche zum Durchsuchen und Darstellen unterschiedlicher Aspekte von Infektionskrankheiten. Langfristig soll es ein umfangreicher Atlas für Infektionen in Maus und Mesch in unterschiedlichen Geweben und unterschiedlichen Pathogenen. Insgesamt stellen wir in dieser Dissertation einen Rahmen zur Untersuchung von Infektionskrankheiten auf Einzelzell-Ebene mit neuen Methoden für die Datenanalyse zur Verfügung und bereiten den Weg für weitere Studien um Wirts-Pathogen Interaktionen auf Einzellzell-Niveau zu untersuchen

Single-cell RNA-sequencing reports growth-condition-specific global transcriptomes of individual bacteria.

Bacteria respond to changes in their environment with specific transcriptional programmes, but even within genetically identical populations these programmes are not homogenously expressed1. Such transcriptional heterogeneity between individual bacteria allows genetically clonal communities to develop a complex array of phenotypes1, examples of which include persisters that resist antibiotic treatment and metabolically specialized cells that emerge under nutrient-limiting conditions2. Fluorescent reporter constructs have played a pivotal role in deciphering heterogeneous gene expression within bacterial populations3 but have been limited to recording the activity of single genes in a few genetically tractable model species, whereas the vast majority of bacteria remain difficult to engineer and/or even to cultivate. Single-cell transcriptomics is revolutionizing the analysis of phenotypic cell-to-cell variation in eukaryotes, but technical hurdles have prevented its robust application to prokaryotes. Here, using an improved poly(A)-independent single-cell RNA-sequencing protocol, we report the faithful capture of growth-dependent gene expression patterns in individual Salmonella and Pseudomonas bacteria across all RNA classes and genomic regions. These transcriptomes provide important reference points for single-cell RNA-sequencing of other bacterial species, mixed microbial communities and host-pathogen interactions

Time-Resolved scRNA-Seq Tracks the Adaptation of a Sensitive MCL Cell Line to Ibrutinib Treatment.

Since the approval of ibrutinib for relapsed/refractory mantle cell lymphoma (MCL), the treatment of this rare mature B-cell neoplasm has taken a great leap forward. Despite promising efficacy of the Bruton tyrosine kinase inhibitor, resistance arises inevitably and the underlying mechanisms remain to be elucidated. Here, we aimed to decipher the response of a sensitive MCL cell line treated with ibrutinib using time-resolved single-cell RNA sequencing. The analysis uncovered five subpopulations and their individual responses to the treatment. The effects on the B cell receptor pathway, cell cycle, surface antigen expression, and metabolism were revealed by the computational analysis and were validated by molecular biological methods. The observed upregulation of B cell receptor signaling, crosstalk with the microenvironment, upregulation of CD52, and metabolic reprogramming towards dependence on oxidative phosphorylation favor resistance to ibrutinib treatment. Targeting these cellular responses provide new therapy options in MCL

Time-Resolved scRNA-Seq Tracks the Adaptation of a Sensitive MCL Cell Line to Ibrutinib Treatment

Since the approval of ibrutinib for relapsed/refractory mantle cell lymphoma (MCL), the treatment of this rare mature B-cell neoplasm has taken a great leap forward. Despite promising efficacy of the Bruton tyrosine kinase inhibitor, resistance arises inevitably and the underlying mechanisms remain to be elucidated. Here, we aimed to decipher the response of a sensitive MCL cell line treated with ibrutinib using time-resolved single-cell RNA sequencing. The analysis uncovered five subpopulations and their individual responses to the treatment. The effects on the B cell receptor pathway, cell cycle, surface antigen expression, and metabolism were revealed by the computational analysis and were validated by molecular biological methods. The observed upregulation of B cell receptor signaling, crosstalk with the microenvironment, upregulation of CD52, and metabolic reprogramming towards dependence on oxidative phosphorylation favor resistance to ibrutinib treatment. Targeting these cellular responses provide new therapy options in MCL

Vector-borne Trypanosoma brucei parasites develop in artificial human skin and persist as skin tissue forms

Abstract Transmission of Trypanosoma brucei by tsetse flies involves the deposition of the cell cycle-arrested metacyclic life cycle stage into mammalian skin at the site of the fly’s bite. We introduce an advanced human skin equivalent and use tsetse flies to naturally infect the skin with trypanosomes. We detail the chronological order of the parasites’ development in the skin by single-cell RNA sequencing and find a rapid activation of metacyclic trypanosomes and differentiation to proliferative parasites. Here we show that after the establishment of a proliferative population, the parasites enter a reversible quiescent state characterized by slow replication and a strongly reduced metabolism. We term these quiescent trypanosomes skin tissue forms, a parasite population that may play an important role in maintaining the infection over long time periods and in asymptomatic infected individuals

Dynamics of monocyte-derived macrophage diversity in experimental myocardial infarction

Aims Macrophages have a critical and dual role in post-ischaemic cardiac repair, as they can foster both tissue healing and damage. Multiple subsets of tissue resident and monocyte-derived macrophages coexist in the infarcted heart, but their precise identity, temporal dynamics, and the mechanisms regulating their acquisition of discrete states are not fully understood. To address this, we used multi-modal single-cell immune profiling, combined with targeted cell depletion and macrophage fate mapping, to precisely map monocyte/macrophage transitions after experimental myocardial infarction. Methods and results We performed single-cell transcriptomic and cell-surface marker profiling of circulating and cardiac immune cells in mice challenged with acute myocardial infarction, and integrated single-cell transcriptomes obtained before and at 1, 3, 5, 7, and 11 days after infarction. Using complementary strategies of CCR2(+) monocyte depletion and fate mapping of tissue resident macrophages, we determined the origin of cardiac macrophage populations. The macrophage landscape of the infarcted heart was dominated by monocyte-derived cells comprising two pro-inflammatory populations defined as Isg15(hi) and MHCII(+)Il1b(+), alongside non-inflammatory Trem2(hi) cells. Trem2(hi) macrophages were observed in the ischaemic area, but not in the remote viable myocardium, and comprised two subpopulations sequentially populating the heart defined as Trem2(hi)Spp1(hi) monocyte-to-macrophage intermediates, and fully differentiated Trem2(hi)Gdf15(hi) macrophages. Cardiac Trem2(hi) macrophages showed similarities to 'lipid-associated macrophages' found in mouse models of metabolic diseases and were observed in the human heart, indicating conserved features of this macrophage state across diseases and species. Ischaemic injury induced a shift of circulating Ly6C(hi) monocytes towards a Chil3(hi) state with granulocyte-like features, but the acquisition of the Trem2(hi) macrophage signature occurred in the ischaemic tissue. In vitro, macrophages acquired features of the Trem2(hi) signature following apoptotic-cell efferocytosis. Conclusion Our work provides a comprehensive map of monocyte/macrophage transitions in the ischaemic heart, constituting a valuable resource for further investigating how these cells may be harnessed and modulated to promote post-ischaemic heart repair

Dynamics of Cardiac Neutrophil Diversity in Murine Myocardial Infarction.

We employed single-cell transcriptomics combined with cell surface epitope detection by sequencing to investigate temporal neutrophil diversity in the blood and heart after murine myocardial infarction. At day 1, 3, and 5 after infarction, cardiac Ly6G+ (lymphocyte antigen 6G) neutrophils could be delineated into 6 distinct clusters with specific time-dependent patterning and proportions. At day 1, neutrophils were characterized by a gene expression profile proximal to bone marrow neutrophils (Cd177, Lcn2, Fpr1), and putative activity of transcriptional regulators involved in hypoxic response (Hif1a) and emergency granulopoiesis (Cebpb). At 3 and 5 days, 2 major subsets of Siglecfhi (enriched for eg, Icam1 and Tnf) and Siglecflow (Slpi, Ifitm1) neutrophils were found. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) analysis in blood and heart revealed that while circulating neutrophils undergo a process of aging characterized by loss of surface CD62L and upregulation of Cxcr4, heart infiltrating neutrophils acquired a unique SiglecFhi signature. SiglecFhi neutrophils were absent from the bone marrow and spleen, indicating local acquisition of the SiglecFhi signature. Reducing the influx of blood neutrophils by anti-Ly6G treatment increased proportions of cardiac SiglecFhi neutrophils, suggesting accumulation of locally aged neutrophils. Computational analysis of ligand/receptor interactions revealed putative pathways mediating neutrophil to macrophage communication in the myocardium. Finally, SiglecFhi neutrophils were also found in atherosclerotic vessels, revealing that they arise across distinct contexts of cardiovascular inflammation

Initial HCV infection of adult hepatocytes triggers a temporally structured transcriptional program containing diverse pro- and anti-viral elements.

Transcriptional profiling provides global snapshots of virus-mediated cellular reprogramming, which can simultaneously encompass pro- and antiviral components. To determine early transcriptional signatures associated with HCV infection of authentic target cells, we performed ex vivo infections of adult primary human hepatocytes (PHHs) from seven donors. Longitudinal sampling identified minimal gene dysregulation at six hours post infection (hpi). In contrast, at 72 hpi, massive increases in the breadth and magnitude of HCV-induced gene dysregulation were apparent, affecting gene classes associated with diverse biological processes. Comparison with HCV-induced transcriptional dysregulation in Huh-7.5 cells identified limited overlap between the two systems. Of note, in PHHs, HCV infection initiated broad upregulation of canonical interferon (IFN)-mediated defense programs, limiting viral RNA replication and abrogating virion release. We further find that constitutive expression of IRF1 in PHHs maintains a steady-state antiviral program in the absence of infection, which can additionally reduce HCV RNA translation and replication. We also detected infection-induced downregulation of ∼90 genes encoding components of the EIF2 translation initiation complex and ribosomal subunits in PHHs, consistent with a signature of translational shutoff. As HCV polyprotein translation occurs independently of the EIF2 complex, this process is likely pro-viral: only translation initiation of host transcripts is arrested. The combination of antiviral intrinsic and inducible immunity, balanced against pro-viral programs, including translational arrest, maintains HCV replication at a low-level in PHHs. This may ultimately keep HCV under the radar of extra-hepatocyte immune surveillance while initial infection is established, promoting tolerance, preventing clearance and facilitating progression to chronicity.IMPORTANCEAcute HCV infections are often asymptomatic and therefore frequently undiagnosed. We endeavored to recreate this understudied phase of HCV infection using explanted PHHs and monitored host responses to initial infection. We detected temporally distinct virus-induced perturbations in the transcriptional landscape, which were initially narrow but massively amplified in breadth and magnitude over time. At 72 hpi, we detected dysregulation of diverse gene programs, concurrently promoting both virus clearance and virus persistence. On the one hand, baseline expression of IRF1 combined with infection-induced upregulation of IFN-mediated effector genes suppresses virus propagation. On the other, we detect transcriptional signatures of host translational inhibition, which likely reduces processing of IFN-regulated gene transcripts and facilitates virus survival. Together, our data provide important insights into constitutive and virus-induced transcriptional programs in PHHs, and identifies simultaneous antagonistic dysregulation of pro-and anti-viral programs which may facilitate host tolerance and promote viral persistence

Atlas of the Immune Cell Repertoire in Mouse Atherosclerosis Defined by Single-Cell RNA-Sequencing and Mass Cytometry.

Using single-cell RNA-sequencing of aortic leukocytes from chow diet- and Western diet-fed Apoe-/- and Ldlr-/- mice, we detected 11 principal leukocyte clusters with distinct phenotypic and spatial characteristics while the cellular repertoire in healthy aortas was less diverse. Gene set enrichment analysis on the single-cell level established that multiple pathways, such as for lipid metabolism, proliferation, and cytokine secretion, were confined to particular leukocyte clusters. Leukocyte populations were differentially regulated in atherosclerotic Apoe-/- and Ldlr-/- mice. We confirmed the phenotypic diversity of these clusters with a novel mass cytometry 35-marker panel with metal-labeled antibodies and conventional flow cytometry. Cell populations retrieved by these protein-based approaches were highly correlated to transcriptionally defined clusters. In an integrated screening strategy of single-cell RNA-sequencing, mass cytometry, and fluorescence-activated cell sorting, we detected 3 principal B-cell subsets with alterations in surface markers, functional pathways, and in vitro cytokine secretion. Leukocyte cluster gene signatures revealed leukocyte frequencies in 126 human plaques by a genetic deconvolution strategy. This approach revealed that human carotid plaques and microdissected mouse plaques were mostly populated by macrophages, T-cells, and monocytes. In addition, the frequency of genetically defined leukocyte populations in carotid plaques predicted cardiovascular events in patients

Neonatally imprinted stromal cell subsets induce tolerogenic dendritic cells in mesenteric lymph nodes

Gut-draining mesenteric lymph nodes (mLNs) are important for inducing peripheral tolerance towards food and commensal antigens by providing an optimal microenvironment for de novo generation of Foxp3(+) regulatory T cells (Tregs). We previously identified microbiota-imprinted mLN stromal cells as a critical component in tolerance induction. Here we show that this imprinting process already takes place in the neonatal phase, and renders the mLN stromal cell compartment resistant to inflammatory perturbations later in life. LN transplantation and single-cell RNA-seq uncover stably imprinted expression signatures in mLN fibroblastic stromal cells. Subsetting common stromal cells across gut-draining mLNs and skin-draining LNs further refine their location-specific immunomodulatory functions, such as subset-specific expression of Aldh1a2/3. Finally, we demonstrate that mLN stromal cells shape resident dendritic cells to attain high Treg-inducing capacity in a Bmp2-dependent manner. Thus, crosstalk between mLN stromal and resident dendritic cells provides a robust regulatory mechanism for the maintenance of intestinal tolerance