The identification of informative genes from multiple datasets with increasing complexity

Background In microarray data analysis, factors such as data quality, biological variation, and the increasingly multi-layered nature of more complex biological systems complicates the modelling of regulatory networks that can represent and capture the interactions among genes. We believe that the use of multiple datasets derived from related biological systems leads to more robust models. Therefore, we developed a novel framework for modelling regulatory networks that involves training and evaluation on independent datasets. Our approach includes the following steps: (1) ordering the datasets based on their level of noise and informativeness; (2) selection of a Bayesian classifier with an appropriate level of complexity by evaluation of predictive performance on independent data sets; (3) comparing the different gene selections and the influence of increasing the model complexity; (4) functional analysis of the informative genes. Results In this paper, we identify the most appropriate model complexity using cross-validation and independent test set validation for predicting gene expression in three published datasets related to myogenesis and muscle differentiation. Furthermore, we demonstrate that models trained on simpler datasets can be used to identify interactions among genes and select the most informative. We also show that these models can explain the myogenesis-related genes (genes of interest) significantly better than others (P < 0.004) since the improvement in their rankings is much more pronounced. Finally, after further evaluating our results on synthetic datasets, we show that our approach outperforms a concordance method by Lai et al. in identifying informative genes from multiple datasets with increasing complexity whilst additionally modelling the interaction between genes. Conclusions We show that Bayesian networks derived from simpler controlled systems have better performance than those trained on datasets from more complex biological systems. Further, we present that highly predictive and consistent genes, from the pool of differentially expressed genes, across independent datasets are more likely to be fundamentally involved in the biological process under study. We conclude that networks trained on simpler controlled systems, such as in vitro experiments, can be used to model and capture interactions among genes in more complex datasets, such as in vivo experiments, where these interactions would otherwise be concealed by a multitude of other ongoing events

Using a 3D virtual muscle model to link gene expression changes during myogenesis to protein spatial location in muscle

Background: Myogenesis is an ordered process whereby mononucleated muscle precursor cells (myoblasts) fuse into multinucleated myotubes that eventually differentiate into myofibres, involving substantial changes in gene expression and the organisation of structural components of the cells. To gain further insight into the orchestration of these structural changes we have overlaid the spatial organisation of the protein components of a muscle cell with their gene expression changes during differentiation using a new 3D visualisation tool: the Virtual Muscle 3D (VMus3D)

Ataxin-3 Plays a Role in Mouse Myogenic Differentiation through Regulation of Integrin Subunit Levels

BACKGROUND: During myogenesis several transcription factors and regulators of protein synthesis and assembly are rapidly degraded by the ubiquitin-proteasome system (UPS). Given the potential role of the deubiquitinating enzyme (DUB) ataxin-3 in the UPS, and the high expression of the murine ataxin-3 homolog in muscle during embryogenesis, we sought to define its role in muscle differentiation. METHODOLOGY/PRINCIPAL FINDINGS: Using immunofluorescence analysis, we found murine ataxin-3 (mATX3) to be highly expressed in the differentiated myotome of E9.5 mouse embryos. C2C12 myoblasts depleted of mATX3 by RNA interference exhibited a round morphology, cell misalignment, and a delay in differentiation following myogenesis induction. Interestingly, these cells showed a down-regulation of alpha5 and alpha7 integrin subunit levels both by immunoblotting and immunofluorescence. Mouse ATX3 was found to interact with alpha5 integrin subunit and to stabilize this protein by repressing its degradation through the UPS. Proteomic analysis of mATX3-depleted C2C12 cells revealed alteration of the levels of several proteins related to integrin signaling. CONCLUSIONS: Ataxin-3 is important for myogenesis through regulation of integrin subunit levels.This work was financed by the Fundacao para a Ciencia e a Tecnologia (FCT) (POCI/SAU-MMO/60412/2002) and by National Institutes of Health/National Institute of Neurological Disorders and Stroke (NIH/NINDS) grant RO1 NS038712 to HLP. MCC, FB, AJR, and RJT were supported by the FCT fellowships (SFRH/BD/9759/2003 and SFRH/BPD/28560/2006), (SFRH/BPD/17368/2004), (SFRH/BD/17066/2004), (SFRH/BD/29947/2006), respectively. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

β1-Syntrophin Modulation by miR-222 in mdx Mice

Background: In mdx mice, the absence of dystrophin leads to the deficiency of other components of the dystrophin-glycoprotein complex (DAPC), making skeletal muscle fibers more susceptible to necrosis. The mechanisms involved in the disappearance of the DAPC are not completely understood. The muscles of mdx mice express normal amounts of mRNA for the DAPC components, thus suggesting post-transcriptional regulation. Methodology/Principal Findings: We investigated the hypothesis that DAPC reduction could be associated with the microRNA system. Among the possible microRNAs (miRs) found to be upregulated in the skeletal muscle tissue of mdx compared to wt mice, we demonstrated that miR-222 specifically binds to the 3′-UTR of β1-syntrophin and participates in the downregulation of β1-syntrophin. In addition, we documented an altered regulation of the 3′-UTR of β1-syntrophin in muscle tissue from dystrophic mice. Conclusion/Significance: These results show the importance of the microRNA system in the regulation of DAPC components in dystrophic muscle, and suggest a potential role of miRs in the pathophysiology of dystrophy. © 2010 De Arcangelis et al

Comparative Analyses of SUV420H1 Isoforms and SUV420H2 Reveal Differences in Their Cellular Localization and Effects on Myogenic Differentiation

Methylation of histone H4 on lysine 20 plays critical roles in chromatin structure and function via mono- (H4K20me1), di- (H4K20me2), and trimethyl (H4K20me3) derivatives. In previous analyses of histone methylation dynamics in mid-gestation mouse embryos, we documented marked changes in H4K20 methylation during cell differentiation. These changes were particularly robust during myogenesis, both in vivo and in cell culture, where we observed a transition from H4K20me1 to H4K20me3. To assess the significance of this change, we used a gain-of-function strategy involving the lysine methyltransferases SUV420H1 and SUV420H2, which catalyze H4K20me2 and H4K20me3. At the same time, we characterized a second isoform of SUV420H1 (designated SUV420H1_i2) and compared the activity of all three SUV420H proteins with regard to localization and H4K20 methylation.Immunofluorescence revealed that exogenous SUV420H1_i2 was distributed throughout the cell, while a substantial portion of SUV420H1_i1 and SUV420H2 displayed the expected association with constitutive heterochromatin. Moreover, SUV420H1_i2 distribution was unaffected by co-expression of heterochromatin protein-1α, which increased the targeting of SUV420H1_i1 and SUV420H2 to regions of pericentromeric heterochromatin. Consistent with their distributions, SUV420H1_i2 caused an increase in H4K20me3 levels throughout the nucleus, whereas SUV420H1_i1 and SUV420H2 facilitated an increase in pericentric H4K20me3. Striking differences continued when the SUV420H proteins were tested in the C2C12 myogenic model system. Specifically, although SUV420H1_i2 induced precocious appearance of the differentiation marker Myogenin in the presence of mitogens, only SUV420H2 maintained a Myogenin-enriched population over the course of differentiation. Paradoxically, SUV420H1_i1 could not be expressed in C2C12 cells, which suggests it is under post-transcriptional or post-translational control.These data indicate that SUV420H proteins differ substantially in their localization and activity. Importantly, SUV420H2 can induce a transition from H4K20me1 to H4K20me3 in regions of constitutive heterochromatin that is sufficient to enhance myogenic differentiation, suggesting it can act an as epigenetic ‘switch’ in this process

Validation of Skeletal Muscle cis-Regulatory Module Predictions Reveals Nucleotide Composition Bias in Functional Enhancers

We performed a genome-wide scan for muscle-specific cis-regulatory modules (CRMs) using three computational prediction programs. Based on the predictions, 339 candidate CRMs were tested in cell culture with NIH3T3 fibroblasts and C2C12 myoblasts for capacity to direct selective reporter gene expression to differentiated C2C12 myotubes. A subset of 19 CRMs validated as functional in the assay. The rate of predictive success reveals striking limitations of computational regulatory sequence analysis methods for CRM discovery. Motif-based methods performed no better than predictions based only on sequence conservation. Analysis of the properties of the functional sequences relative to inactive sequences identifies nucleotide sequence composition can be an important characteristic to incorporate in future methods for improved predictive specificity. Muscle-related TFBSs predicted within the functional sequences display greater sequence conservation than non-TFBS flanking regions. Comparison with recent MyoD and histone modification ChIP-Seq data supports the validity of the functional regions

Magic-factor 1, a partial agonist of Met, induces muscle hypertrophy by protecting myogenic progenitors from apoptosis.

Hepatocyte Growth Factor (HGF) is a pleiotropic cytokine of mesenchymal origin that mediates a characteristic array of biological activities including cell proliferation, survival, motility and morphogenesis. Its high affinity receptor, the tyrosine kinase Met, is expressed by a wide range of tissues and can be activated by either paracrine or autocrine stimulation. Adult myogenic precursor cells, the so called satellite cells, express both HGF and Met. Following muscle injury, autocrine HGF-Met stimulation plays a key role in promoting activation and early division of satellite cells, but is shut off in a second phase to allow myogenic differentiation. In culture, HGF stimulation promotes proliferation of muscle precursors thereby inhibiting their differentiation