Fig. 2

Multiomics profiling of rs75776403. A Forest plot showing the meta-analysis statistics of rs75776403 to human body height across ethnicities (Taiwanese, European, and Japanese). Fixed-effect and Knapp-Hartung adjusted random-effect models were both applied. B–F Volcano plots showing the differential transcriptomic B proteomic C metabolomic D, and phosphoproteomic (trypsin- E or gluC-digested F analyses to identify molecular features (mRNAs, proteins, metabolites and phosphorylation sites) significantly associated with rs75776403. These features were categorized into up-regulated (red) and down-regulated (pale green) according to estimates (i.e., statistical model coefficient; x-axis). Horizontal dashed lines indicate P < 0.01 (y-axis). G Combined analysis of gene-set enrichment analysis (GSEA) results from expression data (i.e., transcriptome and proteome). Dashed lines indicate a scaled weighted statistic (x- and y-axes) of + 2 or − 2. H Bar plot showing genes (x-axis) that are associated with CLEC18A transcript levels (P < 0.05). The strength of the association was quantified by Spearman’s rho statistic (y-axis). Genes that are involved in cellular response to thyroid hormone stimulus and corticosteroid receptor signaling are colored in purple and orange, respectively. I Venn diagram showing overlapping of CLEC18A transcript-associated genes across three tissues (brain, pituitary gland and thyroid). J Bar plot showing three genes (x-axis) associated with CLEC18A transcript levels in all brain, pituitary gland and thyroid tissues (P < 0.05). The strength of the association was quantified by Spearman’s rho statistic (y-axis). Genes that are implicated in cellular responses to thyroid hormone stimulus and corticosteroid receptor signaling are colored in purple and orange, respectively