Multianalyte serology in home-sampled blood enables an unbiased assessment of the immune response against SARS-CoV-2.

Serological testing is essential to curb the consequences of the COVID-19 pandemic. However, most assays are still limited to single analytes and samples collected within healthcare. Thus, we establish a multianalyte and multiplexed approach to reliably profile IgG and IgM levels against several versions of SARS-CoV-2 proteins (S, RBD, N) in home-sampled dried blood spots (DBS). We analyse DBS collected during spring of 2020 from 878 random and undiagnosed individuals from the population in Stockholm, Sweden, and use classification approaches to estimate an accumulated seroprevalence of 12.5% (95% CI: 10.3%-14.7%). This includes 5.4% of the samples being IgG+IgM+ against several SARS-CoV-2 proteins, as well as 2.1% being IgG-IgM+ and 5.0% being IgG+IgM- for the virus' S protein. Subjects classified as IgG+ for several SARS-CoV-2 proteins report influenza-like symptoms more frequently than those being IgG+ for only the S protein (OR = 6.1; p < 0.001). Among all seropositive cases, 30% are asymptomatic. Our strategy enables an accurate individual-level and multiplexed assessment of antibodies in home-sampled blood, assisting our understanding about the undiagnosed seroprevalence and diversity of the immune response against the coronavirus.

FUBAR: a fast, unconstrained bayesian approximation for inferring selection.

Model-based analyses of natural selection often categorize sites into a relatively small number of site classes. Forcing each site to belong to one of these classes places unrealistic constraints on the distribution of selection parameters, which can result in misleading inference due to model misspecification. We present an approximate hierarchical Bayesian method using a Markov chain Monte Carlo (MCMC) routine that ensures robustness against model misspecification by averaging over a large number of predefined site classes. This leaves the distribution of selection parameters essentially unconstrained, and also allows sites experiencing positive and purifying selection to be identified orders of magnitude faster than by existing methods. We demonstrate that popular random effects likelihood methods can produce misleading results when sites assigned to the same site class experience different levels of positive or purifying selection--an unavoidable scenario when using a small number of site classes. Our Fast Unconstrained Bayesian AppRoximation (FUBAR) is unaffected by this problem, while achieving higher power than existing unconstrained (fixed effects likelihood) methods. The speed advantage of FUBAR allows us to analyze larger data sets than other methods: We illustrate this on a large influenza hemagglutinin data set (3,142 sequences). FUBAR is available as a batch file within the latest HyPhy distribution (http://www.hyphy.org), as well as on the Datamonkey web server (http://www.datamonkey.org/).

A random effects branch-site model for detecting episodic diversifying selection.

Adaptive evolution frequently occurs in episodic bursts, localized to a few sites in a gene, and to a small number of lineages in a phylogenetic tree. A popular class of "branch-site" evolutionary models provides a statistical framework to search for evidence of such episodic selection. For computational tractability, current branch-site models unrealistically assume that all branches in the tree can be partitioned a priori into two rigid classes--"foreground" branches that are allowed to undergo diversifying selective bursts and "background" branches that are negatively selected or neutral. We demonstrate that this assumption leads to unacceptably high rates of false positives or false negatives when the evolutionary process along background branches strongly deviates from modeling assumptions. To address this problem, we extend Felsenstein's pruning algorithm to allow efficient likelihood computations for models in which variation over branches (and not just sites) is described in the random effects likelihood framework. This enables us to model the process at every branch-site combination as a mixture of three Markov substitution models--our model treats the selective class of every branch at a particular site as an unobserved state that is chosen independently of that at any other branch. When benchmarked on a previously published set of simulated sequences, our method consistently matched or outperformed existing branch-site tests in terms of power and error rates. Using three empirical data sets, previously analyzed for episodic selection, we discuss how modeling assumptions can influence inference in practical situations.