- Open Access
Identification of target genes in cardiomyopathy with fibrosis and cardiac remodeling
Journal of Biomedical Science volume 25, Article number: 63 (2018)
Identify genes probably associated with chronic heart failure and predict potential target genes for dilated cardiomyopathy using bioinformatics analyses.
Gene expression profiles (series number GSE3585 and GSE42955) of cardiomyopathy patients and healthy controls were downloaded from the Expression Omnibus Gene (GEO) database. Differential expression of genes (DEGS) between the two groups of total 14 cardiomyopathy patients and 10 healthy controls were subsequently identified by limma package of R. Database for Annotation, Visualization, and Integrated Discovery (DAVID Tool), which is an analysis of enriched biological processes. Search Tool for the Retrieval Interacting Genes (STRING) was used as well for the analysis of protein-protein interaction network (PPI). Prediction of the potential drugs was suggested based on the preliminarily identified genes using Connectivity Map (CMap).
Eighty-nine DEGs were identified (57 up-regulated and 32 down-regulated). The most enrichment Gene Ontology (GO) terms (P < 0.05) contain genes involved in extracellular matrix (ECM) and biological adhesion signal pathways (P < 0.05, ES > 1.5) such as ECM-receptors, focal adhesion and transforming growth factor beta (TGF-β), etc. Fifty-one differentially expressed genes were found to encode interacting proteins. Eleven key genes along with related transcription factors were identified including CTGF, POSTN, CORIN, FIGF, etc.
Bioinformatics-based analyses reveal the targeted genes probably associated with cardiomyopathy, which provide clues for pharmacological therapies aiming at the targets.
Patients with dilated cardiomyopathy (DCM) often presents with progressive congested heart failure, arrhythmia and thromboembolic disease in forms of left ventricular mural thrombus and/or strike . DCM is seen as a major cause of heart failure apart from coronary heart disease and hypertension . Heart failure is a progressive chronic disease with a 5-year-suvival rate less than 50% . The histologically confirmed diagnosis indicated that prevalence of DCM is 14 to 52% among population proven previous history of myocarditis . High morbidity and mortality underscore the necessity of deeper investigation of the underlying mechanism responsible for development of heart failure in DCM . B-type natriuretic peptides (BNP) is a commonly used biomarker so far for diagnosis of DCM. However, the biomarkers are not so specific since the increased levels of the biomarkers sometimes indicate variety of cardiovascular diseases caused by heterogenetic etiologies, and cannot be explained by impaired left ventricular function alone . The major goal for the treatment of DCM is to reduce the mortality and morbidity rate, to relief symptoms, and to prevent or, in some extent, reverse ventricular remodeling . Fett et al. reported the necessity of application of polymerase chain reaction (PCR) prior to the immunosuppressive therapy . Hamshere at el provided a therapy by administration of granulocyte colony-stimulating factor (G-CSF) with intracoronary autologous bone marrow-derived cells (BMCs) to improve left ventricular systolic function in patients with DCM . Beta-blocker and pace maker inhibitor ivabradine has been proved to have positive effects in reversing cardiac remodeling. However, reactions to beta-blocker or ivabradine vary based on the cases. In some individuals the left ventricular ejection fraction is improved significantly as a result of reversed of attenuated remodeling while in the others it remains . More studies are needed to focus on treatment that improves the outcome of patients with DCM, to precisely make the diagnosis of the disease based on screening of biomarkers. These studies can improve prognosis of patients by lowering the risk of development of heart failure and relevant complications. So it is crucial to understand the mechanism and find biomarkers with a good specificity and sensitivity.
Gene chip microarray database has an open access. Gene chip technique is a widely used approach in analyzing gene function in post-genome era . It is a High-throughput sequencing technique with optimal specificity and sensitivity. Previous studies have partially demonstrated the mechanism underlying DCM by this approach [12, 13]. In this study, two platforms (GPL96 and GPL6244) are incorporated to analyze differential expression genes, enrichment of GO terms or pathways and protein-protein interactions in DCM and predict potential targets and drugs for a better treatment of the disease.
Microarray gene expression
Gene expression profiles with series number GSE3585 based on platform GPL96 and series number GSE42955 on platform GPL6244 were downloaded from the Expression Omnibus Gene (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). Data of seven cardiomyopathy patients and seven healthy controls were randomly selected from each platform, respectively. The total sample of the two patient groups is 14 and the healthy controls are 10 subjects. The samples from all selected subjects had been hybridized on the Affymetrix Human Genome U133A array on a GPL96/GPL6244 platform (Affymetrix, Santa Clara, CA, USA).
Identification of Differential Expressed Genes (DEG)
Limma of R/bio-conductor was used for screening of the DEGs (settings: P < 0.05, log2 (Fold Change)>/=1). Fifty-seven up-regulated (gene set A) and 32 down -regulated (gene set B) genes were identified. Hierarchical clustering and visualization were made by Heat-map package of R.
Enrichment analysis of significant modules
The Database for Annotation, Visualization and Integrated Discovery (DAVID)  provides a non-line comprehensive set of functional annotation tools for biological interpretation of large gene lists. DAVID was used here to group the functions of DEGs in modules, identify enriched biological processes and cellular components, and identify the pathways associated with the DEGs in the most significant modules. Function and pathway terms were retrieved from the Gene Ontology (p < 0.05, Benjamin< 0.01) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (p < 0.05, Benjamin< 0.01) databases, respectively.
Analysis of protein-protein interaction network
The topological properties of the PPI networks were analyzed using Network Analyzer available in Cytoscape. Search Tool for the Retrieval Interacting Genes 10.0 (http://string-db.org/)  provide online analysis of interactions among DEG-encoding proteins. The network of interactions was then imported into Cytoscape , degree ≥5 was set to filter crucial proteins in the middle of the network.
Analysis of transcription factors
PASTAA  was used for predictive analysis of transcription factors of DCMs, P value calculated from hyper geometric distribution was used to evaluate the correlation between the DCMs and transcription factor. TRAP  was used to predict the correlation. Geneset C,D was uploaded to the database and JASPAR, (version 2016)  was used to predict the DNA binding site.
Acquisition of target genes of drugs
Connectivity map  is an online implement that provides gene transcription-expression profile of thousands of genomes reflecting how cultured mammalian cell react to administration of 1309 kinds of bioactive modules in terms of gene-expression. CMap was used to identify potential unknown effect of existing drugs. Up-regulated gene set (gene set A) and down-regulated gene set (gene set B) were uploaded to CMap, EBAYES was used to calculate P-value, PEBAYES-value< 0.05 was considered significant, and could be a potential molecule in treatment of DCM. The half maximal inhibitory concentration (IC50) is a measure of the effectiveness of a substance in inhibiting a specific biological or biochemical function. In previous studies effectiveness of different drugs were demonstrated and published. Thus, results obtained by CMap were uploaded to NCBI PubChem database for verification.
Identification of DCGs
Compared to the control group, 89 differential expression genes (DEG) are identified in the DCM group (57 up-regulated and 32 down-regulated), Fig. 1 is the hierarchical clustering heat-map.
Functional enrichment analysis
In GO functional and enrichment analysis, the most enriched GO terms associated with up-regulated genes are extracellular matrix concerning CC (P = 3.810− 11), extracellular region (P = 5 × 10-14) and biological adhesion (P = 1.5 × 10-6), the most enriched GO terms associated with down-regulated genes are response to injury concerning MF (P = 1.5 × 10-6) and extracellular region (P = 5 × 10-14) concerning CC.
The results of functional annotation about KEGG pathways are shown in Fig. 2. COL1A1, COL1A2, COMP, THBS2, and THBS4, RELN participate in the pathways associated with ECM-receptor interaction and focal adhesion. COMP, THBS2 and THBS4 are significantly enriched in all three pathways.
Topological analysis of PPI networks
Fifty-one significantly enriched gene were updated to STRING to construct the PPI network, and the PPI network was subsequently imported in Cytoscape to construct sub-networks. In sub-network 1, CTGF, COL1A1 (proteins with the highest degree of 8) and protein adjacent to them (COL1A2, THBS2, THBS4, POSTN, COMP, FIGF, PRELP) are highlighted. In sub-network 2 RPS4Y1 is in the central position. Sub-network three exhibit the three proteins with highest LogFC among all the sub-networks, as is all shown in Fig. 3. Table 1 shows fold changes of 11 DCMs and relevant P value, respectively. Figure 4 shows expression level of selected 11 genes between control and DCM samples. Among the 11 genes, nine of them are up-regulated while two are down-regulated.
Analysis of transcription factor
Transcription factor that modulate gene expression in DCM predicted by PASTAA is shown in Table 2. As shown in the Fig. 5a, Tef-1 and TFIIA transcription factor families play an important role in up-regulation of gene expression. TBP and TFIIA show features of constitutive expression as widely existing in transcription. By contrast, as shown in Fig. 5b, CORIN, FIGF of Myb families and Hnf-4 family are predicted to play an important role in down-regulation of gene expression. Figure 6 shows the transcription factor-binding site predicted by JASPAR.
Acquisition of target genes of drugs
Gene set A and B was uploaded to Connectivity Map to search for potential drugs, Table 3 shows predicted drug molecules that may induce/inhibit DCM-associated gene expression. Predicted molecule in Table 2 was uploaded to NCBI PubChem database to search for target genes of receptors. 19 target genes are significantly correlated with DCM, as shown in Table 4. Numbers in the brackets stand for the count of drugs. CTGF, POSTN, CORIN, FIGH have higher amount of drugs. Collagen family and THBS2, THBS4 are also ranked highly.
Microarray analysis, an optimal approach to identify differentially expressed gene, helps to define an early diagnosis and lower misdiagnosis rate . So far, the application of microarray analysis has revealed considerable bioprocesses associated with DCM. Barth et al  predicted and constructed a gene set including 27 genes differentially expressed, among which 25 were later confirmed. Micaela et al  analyzed ion channel-associated gene differential expression. They reported that ion channel-associated gene such as SCN2B. KCNJ5, CLIC2 play an important role in related bioprocesses. The present study aims at underlining the effects of regulation by transcription factors on differential gene expression. Effectiveness of target gene set is confirmed by referring to existing database.
By comparing differentially expressed genes in DCM samples with that from healthy controls, we predict that CORIN, FIGF, CTGF, COL1A1, COL1A2, THBS4, THBS4, POSTN, COMP, PRELP and RPS4Y1 may play a role in development of DCM. By the way, CORIN (encoding production could convert pro-ANP to bioactive ANP) is one of the known marker genes, the purpose of this study is to explore the mechanism underlying the differential expression level.
Cardiac remodeling characterized by collagen deposition in extracellular matrix  and myocardial fibrosis leading to heart failure [24, 25]. Expression level changes of genes encoding proteins associated with: 1) ECM-receptor interaction, 2) focal adhesion, 3) TGF-beta signaling pathway may play an integral role in development of systolic dysfunction. Up-regulation of CTGF in TGF-beta signaling pathway can lead to an increased deposition of type 3 collagen (COL3A1) and type 1 collagen (COL1A1, COL1A2 and COL1A3) [26, 27]. Furthermore, excessive collagen may cause myocardial fibrosis and heart failure. THBS2 and THBS4, existing in all three pathways that encode fibronectin are also modulated by CTGF. Previous study shows that up-regulation of THBS1 and THBS4 may lead to matri-cellular protein deposition, and subsequently resulting in DCM, heart failure or death. FIGF (also known as vascular endothelial growth factor-D; VEGF-D) in ECM-receptor interaction may also participate in the pathology of DCM. Gong X et al.  reported that VEGF preserve cardiac function after intra-myocardial transplantation in a DCM mouse by reducing cellular apoptosis and myocardial fibrosis in addition to enhanced angiogenesis, indicating that down-regulation of FIGF may improve cardiac function and preserve myocardial cells.
Differential expression gene exhibit up-regulation in GO terms of extracellular matrix and focal adhesion . Mal-regulation of extracellular matrix is associated with progression of cardiac remodeling and heart failure. COMP encodes a non-collagen protein in ECM and previous studies have demonstrated that abnormal expression of COMP result in myocardial cell apoptosis and loss of myofilament . COMP, COL1A1, COL1A2 and PRELP are up-regulated in ECM-receptor interaction and focal adhesion, inducing a component change in extracellular matrix .
POSTN bond with heparin, and is showed in the PPI network co-expressed with COL1A1, COL1A2, THBS4 and CTGF. High impression level of POSTN in myocardial fibroblast contributes to cardiac remodeling . Therefore, highly co-expression of both POSTN and CTGF (regulate heparin in fibroblasts) may be associated with antagonism. However, it is unknown whether POSTN is regulated by CTGF. The co-expression of RPS4Y1, DDX3Y, EIF1AY and RPL13 were discovered. Heidecker et al.  demonstrated that they are associated with idiopathic DCM, but the expression level differentiates significantly and the mechanism remains unclear.
Transcription factors may affect the differential expression genes in DCM. Tef-1 (also known as TEAD1, TCF-13), a transcription enhancer, promotes the expression of troponin, myosin and actin [34, 35]. TFIIA, a member of TF family, exists ubiquitously in global tissues, participating in the formation of transcription initiation complex as an intra-nuclear protein . It may trigger the up-regulation of differentially expressed genes. MYB , known as a transcription activator, encodes transcription factors of myb family. MYB has 1) a N-terminal DNA binding domain, 2) a transcription activating domain in the center  3) a C-terminal domain associated with inhibition of transcription. The C-domain may be responsible for down-regulation of FIGF and CORIN. Hnf-4 transcription factor family is mainly expressed in kidney and liver . Given that it is a transcription activator, the inhibition may down-regulate the expression of CORIN, FIGF and subsequentially the down-regulation of ANP level.
Re-localization of drugs is a cost-effective approach so that the development of CMap is engaged in efficiently predicting potential targets that drugs can aim at. Tridihexethyl, fluoxetine and pirlindole can inactivate G-protein-coupled receptor. Fluoxetine, methyldopa and sulfadiazine can inactivate cystic fibrosis trans-membrane conductance regulator and sequentially down-regulate the expression of Collagen, THBS, COMP and PRELP. Methyldopa, sulfadiazine and khellin may further regulate the gene expression by regulating the translation. Besides, up-regulated CTGF, POSTN and down-regulated CORN, FIGF are seen as the major targets in CMap prediction, indicating that combination of multiple drugs may achieve a better treatment.
However, up or down regulation of specific genes may not be responsible for pathogenesis of the relevant diseases [40, 41]. They may be only a downstream reaction or byproducts, such as s NPPB(logFC> 4.2)or NPPA(logFC> 2.5)in DCM. Age, gender, multi-drug therapy, etiology and individual variation may also cause differential gene expression. That is the limitation of the present study. In addition, we all know now that myofibril gene (protein) mutations are associated with the development of various cardiomyopathies including DCM. It is necessary in the future to carry out more analyses including the above-mentioned factors and myofibril protein mutations. Besides, differentially expressed genes in the same expression profile vary among different filters. These factors have also been taken into consideration in drawing any conclusions.
Bioinformatics-based analyses reveal the targeted genes probably associated with cardiomyopathy, which provide clues for pharmacological therapies aiming at the targets. Further studies considering other factors such as age, gender, multi-drug therapy, etiology and individual variation should be carried out in the future and bench work experiments should be performed to verify the relationship between the DEGs and the development of DCM.
Bone marrow-derived cells
B-type natriuretic peptides
Visualization and integrated discovery
Differential expression of genes
Granulocyte colony-stimulating factor
Expression Omnibus Gene
The Kyoto Encyclopedia of Genes and Genomes
Polymerase chain reaction
Protein-protein interaction network
Search tool for the retrieval interacting genes
Transforming growth factor beta
Xiao HB, Brecker SJ, Gibson DG. Effects of abnormal activation on the time course of the left ventricular pressure pulse in dilated cardiomyopathy. Br Heart J. 1992;68:403–7.
Meder B, Rühle F, Weis T, Homuth G, Keller A, Franke J, Peil B, Bermejo JL, Frese K, Huge A. A genome-wide association study identifies 6p21 as novel risk locus for dilated cardiomyopathy. Eur Heart J. 2014;35:1069–77.
Pinto YM, Elliott PM, Arbustini E, Adler Y, Anastasakis A, Böhm M, Duboc D, Gimeno J, de Groote P, Imazio M. Proposal for a revised definition of dilated cardiomyopathy, hypokinetic non-dilated cardiomyopathy, and its implications for clinical practice: a position statement of the ESC working group on myocardial and pericardial diseases. Eur Heart J. 2016;37(23):1850–8.
d’Ambrosio A, Patti G, Manzoli A, Sinagra G, Di Lenarda A, Silvestri F, Di Sciascio G. The fate of acute myocarditis between spontaneous improvement and evolution to dilated cardiomyopathy: a review. Heart. 2001;85:499–504.
Barth AS, Kuner R, Buness A, Ruschhaupt M, Merk S, Zwermann L, Kääb S, Kreuzer E, Steinbeck G, Mansmann U. Identification of a common gene expression signature in dilated cardiomyopathy across independent microarray studies. J Am Coll Cardiol. 2006;48:1610–7.
Yung CK, Halperin VL, Tomaselli GF, Winslow RL. Gene expression profiles in end-stage human idiopathic dilated cardiomyopathy: altered expression of apoptotic and cytoskeletal genes. Genomics. 2004;83:281–97.
de Groote P, Pinet F, Bauters C. New technologies, new therapies: toward personalized medicine in heart failure patients? Eur Heart J. 2012;
Fett JD. Inflammation and virus in dilated cardiomyopathy as indicated by endomyocardial biopsy. Int J Cardiol. 2006;112:125–6.
Hamshere S, Arnous S, Choudhury T, Choudry F, Mozid A, Yeo C, Barrett C, Saunders N, Gulati A, Knight C. Randomized trial of combination cytokine and adult autologous bone marrow progenitor cell administration in patients with non-ischaemic dilated cardiomyopathy: the REGENERATE-DCM clinical trial. Eur Heart J. 2015;36(44):3061–9.
Xie X, Wu B, Chen Y, Li W, Hu X, Chen J. Peripheral endothelial function may predict the effectiveness of beta-blocker therapy in patients with idiopathic dilated cardiomyopathy. Int J Cardiol. 2016;221:128–33.
Jo B-S, Koh I-U, Bae J-B, Yu H-Y, Jeon E-S, Lee H-Y, Kim J-J, Choi M, Choi SS. Methylome analysis reveals alterations in DNA methylation in the regulatory regions of left ventricle development genes in human dilated cardiomyopathy. Genomics. 2016;108:84–92.
Liu Y, Morley M, Brandimarto J, Hannenhalli S, Hu Y, Ashley EA, Tang WH, Moravec CS, Margulies KB, Cappola TP, Li M. MAGNet consortium. RNA-Seq identifies novel myocardial gene expression signatures of heart failure. Genomics. 2015;105(2):83–9. https://doi.org/10.1016/j.ygeno.2014.12.002.
Barth AS, Kuner R, Buness A, Ruschhaupt M, Merk S, Zwermann L, Kääb S, Kreuzer E, Steinbeck G, Mansmann U, Poustka A, Nabauer M, Sültmann H. Identification of a common gene expression signature in dilated cardiomyopathy across independent microarray studies. J Am Coll Cardiol. 2006;48(8):1610–7.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4:1.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, Von Mering C. STRING v9. 1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41:D808–15.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Roider HG, Manke T, O'Keeffe S, Vingron M, Haas SA. PASTAA: identifying transcription factors associated with sets of co-regulated genes. Bioinformatics. 2009;25:435–42.
Wang L, Zhu L, Luan R, Wang L, Fu J, Wang X, Sui L. Analyzing gene expression profiles in dilated cardiomyopathy via bioinformatics methods. Braz J Med Biol Res. 2016;49
Sandelin A, Wasserman WW. Constrained binding site diversity within families of transcription factors enhances pattern discovery bioinformatics. J Mol Biol. 2004;338:207–15.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet J-P, Subramanian A, Ross KN. The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313:1929–35.
Ivandic BT, Mastitsky SE, Schönsiegel F, Bekeredjian R, Eils R, Frey N, Katus HA, Brors B. Whole-genome analysis of gene expression associates the ubiquitin-proteasome system with the cardiomyopathy phenotype in disease-sensitized congenic mouse strains. Cardiovasc Res. 2012;94:87–95.
Rosello-Lleti E, Navarro MM, Ortega A, Tarazon E, Sanchez D, Lazaro IS, Almenar L, Salvador A, Portoles M, Otero MR. Differential gene expression of cardiac chloride and potassium ion channels in human dilated non-ischemic cardiomyopathy. Eur Heart J. 2013;34:P4194.
Tsoutsman T, Wang X, Garchow K, Riser B, Twigg S, Semsarian C. CCN2 plays a key role in extracellular matrix gene expression in severe hypertrophic cardiomyopathy and heart failure. J Mol Cell Cardiol. 2013;62:164–78.
Touvron M, Escoubet B, Mericskay M, Angelini A, Lamotte L, Santini MP, Rosenthal N, Daegelen D, Tuil D, Decaux J-F. Locally expressed IGF1 propeptide improves mouse heart function in induced dilated cardiomyopathy by blocking myocardial fibrosis and SRF-dependent CTGF induction. Disease Models and Mechanisms. 2012;5:481–91.
Koshman YE, Sternlicht MD, Kim T, O'Hara CP, Koczor CA, Lewis W, Seeley TW, Lipson KE, Samarel AM. Connective tissue growth factor regulates cardiac function and tissue remodeling in a mouse model of dilated cardiomyopathy. J Mol Cell Cardiol. 2015;89:214–22.
Thiele B-J, Doller A, Kähne T, Pregla R, Hetzer R, Regitz-Zagrosek V. RNA-binding proteins heterogeneous nuclear ribonucleoprotein A1, E1, and K are involved in post-transcriptional control of collagen I and III synthesis. Circ Res. 2004;95:1058–66.
Zhao H-p LD, Zhang W, Zhang L, S-m W, Ma C-m, Qin C, L-f Z. Protective action of tetramethylpyrazine phosphate against dilated cardiomyopathy in cTnTR141W transgenic mice. Acta Pharmacol Sin. 2010;31:281–8.
Gong X, Wang P, Wu Q, Wang S, Yu L, Wang G. Human umbilical cord blood derived mesenchymal stem cells improve cardiac function in cTnT R141W transgenic mouse of dilated cardiomyopathy. Eur J Cell Biol. 2016;95:57–67.
Louzao-Martinez L, Vink A, Harakalova M, Asselbergs FW, Verhaar MC, Cheng C. Characteristic adaptations of the extracellular matrix in dilated cardiomyopathy. Int J Cardiol. 2016;220:634–46.
Huang Y, Xia J, Zheng J, Geng B, Liu P, Yu F, Liu B, Zhang H, Xu M, Ye P. Deficiency of cartilage oligomeric matrix protein causes dilated cardiomyopathy. Basic Res Cardiol. 2013;108:1–21.
Li H, Cui Y, Luan J, Zhang X, Li C, Zhou X, Shi L, Wang H, Han J. PRELP (proline/arginine-rich end leucine-rich repeat protein) promotes osteoblastic differentiation of preosteoblastic MC3T3-E1 cells by regulating the β-catenin pathway. Biochem Biophys Res Commun. 2016;470:558–62.
Guan J, Liu W-Q, Xing M-Q, Shi Y, Tan X-Y, Jiang C-Q, Dai H-Y. Elevated expression of periostin in diabetic cardiomyopathy and the effect of valsartan. BMC Cardiovasc Disord. 2015;15:90.
Heidecker B, Lamirault G, Kasper EK, Wittstein IS, Champion HC, Breton E, Russell SD, Hall J, Kittleson MM, Baughman KL. The gene expression profile of patients with new-onset heart failure reveals important gender-specific differences. Eur Heart J. 2009;
Benhaddou A, Keime C, Ye T, Morlon A, Michel I, Jost B, Mengus G, Davidson I. Transcription factor TEAD4 regulates expression of myogenin and the unfolded protein response genes during C2C12 cell differentiation. Cell Death Differ. 2012;19:220–31.
Jiang S-W, Trujillo MA, Sakagashira M, Wilke RA, Eberhardt NL. Novel human TEF-1 isoforms exhibit altered DNA binding and functional properties. Biochemistry. 2000;39:3505–13.
Høiby T, Zhou H, Mitsiou DJ, Stunnenberg HG. A facelift for the general transcription factor TFIIA. Biochimica et Biophysica Acta (BBA)-Gene Structure and Expression. 2007;1769:429–36.
Chen Y, Xu H, Liu J, Zhang C, Leutz A, Mo X. The c-Myb functions as a downstream target of PDGF-mediated survival signal in vascular smooth muscle cells. Biochem Biophys Res Commun. 2007;360:433–6.
Vargova K, Curik N, Burda P, Basova P, Kulvait V, Pospisil V, Savvulidi F, Kokavec J, Necas E, Berkova A. MYB transcriptionally regulates the miR-155 host gene in chronic lymphocytic leukemia. Blood. 2011;117:3816–25.
Chartier FL, Bossu J-P, Laudet V, Fruchart J-C, Laine B. Cloning and sequencing of cDNAs encoding the human hepatocyte nuclear factor 4 indicate the presence of two isoforms in human liver. Gene. 1994;147:269–72.
Fan S, Li X, Tie L, Pan Y, Li X. KIAA0101 is associated with human renal cell carcinoma proliferation and migration induced by erythropoietin. Oncotarget. 2016;7(12):13520–37. https://doi.org/10.18632/oncotarget.5876.
Qiu LL, Xu L, Guo XM, Li ZT, Wan F, Liu XP, Chen GH, Chang GB. Gene expression changes in chicken NLRC5 signal pathway associated with in vitro avian leukosis virus subgroup J infection. Genet Mol Res. 2016;15(1) https://doi.org/10.4238/gmr.15017640.
We would like to give our sincere appreciation to the editor and reviewers of the journal for their helpful comments on this article.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.