Replicative history marks transcriptional and functional disparity in the CD8+ T cell memory pool.

Clonal expansion is a core aspect of T cell immunity. However, little is known with respect to the relationship between replicative history and the formation of distinct CD8+ memory T cell subgroups. To address this issue, we developed a genetic-tracing approach, termed the DivisionRecorder, that reports the extent of past proliferation of cell pools in vivo. Using this system to genetically 'record' the replicative history of different CD8+ T cell populations throughout a pathogen-specific immune response, we demonstrate that the central memory T (TCM) cell pool is marked by a higher number of prior divisions than the effector memory T cell pool, owing to the combination of strong proliferative activity during the acute immune response and selective proliferative activity after pathogen clearance. Furthermore, by combining DivisionRecorder analysis with single-cell transcriptomics and functional experiments, we show that replicative history identifies distinct cell pools within the TCM compartment. Specifically, we demonstrate that lowly divided TCM cells display enriched expression of stem-cell-associated genes, exist in a relatively quiescent state, and are superior in eliciting a proliferative recall response upon activation. These data provide the first evidence that a stem-cell-like memory T cell pool that reconstitutes the CD8+ T cell effector pool upon reinfection is marked by prior quiescence.

Lineage tracing reveals B cell antibody class switching is stochastic, cell-autonomous, and tuneable.

To optimize immunity to pathogens, B lymphocytes generate plasma cells with functionally diverse antibody isotypes. By lineage tracing single cells within differentiating B cell clones, we identified the heritability of discrete fate controlling mechanisms to inform a general mathematical model of B cell fate regulation. Founder cells highly influenced clonal plasma-cell fate, whereas class switch recombination (CSR) was variegated within clones. In turn, these CSR patterns resulted from independent all-or-none expression of both activation-induced cytidine deaminase (AID) and IgH germline transcription (GLT), with the latter being randomly re-expressed after each cell division. A stochastic model premised on these molecular transition rules accurately predicted antibody switching outcomes under varied conditions in vitro and during an immune response in vivo. Thus, the generation of functionally diverse antibody types follows rules of autonomous cellular programming that can be adapted and modeled for the rational control of antibody classes for potential therapeutic benefit.

The Branching Point in Erythro-Myeloid Differentiation.

Development of mature blood cell progenies from hematopoietic stem cells involves the transition through lineage-restricted progenitors. The first branching point along this developmental process is thought to separate the erythro-myeloid and lymphoid lineage fate by yielding two intermediate progenitors, the common myeloid and the common lymphoid progenitors (CMPs and CLPs). Here, we use single-cell lineage tracing to demonstrate that so-called CMPs are highly heterogeneous with respect to cellular output, with most individual CMPs yielding either only erythrocytes or only myeloid cells after transplantation. Furthermore, based on the labeling of earlier progenitors, we show that the divergence between the myeloid and erythroid lineage develops within multipotent progenitors (MPP). These data provide evidence for a model of hematopoietic branching in which multiple distinct lineage commitments occur in parallel within the MPP pool.

Towards developing forensically relevant single-cell pipelines by incorporating direct-to-PCR extraction: compatibility, signal quality, and allele detection.

Current analysis of forensic DNA stains relies on the probabilistic interpretation of bulk-processed samples that represent mixed profiles consisting of an unknown number of potentially partial representations of each contributor. Single-cell methods, in contrast, offer a solution to the forensic DNA mixture problem by incorporating a step that separates cells before extraction. A forensically relevant single-cell pipeline relies on efficient direct-to-PCR extractions that are compatible with standard downstream forensic reagents. Here we demonstrate the feasibility of implementing single-cell pipelines into the forensic process by exploring four metrics of electropherogram (EPG) signal quality-i.e., allele detection rates, peak heights, peak height ratios, and peak height balance across low- to high-molecular-weight short tandem repeat (STR) markers-obtained with four direct-to-PCR extraction treatments and a common post-PCR laboratory procedure. Each treatment was used to extract DNA from 102 single buccal cells, whereupon the amplification reagents were immediately added to the tube and the DNA was amplified/injected using post-PCR conditions known to elicit a limit of detection (LoD) of one DNA molecule. The results show that most cells, regardless of extraction treatment, rendered EPGs with at least a 50% true positive allele detection rate and that allele drop-out was not cell independent. Statistical tests demonstrated that extraction treatments significantly impacted all metrics of EPG quality, where the Arcturus® PicoPure™ extraction method resulted in the lowest median allele drop-out rate, highest median average peak height, highest median average peak height ratio, and least negative median values of EPG sloping for GlobalFiler™ STR loci amplified at half volume. We, therefore, conclude the feasibility of implementing single-cell pipelines for casework purposes and demonstrate that inferential systems assuming cell independence will not be appropriate in the probabilistic interpretation of a collection of single-cell EPGs.

Multiplexed Division Tracking Dyes for Proliferation-Based Clonal Lineage Tracing.

The generation of cellular heterogeneity is an essential feature of immune responses. Understanding the heritability and asymmetry of phenotypic changes throughout this process requires determination of clonal-level contributions to fate selection. Evaluating intraclonal and interclonal heterogeneity and the influence of distinct fate determinants in large numbers of cell lineages, however, is usually laborious, requiring familial tracing and fate mapping. In this study, we introduce a novel, accessible, high-throughput method for measuring familial fate changes with accompanying statistical tools for testing hypotheses. The method combines multiplexing of division tracking dyes with detection of phenotypic markers to reveal clonal lineage properties. We illustrate the method by studying in vitro-activated mouse CD8+ T cell cultures, reporting division and phenotypic changes at the level of families. This approach has broad utility as it is flexible and adaptable to many cell types and to modifications of in vitro, and potentially in vivo, fate monitoring systems.

Sample path properties of the average generation of a Bellman-Harris process.

Motivated by a recently proposed design for a DNA coded randomised algorithm that enables inference of the average generation of a collection of cells descendent from a common progenitor, here we establish strong convergence properties for the average generation of a super-critical Bellman-Harris process. We further extend those results to a two-type Bellman-Harris process where one type can give rise to the other, but not vice versa. These results further affirm the estimation method's potential utility by establishing its long run accuracy on individual sample-paths, and significantly expanding its remit to encompass cellular development that gives rise to differentiated offspring with distinct population dynamics.

Exploring STR signal in the single- and multicopy number regimes: Deductions from an in silico model of the entire DNA laboratory process.

Short tandem repeat (STR) profiling from DNA samples has long been the bedrock of human identification. The laboratory process is composed of multiple procedures that include quantification, sample dilution, PCR, electrophoresis, and fragment analysis. The end product is a short tandem repeat electropherogram comprised of signal from allele, artifacts, and instrument noise. In order to optimize or alter laboratory protocols, a large number of validation samples must be created at significant expense. As a tool to support that process and to enable the exploration of complex scenarios without costly sample creation, a mechanistic stochastic model that incorporates each of the aforementioned processing features is described herein. The model allows rapid in silico simulation of electropherograms from multicontributor samples and enables detailed investigations of involved scenarios. An implementation of the model that is parameterized by extensive laboratory data is publically available. To illustrate its utility, the model was employed in order to evaluate the effects of sample dilutions, injection time, and cycle number on peak height, and the nature of stutter ratios at low template. We verify the model's findings by comparison with experimentally generated data.

Inferring average generation via division-linked labeling.

For proliferating cells subject to both division and death, how can one estimate the average generation number of the living population without continuous observation or a division-diluting dye? In this paper we provide a method for cell systems such that at each division there is an unlikely, heritable one-way label change that has no impact other than to serve as a distinguishing marker. If the probability of label change per cell generation can be determined and the proportion of labeled cells at a given time point can be measured, we establish that the average generation number of living cells can be estimated. Crucially, the estimator does not depend on knowledge of the statistics of cell cycle, death rates or total cell numbers. We explore the estimator's features through comparison with physiologically parameterized stochastic simulations and extrapolations from published data, using it to suggest new experimental designs.

Single-cell investigative genetics: Single-cell data produces genotype distributions concentrated at the true genotype across all mixture complexities.

In the absence of a suspect the forensic aim is investigative, and the focus is one of discerning what genotypes best explain the evidence. In traditional systems, the list of candidate genotypes may become vast if the sample contains DNA from many donors or the information from a minor contributor is swamped by that of major contributors, leading to lower evidential value for a true donor's contribution and, as a result, possibly overlooked or inefficient investigative leads. Recent developments in single-cell analysis offer a way forward, by producing data capable of discriminating genotypes. This is accomplished by first clustering single-cell data by similarity without reference to a known genotype. With good clustering it is reasonable to assume that the scEPGs in a cluster are of a single contributor. With that assumption we determine the probability of a cluster's content given each possible genotype at each locus, which is then used to determine the posterior probability mass distribution for all genotypes by application of Bayes' rule. A decision criterion is then applied such that the sum of the ranked probabilities of all genotypes falling in the set is at least 1-α. This is the credible genotype set and is used to inform database search criteria. Within this work we demonstrate the salience of single-cell analysis by performance testing a set of 630 previously constructed admixtures containing up to 5 donors of balanced and unbalanced contributions. We use scEPGs that were generated by isolating single cells, employing a direct-to-PCR extraction treatment, amplifying STRs that are compliant with existing national databases and applying post-PCR treatments that elicit a detection limit of one DNA copy. We determined that, for these test data, 99.3% of the true genotypes are included in the 99.8% credible set, regardless of the number of donors that comprised the mixture. We also determined that the most probable genotype was the true genotype for 97% of the loci when the number of cells in a cluster was at least two. Since efficient investigative leads will be borne by posterior mass distributions that are narrow and concentrated at the true genotype, we report that, for this test set, 47,900 (86%) loci returned only one credible genotype and of these 47,551 (99%) were the true genotype. When determining the LR for true contributors, 91% of the clusters rendered LR>1018, showing the potential of single-cell data to positively affect investigative reporting.

Optimal age-specific vaccination control for COVID-19: An Irish case study.

The outbreak of a novel coronavirus causing severe acute respiratory syndrome in December 2019 has escalated into a worldwide pandemic. In this work, we propose a compartmental model to describe the dynamics of transmission of infection and use it to obtain the optimal vaccination control. The model accounts for the various stages of the vaccination, and the optimisation is focused on minimising the infections to protect the population and relieve the healthcare system. As a case study, we selected the Republic of Ireland. We use data provided by Ireland's COVID-19 Data-Hub and simulate the evolution of the pandemic with and without the vaccination in place for two different scenarios, one representative of a national lockdown situation and the other indicating looser restrictions in place. One of the main findings of our work is that the optimal approach would involve a vaccination programme where the older population is vaccinated in larger numbers earlier while simultaneously part of the younger population also gets vaccinated to lower the risk of transmission between groups. We compare our simulated results with those of the vaccination policy taken by the Irish government to explore the advantages of our optimisation method. Our comparison suggests that a similar reduction in cases may have been possible even with a reduced set of vaccinations available for use.

Simultaneous Assessment of Kinship, Division Number, and Phenotype via Flow Cytometry for Hematopoietic Stem and Progenitor Cells.

Few techniques can assess phenotype and fate for the same cell simultaneously. Most of the current protocols used to characterize phenotype, although able to generate large datasets, necessitate the destruction of the cell of interest, making it impossible to assess its functional fate. Heterogeneous biological differentiating systems like hematopoiesis are therefore difficult to describe. Building on cell division tracking dyes, we further developed a protocol to simultaneously determine kinship, division number, and differentiation status for many single hematopoietic progenitors. This protocol allows the assessment of the ex vivo differentiation potential of murine and human hematopoietic progenitors, isolated from various biological sources. Moreover, as it is based on flow cytometry and a limited number of reagents, it can quickly generate a large amount of data, at the single-cell level, in a relatively inexpensive manner. We also provide the analytical pipeline for single-cell analysis, combined with a robust statistical framework. As this protocol allows the linking of cell division and differentiation at the single-cell level, it can be used to quantitatively assess symmetric and asymmetric fate commitment, the balance between self-renewal and differentiation, and the number of divisions for a given commitment fate. Altogether, this protocol can be used in experimental designs aiming to unravel the biological differences between hematopoietic progenitors, from a single-cell perspective.

Evidentiary evaluation of single cells renders highly informative forensic comparisons across multifarious admixtures.

The consistency between DNA evidence and person(s) of interest (PoI) is summarized by a likelihood ratio (LR): the probability of the data given the PoI contributed divided by the probability given they did not. It is often the case that there are several PoI who may have individually or jointly contributed to the stain. If there is more than one PoI, or the number of contributors (NoC) cannot easily be determined, then several sets of hypotheses are needed, requiring significant resources to complete the interpretation. Recent technological developments in laboratory systems offer a way forward, by enabling production of single cell data. Though single-cell data may be procured by next generation sequencing or capillary electrophoresis workflows, in this work we focus our attention on assessing the consistency between PoIs and a collection of single cell electropherograms (scEPGs) from diploid cells - i.e., leukocytes and epithelial cells. Specifically, we introduce a framework that: I) clusters scEPGs into collections, each originating from one genetic source; II) for each PoI, determines a LR for each cluster of scEPGs; and III) by averaging the likelihood ratios for each PoI across all clusters provides a whole-sample weight of evidence summary. By using Model Based Clustering (MBC) in step I) and an algorithm, named EESCIt for Evidentiary Evaluation of Single Cells, that computes single-cell LRs in step II), we show that 99% of the comparisons rendered log LR values > 0 for true contributors, and of these all but one gave log LR > 5, regardless of the number of donors or whether the smallest contributor donated less than 20% of the cells, greatly expanding the collection of cases for which DNA forensics provides informative results.

Maintaining naivety of T cells.

Mathematical models encoding biological hypotheses reveal new insight into the dynamics of naive immune cells in mice from birth to old age.

High-quality data from a forensically relevant single-cell pipeline enabled by low PBS and proteinase K concentrations.

Interpreting forensic DNA signal is arduous since the total intensity is a cacophony of signal from noise, artifact, and allele from an unknown number of contributors (NOC). An alternate to traditional bulk-processing pipelines is a single-cell one, where the sample is collected, and each cell is sequestered resulting in n single-source, single-cell EPGs (scEPG) that must be interpreted using applicable strategies. As with all forensic DNA interpretation strategies, high quality electropherograms are required; thus, to enhance the credibility of single-cell forensics, it is necessary to produce an efficient direct-to-PCR treatment that is compatible with prevailing downstream laboratory processes. We incorporated the semi-automated micro-fluidic DEPArray™ technology into the single-cell laboratory and optimized its implementation by testing the effects of four laboratory treatments on single-cell profiles. We focused on testing effects of phosphate buffer saline (PBS) since it is an important reagent that mitigates cell rupture but is also a PCR inhibitor. Specifically, we explored the effect of decreasing PBS concentrations on five electropherogram-quality metrics from 241 leukocytes: profile drop-out, allele drop-out, allele peak heights, peak height ratios, and scEPG sloping. In an effort to improve reagent use, we also assessed two concentrations of proteinase K. The results indicate that decreasing PBS concentrations to 0.5X or 0.25X improves scEPG quality, while modest modifications to proteinase K concentrations did not significantly impact it. We, therefore, conclude that a lower than recommended proteinase K concentration coupled with a lower than recommended PBS concentration results in enhanced scEPGs within the semi-automated single-cell pipeline.

The a posteriori probability of the number of contributors when conditioned on an assumed contributor.

Forensic DNA signal is notoriously challenging to assess, requiring computational tools to support its interpretation. Over-expressions of stutter, allele drop-out, allele drop-in, degradation, differential degradation, and the like, make forensic DNA profiles too complicated to evaluate by manual methods. In response, computational tools that make point estimates on the Number of Contributors (NOC) to a sample have been developed, as have Bayesian methods that evaluate an A Posteriori Probability (APP) distribution on the NOC. In cases where an overly narrow NOC range is assumed, the downstream strength of evidence may be incomplete insofar as the evidence is evaluated with an inadequate set of propositions. In the current paper, we extend previous work on NOCIt, a Bayesian method that determines an APP on the NOC given an electropherogram, by reporting on an implementation where the user can add assumed contributors. NOCIt is a continuous system that incorporates models of peak height (including degradation and differential degradation), forward and reverse stutter, noise, and allelic drop-out, while being cognizant of allele frequencies in a reference population. When conditioned on a known contributor, we found that the mode of the APP distribution can shift to one greater when compared with the circumstance where no known contributor is assumed, and that occurred most often when the assumed contributor was the minor constituent to the mixture. In a development of a result of Slooten and Caliebe (FSI:G, 2018) that, under suitable assumptions, establishes the NOC can be treated as a nuisance variable in the computation of a likelihood ratio between the prosecution and defense hypotheses, we show that this computation must not only use coincident models, but also coincident contextual information. The results reported here, therefore, illustrate the power of modern probabilistic systems to assess full weights-of-evidence, and to provide information on reasonable NOC ranges across multiple contexts.

HSPCs display within-family homogeneity in differentiation and proliferation despite population heterogeneity.

High-throughput single-cell methods have uncovered substantial heterogeneity in the pool of hematopoietic stem and progenitor cells (HSPCs), but how much instruction is inherited by offspring from their heterogeneous ancestors remains unanswered. Using a method that enables simultaneous determination of common ancestor, division number, and differentiation status of a large collection of single cells, our data revealed that murine cells that derived from a common ancestor had significant similarities in their division progression and differentiation outcomes. Although each family diversifies, the overall collection of cell types observed is composed of homogeneous families. Heterogeneity between families could be explained, in part, by differences in ancestral expression of cell surface markers. Our analyses demonstrate that fate decisions of cells are largely inherited from ancestor cells, indicating the importance of common ancestor effects. These results may have ramifications for bone marrow transplantation and leukemia, where substantial heterogeneity in HSPC behavior is observed.

Cyton2: A Model of Immune Cell Population Dynamics That Includes Familial Instructional Inheritance.

Lymphocytes are the central actors in adaptive immune responses. When challenged with antigen, a small number of B and T cells have a cognate receptor capable of recognising and responding to the insult. These cells proliferate, building an exponentially growing, differentiating clone army to fight off the threat, before ceasing to divide and dying over a period of weeks, leaving in their wake memory cells that are primed to rapidly respond to any repeated infection. Due to the non-linearity of lymphocyte population dynamics, mathematical models are needed to interrogate data from experimental studies. Due to lack of evidence to the contrary and appealing to arguments based on Occam's Razor, in these models newly born progeny are typically assumed to behave independently of their predecessors. Recent experimental studies, however, challenge that assumption, making clear that there is substantial inheritance of timed fate changes from each cell by its offspring, calling for a revision to the existing mathematical modelling paradigms used for information extraction. By assessing long-term live-cell imaging of stimulated murine B and T cells in vitro, we distilled the key phenomena of these within-family inheritances and used them to develop a new mathematical model, Cyton2, that encapsulates them. We establish the model's consistency with these newly observed fine-grained features. Two natural concerns for any model that includes familial correlations would be that it is overparameterised or computationally inefficient in data fitting, but neither is the case for Cyton2. We demonstrate Cyton2's utility by challenging it with high-throughput flow cytometry data, which confirms the robustness of its parameter estimation as well as its ability to extract biological meaning from complex mixed stimulation experiments. Cyton2, therefore, offers an alternate mathematical model, one that is, more aligned to experimental observation, for drawing inferences on lymphocyte population dynamics.

A large-scale validation of NOCIt's a posteriori probability of the number of contributors and its integration into forensic interpretation pipelines.

Forensic DNA signal is notoriously challenging to interpret and requires the implementation of computational tools that support its interpretation. While data from high-copy, low-contributor samples result in electropherogram signal that is readily interpreted by probabilistic methods, electropherogram signal from forensic stains is often garnered from low-copy, high-contributor-number samples and is frequently obfuscated by allele sharing, allele drop-out, stutter and noise. Since forensic DNA profiles are too complicated to quantitatively assess by manual methods, continuous, probabilistic frameworks that draw inferences on the Number of Contributors (NOC) and compute the Likelihood Ratio (LR) given the prosecution's and defense's hypotheses have been developed. In the current paper, we validate a new version of the NOCIt inference platform that determines an A Posteriori Probability (APP) distribution of the number of contributors given an electropherogram. NOCIt is a continuous inference system that incorporates models of peak height (including degradation and differential degradation), forward and reverse stutter, noise and allelic drop-out while taking into account allele frequencies in a reference population. We established the algorithm's performance by conducting tests on samples that were representative of types often encountered in practice. In total, we tested NOCIt's performance on 815 degraded, UV-damaged, inhibited, differentially degraded, or uncompromised DNA mixture samples containing up to 5 contributors. We found that the model makes accurate, repeatable and reliable inferences about the NOCs and significantly outperformed methods that rely on signal filtering. By leveraging recent theoretical results of Slooten and Caliebe (FSI:G, 2018) that, under suitable assumptions, establish the NOC can be treated as a nuisance variable, we demonstrated that when NOCIt's APP is used in conjunction with a downstream likelihood ratio (LR) inference system that employs the same probabilistic model, a full evaluation across multiple contributor numbers is rendered. This work, therefore, illustrates the power of modern probabilistic systems to report holistic and interpretable weights-of-evidence to the trier-of-fact without assigning a specified number of contributors or filtering signal.

Manipulating niche composition limits damage to haematopoietic stem cells during Plasmodium infection.

Severe infections are a major stress on haematopoiesis, where the consequences for haematopoietic stem cells (HSCs) have only recently started to emerge. HSC function critically depends on the integrity of complex bone marrow (BM) niches; however, what role the BM microenvironment plays in mediating the effects of infection on HSCs remains an open question. Here, using a murine model of malaria and combining single-cell RNA sequencing, mathematical modelling, transplantation assays and intravital microscopy, we show that haematopoiesis is reprogrammed upon infection, whereby the HSC compartment turns over substantially faster than at steady-state and HSC function is drastically affected. Interferon is found to affect both haematopoietic and mesenchymal BM cells and we specifically identify a dramatic loss of osteoblasts and alterations in endothelial cell function. Osteo-active parathyroid hormone treatment abolishes infection-triggered HSC proliferation and-coupled with reactive oxygen species quenching-enables partial rescuing of HSC function.