PAX8 and endocrine system disorder: As expected, the 81 genes (Table S4) excluded because of divergent expression are enriched for thyroid hormonogenesis and function: thyroglobulin (TG), thyroid peroxidase (TPO), deiodinase type II (DIO2), deiodinase type I (DIO1), dual oxidase 2 (DUOX2), paired box gene 8 (PAX8), thyroid stimulating hormone receptor (TSHR) and thyroid hormone responsive SPOT14 homolog (THRSP) (i.e for GO term endocrine system disorders, the enrichment score is 2; p<0.01 and Benjamini-Hochberg multiple correction p<0.01).