Comprehensive Analysis of circRNA Expression Profiles in Human Brown Adipose Tissue

Figures & data

Figure 1 H&E staining and UCP1 expression in BAT and WAT. (A) BAT and WAT were stained with H&E, scale bar 200μm. (B) UCP1 mRNA expression of BAT and WAT. n=3, **p<0.01.

Figure 2 DECs in human BAT and WAT. (A) The box plot displayed the normalized intensity distribution of all data. (B) Hierarchical clustering maps of DECs were shown in BAT and WAT. Red represents relatively high expression; green represents relatively low expression. (C) Scatter plots were used to assess changes in circRNAs expression between BAT and WAT. CircRNAs above the top green line and below the bottom green line indicate more than a 1.5-fold variation in circRNAs between the two comparison samples. (D) The volcano plot showed statistically significant DECs between BAT and WAT (fold change >2.0, P<0.05). (E) Classification and distribution of DECs in human chromosomes.

Table 1 Upregulated circRNAs in Human BAT Vs WAT

circRNA	P-value	FDR	FC	Best Transcript	Gene Symbol
hsa_circ_0006168	0.20031454	0.484721846	11.5191529	NM_144571	CNOT6L
hsa_circ_0004789	0.101221236	0.476643198	4.1624141	NM_022739	SMURF2
hsa_circ_0026337	0.145820553	0.476643198	3.922333	NM_014191	SCN8A
hsa_circ_0007507	0.014112967	0.460528239	3.7781622	NM_002890	RASA1
hsa_circ_0000542	0.04078890653	0.476643198	2.0989004	NM_002892	ARID4A
hsa_circ_0005606	0.037371486	0.476643198	1.9468111	NM_020825	CRAMP1L
hsa_circ_0005916	0.004717319	0.412506315	1.90146	NM_153348	FBXW8
hsa_circ_0008417	0.039665096	0.476643198	1.7857437	NM_024717	MCTP1
hsa_circ_0072568	0.01104448507	0.450494814	1.6449448	NM_001197220	PDE4D
hsa_circ_0022382	0.03658574027	0.476643198	1.6236087	NM_004265	FADS2

Abbreviations: FDR, false discovery rate; FC, fold change.

Table 2 Downregulated circRNAs in Human BAT Vs WAT

circRNA	P-value	FDR	FC	Best Transcript	Gene Symbol
hsa_circ_0030162	0.175399345	0.476797797	10.4207494	NM_003295	TPT1
hsa_circ_0059665	0.122011263	0.476643198	9.9068558	NM_001042472	ABHD12
hsa_circ_0081881	0.045191535	0.476643198	5.8005395	NM_002736	PRKAR2B
hsa_circ_0034510	0.008072181	0.442202588	3.9667945	NM_003246	THBS1
hsa_circ_0003416	0.049459158	0.476643198	3.7847447	NM_021109	TMSB4X
hsa_circ_0017693	0.014205109	0.460528239	3.7760767	NM_024693	ECHDC3
hsa_circ_0000102	0.041575465	0.476643198	3.6706236	NM_001048210	CLCC1
hsa_circ_0005778	0.029609067	0.476643198	3.4702494	NM_001746	CANX
hsa_circ_0003183	0.048432575	0.476643198	3.3750592	NM_002865	RAB2A
hsa_circ_0068465	0.012177718	0.454378693	3.3168519	NM_001967	EIF4A2

Abbreviations: FDR, false discovery rate; FC, fold change.

Figure 3 GO and pathway analysis of DECs in BAT and WAT. GO analysis of upregulated (A) and downregulated (B) DECs in BAT and WAT. Pathway analysis of upregulated (C) and downregulated (D) DECs in BAT and WAT.

Figure 4 Validated of DECs by qRT-PCR. The expression of hsa_circ_0006168 (A), hsa_circ_0026337 (B), hsa_circ_000750 (C), hsa_circ_0030162 (D), hsa_circ_0004789 (E) and hsa_circ_0059665 (F) in BAT and WAT. n=3, *p<0.05. Results are expressed as the means.

Figure 5 Prediction of the circRNA–microRNA–mRNA interaction network. The target microRNAs of hsa_circ_0006168 (A), hsa_circ_0026337 (B), hsa_circ_000750 (C), and hsa_circ_0030162 (D) were predicted by Arraystar’s home-made software and target mRNA of microRNA was predicted by DIANA TOOLS software. circRNA–microRNA–mRNA interaction network was visualized and constructed by Cytoscape 2.8.2.