Controlling organoid symmetry breaking uncovers an excitable system underlying human axial elongation.

The human embryo breaks symmetry to form the anterior-posterior axis of the body. As the embryo elongates along this axis, progenitors in the tail bud give rise to tissues that generate spinal cord, skeleton, and musculature. This raises the question of how the embryo achieves axial elongation and patterning. While ethics necessitate in vitro studies, the variability of organoid systems has hindered mechanistic insights. Here, we developed a bioengineering and machine learning framework that optimizes organoid symmetry breaking by tuning their spatial coupling. This framework enabled reproducible generation of axially elongating organoids, each possessing a tail bud and neural tube. We discovered that an excitable system composed of WNT/FGF signaling drives elongation by inducing a neuromesodermal progenitor-like signaling center. We discovered that instabilities in the excitable system are suppressed by secreted WNT inhibitors. Absence of these inhibitors led to ectopic tail buds and branches. Our results identify mechanisms governing stable human axial elongation.

Cross-prediction-powered inference.

While reliable data-driven decision-making hinges on high-quality labeled data, the acquisition of quality labels often involves laborious human annotations or slow and expensive scientific measurements. Machine learning is becoming an appealing alternative as sophisticated predictive techniques are being used to quickly and cheaply produce large amounts of predicted labels; e.g., predicted protein structures are used to supplement experimentally derived structures, predictions of socioeconomic indicators from satellite imagery are used to supplement accurate survey data, and so on. Since predictions are imperfect and potentially biased, this practice brings into question the validity of downstream inferences. We introduce cross-prediction: a method for valid inference powered by machine learning. With a small labeled dataset and a large unlabeled dataset, cross-prediction imputes the missing labels via machine learning and applies a form of debiasing to remedy the prediction inaccuracies. The resulting inferences achieve the desired error probability and are more powerful than those that only leverage the labeled data. Closely related is the recent proposal of prediction-powered inference [A. N. Angelopoulos, S. Bates, C. Fannjiang, M. I. Jordan, T. Zrnic, Science 382, 669-674 (2023)], which assumes that a good pretrained model is already available. We show that cross-prediction is consistently more powerful than an adaptation of prediction-powered inference in which a fraction of the labeled data is split off and used to train the model. Finally, we observe that cross-prediction gives more stable conclusions than its competitors; its CIs typically have significantly lower variability.

Disentangling the causes of mumps reemergence in the United States.

Over the past two decades, multiple countries with high vaccine coverage have experienced resurgent outbreaks of mumps. Worryingly, in these countries, a high proportion of cases have been among those who have completed the recommended vaccination schedule, raising alarm about the effectiveness of existing vaccines. Two putative mechanisms of vaccine failure have been proposed as driving observed trends: 1) gradual waning of vaccine-derived immunity (necessitating additional booster doses) and 2) the introduction of novel viral genotypes capable of evading vaccinal immunity. Focusing on the United States, we conduct statistical likelihood-based hypothesis testing using a mechanistic transmission model on age-structured epidemiological, demographic, and vaccine uptake time series data. We find that the data are most consistent with the waning hypothesis and estimate that 32.8% (32%, 33.5%) of individuals lose vaccine-derived immunity by age 18 y. Furthermore, we show using our transmission model how waning vaccine immunity reproduces qualitative and quantitatively consistent features of epidemiological data, namely 1) the shift in mumps incidence toward older individuals, 2) the recent recurrence of mumps outbreaks, and 3) the high proportion of mumps cases among previously vaccinated individuals.

An inference model gives insights into innate immune adaptation and repertoire diversity.

The innate immune system is the body's first line of defense against infection. Natural killer (NK) cells, a vital part of the innate immune system, help to control infection and eliminate cancer. Studies have identified a vast array of receptors that NK cells use to discriminate between healthy and unhealthy cells. However, at present, it is difficult to explain how NK cells will respond to novel stimuli in different environments. In addition, the expression of different receptors on individual NK cells is highly stochastic, but the reason for these variegated expression patterns is unclear. Here, we studied the recognition of unhealthy target cells as an inference problem, where NK cells must distinguish between healthy targets with normal variability in ligand expression and ones that are clear "outliers." Our mathematical model fits well with experimental data, including NK cells' adaptation to changing environments and responses to different target cells. Furthermore, we find that stochastic, "sparse" receptor expression profiles are best able to detect a variety of possible threats, in agreement with experimental studies of the NK cell repertoire. While our study was specifically motivated by NK cells, our model is general and could also apply more broadly to explain principles of target recognition for other immune cell types.

Optimal inference of molecular interaction dynamics in FRET microscopy.

Intensity-based time-lapse fluorescence resonance energy transfer (FRET) microscopy has been a major tool for investigating cellular processes, converting otherwise unobservable molecular interactions into fluorescence time series. However, inferring the molecular interaction dynamics from the observables remains a challenging inverse problem, particularly when measurement noise and photobleaching are nonnegligible-a common situation in single-cell analysis. The conventional approach is to process the time-series data algebraically, but such methods inevitably accumulate the measurement noise and reduce the signal-to-noise ratio (SNR), limiting the scope of FRET microscopy. Here, we introduce an alternative probabilistic approach, B-FRET, generally applicable to standard 3-cube FRET-imaging data. Based on Bayesian filtering theory, B-FRET implements a statistically optimal way to infer molecular interactions and thus drastically improves the SNR. We validate B-FRET using simulated data and then apply it to real data, including the notoriously noisy in vivo FRET time series from individual bacterial cells to reveal signaling dynamics otherwise hidden in the noise.

Inferring the T cell repertoire dynamics of healthy individuals.

The adaptive immune system is a diverse ecosystem that responds to pathogens by selecting cells with specific receptors. While clonal expansion in response to particular immune challenges has been extensively studied, we do not know the neutral dynamics that drive the immune system in the absence of strong stimuli. Here, we learn the parameters that underlie the clonal dynamics of the T cell repertoire in healthy individuals of different ages, by applying Bayesian inference to longitudinal immune repertoire sequencing (RepSeq) data. Quantifying the experimental noise accurately for a given RepSeq technique allows us to disentangle real changes in clonal frequencies from noise. We find that the data are consistent with clone sizes following a geometric Brownian motion and show that its predicted steady state is in quantitative agreement with the observed power-law behavior of the clone-size distribution. The inferred turnover time scale of the repertoire increases with patient age and depends on the clone size in some individuals.

Phylogenetic inference of where species spread or split across barriers.

SignificanceGeography molds how species evolve in space. Strong geographical barriers to movement, for instance, both inhibit dispersal between regions and allow isolated populations to diverge as new species. Weak barriers, by contrast, permit species range expansion and persistence. These factors present a conundrum: How strong must a barrier be before between-region speciation outpaces dispersal? We designed a phylogenetic model of dispersal, extinction, and speciation that allows regional features to influence rates of biogeographic change and applied it to the neotropical radiation of Anolis lizards. Separation by water induces a threefold steeper barrier to movement than equivalent distances over land. Our model will help biologists detect relationships between evolutionary processes and the spatial contexts in which they operate.

The value of genetic data from 665,460 individuals in managing iron deficiency anaemia and suitability to donate blood.

BACKGROUND AND OBJECTIVES: Although the genetic determinants of haemoglobin and ferritin have been widely studied, those of the clinically and globally relevant iron deficiency anaemia (IDA) and deferral due to hypohaemoglobinemia (Hb-deferral) are unclear. In this investigation, we aimed to quantify the value of genetic information in predicting IDA and Hb-deferral. MATERIALS AND METHODS: We analysed genetic data from up to 665,460 participants of the FinnGen, Blood Service Biobank and UK Biobank, and used INTERVAL (N = 39,979) for validation. We performed genome-wide association studies (GWASs) of IDA and Hb-deferral and utilized publicly available genetic associations to compute polygenic scores for IDA, ferritin and Hb. We fitted models to estimate the effect sizes of these polygenic risk scores (PRSs) on IDA and Hb-deferral risk while accounting for the individual's age, sex, weight, height, smoking status and blood donation history. RESULTS: Significant variants in GWASs of IDA and Hb-deferral appear to be a small subset of variants associated with ferritin and Hb. Effect sizes of genetic predictors of IDA and Hb-deferral are similar to those of age and weight which are typically used in blood donor management. A total genetic score for Hb-deferral was estimated for each individual. The odds ratio estimate between first decile against that at ninth decile of total genetic score distribution ranged from 1.4 to 2.2. CONCLUSION: The value of genetic data in predicting IDA or suitability to donate blood appears to be on a practically useful level.

Robust inference methods for meta-analysis involving influential outlying studies.

Meta-analysis is an essential tool to comprehensively synthesize and quantitatively evaluate results of multiple clinical studies in evidence-based medicine. In many meta-analyses, the characteristics of some studies might markedly differ from those of the others, and these outlying studies can generate biases and potentially yield misleading results. In this article, we provide effective robust statistical inference methods using generalized likelihoods based on the density power divergence. The robust inference methods are designed to adjust the influences of outliers through the use of modified estimating equations based on a robust criterion, even when multiple and serious influential outliers are present. We provide the robust estimators, statistical tests, and confidence intervals via the generalized likelihoods for the fixed-effect and random-effects models of meta-analysis. We also assess the contribution rates of individual studies to the robust overall estimators that indicate how the influences of outlying studies are adjusted. Through simulations and applications to two recently published systematic reviews, we demonstrate that the overall conclusions and interpretations of meta-analyses can be markedly changed if the robust inference methods are applied and that only the conventional inference methods might produce misleading evidence. These methods would be recommended to be used at least as a sensitivity analysis method in the practice of meta-analysis. We have also developed an R package, robustmeta, that implements the robust inference methods.

Confidence sets for a level set in linear regression.

Regression modeling is the workhorse of statistics and there is a vast literature on estimation of the regression function. It has been realized in recent years that in regression analysis the ultimate aim may be the estimation of a level set of the regression function, ie, the set of covariate values for which the regression function exceeds a predefined level, instead of the estimation of the regression function itself. The published work on estimation of the level set has thus far focused mainly on nonparametric regression, especially on point estimation. In this article, the construction of confidence sets for the level set of linear regression is considered. In particular, 1 - α $$ 1-\alpha $$ level upper, lower and two-sided confidence sets are constructed for the normal-error linear regression. It is shown that these confidence sets can be easily constructed from the corresponding 1 - α $$ 1-\alpha $$ level simultaneous confidence bands. It is also pointed out that the construction method is readily applicable to other parametric regression models where the mean response depends on a linear predictor through a monotonic link function, which include generalized linear models, linear mixed models and generalized linear mixed models. Therefore, the method proposed in this article is widely applicable. Simulation studies with both linear and generalized linear models are conducted to assess the method and real examples are used to illustrate the method.

Extracorporeal cardiopulmonary resuscitation versus standard treatment for refractory out-of-hospital cardiac arrest: a Bayesian meta-analysis.

BACKGROUND: The outcomes of several randomized trials on extracorporeal cardiopulmonary resuscitation (ECPR) in patients with refractory out-of-hospital cardiac arrest were examined using frequentist methods, resulting in a dichotomous interpretation of results based on p-values rather than in the probability of clinically relevant treatment effects. To determine such a probability of a clinically relevant ECPR-based treatment effect on neurological outcomes, the authors of these trials performed a Bayesian meta-analysis of the totality of randomized ECPR evidence. METHODS: A systematic search was applied to three electronic databases. Randomized trials that compared ECPR-based treatment with conventional CPR for refractory out-of-hospital cardiac arrest were included. The study was preregistered in INPLASY (INPLASY2023120060). The primary Bayesian hierarchical meta-analysis estimated the difference in 6-month neurologically favorable survival in patients with all rhythms, and a secondary analysis assessed this difference in patients with shockable rhythms (Bayesian hierarchical random-effects model). Primary Bayesian analyses were performed under vague priors. Outcomes were formulated as estimated median relative risks, mean absolute risk differences, and numbers needed to treat with corresponding 95% credible intervals (CrIs). The posterior probabilities of various clinically relevant absolute risk difference thresholds were estimated. RESULTS: Three randomized trials were included in the analysis (ECPR, n = 209 patients; conventional CPR, n = 211 patients). The estimated median relative risk of ECPR for 6-month neurologically favorable survival was 1.47 (95%CrI 0.73-3.32) with a mean absolute risk difference of 8.7% (- 5.0; 42.7%) in patients with all rhythms, and the median relative risk was 1.54 (95%CrI 0.79-3.71) with a mean absolute risk difference of 10.8% (95%CrI - 4.2; 73.9%) in patients with shockable rhythms. The posterior probabilities of an absolute risk difference > 0% and > 5% were 91.0% and 71.1% in patients with all rhythms and 92.4% and 75.8% in patients with shockable rhythms, respectively. CONCLUSION: The current Bayesian meta-analysis found a 71.1% and 75.8% posterior probability of a clinically relevant ECPR-based treatment effect on 6-month neurologically favorable survival in patients with all rhythms and shockable rhythms. These results must be interpreted within the context of the reported credible intervals and varying designs of the randomized trials. REGISTRATION: INPLASY (INPLASY2023120060, December 14th, 2023, https://doi.org/10.37766/inplasy2023.12.0060 ).

Valence and interactions in judicial voting.

The collective statistics of voting on judicial courts present hints about their inner workings. Many approaches for studying these statistics, however, assume that judges' decisions are conditionally independent: a judge reaches a decision based on the case at hand and his or her personal views. In reality, judges interact. We develop a minimal model that accounts for judge bias, depending on the context of the case, and peer interaction. We apply the model to voting data from the US Supreme Court. We find strong evidence that interaction is an important factor across natural courts from 1946 to 2021. We also find that, after accounting for interaction, the recovered biases differ from highly cited ideological scores. Our method exemplifies how physics and complexity-inspired modelling can drive the development of theoretical models and improved measures for political voting. This article is part of the theme issue 'A complexity science approach to law and governance'.

Anomalous influenza seasonality in the United States and the emergence of novel influenza B viruses.

The 2019/2020 influenza season in the United States began earlier than any season since the 2009 H1N1 pandemic, with an increase in influenza-like illnesses observed as early as August. Also noteworthy was the numerical domination of influenza B cases early in this influenza season, in contrast to their typically later peak in the past. Here, we dissect the 2019/2020 influenza season not only with regard to its unusually early activity, but also with regard to the relative dynamics of type A and type B cases. We propose that the recent expansion of a novel influenza B/Victoria clade may be associated with this shift in the composition and kinetics of the influenza season in the United States. We use epidemiological transmission models to explore whether changes in the effective reproduction number or short-term cross-immunity between these viruses can explain the dynamics of influenza A and B seasonality. We find support for an increase in the effective reproduction number of influenza B, rather than support for cross-type immunity-driven dynamics. Our findings have clear implications for optimal vaccination strategies.

Deep generative selection models of T and B cell receptor repertoires with soNNia.

Subclasses of lymphocytes carry different functional roles to work together and produce an immune response and lasting immunity. Additionally to these functional roles, T and B cell lymphocytes rely on the diversity of their receptor chains to recognize different pathogens. The lymphocyte subclasses emerge from common ancestors generated with the same diversity of receptors during selection processes. Here, we leverage biophysical models of receptor generation with machine learning models of selection to identify specific sequence features characteristic of functional lymphocyte repertoires and subrepertoires. Specifically, using only repertoire-level sequence information, we classify CD4+ and CD8+ T cells, find correlations between receptor chains arising during selection, and identify T cell subsets that are targets of pathogenic epitopes. We also show examples of when simple linear classifiers do as well as more complex machine learning methods.

Emergence of hierarchy in networked endorsement dynamics.

Many social and biological systems are characterized by enduring hierarchies, including those organized around prestige in academia, dominance in animal groups, and desirability in online dating. Despite their ubiquity, the general mechanisms that explain the creation and endurance of such hierarchies are not well understood. We introduce a generative model for the dynamics of hierarchies using time-varying networks, in which new links are formed based on the preferences of nodes in the current network and old links are forgotten over time. The model produces a range of hierarchical structures, ranging from egalitarianism to bistable hierarchies, and we derive critical points that separate these regimes in the limit of long system memory. Importantly, our model supports statistical inference, allowing for a principled comparison of generative mechanisms using data. We apply the model to study hierarchical structures in empirical data on hiring patterns among mathematicians, dominance relations among parakeets, and friendships among members of a fraternity, observing several persistent patterns as well as interpretable differences in the generative mechanisms favored by each. Our work contributes to the growing literature on statistically grounded models of time-varying networks.

Multimodel Approaches Are Not the Best Way to Understand Multifactorial Systems.

Information-theoretic (IT) and multi-model averaging (MMA) statistical approaches are widely used but suboptimal tools for pursuing a multifactorial approach (also known as the method of multiple working hypotheses) in ecology. (1) Conceptually, IT encourages ecologists to perform tests on sets of artificially simplified models. (2) MMA improves on IT model selection by implementing a simple form of shrinkage estimation (a way to make accurate predictions from a model with many parameters relative to the amount of data, by "shrinking" parameter estimates toward zero). However, other shrinkage estimators such as penalized regression or Bayesian hierarchical models with regularizing priors are more computationally efficient and better supported theoretically. (3) In general, the procedures for extracting confidence intervals from MMA are overconfident, providing overly narrow intervals. If researchers want to use limited data sets to accurately estimate the strength of multiple competing ecological processes along with reliable confidence intervals, the current best approach is to use full (maximal) statistical models (possibly with Bayesian priors) after making principled, a priori decisions about model complexity.

Research on Active Safety Situation of Road Passenger Transportation Enterprises: Evaluation, Prediction, and Analysis.

The road passenger transportation enterprise is a complex system, requiring a clear understanding of their active safety situation (ASS), trends, and influencing factors. This facilitates transportation authorities to promptly receive signals and take effective measures. Through exploratory factor analysis and confirmatory factor analysis, we delved into potential factors for evaluating ASS and extracted an ASS index. To predict obtaining a higher ASS information rate, we compared multiple time series models, including GRU (gated recurrent unit), LSTM (long short-term memory), ARIMA, Prophet, Conv_LSTM, and TCN (temporal convolutional network). This paper proposed the WDA-DBN (water drop algorithm-Deep Belief Network) model and employed DEEPSHAP to identify factors with higher ASS information content. TCN and GRU performed well in the prediction. Compared to the other models, WDA-DBN exhibited the best performance in terms of MSE and MAE. Overall, deep learning models outperform econometric models in terms of information processing. The total time spent processing alarms positively influences ASS, while variables such as fatigue driving occurrences, abnormal driving occurrences, and nighttime driving alarm occurrences have a negative impact on ASS.

MSstats Version 4.0: Statistical Analyses of Quantitative Mass Spectrometry-Based Proteomic Experiments with Chromatography-Based Quantification at Scale.

The MSstats R-Bioconductor family of packages is widely used for statistical analyses of quantitative bottom-up mass spectrometry-based proteomic experiments to detect differentially abundant proteins. It is applicable to a variety of experimental designs and data acquisition strategies and is compatible with many data processing tools used to identify and quantify spectral features. In the face of ever-increasing complexities of experiments and data processing strategies, the core package of the family, with the same name MSstats, has undergone a series of substantial updates. Its new version MSstats v4.0 improves the usability, versatility, and accuracy of statistical methodology, and the usage of computational resources. New converters integrate the output of upstream processing tools directly with MSstats, requiring less manual work by the user. The package's statistical models have been updated to a more robust workflow. Finally, MSstats' code has been substantially refactored to improve memory use and computation speed. Here we detail these updates, highlighting methodological differences between the new and old versions. An empirical comparison of MSstats v4.0 to its previous implementations, as well as to the packages MSqRob and DEqMS, on controlled mixtures and biological experiments demonstrated a stronger performance and better usability of MSstats v4.0 as compared to existing methods.

Neural correlates of valence-dependent belief and value updating during uncertainty reduction: An fNIRS study.

Selective use of new information is crucial for adaptive decision-making. Combining a gamble bidding task with assessing cortical responses using functional near-infrared spectroscopy (fNIRS), we investigated potential effects of information valence on behavioral and neural processes of belief and value updating during uncertainty reduction in young adults. By modeling changes in the participants' expressed subjective values using a Bayesian model, we dissociated processes of (i) updating beliefs about statistical properties of the gamble, (ii) updating values of a gamble based on new information about its winning probabilities, as well as (iii) expectancy violation. The results showed that participants used new information to update their beliefs and values about the gambles in a quasi-optimal manner, as reflected in the selective updating only in situations with reducible uncertainty. Furthermore, their updating was valence-dependent: information indicating an increase in winning probability was underweighted, whereas information about a decrease in winning probability was updated in good agreement with predictions of the Bayesian decision theory. Results of model-based and moderation analyses showed that this valence-dependent asymmetry was associated with a distinct contribution of expectancy violation, besides belief updating, to value updating after experiencing new positive information regarding winning probabilities. In line with the behavioral results, we replicated previous findings showing involvements of frontoparietal brain regions in the different components of updating. Furthermore, this study provided novel results suggesting a valence-dependent recruitment of brain regions. Individuals with stronger oxyhemoglobin responses during value updating was more in line with predictions of the Bayesian model while integrating new information that indicates an increase in winning probability. Taken together, this study provides first results showing expectancy violation as a contributing factor to sub-optimal valence-dependent updating during uncertainty reduction and suggests limitations of normative Bayesian decision theory.

Dissecting Genomic Determinants of Positive Selection with an Evolution-Guided Regression Model.

In evolutionary genomics, it is fundamentally important to understand how characteristics of genomic sequences, such as gene expression level, determine the rate of adaptive evolution. While numerous statistical methods, such as the McDonald-Kreitman (MK) test, are available to examine the association between genomic features and the rate of adaptation, we currently lack a statistical approach to disentangle the independent effect of a genomic feature from the effects of other correlated genomic features. To address this problem, I present a novel statistical model, the MK regression, which augments the MK test with a generalized linear model. Analogous to the classical multiple regression model, the MK regression can analyze multiple genomic features simultaneously to infer the independent effect of a genomic feature, holding constant all other genomic features. Using the MK regression, I identify numerous genomic features driving positive selection in chimpanzees. These features include well-known ones, such as local mutation rate, residue exposure level, tissue specificity, and immune genes, as well as new features not previously reported, such as gene expression level and metabolic genes. In particular, I show that highly expressed genes may have a higher adaptation rate than their weakly expressed counterparts, even though a higher expression level may impose stronger negative selection. Also, I show that metabolic genes may have a higher adaptation rate than their nonmetabolic counterparts, possibly due to recent changes in diet in primate evolution. Overall, the MK regression is a powerful approach to elucidate the genomic basis of adaptation.