RESUMEN
A key issue for research on COVID-19 pathogenesis is the lack of biopsies from patients and of samples at the onset of infection. To overcome these hurdles, hamsters were shown to be useful models for studying this disease. Here, we further leverage the model to molecularly survey the disease progression from time-resolved single-cell RNA sequencing data collected from healthy and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected Syrian and Roborovski hamster lungs. We compare our data to human COVID-19 studies, including bronchoalveolar lavage, nasal swab, and postmortem lung tissue, and identify a shared axis of inflammation dominated by macrophages, neutrophils, and endothelial cells, which we show to be transient in Syrian and terminal in Roborovski hamsters. Our data suggest that, following SARS-CoV-2 infection, commitment to a type 1- or type 3-biased immunity determines moderate versus severe COVID-19 outcomes, respectively.
Asunto(s)
COVID-19 , Células Endoteliales , Pulmón , Neutrófilos , SARS-CoV-2 , Análisis de la Célula Individual , COVID-19/inmunología , COVID-19/virología , COVID-19/patología , Animales , Humanos , Neutrófilos/inmunología , SARS-CoV-2/inmunología , Pulmón/patología , Pulmón/virología , Pulmón/inmunología , Cricetinae , Células Endoteliales/virología , Células Endoteliales/patología , Inflamación/patología , Mesocricetus , Modelos Animales de Enfermedad , Masculino , Especificidad de la EspecieRESUMEN
Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.
Asunto(s)
ADN Mitocondrial , Enfermedades Mitocondriales , Humanos , ADN Mitocondrial/genética , Multiómica , Enfermedades Mitocondriales/genética , Mitocondrias/genética , MutaciónRESUMEN
Vaccines play a critical role in combating the COVID-19 pandemic. Future control of the pandemic requires improved vaccines with high efficacy against newly emerging SARS-CoV-2 variants and the ability to reduce virus transmission. Here we compare immune responses and preclinical efficacy of the mRNA vaccine BNT162b2, the adenovirus-vectored spike vaccine Ad2-spike and the live-attenuated virus vaccine candidate sCPD9 in Syrian hamsters, using both homogeneous and heterologous vaccination regimens. Comparative vaccine efficacy was assessed by employing readouts from virus titrations to single-cell RNA sequencing. Our results show that sCPD9 vaccination elicited the most robust immunity, including rapid viral clearance, reduced tissue damage, fast differentiation of pre-plasmablasts, strong systemic and mucosal humoral responses, and rapid recall of memory T cells from lung tissue after challenge with heterologous SARS-CoV-2. Overall, our results demonstrate that live-attenuated vaccines offer advantages over currently available COVID-19 vaccines.
Asunto(s)
COVID-19 , SARS-CoV-2 , Animales , Cricetinae , Humanos , Vacunas Atenuadas , COVID-19/prevención & control , Vacunas contra la COVID-19 , Vacuna BNT162 , Pandemias , MesocricetusRESUMEN
Natural sequence variation within mitochondrial DNA (mtDNA) contributes to human phenotypes and may serve as natural genetic markers in human cells for clonal and lineage tracing. We recently developed a single-cell multi-omic approach, called 'mitochondrial single-cell assay for transposase-accessible chromatin with sequencing' (mtscATAC-seq), enabling concomitant high-throughput mtDNA genotyping and accessible chromatin profiling. Specifically, our technique allows the mitochondrial genome-wide inference of mtDNA variant heteroplasmy along with information on cell state and accessible chromatin variation in individual cells. Leveraging somatic mtDNA mutations, our method further enables inference of clonal relationships among native ex vivo-derived human cells not amenable to genetic engineering-based clonal tracing approaches. Here, we provide a step-by-step protocol for the use of mtscATAC-seq, including various cell-processing and flow cytometry workflows, by using primary hematopoietic cells, subsequent single-cell genomic library preparation and sequencing that collectively take ~3-4 days to complete. We discuss experimental and computational data quality control metrics and considerations for the extension to other mammalian tissues. Overall, mtscATAC-seq provides a broadly applicable platform to map clonal relationships between cells in human tissues, investigate fundamental aspects of mitochondrial genetics and enable additional modes of multi-omic discovery.
Asunto(s)
Secuenciación de Inmunoprecipitación de Cromatina , Cromatina , Animales , Humanos , Cromatina/genética , Multiómica , Análisis de Secuencia de ADN/métodos , Secuenciación de Nucleótidos de Alto Rendimiento/métodos , ADN Mitocondrial/genética , Genotipo , Mamíferos/genéticaRESUMEN
Detailed knowledge of the molecular biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is crucial for understanding of viral replication, host responses, and disease progression. Here, we report gene expression profiles of three SARS-CoV- and SARS-CoV-2-infected human cell lines. SARS-CoV-2 elicited an approximately two-fold higher stimulation of the innate immune response compared to SARS-CoV in the human epithelial cell line Calu-3, including induction of miRNA-155. Single-cell RNA sequencing of infected cells showed that genes induced by virus infections were broadly upregulated, whereas interferon beta/lambda genes, a pro-inflammatory cytokines such as IL-6, were expressed only in small subsets of infected cells. Temporal analysis suggested that transcriptional activities of interferon regulatory factors precede those of nuclear factor κB. Lastly, we identified heat shock protein 90 (HSP90) as a protein relevant for the infection. Inhibition of the HSP90 activity resulted in a reduction of viral replication and pro-inflammatory cytokine expression in primary human airway epithelial cells.
RESUMEN
Characterizing the complex composition of solid tumors is fundamental for understanding tumor initiation, progression and metastasis. While patient-derived samples provide valuable insight, they are heterogeneous on multiple molecular levels, and often originate from advanced tumor stages. Here, we use single-cell transcriptome and epitope profiling together with pathway and lineage analyses to study tumorigenesis from a developmental perspective in a mouse model of salivary gland squamous cell carcinoma. We provide a comprehensive cell atlas and characterize tumor-specific cells. We find that these cells are connected along a reproducible developmental trajectory: initiated in basal cells exhibiting an epithelial-to-mesenchymal transition signature, tumorigenesis proceeds through Wnt-differential cancer stem cell-like subpopulations before differentiating into luminal-like cells. Our work provides unbiased insights into tumor-specific cellular identities in a whole tissue environment, and emphasizes the power of using defined genetic model systems.
Asunto(s)
Carcinogénesis/genética , Carcinogénesis/patología , Animales , Carcinogénesis/inmunología , Carcinoma de Células Escamosas/genética , Carcinoma de Células Escamosas/inmunología , Carcinoma de Células Escamosas/patología , Diferenciación Celular/genética , Diferenciación Celular/inmunología , Linaje de la Célula/genética , Linaje de la Célula/inmunología , Transformación Celular Neoplásica/genética , Transformación Celular Neoplásica/inmunología , Transformación Celular Neoplásica/patología , Transición Epitelial-Mesenquimal/genética , Transición Epitelial-Mesenquimal/inmunología , Perfilación de la Expresión Génica , Ratones , Ratones Endogámicos C57BL , Ratones Transgénicos , Células Madre Neoplásicas/clasificación , Células Madre Neoplásicas/inmunología , Células Madre Neoplásicas/patología , ARN Mensajero/genética , RNA-Seq , Neoplasias de las Glándulas Salivales/genética , Neoplasias de las Glándulas Salivales/inmunología , Neoplasias de las Glándulas Salivales/patología , Análisis de la Célula Individual , Microambiente Tumoral/genética , Microambiente Tumoral/inmunologíaRESUMEN
Herpesvirus infection initiates a range of perturbations in the host cell, which remain poorly understood at the level of individual cells. Here, we quantify the transcriptome of single human primary fibroblasts during the first hours of lytic infection with HSV-1. By applying a generalizable analysis scheme, we define a precise temporal order of early viral gene expression and propose a set-wise emergence of viral genes. We identify host cell genes and pathways relevant for infection by combining three different computational approaches: gene and pathway overdispersion analysis, prediction of cell-state transition probabilities, as well as future cell states. One transcriptional program, which correlates with increased resistance to infection, implicates the transcription factor NRF2. Consequently, Bardoxolone methyl and Sulforaphane, two known NRF2 agonists, impair virus production, suggesting that NRF2 activation restricts viral infection. Our study provides insights into early stages of HSV-1 infection and serves as a general blueprint for the investigation of heterogeneous cell states in virus infection.
Asunto(s)
Fibroblastos/metabolismo , Herpes Simple/genética , Herpesvirus Humano 1 , Interacciones Huésped-Patógeno/genética , Factor 2 Relacionado con NF-E2/genética , Fibroblastos/virología , Perfilación de la Expresión Génica , Redes Reguladoras de Genes , Humanos , Isotiocianatos/farmacología , Factor 2 Relacionado con NF-E2/agonistas , Ácido Oleanólico/análogos & derivados , Ácido Oleanólico/farmacología , Cultivo Primario de Células , Análisis de Secuencia de ARN , Análisis de la Célula Individual , Sulfóxidos , Replicación Viral/efectos de los fármacosRESUMEN
Hedgehog (Hh) signaling plays a key role in the development and maintenance of animal tissues. This signaling is mediated by the atypical G protein-coupled receptor (GPCR) Smoothened (Smo). Smo activation leads to signaling through several well-characterized effectors to activate Hh target gene expression. Recent studies have implicated activation of the heterotrimeric G protein subunit Gαi and the subsequent decrease in cellular cAMP levels in promoting the Hh response in flies and mammals. Although Hh stimulation decreases cAMP levels in some insect cell lines, here using a bioluminescence resonance energy transfer (BRET)-based assay we found that this stimulation had no detectable effect in Drosophila S2-R+ cells. However, we observed an unexpected and significant Gαs-dependent increase in cAMP levels in response to strong Smo activation in Smo-transfected cells. This effect was mediated by Smo's broadly conserved core, and was specifically activated in response to phosphorylation of the Smo C-terminus by GPCR kinase 2 (Gprk2). Genetic analysis of heterotrimeric G protein function in the developing Drosophila wing revealed a positive role for cAMP in the endogenous Hh response. Specifically, we found that mutation or depletion of Gαs diminished low-threshold Hh responses in Drosophila, whereas depletion of Gαi potentiated them (in contrast to previous findings). Our analysis suggested that regulated cAMP production is important for controlling the sensitivity of cellular responses to Hh in Drosophila.
Asunto(s)
AMP Cíclico/metabolismo , Proteínas de Drosophila/metabolismo , Drosophila melanogaster/metabolismo , Subunidades alfa de la Proteína de Unión al GTP Gs/metabolismo , Proteínas Hedgehog/metabolismo , Transducción de Señal , Receptor Smoothened/metabolismo , Animales , FosforilaciónRESUMEN
BACKGROUND: Recent developments in droplet-based microfluidics allow the transcriptional profiling of thousands of individual cells in a quantitative, highly parallel and cost-effective way. A critical, often limiting step is the preparation of cells in an unperturbed state, not altered by stress or ageing. Other challenges are rare cells that need to be collected over several days or samples prepared at different times or locations. METHODS: Here, we used chemical fixation to address these problems. Methanol fixation allowed us to stabilise and preserve dissociated cells for weeks without compromising single-cell RNA sequencing data. RESULTS: By using mixtures of fixed, cultured human and mouse cells, we first showed that individual transcriptomes could be confidently assigned to one of the two species. Single-cell gene expression from live and fixed samples correlated well with bulk mRNA-seq data. We then applied methanol fixation to transcriptionally profile primary cells from dissociated, complex tissues. Low RNA content cells from Drosophila embryos, as well as mouse hindbrain and cerebellum cells prepared by fluorescence-activated cell sorting, were successfully analysed after fixation, storage and single-cell droplet RNA-seq. We were able to identify diverse cell populations, including neuronal subtypes. As an additional resource, we provide 'dropbead', an R package for exploratory data analysis, visualization and filtering of Drop-seq data. CONCLUSIONS: We expect that the availability of a simple cell fixation method will open up many new opportunities in diverse biological contexts to analyse transcriptional dynamics at single-cell resolution.
Asunto(s)
Células Cultivadas/citología , Citometría de Flujo/métodos , Perfilación de la Expresión Génica/métodos , Análisis de la Célula Individual/métodos , Animales , Cerebelo/citología , Drosophila/citología , Embrión no Mamífero/citología , Citometría de Flujo/instrumentación , Perfilación de la Expresión Génica/instrumentación , Humanos , Metanol/química , Ratones , ARN Mensajero/análisis , Rombencéfalo/citología , Análisis de Secuencia de ARN , Análisis de la Célula Individual/instrumentación , Programas InformáticosRESUMEN
BACKGROUND: The C57BL/6J.YA/J mouse strain is a chromosome-substituted line where the original male-specific portion of chromosome Y (MSY) from C57BL/6J mice was substituted for that from A/J mice. In hearts from male C57BL/6J.YA/J and C57BL/6J mice, orchidectomy (ORX) affected in a strictly strain-specific fashion the expression a subset of genes showing enrichment for functional categories, including that of circadian rhythms and cardiac contractility. We further tested whether: (1) there were strain-specific differences in cardiac circadian rhythms; (2) strain-dependent differences in the effects of ORX on contractility genes translated into differences in cardiac functions; and (3) differential contractility responses occurred preferentially at times when circadian rhythms also showed strain-specific differences. METHODS: In hearts from the two above strains, we (1) profiled the expression levels of 15 circadian genes at 4-h intervals across a 24 h period; (2) tested the effects of either ORX or androgen replacement on expression of cardiac contractility genes, and that of ORX on myocardial functional reserve; and (3) verified whether the effects of MSY variants on cardiac contractility-related responses showed synchronicity with differences in circadian rhythms. RESULTS: Among the 15 tested circadian genes, a subset of them were affected by strain (and thus the genetic origin of MSY), which interacted with the amplitude of their peak of maximal expression at 2:00 PM. At that same time-point, ORX decreased (and androgen supplementation increased) the expression of three contractility-related genes, and decreased myocardial relaxation reserve in C57BL/6J.YA/J, but not in C57BL/6J mice. These effects were not detected at 10:00 AM, i.e., at another time-point when circadian genes showed no strain-specific differences. CONCLUSIONS: The results indicate that in mice, androgens have activational effects on cardiac circadian rhythms, contractile gene expression, and myocardial functional reserve. All effects occurred preferentially at the same time of the day, but varied as a function of the genetic origin of MSY. Androgens may therefore be necessary but not sufficient to impart male-specific characteristics to some particular cardiac functions, with genetic material from MSY being one other necessary factor to fully define their range of actions.
RESUMEN
Despite indications that hearts from the C57BL/6N and C57BL/6J mouse substrains differ in terms of their contractility and their responses to stress-induced overload, no information is available about the underlying molecular and cellular mechanisms. We tested whether subacute (48 hours) and chronic (14 days) administration of angiotensin II (500 ng/kg per day) had different effects on the left ventricles of male C57BL/6J and C57BL/6N mice. Despite higher blood pressure in C57BL/6J mice, chronic angiotensin II induced fibrosis and increased the left ventricular weight/body weight ratio and cardiac expression of markers of left ventricular hypertrophy to a greater extent in C57BL/6N mice. Subacute angiotensin II affected a greater number of cardiac genes in C57BL/6N than in C57BL/6J mice. Some of the most prominent differences were observed for markers of (1) macrophage activation and M2 polarization, including 2 genes (osteopontin and galectin-3) whose inactivation was reported as sufficient to prevent angiotensin II-induced myocardial fibrosis; and (2) fibroblast activation. These differences were confirmed in macrophage- and fibroblast-enriched populations of cells isolated from the hearts of experimental mice. When testing F2 animals, the amount of connective tissue present after chronic angiotensin II administration did not cosegregate with the inactivation mutation of the nicotinamide nucleotide transhydrogenase gene from C57BL/6J mice, thus discounting its possible contribution to differences in cardiac remodeling. However, expression levels of osteopontin and galectin-3 were cosegregated in hearts from angiotensin II-treated F2 animals and may represent endophenotypes that could facilitate the identification of genetic regulators of the cardiac fibrogenic response to angiotensin II.
Asunto(s)
Angiotensina II/farmacología , Ventrículos Cardíacos/efectos de los fármacos , Hipertrofia Ventricular Izquierda/patología , Ratones Endogámicos C57BL/genética , Remodelación Ventricular/efectos de los fármacos , Angiotensina II/efectos adversos , Animales , Modelos Animales de Enfermedad , Relación Dosis-Respuesta a Droga , Fibroblastos/efectos de los fármacos , Fibroblastos/patología , Fibrosis/inducido químicamente , Fibrosis/metabolismo , Fibrosis/patología , Galectina 3/metabolismo , Corazón/efectos de los fármacos , Ventrículos Cardíacos/metabolismo , Ventrículos Cardíacos/patología , Hipertrofia Ventricular Izquierda/inducido químicamente , Hipertrofia Ventricular Izquierda/metabolismo , Técnicas In Vitro , Macrófagos/efectos de los fármacos , Macrófagos/patología , Masculino , Ratones , Miocardio/metabolismo , Miocardio/patología , Osteopontina/metabolismo , Remodelación Ventricular/fisiologíaRESUMEN
We have previously reported Lvm1 as a quantitative trait locus (QTL) on chromosome 13 that links to cardiac left ventricular mass (LVM) in a panel of AxB/BxA mouse recombinant inbred strains (RIS). When performing a gene expression QTL (eQTL) analysis, we detected 33 cis-eQTLs that correlated with LVM. Among the latter, a group of eight cis-eQTLs clustered in a genomic region smaller than 6 Mb and surrounding the Lvm1 peak on chr13. Co-variant analysis indicated that all eight genes correlated with the phenotype in a causal rather than a reactive fashion, a finding that (despite its functional interest) did not provide grounds to prioritize any of these candidate genes. As a complementary approach, we performed weighted gene co-expression network analysis, which allowed us to detect 49 modules of highly connected genes. The module that correlated best with LVM: (1) showed linkage to a module QTL whose boundaries matched closely those of the phenotypic Lvm1 QTL on chr13; (2) harbored a disproportionately high proportion of genes originating from a small genomic region on chromosome 13 (including the 8 previously detected cis-eQTL genes); (3) contained genes that, beyond their individual level of expression, correlated with LVM as a function of their inter-connectivity; and (4) showed increased abundance of polymorphic insertion-deletion elements in the same region. Taken together, these data suggest that a domain on chromosome 13 constitutes the biologic principle responsible for the organization and linkage of the gene co-expression module, and indicate a mechanism whereby genetic variants within chromosome domains may associate to phenotypic changes via coordinate changes in the expression of several genes. One other possible implication of these findings is that candidate genes to consider as contributors to a particular phenotype should extend further than those that are closest to the QTL peak.
RESUMEN
Little is known about the functions of chromosome Y (chrY) genes beyond their effects on sex and reproduction. In hearts, postpubertal testosterone affects the size of cells and the expression of genes differently in male C57BL/6J than in their C57.Y(A) counterparts, where the original chrY has been substituted with that from A/J mice. We further compared the 2 strains to better understand how chrY polymorphisms may affect cardiac properties, the latter being sexually dimorphic but unrelated to sex and reproduction. Genomic regions showing occupancy with androgen receptors (ARs) were identified in adult male hearts from both strains by chromatin immunoprecipitation. AR chromatin immunoprecipitation peaks (showing significant enrichment for consensus AR binding sites) were mostly strain specific. Measurements of anogenital distances in male pups showed that the biologic effects of perinatal androgens were greater in C57BL/6J than in C57.Y(A). Although perinatal endocrine manipulations showed that these differences contributed to the strain-specific differences in the response of adult cardiac cells to testosterone, the amounts of androgens produced by fetal testes were not different in each strain. Nonetheless, chrY polymorphisms associated in newborn pups' hearts with strain-specific differences in genomic regions showing either AR occupancy, accessible chromatin sites, or trimethylation of histone H3 Lysine 4 marks, as well as with differential expression of 2 chrY-encoded histone demethylases. In conclusion, the effects of chrY on adult cardiac phenotypes appeared to result from an interaction of this chromosome with the organizational programming effects exerted by the neonatal testosterone surge and show several characteristics of being mediated by an epigenetic remodeling of chromatin.