Lymphocytic Thyroiditis Transcriptomic Profiles Support the Role of Checkpoint Pathways and B Cells in Pathogenesis

Autoimmune thyroid diseases are the most common types of autoimmune diseases, but their physiopathology is still relatively unexplored. Genotype-tissue expression (GTEx) is a publicly available repository containing RNAseq data, including profiles from thyroid. Approximately 14.8% of these glands were affected by focal lymphocytic thyroiditis and 6.3% were annotated as Hashimoto. We interrogated these data to improve the characterization of infiltrating cells and to identify new molecular pathways active in autoimmune thyroiditis. Histological GTEx images of 336 thyroid samples were classified into three categories, that is, non-infiltrated thyroid, small focal infiltrated thyroid, and extensive lymphoid infiltrated thyroid. Differentially expressed genes among these categories were identified and subjected to in silico pathway enrichment analysis accordingly. CIBERSORTx deconvolution was used to characterize infiltrating cells. As expected, most of the transcriptional changes were dependent on tissue infiltration. Upregulated genes in tissues include-in addition to lineage-specific B and T cell genes-a broad representation of inhibitory immune checkpoint receptors expressed by B and T lymphocytes. CIBERSORTx analysis identified 22 types of infiltrating cells showed that T cells predominate 3:1 over B cells in glands with small infiltrates, only by 1.7:1 in those with large infiltrates. Follicular helper and memory CD4 T cells were significantly more abundant in glands with large infiltrates (p < 0.0001), but the most prominent finding in these glands was an almost sixfold increase in the number of naive B cells (p < 0.0001). A predominance of M2 macrophages over M1 and M0 macrophages was observed in the three gland categories (p < 0.001). Analysis of transcriptomic RNA-seq profiles constitutes a rich source of information for the analysis of autoimmune tissues. High-resolution transcriptomic data analysis of thyroid glands indicates the following: (a) in all infiltrated glands, active autoimmune response coexists with suppressor counteracting mechanisms involving several inhibitory checkpoint receptor pairs, (b) glands with small infiltrates contain an unexpected relatively high proportion of B lymphocytes, and (c) in highly infiltrated glands, there is a distinct transcriptomic signature of active tertiary lymphoid organs. These results support the concept that the autoimmune response is amplified in the thyroid tissue.