Our immune system contains multiple tolerance checkpoints to prevent the activation of self-reactive lymphocytes. How some lymphocytes escape these constraints to cause autoimmune disease remains poorly understood. A long-standing hypothesis in autoimmunity posits that somatic mutations in immune regulatory genes may enable self-reactive lymphocytes to bypass tolerance checkpoints. However, testing this hypothesis has proved challenging due to technical limitations in detecting somatic mutations in polyclonal cell populations. Here, we use deep whole-exome and targeted NanoSeq, a highly accurate single-molecule DNA sequencing protocol, to comprehensively search for driver mutations in autoimmune thyroid disease. This revealed a remarkably high number of B cell clones convergently acquiring loss-of-function somatic mutations in key immune checkpoint genes TNFRSF14 (also known as HVEM) and CD274 (encoding PD-L1), as well as less frequent driver mutations in a diversity of other immune genes. In highly inflamed biopsies, we detected tens to hundreds of independent immune-checkpoint mutant clones. Laser capture microdissection, methylation sequencing, spatial transcriptomics, immunohistochemistry and single-nucleus DNA sequencing localised these mutations to non-naive B cells and revealed frequent co-occurrence of drivers in single cells. We found widespread TNFRSF14 biallelic loss in mutant B cells, and several clones with as many as 4-6 driver mutations. Whilst each clone accounts for a small proportion of cells (typically <1%), the myriad mutant clones in each donor collectively amounted to a substantial fraction of B cells harbouring one or more driver mutations. Our results support the hypothesis that somatic mutations in autoimmune lymphocytes may allow them to escape tolerance constraints through a polyclonal cascade of somatic evolution. These findings provide new insights into the molecular basis of autoimmune disease and suggest novel diagnostic and therapeutic avenues.
