Possible New Candidates Involved to Thalidomide-Related Limbs and Cardiac Defects: A Systems Biology Approach

Thalidomide is a known teratogen that causes malformations especially in heart and limbs. Its mechanism of teratogenicity is still not fully elucidated. Recently, a new target of thalidomide was described, TBX5, and was observed a new interaction between HAND2 and TBX5 that is disrupted in the presence of thalidomide. Therefore, our study aimed to raise potential candidates for thalidomide teratogenesis, through systems biology, evaluating HAND2 and TBX5 interaction and heart and limbs malformations of thalidomide. Genes and proteins related to TBX5 and HAND2 were selected through TF2DNA, REACTOME, Human Phenotype Ontology, and InterPro databases. Networks were assembled using STRING © database. Network analysis were performed in Cytoscape © and R v3.6.2. Differential gene expression (DGE) analysis was performed through gene expression omnibus. We constructed a network for HAND2 and TBX5 interaction; a network for heart and limbs malformations of TE; and the two joined networks. We observed that EP300 protein seemed to be important in all networks. We also looked for proteins containing C2H2 domain in the assembled networks. ZIC3, GLI1, GLI3, ZNF148, and PRDM16 were the ones present in both heart and limbs malformations of TE networks. Furthermore, in the DGE analysis after treatment with thalidomide, we observed that FANCB, ESCO2, and XRCC2 were downregulated and present both in heart and limbs networks. Through systems biology, we were able to point to different new proteins and genes, and selected specially EP300, which was important in all the analyzed networks, to be further evaluated in the TE teratogenicity.