Aging is closely associated with increased susceptibility to breast cancer, yet there have been limited systematic studies of aging-induced alterations in the mammary gland. Here, we leverage high-throughput single-cell RNA sequencing to generate a detailed transcriptomic atlas of young and aged murine mammary tissues. By analyzing epithelial, stromal, and immune cells, we identify age-dependent alterations in cell proportions and gene expression, providing evidence that suggests alveolar maturation and physiological decline. The analysis also uncovers potential pro-tumorigenic mechanisms coupled to the age-associated loss of tumor suppressor function and change in microenvironment. In addition, we identify a rare, age-dependent luminal population co-expressing hormone-sensing and secretory-alveolar lineage markers, as well as two macrophage populations expressing distinct gene signatures, underscoring the complex heterogeneity of the mammary epithelia and stroma. Collectively, this rich single-cell atlas reveals the effects of aging on mammary physiology and can serve as a useful resource for understanding aging-associated cancer risk.

Lineage tracing provides key insights into the fate of individual cells in complex organisms. Although effective genetic labeling approaches are available in model systems, in humans, most approaches require detection of nuclear somatic mutations, which have high error rates, limited scale, and do not capture cell state information. Here, we show that somatic mutations in mtDNA can be tracked by single-cell RNA or assay for transposase accessible chromatin (ATAC) sequencing. We leverage somatic mtDNA mutations as natural genetic barcodes and demonstrate their utility as highly accurate clonal markers to infer cellular relationships. We track native human cells both in vitro and in vivo and relate clonal dynamics to gene expression and chromatin accessibility. Our approach should allow clonal tracking at a 1,000-fold greater scale than with nuclear genome sequencing, with simultaneous information on cell state, opening the way to chart cellular dynamics in human health and disease.

The use of whole-genome phylogenetic analysis has revolutionized our understanding of the evolution and spread of many important bacterial pathogens due to the high resolution view it provides. However, the majority of such analyses do not consider the potential role of accessory genes when inferring evolutionary trajectories. Moreover, the recently discovered importance of the switching of gene regulatory elements suggests that an exhaustive analysis, combining information from core and accessory genes with regulatory elements could provide unparalleled detail of the evolution of a bacterial population. Here we demonstrate this principle by applying it to a worldwide multi-host sample of the important pathogenic E. coli lineage ST131. Our approach reveals the existence of multiple circulating subtypes of the major drug-resistant clade of ST131 and provides the first ever population level evidence of core genome substitutions in gene regulatory regions associated with the acquisition and maintenance of different accessory genome elements.

Human cancer cell lines are the workhorse of cancer research. Although cell lines are known to evolve in culture, the extent of the resultant genetic and transcriptional heterogeneity and its functional consequences remain understudied. Here we use genomic analyses of 106 human cell lines grown in two laboratories to show extensive clonal diversity. Further comprehensive genomic characterization of 27 strains of the common breast cancer cell line MCF7 uncovered rapid genetic diversification. Similar results were obtained with multiple strains of 13 additional cell lines. Notably, genetic changes were associated with differential activation of gene expression programs and marked differences in cell morphology and proliferation. Barcoding experiments showed that cell line evolution occurs as a result of positive clonal selection that is highly sensitive to culture conditions. Analyses of single-cell-derived clones demonstrated that continuous instability quickly translates into heterogeneity of the cell line. When the 27 MCF7 strains were tested against 321 anti-cancer compounds, we uncovered considerably different drug responses: at least 75% of compounds that strongly inhibited some strains were completely inactive in others. This study documents the extent, origins and consequences of genetic variation within cell lines, and provides a framework for researchers to measure such variation in efforts to support maximally reproducible cancer research.

Here we present an extensive genomic and genetic analysis of Escherichia coli strains of serotype O78 that represent the major cause of avian colisepticemia, an invasive infection caused by avian pathogenic Escherichia coli (APEC) strains. It is associated with high mortality and morbidity, resulting in significant economic consequences for the poultry industry. To understand the genetic basis of the virulence of avian septicemic E. coli, we sequenced the entire genome of a clinical isolate of serotype O78-O78:H19 ST88 isolate 789 (O78-9)-and compared it with three publicly available APEC O78 sequences and one complete genome of APEC serotype O1 strain. Although there was a large variability in genome content between the APEC strains, several genes were conserved, which are potentially critical for colisepticemia. Some of these genes are present in multiple copies per genome or code for gene products with overlapping function, signifying their importance. A systematic deletion of each of these virulence-related genes identified three systems that are conserved in all septicemic strains examined and are critical for serum survival, a prerequisite for septicemia. These are the plasmid-encoded protein, the defective ETT2 (E. coli type 3 secretion system 2) type 3 secretion system ETT2sepsis, and iron uptake systems. Strain O78-9 is the only APEC O78 strain that also carried the regulon coding for yersiniabactin, the iron binding system of the Yersinia high-pathogenicity island. Interestingly, this system is the only one that cannot be complemented by other iron uptake systems under iron limitation and in serum.

Anticancer therapies have been limited by the emergence of mutations and other adaptations. In bacteria, antibiotics activate the SOS response, which mobilizes error-prone factors that allow for continuous replication at the cost of mutagenesis. We investigated whether the treatment of lung cancer with EGFR inhibitors (EGFRi) similarly engages hypermutators. In cycling drug-tolerant persister (DTP) cells and in EGFRi-treated patients presenting residual disease, we observed upregulation of GAS6, whereas ablation of GAS6's receptor, AXL, eradicated resistance. Reciprocally, AXL overexpression enhanced DTP survival and accelerated the emergence of T790M, an EGFR mutation typical to resistant cells. Mechanistically, AXL induces low-fidelity DNA polymerases and activates their organizer, RAD18, by promoting neddylation. Metabolomics uncovered another hypermutator, AXL-driven activation of MYC, and increased purine synthesis that is unbalanced by pyrimidines. Aligning anti-AXL combination treatments with the transition from DTPs to resistant cells cured patient-derived xenografts. Hence, similar to bacteria, tumors tolerate therapy by engaging pharmacologically targetable endogenous mutators.

Drug-resistant extraintestinal pathogenic Escherichia coli (ExPEC) strains are the major cause of colisepticemia (colibacillosis), a condition that has become an increasing public health problem in recent years. ExPEC strains are characterized by high resistance to serum, which is otherwise highly toxic to most bacteria. To understand how these bacteria survive and grow in serum, we performed system-wide analyses of their response to serum, making a clear distinction between the responses to nutritional immunity and innate immunity. Thus, mild heat inactivation of serum destroys the immune complement and abolishes the bactericidal effect of serum (inactive serum), making it possible to examine nutritional immunity. We used a combination of deep RNA sequencing and proteomics in order to characterize ExPEC genes whose expression is affected by the nutritional stress of serum and by the immune complement. The major change in gene expression induced by serum-active and inactive-involved metabolic genes. In particular, the serum metabolic response is coordinated by three transcriptional regulators, Fur, BasR, and CysB. Fur alone was responsible for more than 80% of the serum-induced transcriptional response. Consistent with its role as a major serum response regulator, deletion of Fur renders the bacteria completely serum sensitive. These results highlight the role of metabolic adaptation in colisepticemia and virulence.