Expression Analysis of Wound Healing Genes in Human Periapical Granulomas of Progressive and Stable Nature

Introduction: Wound healing process involves the activation of extracellular matrix components, remodeling enzymes, cellular adhesion molecules, growth factors, cytokines and chemokines genes. However, the molecular patterns underlying the healing process periapical environment remain unclear. Here we hypothesized that endodontic infection might result in an imbalance in the expression of wound healing genes involved in the pathogenesis of periapical lesions. Furthermore, we suggest that differential expression of wound healing markers in active and latent granulomas could account for different clinical outcomes for such lesions. Methods: Study samples consisted of 93 periapical granulomas collected after endodontic surgeries and 24 healthy periodontal ligament tissues collected from premolars extracted for orthodontic purposes as control samples. Of these, 10 periapical granulomas and 5 healthy periapical tissues were used for expression analysis of 84 wound healing genes by using a pathway-specific real-time polymerase chain reaction array. The remaining 83 granulomas and all 24 control specimens were used to validate the obtained array data by real-time polymerase chain reaction. Observed variations in expression of wound healing genes were analyzed according to the classification of periapical granulomas as active/progressive versus inactive/stable (as determined by receptor activator for nuclear factor kappa B ligand/osteoprotegerin expression ratio). Results: We observed a marked increase of 5-fold or greater in SERPINE1, TIMP1, COL1A1, COL5A1, VTN, CTGF, FGF7, TGFB1, TNF, CXCL11, ITGA4, and ITGA5 genes in the periapical granulomas when compared with control samples. SERPINE1, TIMP1, COL1A1, TGFB1, and ITGA4 mRNA expression was significantly higher in inactive compared with active periapical granulomas (P < .001), whereas TNF and CXCL11 mRNA expression was higher in active lesions (P < .001). Conclusions: The identification of novel gene targets that curb the progression status of periapical lesions might contribute to a more accurate diagnosis and lead to treatment modalities more conducive to endodontic success. (J Endod 2012;38:185-190)

A feature selection approach for identification of signature genes from SAGE data

Abstract Background One goal of gene expression profiling is to identify signature genes that robustly distinguish different types or grades of tumors. Several tumor classifiers based on expression profiling have been proposed using microarray technique. Due to important differences in the probabilistic models of microarray and SAGE technologies, it is important to develop suitable techniques to select specific genes from SAGE measurements. Results A new framework to select specific genes that distinguish different biological states based on the analysis of SAGE data is proposed. The new framework applies the bolstered error for the identification of strong genes that separate the biological states in a feature space defined by the gene expression of a training set. Credibility intervals defined from a probabilistic model of SAGE measurements are used to identify the genes that distinguish the different states with more reliability among all gene groups selected by the strong genes method. A score taking into account the credibility and the bolstered error values in order to rank the groups of considered genes is proposed. Results obtained using SAGE data from gliomas are presented, thus corroborating the introduced methodology. Conclusion The model representing counting data, such as SAGE, provides additional statistical information that allows a more robust analysis. The additional statistical information provided by the probabilistic model is incorporated in the methodology described in the paper. The introduced method is suitable to identify signature genes that lead to a good separation of the biological states using SAGE and may be adapted for other counting methods such as Massive Parallel Signature Sequencing (MPSS) or the recent Sequencing-By-Synthesis (SBS) technique. Some of such genes identified by the proposed method may be useful to generate classifiers.

Different patterns of evolution for duplicated DNA repair genes in bacteria of the Xanthomonadales group

Abstract Background DNA repair genes encode proteins that protect organisms against genetic damage generated by environmental agents and by-products of cell metabolism. The importance of these genes in life maintenance is supported by their high conservation, and the presence of duplications of such genes may be easily traced, especially in prokaryotic genomes. Results The genome sequences of two Xanthomonas species were used as the basis for phylogenetic analyses of genes related to DNA repair that were found duplicated. Although 16S rRNA phylogenetic analyses confirm their classification at the basis of the gamma proteobacteria subdivision, differences were found in the origin of the various genes investigated. Except for lexA, detected as a recent duplication, most of the genes in more than one copy are represented by two highly divergent orthologs. Basically, one of such duplications is frequently positioned close to other gamma proteobacteria, but the second is often positioned close to unrelated bacteria. These orthologs may have occurred from old duplication events, followed by extensive gene loss, or were originated from lateral gene transfer (LGT), as is the case of the uvrD homolog. Conclusions Duplications of DNA repair related genes may result in redundancy and also improve the organisms' responses to environmental challenges. Most of such duplications, in Xanthomonas, seem to have arisen from old events and possibly enlarge both functional and evolutionary genome potentiality.

RNA interference-mediated knockdown of CD49e (α5 integrin chain) in human thymic epithelial cells modulates the expression of multiple genes and decreases thymocyte adhesion

Abstract Background The thymus is a central lymphoid organ, in which bone marrow-derived T cell precursors undergo a complex process of maturation. Developing thymocytes interact with thymic microenvironment in a defined spatial order. A component of thymic microenvironment, the thymic epithelial cells, is crucial for the maturation of T-lymphocytes through cell-cell contact, cell matrix interactions and secretory of cytokines/chemokines. There is evidence that extracellular matrix molecules play a fundamental role in guiding differentiating thymocytes in both cortical and medullary regions of the thymic lobules. The interaction between the integrin α5β1 (CD49e/CD29; VLA-5) and fibronectin is relevant for thymocyte adhesion and migration within the thymic tissue. Our previous results have shown that adhesion of thymocytes to cultured TEC line is enhanced in the presence of fibronectin, and can be blocked with anti-VLA-5 antibody. Results Herein, we studied the role of CD49e expressed by the human thymic epithelium. For this purpose we knocked down the CD49e by means of RNA interference. This procedure resulted in the modulation of more than 100 genes, some of them coding for other proteins also involved in adhesion of thymocytes; others related to signaling pathways triggered after integrin activation, or even involved in the control of F-actin stress fiber formation. Functionally, we demonstrated that disruption of VLA-5 in human TEC by CD49e-siRNA-induced gene knockdown decreased the ability of TEC to promote thymocyte adhesion. Such a decrease comprised all CD4/CD8-defined thymocyte subsets. Conclusion Conceptually, our findings unravel the complexity of gene regulation, as regards key genes involved in the heterocellular cell adhesion between developing thymocytes and the major component of the thymic microenvironment, an interaction that is a mandatory event for proper intrathymic T cell differentiation.

Cell organisation, sulphur metabolism and ion transport-related genes are differentially expressed in Paracoccidioides brasiliensis mycelium and yeast cells

Abstract Background Mycelium-to-yeast transition in the human host is essential for pathogenicity by the fungus Paracoccidioides brasiliensis and both cell types are therefore critical to the establishment of paracoccidioidomycosis (PCM), a systemic mycosis endemic to Latin America. The infected population is of about 10 million individuals, 2% of whom will eventually develop the disease. Previously, transcriptome analysis of mycelium and yeast cells resulted in the assembly of 6,022 sequence groups. Gene expression analysis, using both in silico EST subtraction and cDNA microarray, revealed genes that were differential to yeast or mycelium, and we discussed those involved in sugar metabolism. To advance our understanding of molecular mechanisms of dimorphic transition, we performed an extended analysis of gene expression profiles using the methods mentioned above. Results In this work, continuous data mining revealed 66 new differentially expressed sequences that were MIPS(Munich Information Center for Protein Sequences)-categorised according to the cellular process in which they are presumably involved. Two well represented classes were chosen for further analysis: (i) control of cell organisation – cell wall, membrane and cytoskeleton, whose representatives were hex (encoding for a hexagonal peroxisome protein), bgl (encoding for a 1,3-β-glucosidase) in mycelium cells; and ags (an α-1,3-glucan synthase), cda (a chitin deacetylase) and vrp (a verprolin) in yeast cells; (ii) ion metabolism and transport – two genes putatively implicated in ion transport were confirmed to be highly expressed in mycelium cells – isc and ktp, respectively an iron-sulphur cluster-like protein and a cation transporter; and a putative P-type cation pump (pct) in yeast. Also, several enzymes from the cysteine de novo biosynthesis pathway were shown to be up regulated in the yeast form, including ATP sulphurylase, APS kinase and also PAPS reductase. Conclusion Taken together, these data show that several genes involved in cell organisation and ion metabolism/transport are expressed differentially along dimorphic transition. Hyper expression in yeast of the enzymes of sulphur metabolism reinforced that this metabolic pathway could be important for this process. Understanding these changes by functional analysis of such genes may lead to a better understanding of the infective process, thus providing new targets and strategies to control PCM.

Analysis of Schistosoma mansoni genes shared with Deuterostomia and with possible roles in host interactions

Abstract Background: Schistosoma mansoni is a blood helminth parasite that causes schistosomiasis, a disease that affects 200 million people in the world. Many orthologs of known mammalian genes have been discovered in this parasite and evidence is accumulating that some of these genes encode proteins linked to signaling pathways in the parasite that appear to be involved with growth or development, suggesting a complex co-evolutionary process. Results: In this work we found 427 genes conserved in the Deuterostomia group that have orthologs in S. mansoni and no members in any nematodes and insects so far sequenced. Among these genes we have identified Insulin Induced Gene (INSIG), Interferon Regulatory Factor (IRF) and vasohibin orthologs, known to be involved in mammals in mevalonate metabolism, immune response and angiogenesis control, respectively. We have chosen these three genes for a more detailed characterization, which included extension of their cloned messages to obtain full-length sequences. Interestingly, SmINSIG showed a 10-fold higher expression in adult females as opposed to males, in accordance with its possible role in regulating egg production. SmIRF has a DNA binding domain, a tryptophan-rich N-terminal region and several predicted phosphorylation sites, usually important for IRF activity. Fourteen different alternatively spliced forms of the S. mansoni vasohibin (SmVASL) gene were detected that encode seven different protein isoforms including one with a complete C-terminal end, and other isoforms with shorter C-terminal portions. Using S. mansoni homologs, we have employed a parsimonious rationale to compute the total gene losses/gains in nematodes, arthropods and deuterostomes under either the Coelomata or the Ecdysozoa evolutionary hypotheses; our results show a lower losses/gains number under the latter hypothesis. Conclusion: The genes discussed which are conserved between S. mansoni and deuterostomes, probably have an ancient origin and were lost in Ecdysozoa, being still present in Lophotrochozoa. Given their known functions in Deuterostomia, it is possible that some of them have been co-opted to perform functions related (directly or indirectly) to host adaptation or interaction with host signaling processes.

Genome analysis of DNA repair genes in the alpha proteobacterium Caulobacter crescentus

Abstract Background The integrity of DNA molecules is fundamental for maintaining life. The DNA repair proteins protect organisms against genetic damage, by removal of DNA lesions or helping to tolerate them. DNA repair genes are best known from the gamma-proteobacterium Escherichia coli, which is the most understood bacterial model. However, genome sequencing raises questions regarding uniformity and ubiquity of these DNA repair genes and pathways, reinforcing the need for identifying genes and proteins, which may respond to DNA damage in other bacteria. Results In this study, we employed a bioinformatic approach, to analyse and describe the open reading frames potentially related to DNA repair from the genome of the alpha-proteobacterium Caulobacter crescentus. This was performed by comparison with known DNA repair related genes found in public databases. As expected, although C. crescentus and E. coli bacteria belong to separate phylogenetic groups, many of their DNA repair genes are very similar. However, some important DNA repair genes are absent in the C. crescentus genome and other interesting functionally related gene duplications are present, which do not occur in E. coli. These include DNA ligases, exonuclease III (xthA), endonuclease III (nth), O6-methylguanine-DNA methyltransferase (ada gene), photolyase-like genes, and uracil-DNA-glycosylases. On the other hand, the genes imuA and imuB, which are involved in DNA damage induced mutagenesis, have recently been described in C. crescentus, but are absent in E. coli. Particularly interesting are the potential atypical phylogeny of one of the photolyase genes in alpha-proteobacteria, indicating an origin by horizontal transfer, and the duplication of the Ada orthologs, which have diverse structural configurations, including one that is still unique for C. crescentus. Conclusion The absence and the presence of certain genes are discussed and predictions are made considering the particular aspects of the C. crescentus among other known DNA repair pathways. The observed differences enlarge what is known for DNA repair in the Bacterial world, and provide a useful framework for further experimental studies in this organism.

Gene expression down-regulation in CD90+ prostate tumor-associated stromal cells involves potential organ-specific genes

Abstract Background The prostate stroma is a key mediator of epithelial differentiation and development, and potentially plays a role in the initiation and progression of prostate cancer. The tumor-associated stroma is marked by increased expression of CD90/THY1. Isolation and characterization of these stromal cells could provide valuable insight into the biology of the tumor microenvironment. Methods Prostate CD90+ stromal fibromuscular cells from tumor specimens were isolated by cell-sorting and analyzed by DNA microarray. Dataset analysis was used to compare gene expression between histologically normal and tumor-associated stromal cells. For comparison, stromal cells were also isolated and analyzed from the urinary bladder. Results The tumor-associated stromal cells were found to have decreased expression of genes involved in smooth muscle differentiation, and those detected in prostate but not bladder. Other differential expression between the stromal cell types included that of the CXC-chemokine genes. Conclusion CD90+ prostate tumor-associated stromal cells differed from their normal counterpart in expression of multiple genes, some of which are potentially involved in organ development.

Association of LEC and tnpA Helicobacter pylori genes with gastric cancer in a Brazilian population

Abstract Background H. pylori seroprevalence in Brazilians varies and is dependent on socioeconomic status, sanitation conditions and ethnicity; furthermore, H. pylori is not always associated with the incidence of gastric cancer, suggesting the role of more virulent strains. The purpose of this study was to analyze the association of more virulent H. pylori strains with gastric cancer. Methods DNA was extracted from gastric biopsies of thirty-four cases of gastric cancer (11 intestinal-type, 23 diffuse-type), and thirty-four of patients with endoscopic gastritis. The presence of cagPAI genes (cagA, cagA promoter, cagE, cagM, tnpB, tnpA, cagT and the left end of the cagII (LEC)) and babA were analyzed by PCR. Results Comparison of H. pylori isolates from gastric cancer and gastritis patients showed significant associations of tnpA and LEC with gastric cancer (73.5% [OR, 6.66; 95% CI, 2.30-19.25] and 58.8% [OR, 10.71; 95% CI, 3.07-37.28] of cases, respectively). Other cagPAI genes were detected in both groups at similar frequencies. Conclusions tnpA and LEC of H. pylori cagPAI were associated with gastric cancer; nonetheless, these results were restricted within this group of patients and further studies are needed to confirm these results in a larger sample and determine their role in gastric carcinogenesis.

Body weight, metabolism and clock genes

Biological rhythms are present in the lives of almost all organisms ranging from plants to more evolved creatures. These oscillations allow the anticipation of many physiological and behavioral mechanisms thus enabling coordination of rhythms in a timely manner, adaption to environmental changes and more efficient organization of the cellular processes responsible for survival of both the individual and the species. Many components of energy homeostasis exhibit circadian rhythms, which are regulated by central (suprachiasmatic nucleus) and peripheral (located in other tissues) circadian clocks. Adipocyte plays an important role in the regulation of energy homeostasis, the signaling of satiety and cellular differentiation and proliferation. Also, the adipocyte circadian clock is probably involved in the control of many of these functions. Thus, circadian clocks are implicated in the control of energy balance, feeding behavior and consequently in the regulation of body weight. In this regard, alterations in clock genes and rhythms can interfere with the complex mechanism of metabolic and hormonal anticipation, contributing to multifactorial diseases such as obesity and diabetes. The aim of this review was to define circadian clocks by describing their functioning and role in the whole body and in adipocyte metabolism, as well as their influence on body weight control and the development of obesity.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

Abstract Background RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known. Results A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages. Conclusion We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Identification of differentially expressed genes from Trichoderma harzianum during growth on cell wall of Fusarium solani as a tool for biotechnological application

Background: The species of T. harzianum are well known for their biocontrol activity against many plant pathogens. However, there is a lack of studies concerning its use as a biological control agent against F. solani, a pathogen involved in several crop diseases. In this study, we have used subtractive library hybridization (SSH) and quantitative real-time PCR (RT-qPCR) techniques in order to explore changes in T. harzianum genes expression during growth on cell wall of F. solani (FSCW) or glucose. RT-qPCR was also used to examine the regulation of 18 genes, potentially involved in biocontrol, during confrontation between T. harzianum and F. solani. Results: Data obtained from two subtractive libraries were compared after annotation using the Blast2GO suite. A total of 417 and 78 readable EST sequence were annotated in the FSCW and glucose libraries, respectively. Functional annotation of these genes identified diverse biological processes and molecular functions required during T. harzianum growth on FSCW or glucose. We identified various genes of biotechnological value encoding to proteins which function such as transporters, hydrolytic activity, adherence, appressorium development and pathogenesis. Fifteen genes were up-regulated and sixteen were down-regulated at least at one-time point during growth of T. harzianum in FSCW. During the confrontation assay most of the genes were up-regulated, mainly after contact, when the interaction has been established. Conclusions: This study demonstrates that T. harzianum expressed different genes when grown on FSCW compared to glucose. It provides insights into the mechanisms of gene expression involved in mycoparasitism of T. harzianum against F. solani. The identification and evaluation of these genes may contribute to the development of an efficient biological control agent.

Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera

Background: The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton. Results: Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5′-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation. Conclusions: These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.

Association of genetic variants in the promoter region of genes encoding p22phox (CYBA) and glutamate cysteine ligase catalytic subunit (GCLC) and renal disease in patients with type 1 diabetes mellitus

Background Oxidative stress is recognized as a major pathogenic factor of cellular damage caused by hyperglycemia. NOX/NADPH oxidases generate reactive oxygen species and NOX1, NOX2 and NOX4 isoforms are expressed in kidney and require association with subunit p22phox (encoded by the CYBA gene). Increased expression of p22phox was described in animal models of diabetic nephropathy. In the opposite direction, glutathione is one of the main endogenous antioxidants whose plasmatic concentrations were reported to be reduced in diabetes patients. The aim of the present investigation was to test whether functional single nucleotide polymorphisms (SNPs) in genes involved in the generation of NADPH-dependent O2•- (-675 T → A in CYBA, unregistered) and in glutathione metabolism (-129 C → T in GCLC [rs17883901] and -65 T → C in GPX3 [rs8177412]) confer susceptibility to renal disease in type 1 diabetes patients. Methods 401 patients were sorted into two groups according to the presence (n = 104) or absence (n = 196) of overt diabetic nephropathy or according to glomerular filtration rate (GFR) estimated by Modification of Diet in Renal Disease (MDRD) equation: ≥ 60 mL (n = 265) or < 60 mL/min/1.73 m2 (n = 136) and were genotyped. Results No differences were found in the frequency of genotypes between diabetic and non-diabetic subjects. The frequency of GFR < 60 mL/min was significantly lower in the group of patients carrying CYBA genotypes T/A+A/A (18.7%) than in the group carrying the T/T genotype (35.3%) (P = 0.0143) and the frequency of GFR < 60 mL/min was significantly higher in the group of patients carrying GCLC genotypes C/T+T/T (47.1%) than in the group carrying the C/C genotype (31.1%) (p = 0.0082). Logistic regression analysis identified the presence of at least one A allele of the CYBA SNP as an independent protection factor against decreased GFR (OR = 0.38, CI95% 0.14-0.88, p = 0.0354) and the presence of at least one T allele of the GCLC rs17883901 SNP as an independent risk factor for decreased GFR (OR = 2.40, CI95% 1.27-4.56, p = 0.0068). Conclusions The functional SNPs CYBA -675 T → A and GCLC rs17883901, probably associated with cellular redox imbalances, modulate the risk for renal disease in the studied population of type 1 diabetes patients and require validation in additional cohorts.

Unveiling novel genes upregulated by both rhBMP2 and rhBMP7 during early osteoblastic transdifferentiation of C2C12 cells

Abstract Findings We set out to analyse the gene expression profile of pre-osteoblastic C2C12 cells during osteodifferentiation induced by both rhBMP2 and rhBMP7 using DNA microarrays. Induced and repressed genes were intercepted, resulting in 1,318 induced genes and 704 repressed genes by both rhBMP2 and rhBMP7. We selected and validated, by RT-qPCR, 24 genes which were upregulated by rhBMP2 and rhBMP7; of these, 13 are related to transcription (Runx2, Dlx1, Dlx2, Dlx5, Id1, Id2, Id3, Fkhr1, Osx, Hoxc8, Glis1, Glis3 and Cfdp1), four are associated with cell signalling pathways (Lrp6, Dvl1, Ecsit and PKCδ) and seven are associated with the extracellular matrix (Ltbp2, Grn, Postn, Plod1, BMP1, Htra1 and IGFBP-rP10). The novel identified genes include: Hoxc8, Glis1, Glis3, Ecsit, PKCδ, LrP6, Dvl1, Grn, BMP1, Ltbp2, Plod1, Htra1 and IGFBP-rP10. Background BMPs (bone morphogenetic proteins) are members of the TGFβ (transforming growth factor-β) super-family of proteins, which regulate growth and differentiation of different cell types in various tissues, and play a critical role in the differentiation of mesenchymal cells into osteoblasts. In particular, rhBMP2 and rhBMP7 promote osteoinduction in vitro and in vivo, and both proteins are therapeutically applied in orthopaedics and dentistry. Conclusion Using DNA microarrays and RT-qPCR, we identified both previously known and novel genes which are upregulated by rhBMP2 and rhBMP7 during the onset of osteoblastic transdifferentiation of pre-myoblastic C2C12 cells. Subsequent studies of these genes in C2C12 and mesenchymal or pre-osteoblastic cells should reveal more details about their role during this type of cellular differentiation induced by BMP2 or BMP7. These studies are relevant to better understanding the molecular mechanisms underlying osteoblastic differentiation and bone repair.