Hashimoto’s thyroiditis (HT) is the most common autoimmune disease characterized by lymphocytic infiltration and thyrocyte destruction. Dissection of the interaction between the thyroidal stromal microenvironment and the infiltrating immune cells might lead to a better understanding of HT pathogenesis. Here we show, using single-cell RNA-sequencing, that three thyroidal stromal cell subsets, ACKR1+ endothelial cells and CCL21+ myofibroblasts and CCL21+ fibroblasts, contribute to the thyroidal tissue microenvironment in HT. These cell types occupy distinct histological locations within the thyroid gland. Our experiments suggest that they might facilitate lymphocyte trafficking from the blood to thyroid tissues, and T cell zone CCL21+ fibroblasts may also promote the formation of tertiary lymphoid organs characteristic to HT. Our study also demonstrates the presence of inflammatory macrophages and dendritic cells expressing high levels of IL-1β in the thyroid, which may contribute to thyrocyte destruction in HT patients. Our findings thus provide a deeper insight into the cellular interactions that might prompt the pathogenesis of HT.

Rui-jia Zhang, Liu Yang and Feng Sun are co-first authors. Huai-dong Song and Mei Dong are cocorresponding authors.

Background. BRAF exon 15 p.V600E (BRAF V600E) mutation has been established as an important molecular marker for papillary thyroid carcinoma diagnosis by ultrasound-guided fine-needle aspiration biopsy (FNAB). Sanger sequencing is the gold standard for detecting BRAF V600E mutations but fails to identify low-frequency mutations. However, droplet digital PCR (ddPCR) is a popular new method for detecting low-frequency mutations. Here, we compare the efficiency of droplet digital PCR (ddPCR) and Sanger sequencing for detection of the BRAF V600E mutation in thyroid fine-needle aspiration (FNA) samples. Methods. Thyroid fine-needle aspiration samples from 278 patients with 310 thyroid nodules were collected. Sanger sequencing and ddPCR were conducted to detect the BRAF V600E mutation. Results. The BRAF V600E mutation was found in 94 nodules (30.32%) by ddPCR and 40 nodules (12.90%) by Sanger sequencing in 310 FNA samples. A total of 119 nodules were confirmed PTC by postsurgical pathology. Among which the BRAF mutation was found in 80 (67.23%) nodules by ddPCR and 31 (26.05%) by Sanger sequencing. All nodules carrying the mutation detected by Sanger sequencing (SS+) were verified by ddPCR (ddPCR+). Also, all nodules with no mutation detected by ddPCR were interpreted as wild-type by Sanger sequencing (SS−). In addition. Almost all SS+/ddPCR + nodules (95.00%; 38/40) and SS−/ddPCR + nodules (100.00%; 54/54) displayed a BRAF mutation rate of >5% and <15%, respectively, indicating easy misdetection by Sanger sequencing when the mutation rate is between 5 and 15%. Conclusion. ddPCR has higher sensitivity than Sanger sequencing and we propose ddPCR as a supplement to Sanger sequencing in molecular testing of BRAF using FNAB samples.