Challenges in quantifying functional redundancy and selection in microbial communities.

Microbiomes can exhibit large variations in species abundances but high reproducibility in abundances of functional units, an observation often considered evidence for functional redundancy. Based on such reduction in functional variability, selection is hypothesized to act on functional units in these ecosystems. However, the link between functional redundancy and selection remains unclear. Here, we show that reduction in functional variability does not always imply selection on functional profiles. We propose empirical null models to account for the confounding effects of statistical averaging and bias toward environment-independent beneficial functions. We apply our models to existing data sets, and find that the abundances of metabolic groups within microbial communities from bromeliad foliage do not exhibit any evidence of the previously hypothesized functional selection. By contrast, communities of soil bacteria or human gut commensals grown in vitro are selected for metabolic capabilities. By separating the effects of averaging and functional bias on functional variability, we find that the appearance of functional selection in gut microbiome samples from the Human Microbiome Project is artifactual, and that there is no evidence of selection for any molecular function represented by KEGG orthology. These concepts articulate a basic framework for quantifying functional redundancy and selection, advancing our understanding of the mapping between microbiome taxonomy and function.

Resource competition predicts assembly of gut bacterial communities in vitro.

Microbial community dynamics arise through interspecies interactions, including resource competition, cross-feeding and pH modulation. The individual contributions of these mechanisms to community structure are challenging to untangle. Here we develop a framework to estimate multispecies niche overlaps by combining metabolomics data of individual species, growth measurements in spent media and mathematical models. We applied our framework to an in vitro model system comprising 15 human gut commensals in complex media and showed that a simple model of resource competition accounted for most pairwise interactions. Next, we built a coarse-grained consumer-resource model by grouping metabolomic features depleted by the same set of species and showed that this model predicted the composition of 2-member to 15-member communities with reasonable accuracy. Furthermore, we found that incorporation of cross-feeding and pH-mediated interactions improved model predictions of species coexistence. Our theoretical model and experimental framework can be applied to characterize interspecies interactions in bacterial communities in vitro.

Modulation of antibiotic effects on microbial communities by resource competition.

Antibiotic treatment significantly impacts the human gut microbiota, but quantitative understanding of how antibiotics affect community diversity is lacking. Here, we build on classical ecological models of resource competition to investigate community responses to species-specific death rates, as induced by antibiotic activity or other growth-inhibiting factors such as bacteriophages. Our analyses highlight the complex dependence of species coexistence that can arise from the interplay of resource competition and antibiotic activity, independent of other biological mechanisms. In particular, we identify resource competition structures that cause richness to depend on the order of sequential application of antibiotics (non-transitivity), and the emergence of synergistic and antagonistic effects under simultaneous application of multiple antibiotics (non-additivity). These complex behaviors can be prevalent, especially when generalist consumers are targeted. Communities can be prone to either synergism or antagonism, but typically not both, and antagonism is more common. Furthermore, we identify a striking overlap in competition structures that lead to non-transitivity during antibiotic sequences and those that lead to non-additivity during antibiotic combination. In sum, our results establish a broadly applicable framework for predicting microbial community dynamics under deleterious perturbations.

Assembly of gut-derived bacterial communities follows "early-bird" resource utilization dynamics.

Diet can impact host health through changes to the gut microbiota, yet we lack mechanistic understanding linking nutrient availability and microbiota composition. Here, we use thousands of microbial communities cultured in vitro from human feces to uncover simple assembly rules and develop a predictive model of community composition upon addition of single nutrients from central carbon metabolism to a complex medium. Community membership was largely determined by the donor feces, whereas relative abundances were determined by the supplemental carbon source. The absolute abundance of most taxa was independent of the supplementing nutrient, due to the ability of fast-growing organisms to quickly exhaust their niche in the complex medium and then exploit and monopolize the supplemental carbon source. Relative abundances of dominant taxa could be predicted from the nutritional preferences and growth dynamics of species in isolation, and exceptions were consistent with strain-level variation in growth capabilities. Our study reveals that community assembly follows simple rules of nutrient utilization dynamics and provides a predictive framework for manipulating gut commensal communities through nutritional perturbations.

Design, construction, and in vivo augmentation of a complex gut microbiome.

Efforts to model the human gut microbiome in mice have led to important insights into the mechanisms of host-microbe interactions. However, the model communities studied to date have been defined or complex, but not both, limiting their utility. Here, we construct and characterize in vitro a defined community of 104 bacterial species composed of the most common taxa from the human gut microbiota (hCom1). We then used an iterative experimental process to fill open niches: germ-free mice were colonized with hCom1 and then challenged with a human fecal sample. We identified new species that engrafted following fecal challenge and added them to hCom1, yielding hCom2. In gnotobiotic mice, hCom2 exhibited increased stability to fecal challenge and robust colonization resistance against pathogenic Escherichia coli. Mice colonized by either hCom2 or a human fecal community are phenotypically similar, suggesting that this consortium will enable a mechanistic interrogation of species and genes on microbiome-associated phenotypes.

Competition for fluctuating resources reproduces statistics of species abundance over time across wide-ranging microbiotas.

Across diverse microbiotas, species abundances vary in time with distinctive statistical behaviors that appear to generalize across hosts, but the origins and implications of these patterns remain unclear. Here, we show that many of these macroecological patterns can be quantitatively recapitulated by a simple class of consumer-resource models, in which the metabolic capabilities of different species are randomly drawn from a common statistical distribution. Our model parametrizes the consumer-resource properties of a community using only a small number of global parameters, including the total number of resources, typical resource fluctuations over time, and the average overlap in resource-consumption profiles across species. We show that variation in these macroscopic parameters strongly affects the time series statistics generated by the model, and we identify specific sets of global parameters that can recapitulate macroecological patterns across wide-ranging microbiotas, including the human gut, saliva, and vagina, as well as mouse gut and rice, without needing to specify microscopic details of resource consumption. These findings suggest that resource competition may be a dominant driver of community dynamics. Our work unifies numerous time series patterns under a simple model, and provides an accessible framework to infer macroscopic parameters of effective resource competition from longitudinal studies of microbial communities.

Combining CD47 blockade with trastuzumab eliminates HER2-positive breast cancer cells and overcomes trastuzumab tolerance.

Trastuzumab, a targeted anti-human epidermal-growth-factor receptor-2 (HER2) monoclonal antibody, represents a mainstay in the treatment of HER2-positive (HER2+) breast cancer. Although trastuzumab treatment is highly efficacious for early-stage HER2+ breast cancer, the majority of advanced-stage HER2+ breast cancer patients who initially respond to trastuzumab acquire resistance to treatment and relapse, despite persistence of HER2 gene amplification/overexpression. Here, we sought to leverage HER2 overexpression to engage antibody-dependent cellular phagocytosis (ADCP) through a combination of trastuzumab and anti-CD47 macrophage checkpoint immunotherapy. We have previously shown that blockade of CD47, a surface protein expressed by many malignancies (including HER2+ breast cancer), is an effective anticancer therapy. CD47 functions as a "don't eat me" signal through its interaction with signal regulatory protein-α (SIRPα) on macrophages to inhibit phagocytosis. Hu5F9-G4 (magrolimab), a humanized monoclonal antibody against CD47, blocks CD47's "don't eat me" signal, thereby facilitating macrophage-mediated phagocytosis. Preclinical studies have shown that combining Hu5F9-G4 with tumor-targeting antibodies, such as rituximab, further enhances Hu5F9-G4's anticancer effects via ADCP. Clinical trials have additionally demonstrated that Hu5F9-G4, in combination with rituximab, produced objective responses in patients whose diffuse large B cell lymphomas had developed resistance to rituximab and chemotherapy. These studies led us to hypothesize that combining Hu5F9-G4 with trastuzumab would produce an anticancer effect in antibody-dependent cellular cytotoxicity (ADCC)-tolerant HER2+ breast cancer. This combination significantly suppressed the growth of ADCC-tolerant HER2+ breast cancers via Fc-dependent ADCP. Our study demonstrates that combining trastuzumab and Hu5F9-G4 represents a potential new treatment option for HER2+ breast cancer patients, even for patients whose tumors have progressed after trastuzumab.

Hepatitis A Virus Infections Among Men Who Have Sex with Men - Eight U.S. States, 2017-2018.

During 1995-2011, the overall incidence of hepatitis A decreased by 95% in the United States from 12 cases per 100,000 population during 1995 to 0.4 cases per 100,000 population during 2011, and then plateaued during 2012─2015. The incidence increased by 294% during 2016-2018 compared with the incidence during 2013-2015, with most cases occurring among populations at high risk for hepatitis A infection, including persons who use illicit drugs (injection and noninjection), persons who experience homelessness, and men who have sex with men (MSM) (1-3). Previous outbreaks among persons who use illicit drugs and MSM led to recommendations issued in 1996 by the Advisory Committee on Immunization Practices (ACIP) for routine hepatitis A vaccination of persons in these populations (4). Despite these long-standing recommendations, vaccination coverage rates among MSM remain low (5). In 2017, the New York City Department of Health and Mental Hygiene contacted CDC after public health officials noted an increase in hepatitis A infections among MSM. Laboratory testing* of clinical specimens identified strains of the hepatitis A virus (HAV) that subsequently matched strains recovered from MSM in other states. During January 1, 2017-October 31, 2018, CDC received reports of 260 cases of hepatitis A among MSM from health departments in eight states, a substantial increase from the 16 cases reported from all 50 states during 2013-2015. Forty-eight percent (124 of 258) of MSM patients were hospitalized for a median of 3 days. No deaths were reported. In response to these cases, CDC supported state and local health departments with public health intervention efforts to decrease HAV transmission among MSM populations. These efforts included organizing multistate calls among health departments to share information, providing guidance on developing targeted outreach and managing supplies for vaccine campaigns, and conducting laboratory testing of clinical specimens. Targeted outreach for MSM to increase awareness about hepatitis A infection and improve access to vaccination services, such as providing convenient locations for vaccination, are needed to prevent outbreaks among MSM.

A Mechanistic Model of the Regulation of Division Timing by the Circadian Clock in Cyanobacteria.

The cyanobacterium Synechococcus elongatus possesses a circadian clock in the form of a group of proteins whose concentrations and phosphorylation states oscillate with daily periodicity under constant conditions. The circadian clock regulates the cell cycle such that the timing of the cell divisions is biased toward certain times during the circadian period, but the mechanism underlying this phenomenon remains unclear. Here, we propose a mechanism in which a protein limiting for division accumulates at a rate proportional to the cell volume growth and is modulated by the clock. This "modulated rate" model, in which the clock signal is integrated over time to affect division timing, differs fundamentally from the previously proposed "gating" concept, in which the clock is assumed to suppress divisions during a specific time window. We found that although both models can capture the single-cell statistics of division timing in S. elongatus, only the modulated rate model robustly places divisions away from darkness during changes in the environment. Moreover, within the framework of the modulated rate model, existing experiments on S. elongatus are consistent with the simple mechanism that division timing is regulated by the accumulation of a division limiting protein in a phase with genes whose activity peaks at dusk.

Deferiprone: Pan-selective Histone Lysine Demethylase Inhibition Activity and Structure Activity Relationship Study.

Deferiprone (DFP) is a hydroxypyridinone-derived iron chelator currently in clinical use for iron chelation therapy. DFP has also been known to elicit antiproliferative activities, yet the mechanism of this effect has remained elusive. We herein report that DFP chelates the Fe2+ ion at the active sites of selected iron-dependent histone lysine demethylases (KDMs), resulting in pan inhibition of a subfamily of KDMs. Specifically, DFP inhibits the demethylase activities of six KDMs - 2A, 2B, 5C, 6A, 7A and 7B - with low micromolar IC50s while considerably less active or inactive against eleven KDMs - 1A, 3A, 3B, 4A-E, 5A, 5B and 6B. The KDM that is most sensitive to DFP, KDM6A, has an IC50 that is between 7- and 70-fold lower than the iron binding equivalence concentrations at which DFP inhibits ribonucleotide reductase (RNR) activities and/or reduces the labile intracellular zinc ion pool. In breast cancer cell lines, DFP potently inhibits the demethylation of H3K4me3 and H3K27me3, two chromatin posttranslational marks that are subject to removal by several KDM subfamilies which are inhibited by DFP in cell-free assay. These data strongly suggest that DFP derives its anti-proliferative activity largely from the inhibition of a sub-set of KDMs. The docked poses adopted by DFP at the KDM active sites enabled identification of new DFP-based KDM inhibitors which are more cytotoxic to cancer cell lines. We also found that a cohort of these agents inhibited HP1-mediated gene silencing and one lead compound potently inhibited breast tumor growth in murine xenograft models. Overall, this study identified a new chemical scaffold capable of inhibiting KDM enzymes, globally changing histone modification profiles, and with specific anti-tumor activities.

TiO2 nanoparticles generate superoxide and alter gene expression in human lung cells.

TiO2 nanoparticles are widely used in consumer products and industrial applications, yet little is understood regarding how the inhalation of these nanoparticles impacts long-term health. This is especially important for the occupational safety of workers who process these materials. We used RNA sequencing to probe changes in gene expression and fluorescence microscopy to image intracellular reactive oxygen species (ROS) in human lung cells incubated with low, non-cytotoxic, concentrations of TiO2 nanoparticles. Experiments were designed to measure changes in gene expression following an acute exposure to TiO2 nanoparticles and changes inherited by progeny cells. We observe that TiO2 nanoparticles lead to significant (>2000 differentially expressed genes) changes in gene expression following a 24 hour incubation. Following this acute exposure, the response dissipates with only 34 differentially expressed genes in progeny cells. The progeny cells adapt to this initial exposure, observed when re-challenged with a second acute TiO2 nanoparticle exposure. Accompanying these changes in gene expression is the production of intracellular ROS, specifically superoxide, along with changes in oxidative stress-related genes. These experiments suggest that TiO2 nanoparticles adapt to oxidative stress through transcriptional changes over multiple generations of cells.

Modeling Cell Size Regulation: From Single-Cell-Level Statistics to Molecular Mechanisms and Population-Level Effects.

Most microorganisms regulate their cell size. In this article, we review some of the mathematical formulations of the problem of cell size regulation. We focus on coarse-grained stochastic models and the statistics that they generate. We review the biologically relevant insights obtained from these models. We then describe cell cycle regulation and its molecular implementations, protein number regulation, and population growth, all in relation to size regulation. Finally, we discuss several future directions for developing understanding beyond phenomenological models of cell size regulation.

Engagement of MHC class I by the inhibitory receptor LILRB1 suppresses macrophages and is a target of cancer immunotherapy.

Exciting progress in the field of cancer immunotherapy has renewed the urgency of the need for basic studies of immunoregulation in both adaptive cell lineages and innate cell lineages. Here we found a central role for major histocompatibility complex (MHC) class I in controlling the phagocytic function of macrophages. Our results demonstrated that expression of the common MHC class I component ß2-microglobulin (ß2M) by cancer cells directly protected them from phagocytosis. We further showed that this protection was mediated by the inhibitory receptor LILRB1, whose expression was upregulated on the surface of macrophages, including tumor-associated macrophages. Disruption of either MHC class I or LILRB1 potentiated phagocytosis of tumor cells both in vitro and in vivo, which defines the MHC class I-LILRB1 signaling axis as an important regulator of the effector function of innate immune cells, a potential biomarker for therapeutic response to agents directed against the signal-regulatory protein CD47 and a potential target of anti-cancer immunotherapy.

Archaeal cells share common size control with bacteria despite noisier growth and division.

In nature, microorganisms exhibit different volumes spanning six orders of magnitude 1 . Despite their capability to create different sizes, a clonal population in a given environment maintains a uniform size across individual cells. Recent studies in eukaryotic and bacterial organisms showed that this homogeneity in cell size can be accomplished by growing a constant size between two cell cycle events (that is, the adder model 2-6 ). Demonstration of the adder model led to the hypothesis that this phenomenon is a consequence of convergent evolution. Given that archaeal cells share characteristics with both bacteria and eukaryotes, we investigated whether and how archaeal cells exhibit control over cell size. To this end, we developed a soft-lithography method of growing the archaeal cells to enable quantitative time-lapse imaging and single-cell analysis, which would be useful for other microorganisms. Using this method, we demonstrated that Halobacterium salinarum, a hypersaline-adapted archaeal organism, grows exponentially at the single-cell level and maintains a narrow-size distribution by adding a constant length between cell division events. Interestingly, the archaeal cells exhibited greater variability in cell division placement and exponential growth rate across individual cells in a population relative to those observed in Escherichia coli 6-9 . Here, we present a theoretical framework that explains how these larger fluctuations in archaeal cell cycle events contribute to cell size variability and control.

A Parallel Adder Coordinates Mycobacterial Cell-Cycle Progression and Cell-Size Homeostasis in the Context of Asymmetric Growth and Organization.

In model bacteria, such as E. coli and B. subtilis, regulation of cell-cycle progression and cellular organization achieves consistency in cell size, replication dynamics, and chromosome positioning [1-3]. Mycobacteria elongate and divide asymmetrically, giving rise to significant variation in cell size and elongation rate among closely related cells [4, 5]. Given the physical asymmetry of mycobacteria, the models that describe coordination of cellular organization and cell-cycle progression in model bacteria are not directly translatable [1, 2, 6-8]. Here, we used time-lapse microscopy and fluorescent reporters of DNA replication and chromosome positioning to examine the coordination of growth, division, and chromosome dynamics at a single-cell level in Mycobacterium smegmatis (M. smegmatis) and Mycobacterium bovis Bacillus Calmette-Guérin (BCG). By analyzing chromosome and replisome localization, we demonstrated that chromosome positioning is asymmetric and proportional to cell size. Furthermore, we found that cellular asymmetry is maintained throughout the cell cycle and is not established at division. Using measurements and stochastic modeling of mycobacterial cell size and cell-cycle timing in both slow and fast growth conditions, we found that well-studied models of cell-size control are insufficient to explain the mycobacterial cell cycle. Instead, we showed that mycobacterial cell-cycle progression is regulated by an unprecedented mechanism involving parallel adders (i.e., constant growth increments) that start at replication initiation. Together, these adders enable mycobacterial populations to regulate cell size, growth, and heterogeneity in the face of varying environments.

Details Matter: Noise and Model Structure Set the Relationship between Cell Size and Cell Cycle Timing.

Organisms across all domains of life regulate the size of their cells. However, the means by which this is done is poorly understood. We study two abstracted "molecular" models for size regulation: inhibitor dilution and initiator accumulation. We apply the models to two settings: bacteria like Escherichia coli, that grow fully before they set a division plane and divide into two equally sized cells, and cells that form a bud early in the cell division cycle, confine new growth to that bud, and divide at the connection between that bud and the mother cell, like the budding yeast Saccharomyces cerevisiae. In budding cells, delaying cell division until buds reach the same size as their mother leads to very weak size control, with average cell size and standard deviation of cell size increasing over time and saturating up to 100-fold higher than those values for cells that divide when the bud is still substantially smaller than its mother. In budding yeast, both inhibitor dilution or initiator accumulation models are consistent with the observation that the daughters of diploid cells add a constant volume before they divide. This "adder" behavior has also been observed in bacteria. We find that in bacteria an inhibitor dilution model produces adder correlations that are not robust to noise in the timing of DNA replication initiation or in the timing from initiation of DNA replication to cell division (the C+D period). In contrast, in bacteria an initiator accumulation model yields robust adder correlations in the regime where noise in the timing of DNA replication initiation is much greater than noise in the C + D period, as reported previously (Ho and Amir, 2015). In bacteria, division into two equally sized cells does not broaden the size distribution.

Interrogating the Escherichia coli cell cycle by cell dimension perturbations.

Bacteria tightly regulate and coordinate the various events in their cell cycles to duplicate themselves accurately and to control their cell sizes. Growth of Escherichia coli, in particular, follows a relation known as Schaechter's growth law. This law says that the average cell volume scales exponentially with growth rate, with a scaling exponent equal to the time from initiation of a round of DNA replication to the cell division at which the corresponding sister chromosomes segregate. Here, we sought to test the robustness of the growth law to systematic perturbations in cell dimensions achieved by varying the expression levels of mreB and ftsZ We found that decreasing the mreB level resulted in increased cell width, with little change in cell length, whereas decreasing the ftsZ level resulted in increased cell length. Furthermore, the time from replication termination to cell division increased with the perturbed dimension in both cases. Moreover, the growth law remained valid over a range of growth conditions and dimension perturbations. The growth law can be quantitatively interpreted as a consequence of a tight coupling of cell division to replication initiation. Thus, its robustness to perturbations in cell dimensions strongly supports models in which the timing of replication initiation governs that of cell division, and cell volume is the key phenomenological variable governing the timing of replication initiation. These conclusions are discussed in the context of our recently proposed "adder-per-origin" model, in which cells add a constant volume per origin between initiations and divide a constant time after initiation.

Hematopoietic stem cell transplantation in immunocompetent hosts without radiation or chemotherapy.

Hematopoietic stem cell (HSC) transplantation can cure diverse diseases of the blood system, including hematologic malignancies, anemias, and autoimmune disorders. However, patients must undergo toxic conditioning regimens that use chemotherapy and/or radiation to eliminate host HSCs and enable donor HSC engraftment. Previous studies have shown that anti-c-Kit monoclonal antibodies deplete HSCs from bone marrow niches, allowing donor HSC engraftment in immunodeficient mice. We show that host HSC clearance is dependent on Fc-mediated antibody effector functions, and enhancing effector activity through blockade of CD47, a myeloid-specific immune checkpoint, extends anti-c-Kit conditioning to fully immunocompetent mice. The combined treatment leads to elimination of >99% of host HSCs and robust multilineage blood reconstitution after HSC transplantation. This targeted conditioning regimen that uses only biologic agents has the potential to transform the practice of HSC transplantation and enable its use in a wider spectrum of patients.

Simultaneous regulation of cell size and chromosome replication in bacteria.

Bacteria are able to maintain a narrow distribution of cell sizes by regulating the timing of cell divisions. In rich nutrient conditions, cells divide much faster than their chromosomes replicate. This implies that cells maintain multiple rounds of chromosome replication per cell division by regulating the timing of chromosome replications. Here, we show that both cell size and chromosome replication may be simultaneously regulated by the long-standing initiator accumulation strategy. The strategy proposes that initiators are produced in proportion to the volume increase and is accumulated at each origin of replication, and chromosome replication is initiated when a critical amount per origin has accumulated. We show that this model maps to the incremental model of size control, which was previously shown to reproduce experimentally observed correlations between various events in the cell cycle and explains the exponential dependence of cell size on the growth rate of the cell. Furthermore, we show that this model also leads to the efficient regulation of the timing of initiation and the number of origins consistent with existing experimental results.

Macrophages eat cancer cells using their own calreticulin as a guide: roles of TLR and Btk.

Macrophage-mediated programmed cell removal (PrCR) is an important mechanism of eliminating diseased and damaged cells before programmed cell death. The induction of PrCR by eat-me signals on tumor cells is countered by don't-eat-me signals such as CD47, which binds macrophage signal-regulatory protein α to inhibit phagocytosis. Blockade of CD47 on tumor cells leads to phagocytosis by macrophages. Here we demonstrate that the activation of Toll-like receptor (TLR) signaling pathways in macrophages synergizes with blocking CD47 on tumor cells to enhance PrCR. Bruton's tyrosine kinase (Btk) mediates TLR signaling in macrophages. Calreticulin, previously shown to be an eat-me signal on cancer cells, is activated in macrophages for secretion and cell-surface exposure by TLR and Btk to target cancer cells for phagocytosis, even if the cancer cells themselves do not express calreticulin.