Cellular adaptation to cancer therapy along a resistance continuum.

Advancements in precision oncology over the past decades have led to new therapeutic interventions, but the efficacy of such treatments is generally limited by an adaptive process that fosters drug resistance1. In addition to genetic mutations2, recent research has identified a role for non-genetic plasticity in transient drug tolerance3 and the acquisition of stable resistance4,5. However, the dynamics of cell-state transitions that occur in the adaptation to cancer therapies remain unknown and require a systems-level longitudinal framework. Here we demonstrate that resistance develops through trajectories of cell-state transitions accompanied by a progressive increase in cell fitness, which we denote as the 'resistance continuum'. This cellular adaptation involves a stepwise assembly of gene expression programmes and epigenetically reinforced cell states underpinned by phenotypic plasticity, adaptation to stress and metabolic reprogramming. Our results support the notion that epithelial-to-mesenchymal transition or stemness programmes-often considered a proxy for phenotypic plasticity-enable adaptation, rather than a full resistance mechanism. Through systematic genetic perturbations, we identify the acquisition of metabolic dependencies, exposing vulnerabilities that can potentially be exploited therapeutically. The concept of the resistance continuum highlights the dynamic nature of cellular adaptation and calls for complementary therapies directed at the mechanisms underlying adaptive cell-state transitions.

Host inflammatory dynamics reveal placental immune modulation by Group B Streptococcus during pregnancy.

Group B Streptococcus (GBS) is a pathobiont that can ascend to the placenta and cause adverse pregnancy outcomes, in part through production of the toxin ß-hemolysin/cytolysin (ß-h/c). Innate immune cells have been implicated in the response to GBS infection, but the impact of ß-h/c on their response is poorly defined. We show that GBS modulates innate immune cell states by subversion of host inflammation through ß-h/c, allowing worse outcomes. We used an ascending mouse model of GBS infection to measure placental cell state changes over time following infection with a ß-h/c-deficient and isogenic wild type GBS strain. Transcriptomic analysis suggests that ß-h/c-producing GBS elicit a worse phenotype through suppression of host inflammatory signaling in placental macrophages and neutrophils, and comparison of human placental macrophages infected with the same strains recapitulates these results. Our findings have implications for identification of new targets in GBS disease to support host defense against pathogenic challenge.

The tempo and mode of gene regulatory programs during bacterial infection.

Innate immune recognition of bacterial pathogens is a key determinant of the ensuing systemic response, and host or pathogen heterogeneity in this early interaction can impact the course of infection. To gain insight into host response heterogeneity, we investigate macrophage inflammatory dynamics using primary human macrophages infected with Group B Streptococcus. Transcriptomic analysis reveals discrete cellular states within responding macrophages, one of which consists of four sub-states, reflecting inflammatory activation. Infection with six additional bacterial species-Staphylococcus aureus, Listeria monocytogenes, Enterococcus faecalis, Yersinia pseudotuberculosis, Shigella flexneri, and Salmonella enterica-recapitulates these states, though at different frequencies. We show that modulating the duration of infection and the presence of a toxin impacts inflammatory trajectory dynamics. We provide evidence for this trajectory in infected macrophages in an in vivo model of Staphylococcus aureus infection. Our cell-state analysis defines a framework for understanding inflammatory activation dynamics in response to bacterial infection.

Bimodal Evolutionary Developmental miRNA Program in Animal Embryogenesis.

miRNAs play essential roles in the mechanics of gene regulation, however, on an organismal-scale, the processes in which they are deployed are not well understood. Here, we adopt an evolutionary developmental approach to study miRNA function by examining their expression throughout embryogenesis in both Caenorhabditis elegans and Drosophila melanogaster. We find that, in both species, miRNA transcriptomic shifts in a punctuated fashion during the mid-developmental transition, specifying two dominant modes of early and late expression profiles. Strikingly, late-expressed miRNAs are enriched for phylogenetic conservation and function by fine-tuning the expression of their targets, implicating a role in the canalization of cell types during differentiation. In contrast, early expressed miRNAs are inversely expressed with their targets suggesting strong target-inhibition. Taken together, our work exposes a bimodal role for miRNA function during animal development, involving late-expressed physiological roles and early expressed repressive roles.

scDual-Seq: mapping the gene regulatory program of Salmonella infection by host and pathogen single-cell RNA-sequencing.

The interaction between a pathogen and a host is a highly dynamic process in which both agents activate complex programs. Here, we introduce a single-cell RNA-sequencing method, scDual-Seq, that simultaneously captures both host and pathogen transcriptomes. We use it to study the process of infection of individual mouse macrophages with the intracellular pathogen Salmonella typhimurium. Among the infected macrophages, we find three subpopulations and we show evidence for a linear progression through these subpopulations, supporting a model in which these three states correspond to consecutive stages of infection.

CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq.

Single-cell transcriptomics requires a method that is sensitive, accurate, and reproducible. Here, we present CEL-Seq2, a modified version of our CEL-Seq method, with threefold higher sensitivity, lower costs, and less hands-on time. We implemented CEL-Seq2 on Fluidigm's C1 system, providing its first single-cell, on-chip barcoding method, and we detected gene expression changes accompanying the progression through the cell cycle in mouse fibroblast cells. We also compare with Smart-Seq to demonstrate CEL-Seq2's increased sensitivity relative to other available methods. Collectively, the improvements make CEL-Seq2 uniquely suited to single-cell RNA-Seq analysis in terms of economics, resolution, and ease of use.

Seeing is believing: new methods for in situ single-cell transcriptomics.

New methods employ RNA-seq to study single cells within complex tissues by in situ sequencing or mRNA capture from single photoactivated cells.

RNA-DNA differences in human mitochondria restore ancestral form of 16S ribosomal RNA.

RNA transcripts are generally identical to the underlying DNA sequences. Nevertheless, RNA-DNA differences (RDDs) were found in the nuclear human genome and in plants and animals but not in human mitochondria. Here, by deep sequencing of human mitochondrial DNA (mtDNA) and RNA, we identified three RDD sites at mtDNA positions 295 (C-to-U), 13710 (A-to-U, A-to-G), and 2617 (A-to-U, A-to-G). Position 2617, within the 16S rRNA, harbored the most prevalent RDDs (>30% A-to-U and ∼15% A-to-G of the reads in all tested samples). The 2617 RDDs appeared already at the precursor polycistrone mitochondrial transcript. By using traditional Sanger sequencing, we identified the A-to-U RDD in six different cell lines and representative primates (Gorilla gorilla, Pongo pigmaeus, and Macaca mulatta), suggesting conservation of the mechanism generating such RDD. Phylogenetic analysis of more than 1700 vertebrate mtDNA sequences supported a thymine as the primate ancestral allele at position 2617, suggesting that the 2617 RDD recapitulates the ancestral 16S rRNA. Modeling U or G (the RDDs) at position 2617 stabilized the large ribosomal subunit structure in contrast to destabilization by an A (the pre-RDDs). Hence, these mitochondrial RDDs are likely functional.

Mitochondrial DNA heteroplasmy in diabetes and normal adults: role of acquired and inherited mutational patterns in twins.

Heteroplasmy, the mixture of mitochondrial genomes (mtDNA), varies among individuals and cells. Heteroplasmy levels alter the penetrance of pathological mtDNA mutations, and the susceptibility to age-related diseases such as Parkinson's disease. Although mitochondrial dysfunction occurs in age-related type 2 diabetes mellitus (T2DM), the involvement of heteroplasmy in diabetes is unclear. We hypothesized that the heteroplasmic mutational (HM) pattern may change in T2DM. To test this, we used next-generation sequencing, i.e. massive parallel sequencing (MPS), along with PCR-cloning-Sanger sequencing to analyze HM in blood and skeletal muscle DNA samples from monozygotic (MZ) twins either concordant or discordant for T2DM. Great variability was identified in the repertoires and amounts of HMs among individuals, with a tendency towards more mutations in skeletal muscle than in blood. Whereas many HMs were unique, many were either shared among twin pairs or among tissues of the same individual, regardless of their prevalence. This suggested a heritable influence on even low abundance HMs. We found no clear differences between T2DM and controls. However, we found ~5-fold increase of HMs in non-coding sequences implying the influence of negative selection (P < 0.001). This negative selection was evident both in moderate to highly abundant heteroplasmy (>5% of the molecules per sample) and in low abundance heteroplasmy (<5% of the molecules). Although our study found no evidence supporting the involvement of HMs in the etiology of T2DM, the twin study found clear evidence of a heritable influence on the accumulation of HMs as well as the signatures of selection in heteroplasmic mutations.