Peripheral mechanisms preventing autoimmunity and maintaining tolerance to commensal microbiota involve CD4+Foxp3+ regulatory T cells1,2 generated in the thymus (tTregs) or extrathymically by induction of naive CD4+Foxp3− T cells (iTregs). Prior studies suggested that the T cell receptor (TCR) repertoires of tTregs and iTregs are biased towards self and non-self antigens, respectively 3–6 but their relative contribution in controlling immunopathology, e.g. colitis and other untoward inflammatory responses triggered by different types of antigens, remains unresolved 7. The intestine, and especially the colon, is a particularly suitable organ to study this question, given the variety of self-, microbiota- and food-derived antigens to which Tregs and other T cell populations are exposed. Intestinal environments can enhance conversion to a regulatory lineage 8,9 and favor tolerogenic presentation of antigens to naive CD4+ T cells 10,11, suggesting that intestinal homeostasis depends on microbiota-specific iTregs 12–15. Here, to identify the origin and antigen-specificity of intestinal Tregs, we performed single cell as well as high-throughput (HT) sequencing of the TCR repertoires of CD4+Foxp3+ and CD4+Foxp3− T cells and analyzed their reactivity against specific commensal species. We show that tTregs constitute the majority of Tregs in all lymphoid and intestinal organs, including colon, where their repertoire is heavily influenced by the composition of the microbiota. Our results suggest that tTregs, and not iTregs, dominantly mediate tolerance to antigens produced by intestinal commensals.

Helios transcription factor and semaphorin receptor Nrp-1 were originally described as constitutively expressed at high levels on CD4+Foxp3+ T regulatory cells of intrathymic origin (tTregs). On the other hand, CD4+Foxp3+ Tregs generated in the periphery (pTregs) or induced ex vivo (iTregs) were reported to express low levels of Helios and Nrp-1. Soon afterwards the reliability of Nrp-1 and Helios as markers discriminating between tTregs and pTregs was questioned and until now no consensus has been reached. Here, we used several genetically modified mouse strains that favor pTregs or tTregs formation and analyzed the TCR repertoire of these cells. We found that Tregs with variable levels of Nrp-1 and Helios were abundant in mice with compromised ability to support natural differentiation of tTregs or pTregs. We also report that TCR repertoires of Treg clones expressing high or low levels of Nrp-1 or Helios are similar and more alike repertoire of CD4+Foxp3+ than repertoire of CD4+Foxp3- thymocytes. These results show that high vs. low expression of Nrp-1 or Helios does not unequivocally identify Treg clones of thymic or peripheral origin.

The paper presents some novel approaches to the empirical analysis of diversity and similarity (overlap) in biological or ecological systems. The analysis is motivated by the molecular studies of highly diverse mammalian T-cell receptor (TCR) populations, and is related to the classical statistical problem of analyzing two-way contingency tables with missing cells and low cell counts. The new measures of diversity and overlap are proposed, based on the information-theoretic as well as geometric considerations, with the capacity to naturally up-weight or down-weight the rare and abundant population species. The consistent estimates are derived by applying the Good-Turing sample-coverage correction. In particular, novel consistent estimates of the Shannon entropy function and the Morisita-Horn index are provided. Data from TCR populations in mice are used to illustrate the empirical performance of the proposed methods vis a vis the existing alternatives.

BackgroundRibosomal deficits are documented in mild cognitive impairment (MCI), which often represents an early stage Alzheimer's disease (AD), as well as in advanced AD. The nucleolar rRNA genes (rDNA), transcription of which is critical for ribosomal biogenesis, are regulated by epigenetic silencing including promoter CpG methylation.Methodology/Principal FindingsTo assess whether CpG methylation of the rDNA promoter was dysregulated across the AD spectrum, we analyzed brain samples from 10 MCI-, 23 AD-, and, 24 age-matched control individuals using bisulfite mapping. The rDNA promoter became hypermethylated in cerebro-cortical samples from MCI and AD groups. In parietal cortex, the rDNA promoter was hypermethylated more in MCI than in advanced AD. The cytosine methylation of total genomic DNA was similar in AD, MCI, and control samples. Consistent with a notion that hypermethylation-mediated silencing of the nucleolar chromatin stabilizes rDNA loci, preventing their senescence-associated loss, genomic rDNA content was elevated in cerebrocortical samples from MCI and AD groups.Conclusions/SignificanceIn conclusion, rDNA hypermethylation could be a new epigenetic marker of AD. Moreover, silencing of nucleolar chromatin may occur during early stages of AD pathology and play a role in AD-related ribosomal deficits and, ultimately, dementia.

Examining enzyme kinetics is critical for understanding cellular systems and for using enzymes in industry. The Michaelis-Menten equation has been widely used for over a century to estimate the enzyme kinetic parameters from reaction progress curves of substrates, which is known as the progress curve assay. However, this canonical approach works in limited conditions, such as when there is a large excess of substrate over enzyme. Even when this condition is satisfied, the identifiability of parameters is not always guaranteed, and often not verifiable in practice. To overcome such limitations of the canonical approach for the progress curve assay, here we propose a Bayesian approach based on an equation derived with the total quasi-steady-state approximation. In contrast to the canonical approach, estimates obtained with this proposed approach exhibit little bias for any combination of enzyme and substrate concentrations. Importantly, unlike the canonical approach, an optimal experiment to identify parameters with certainty can be easily designed without any prior information. Indeed, with this proposed design, the kinetic parameters of diverse enzymes with disparate catalytic efficiencies, such as chymotrypsin, fumarase, and urease, can be accurately and precisely estimated from a minimal amount of timecourse data. A publicly accessible computational package performing such accurate and efficient Bayesian inference for enzyme kinetics is provided.

In this paper, we show that solutions to ordinary differential equations describing the large-population limits of Markovian stochastic epidemic models can be interpreted as survival or cumulative hazard functions when analysing data on individuals sampled from the population. We refer to the individual-level survival and hazard functions derived from population-level equations as a survival dynamical system (SDS). To illustrate how population-level dynamics imply probability laws for individual-level infection and recovery times that can be used for statistical inference, we show numerical examples based on synthetic data. In these examples, we show that an SDS analysis compares favourably with a complete-data maximum-likelihood analysis. Finally, we use the SDS approach to analyse data from a 2009 influenza A(H1N1) outbreak at Washington State University.

In modern molecular biology one of the standard ways of analyzing a vertebrate immune system is to sequence and compare the counts of specific antigen receptor clones (either immunoglobulins or T-cell receptors) derived from various tissues under different experimental or clinical conditions. The resulting statistical challenges are difficult and do not fit readily into the standard statistical framework of contingency tables primarily due to the serious under-sampling of the receptor populations. This under-sampling is caused, on one hand, by the extreme diversity of antigen receptor repertoires maintained by the immune system and, on the other, by the high cost and labor intensity of the receptor data collection process. In most of the recent immunological literature the differences across antigen receptor populations are examined via non-parametric statistical measures of the species overlap and diversity borrowed from ecological studies. While this approach is robust in a wide range of situations, it seems to provide little insight into the underlying clonal size distribution and the overall mechanism differentiating the receptor populations. As a possible alternative, the current paper presents a parametric method that adjusts for the data under-sampling as well as provides a unifying approach to a simultaneous comparison of multiple receptor groups by means of the modern statistical tools of unsupervised learning. The parametric model is based on a flexible multivariate Poisson-lognormal distribution and is seen to be a natural generalization of the univariate Poisson-lognormal models used in the ecological studies of biodiversity patterns. The procedure for evaluating a model's fit is described along with the public domain software developed to perform the necessary diagnostics. The model-driven analysis is seen to compare favorably vis a vis traditional methods when applied to the data from T-cell receptors in transgenic mice populations.

We consider a Markovian SIR-type (Susceptible → Infected → Recovered) stochastic epidemic process with multiple modes of transmission on a contact network. The network is given by a random graph following a multilayer configuration model where edges in different layers correspond to potentially infectious contacts of different types. We assume that the graph structure evolves in response to the epidemic via activation or deactivation of edges of infectious nodes. We derive a large graph limit theorem that gives a system of ordinary differential equations (ODEs) describing the evolution of quantities of interest, such as the proportions of infected and susceptible vertices, as the number of nodes tends to infinity. Analysis of the limiting system elucidates how the coupling of edge activation and deactivation to infection status affects disease dynamics, as illustrated by a two-layer network example with edge types corresponding to community and healthcare contacts. Our theorem extends some earlier results describing the deterministic limit of stochastic SIR processes on static, single-layer configuration model graphs. We also describe precisely the conditions for equivalence between our limiting ODEs and the systems obtained via pair approximation, which are widely used in the epidemiological and ecological literature to approximate disease dynamics on networks. The flexible modeling framework and asymptotic results have potential application to many disease settings including Ebola dynamics in West Africa, which was the original motivation for this study.