enclone: precision clonotyping and analysis of immune receptors

Abstract

Half a billion years of evolutionary battle forged the vertebrate adaptive immune system, an astonishingly versatile factory for molecules that can adapt to arbitrary attacks. The history of an individual encounter is chronicled within a clonotype: the descendants of a single fully rearranged adaptive immune cell. For B cells, reading this immune history for an individual remains a fundamental challenge of modern immunology. Identification of such clonotypes is a magnificently challenging problem for three reasons: O_LIThe cell history is inferred rather than directly observed: the only available data are the sequences of V(D)J molecules occurring in a sample of cells. C_LIO_LIEach immune receptor is a pair of V(D)J molecules. Identifying these pairs at scale is a technological challenge and cannot be done with perfect accuracy--real samples are mixtures of cells and fragments thereof. C_LIO_LIThese molecules can be intensely mutated during the optimization of the response to particular antigens, blurring distinctions between kindred molecules. C_LI It is thus impossible to determine clonotypes exactly. All solutions to this problem make a trade-off between sensitivity and specificity; useful solutions must address actual artifacts found in real data. We present enclone1, a system for computing approximate clonotypes from single cell data, and demonstrate its use and value with the 10x Genomics Immune Profiling Solution. To test it, we generate data for 1.6 million individual B cells, from four humans, including deliberately enriched memory cells, to tax the algorithm and provide a resource for the community. We analytically determine the specificity of enclones clonotyping algorithm, showing that on this dataset the probability of co-clonotyping two unrelated B cells is around 10-9. We prove that using only heavy chains increases the error rate by two orders of magnitude. enclone comprises a comprehensive toolkit for the analysis and display of immune receptor data. It is ultra-fast, easy to install, has public source code, comes with public data, and is documented at bit.ly/enclone. It has three "flavors" of use: (1) as a command-line tool run from a terminal window, that yields visual output; (2) as a command-line tool that yields parseable output that can be fed to other programs; and (3) as a graphical version (GUI).