Search

To investigate the phosphorylation-based signaling and protein changes occurring early in epileptogenesis, the hippocampi of mice treated with pilocarpine were examined by quantitative mass spectrometry at 4 and 24 h post-status epilepticus at vast depth. Hundreds of posttranscriptional regulatory proteins were the major early targets of increased phosphorylation. At 24 h, many protein level changes were detected and the phosphoproteome continued to be perturbed. The major targets of decreased phosphorylation at 4 and 24 h were a subset of postsynaptic density scaffold proteins, ion channels, and neurotransmitter receptors. Many proteins targeted by dephosphorylation at 4 h also had decreased protein abundance at 24 h, indicating a phosphatase-mediated weakening of synapses. Increased translation was indicated by protein changes at 24 h. These observations, and many additional indicators within this multiomic resource, suggest that early epileptogenesis is characterized by signaling that stimulates both growth and a homeostatic response that weakens excitability., Competing Interests: The authors declare no competing interests., (© 2024 The Author(s).)

Transient brain insults including status epilepticus (SE) can initiate a process termed 'epileptogenesis' that results in chronic temporal lobe epilepsy. As a consequence, the entire tri-synaptic circuit of the hippocampus is fundamentally impaired. A key role in epileptogenesis has been attributed to the CA1 region as the last relay station in the hippocampal circuit and as site of aberrant plasticity, e.g. mediated by acquired channelopathies. The transcriptional profiles of the distinct hippocampal neurons are highly dynamic during epileptogenesis. Here, we aimed to elucidate the early SE-elicited mRNA signature changes and the respective upstream regulatory cascades in CA1. RNA sequencing of CA1 was performed in the mouse pilocarpine-induced SE model at multiple time points ranging from 6 to 72 h after the initial insult. Bioinformatics was used to decipher altered gene expression, signalling cascades and their corresponding cell type profiles. Robust transcriptomic changes were detected at 6 h after SE and at subsequent time points during early epileptogenesis. Major differentially expressed mRNAs encoded primarily immediate early and excitability-related gene products, as well as genes encoding immune signalling factors. Binding sites for the transcription factors Nfkb1, Spi1, Irf8, and two Runx family members, were enriched within promoters of differentially expressed genes related to major inflammatory processes, whereas the transcriptional repressors Suz12, Nfe2l2 and Rest were associated with hyperexcitability and GABA / glutamate receptor activity. CA1 quickly responds to SE by inducing transcription of genes linked to inflammation and excitation stress. Transcription factors mediating this transcriptomic switch represent targets for new highly selected, cell type and time window-specific anti-epileptogenic strategies., (© 2023. The Author(s).)

Neuronal communication relies on the release of neurotransmitters from various populations of synaptic vesicles. Despite displaying vastly different release probabilities and mobilities, the reserve and recycling pool of vesicles co-exist within a single cluster suggesting that small synaptic biomolecular condensates could regulate their nanoscale distribution. Here, we performed a large-scale activity-dependent phosphoproteome analysis of hippocampal neurons in vitro and identified Tau as a highly phosphorylated and disordered candidate protein. Single-molecule super-resolution microscopy revealed that Tau undergoes liquid-liquid phase separation to generate presynaptic nanoclusters whose density and number are regulated by activity. This activity-dependent diffusion process allows Tau to translocate into the presynapse where it forms biomolecular condensates, to selectively control the mobility of recycling vesicles. Tau, therefore, forms presynaptic nano-biomolecular condensates that regulate the nanoscale organization of synaptic vesicles in an activity-dependent manner., (© 2023. The Author(s).)

Stable function of networks requires that synapses adapt their strength to levels of neuronal activity, and failure to do so results in cognitive disorders. How such homeostatic regulation may be implemented in mammalian synapses remains poorly understood. Here we show that the phosphorylation status of several positions of the active-zone (AZ) protein RIM1 are relevant for synaptic glutamate release. Position RIMS1045 is necessary and sufficient for expression of silencing-induced homeostatic plasticity and is kept phosphorylated by serine arginine protein kinase 2 (SRPK2). SRPK2-induced upscaling of synaptic release leads to additional RIM1 nanoclusters and docked vesicles at the AZ and is not observed in the absence of RIM1 and occluded by RIM S1045E . Our data suggest that SRPK2 and RIM1 represent a presynaptic phosphosignaling hub that is involved in the homeostatic balance of synaptic coupling of neuronal networks., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

The extensive array of basic helix-loop-helix (bHLH) transcription factors and their combinations as dimers underpin the diversity of molecular function required for cell type specification during embryogenesis. The bHLH factor TWIST1 plays pleiotropic roles during development. However, which combinations of TWIST1 dimers are involved and what impact each dimer imposes on the gene regulation network controlled by TWIST1 remain elusive. In this work, proteomic profiling of human TWIST1-expressing cell lines and transcriptome analysis of mouse cranial mesenchyme have revealed that TWIST1 homodimers and heterodimers with TCF3, TCF4, and TCF12 E-proteins are the predominant dimer combinations. Disease-causing mutations in TWIST1 can impact dimer formation or shift the balance of different types of TWIST1 dimers in the cell, which may underpin the defective differentiation of the craniofacial mesenchyme. Functional analyses of the loss and gain of TWIST1-E-protein dimer activity have revealed previously unappreciated roles in guiding lineage differentiation of embryonic stem cells: TWIST1-E-protein heterodimers activate the differentiation of mesoderm and neural crest cells, which is accompanied by the epithelial-to-mesenchymal transition. At the same time, TWIST1 homodimers maintain the stem cells in a progenitor state and block entry to the endoderm lineage., (Copyright © 2020 American Society for Microbiology.)

Reliable extraction and sensitive detection of RNA from human peripheral blood mononuclear cells (PBMCs) is critical for a broad spectrum of immunology research and clinical diagnostics. RNA analysis platforms are dependent upon high-quality and high-quantity RNA; however, sensitive detection of specific responses associated with high-quality RNA extractions from human samples with limited PBMCs can be challenging. Furthermore, the comparative sensitivity between RNA quantification and best-practice protein quantification is poorly defined. Therefore, we provide herein a critical evaluation of the wide variety of current generation of RNA-based kits for PBMCs, representative of several strategies designed to maximize sensitivity. We assess these kits with a reverse transcription quantitative PCR (RT-qPCR) assay optimized for both analytically and diagnostically sensitive cell-based RNA-based applications. Specifically, three RNA extraction kits, one post-extraction RNA purification/concentration kit, four SYBR master-mix kits, and four reverse transcription kits were tested. RNA extraction and RT-qPCR reaction efficiency were evaluated with commonly used reference and cytokine genes. Significant variation in RNA expression of reference genes was apparent, and absolute quantification based on cell number was established as an effective RT-qPCR normalization strategy. We defined an optimized RNA extraction and RT-qPCR protocol with an analytical sensitivity capable of single cell RNA detection. The diagnostic sensitivity of this assay was sufficient to show a CD8 + T cell peptide epitope hierarchy with as few as 1 × 10 4 cells. Finally, we compared our optimized RNA extraction and RT-qPCR protocol with current best-practice immune assays and demonstrated that our assay is a sensitive alternative to protein-based assays for peptide-specific responses, especially with limited PBMCs number. This protocol with high analytical and diagnostic sensitivity has broad applicability for both primary research and clinical practice., (Copyright © 2020 Browne, Brady, Waardenberg, Loiseau and Doolan.)

Tuberculosis (TB) is the deadliest infectious disease worldwide. Bacille-Calmette-Guérin (BCG), the only licensed TB vaccine, affords variable protection against TB but remains the gold standard. BCG improvement is focused around three strategies: recombinant BCG strains, heterologous routes of administration, and booster vaccination. It is currently unknown whether combining these strategies is beneficial. The preclinical evaluation for new TB vaccines is heavily skewed toward immunogenicity and efficacy; however, safety and efficacy are the dominant considerations in human use. To facilitate stage gating of TB vaccines, we developed a simple empirical model to systematically rank vaccination strategies by integrating multiple measurements of safety, immunogenicity, and efficacy. We assessed 24 vaccination regimens, composed of three BCG strains and eight combinations of delivery. The model presented here highlights that mucosal booster vaccination may cause adverse outcomes and provides a much needed strategy to evaluate and rank data obtained from TB vaccine studies using different routes, strains, or animal models., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

Extensive evaluation of RNA-seq methods have demonstrated that no single algorithm consistently outperforms all others. Removal of unwanted variation (RUV) has also been proposed as a method for stabilizing differential expression (DE) results. Despite this, it remains a challenge to run multiple RNA-seq algorithms to identify significant differences common to multiple algorithms, whilst also integrating and assessing the impact of RUV into all algorithms. consensusDE was developed to automate the process of identifying significant DE by combining the results from multiple algorithms with minimal user input and with the option to automatically integrate RUV. consensusDE only requires a table describing the sample groups, a directory containing BAM files or preprocessed count tables and an optional transcript database for annotation. It supports merging of technical replicates, paired analyses and outputs a compendium of plots to guide the user in subsequent analyses. Herein, we assess the ability of RUV to improve DE stability when combined with multiple algorithms and between algorithms, through application to real and simulated data. We find that, although RUV increased fold change stability between algorithms, it demonstrated improved FDR in a setting of low replication for the intersect, the effect was algorithm specific and diminished with increased replication, reinforcing increased replication for recovery of true DE genes. We finish by offering some rules and considerations for the application of RUV in a consensus-based setting. consensusDE is freely available, implemented in R and available as a Bioconductor package, under the GPL-3 license, along with a comprehensive vignette describing functionality: http://bioconductor.org/packages/consensusDE/., Competing Interests: The authors declare there are no competing interests., (©2019 Waardenberg and Field.)

Depolarization of presynaptic terminals stimulates calcium influx, which evokes neurotransmitter release and activates phosphorylation-based signalling. Here, we present the first global temporal profile of presynaptic activity-dependent phospho-signalling, which includes two KCl stimulation levels and analysis of the poststimulus period. We profiled 1,917 regulated phosphopeptides and bioinformatically identified six temporal patterns of co-regulated proteins. The presynaptic proteins with large changes in phospho-status were again prominently regulated in the analysis of 7,070 activity-dependent phosphopeptides from KCl-stimulated cultured hippocampal neurons. Active zone scaffold proteins showed a high level of activity-dependent phospho-regulation that far exceeded the response from postsynaptic density scaffold proteins. Accordingly, bassoon was identified as the major target of neuronal phospho-signalling. We developed a probabilistic computational method, KinSwing, which matched protein kinase substrate motifs to regulated phosphorylation sites to reveal underlying protein kinase activity. This approach allowed us to link protein kinases to profiles of co-regulated presynaptic protein networks. Ca2+- and calmodulin-dependent protein kinase IIα (CaMKIIα) responded rapidly, scaled with stimulus strength, and had long-lasting activity. Mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) was the main protein kinase predicted to control a distinct and significant pattern of poststimulus up-regulation of phosphorylation. This work provides a unique resource of activity-dependent phosphorylation sites of synaptosomes and neurons, the vast majority of which have not been investigated with regard to their functional impact. This resource will enable detailed characterization of the phospho-regulated mechanisms impacting the plasticity of neurotransmitter release., Competing Interests: The authors have declared that no competing interests exist.

DNA adenine methyltransferase identification (DamID) is an enzymatic technology for detecting DNA regions targeted by chromatin-associated proteins. Proteins are fused to bacterial DNA adenine methyltransferase (Dam) and expressed in cultured cells or whole organisms. Here, we used DamID to detect DNA regions bound by the cardiac-restricted transcription factors (TFs) NKX2-5 and SRF, and ubiquitously-expressed co-factors ELK1 and ELK4. We compared targets bound by these TFs as N- and C-terminal fusions with Dam, for both wild type (WT) NKX2-5 and mutant proteins mimicking those found in congenital heart disease. Overall, DamID is highly robust: while the orientation of WT Dam fusions can affect the size of the target sets, their signatures remained largely reproducible. Furthermore, a severe NKX2-5 mutant lacking the homeodomain showed strong steric effects negatively impacting target discovery. The extent of steric effect is likely to be dependent on the protein in question and the orientation of Dam fusion., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Background: We contrast the pectoralis muscle transcriptomes of broilers selected from within a single genetic line expressing divergent feed efficiency (FE) in an effort to improve our understanding of the mechanistic basis of FE., Results: Application of a virtual muscle model to gene expression data pointed to a coordinated reduction in slow twitch muscle isoforms of the contractile apparatus (MYH15, TPM3, MYOZ2, TNNI1, MYL2, MYOM3, CSRP3, TNNT2), consistent with diminishment in associated slow machinery (myoglobin and phospholamban) in the high FE animals. These data are in line with the repeated transition from red slow to white fast muscle fibres observed in agricultural species selected on mass and FE. Surprisingly, we found that the expression of 699 genes encoding the broiler mitoproteome is modestly-but significantly-biased towards the high FE group, suggesting a slightly elevated mitochondrial content. This is contrary to expectation based on the slow muscle isoform data and theoretical physiological capacity arguments. Reassuringly, the extreme 40 most DE genes can successfully cluster the 12 individuals into the appropriate FE treatment group. Functional groups contained in this DE gene list include metabolic proteins (including opposing patterns of CA3 and CA4), mitochondrial proteins (CKMT1A), oxidative status (SEPP1, HIG2A) and cholesterol homeostasis (APOA1, INSIG1). We applied a differential network method (Regulatory Impact Factors) whose aim is to use patterns of differential co-expression to detect regulatory molecules transcriptionally rewired between the groups. This analysis clearly points to alterations in progesterone signalling (via the receptor PGR) as the major driver. We show the progesterone receptor localises to the mitochondria in a quail muscle cell line., Conclusions: Progesterone is sometimes used in the cattle industry in exogenous hormone mixes that lead to a ~20% increase in FE. Because the progesterone receptor can localise to avian mitochondria, our data continue to point to muscle mitochondrial metabolism as an important component of the phenotypic expression of variation in broiler FE.

The ability to accurately predict the DNA targets and interacting cofactors of transcriptional regulators from genome-wide data can significantly advance our understanding of gene regulatory networks. NKX2-5 is a homeodomain transcription factor that sits high in the cardiac gene regulatory network and is essential for normal heart development. We previously identified genomic targets for NKX2-5 in mouse HL-1 atrial cardiomyocytes using DNA-adenine methyltransferase identification (DamID). Here, we apply machine learning algorithms and propose a knowledge-based feature selection method for predicting NKX2-5 protein : protein interactions based on motif grammar in genome-wide DNA-binding data. We assessed model performance using leave-one-out cross-validation and a completely independent DamID experiment performed with replicates. In addition to identifying previously described NKX2-5-interacting proteins, including GATA, HAND and TBX family members, a number of novel interactors were identified, with direct protein : protein interactions between NKX2-5 and retinoid X receptor (RXR), paired-related homeobox (PRRX) and Ikaros zinc fingers (IKZF) validated using the yeast two-hybrid assay. We also found that the interaction of RXRα with NKX2-5 mutations found in congenital heart disease (Q187H, R189G and R190H) was altered. These findings highlight an intuitive approach to accessing protein-protein interaction information of transcription factors in DNA-binding experiments., (© 2016 The Authors.)

Ataxia-telangiectasia, mutated (ATM) protein plays a central role in phosphorylating a network of proteins in response to DNA damage. These proteins function in signaling pathways designed to maintain the stability of the genome and minimize the risk of disease by controlling cell cycle checkpoints, initiating DNA repair, and regulating gene expression. ATM kinase can be activated by a variety of stimuli, including oxidative stress. Here, we confirmed activation of cytoplasmic ATM by autophosphorylation at multiple sites. Then we employed a global quantitative phosphoproteomics approach to identify cytoplasmic proteins altered in their phosphorylation state in control and ataxia-telangiectasia (A-T) cells in response to oxidative damage. We demonstrated that ATM was activated by oxidative damage in the cytoplasm as well as in the nucleus and identified a total of 9,833 phosphorylation sites, including 6,686 high-confidence sites mapping to 2,536 unique proteins. A total of 62 differentially phosphorylated peptides were identified; of these, 43 were phosphorylated in control but not in A-T cells, and 19 varied in their level of phosphorylation. Motif enrichment analysis of phosphopeptides revealed that consensus ATM serine glutamine sites were overrepresented. When considering phosphorylation events, only observed in control cells (not observed in A-T cells), with predicted ATM sites phosphoSerine/phosphoThreonine glutamine, we narrowed this list to 11 candidate ATM-dependent cytoplasmic proteins. Two of these 11 were previously described as ATM substrates (HMGA1 and UIMCI/RAP80), another five were identified in a whole cell extract phosphoproteomic screens, and the remaining four proteins had not been identified previously in DNA damage response screens. We validated the phosphorylation of three of these proteins (oxidative stress responsive 1 (OSR1), HDGF, and ccdc82) as ATM dependent after H2O2 exposure, and another protein (S100A11) demonstrated ATM-dependence for translocation from the cytoplasm to the nucleus. These data provide new insights into the activation of ATM by oxidative stress through identification of novel substrates for ATM in the cytoplasm., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

Background: Gene ontology (GO) enrichment is commonly used for inferring biological meaning from systems biology experiments. However, determining differential GO and pathway enrichment between DNA-binding experiments or using the GO structure to classify experiments has received little attention., Results: Herein, we present a bioinformatics tool, CompGO, for identifying Differentially Enriched Gene Ontologies, called DiEGOs, and pathways, through the use of a z-score derivation of log odds ratios, and visualizing these differences at GO and pathway level. Through public experimental data focused on the cardiac transcription factor NKX2-5, we illustrate the problems associated with comparing GO enrichments between experiments using a simple overlap approach., Conclusions: We have developed an R/Bioconductor package, CompGO, which implements a new statistic normally used in epidemiological studies for performing comparative GO analyses and visualizing comparisons from . BED data containing genomic coordinates as well as gene lists as inputs. We justify the statistic through inclusion of experimental data and compare to the commonly used overlap method. CompGO is freely available as a R/Bioconductor package enabling easy integration into existing pipelines and is available at: http://www.bioconductor.org/packages/release/bioc/html/CompGO.html packages/release/bioc/html/CompGO.html.

We take a functional genomics approach to congenital heart disease mechanism. We used DamID to establish a robust set of target genes for NKX2-5 wild type and disease associated NKX2-5 mutations to model loss-of-function in gene regulatory networks. NKX2-5 mutants, including those with a crippled homeodomain, bound hundreds of targets including NKX2-5 wild type targets and a unique set of "off-targets", and retained partial functionality. NKXΔHD, which lacks the homeodomain completely, could heterodimerize with NKX2-5 wild type and its cofactors, including E26 transformation-specific (ETS) family members, through a tyrosine-rich homophilic interaction domain (YRD). Off-targets of NKX2-5 mutants, but not those of an NKX2-5 YRD mutant, showed overrepresentation of ETS binding sites and were occupied by ETS proteins, as determined by DamID. Analysis of kernel transcription factor and ETS targets show that ETS proteins are highly embedded within the cardiac gene regulatory network. Our study reveals binding and activities of NKX2-5 mutations on WT target and off-targets, guided by interactions with their normal cardiac and general cofactors, and suggest a novel type of gain-of-function in congenital heart disease.

Background: Hypertrophic cardiomyopathy (HCM) is caused by mutations in sarcomere protein genes, and left ventricular hypertrophy (LVH) develops as an adaptive response to sarcomere dysfunction. It remains unclear whether persistent expression of the mutant gene is required for LVH or whether early gene expression acts as an immutable inductive trigger., Objectives: The aim of this study was to use a regulatable murine model of HCM to study the reversibility of pathological LVH., Methods: The authors generated a double-transgenic mouse model, tTAxαMHCR403Q, in which expression of the HCM-causing Arg403Gln mutation in the α-myosin heavy chain (MHC) gene is inhibited by doxycycline administration. Cardiac structure and function were evaluated in groups of mice that received doxycycline for varying periods from 0 to 40 weeks of age., Results: Untreated tTAxαMHCR403Q mice showed increased left ventricular (LV) mass, contractile dysfunction, myofibrillar disarray, and fibrosis. In contrast, mice treated with doxycycline from conception to 6 weeks had markedly less LVH and fibrosis at 40 weeks. Transgene inhibition from 6 weeks reduced fibrosis but did not prevent LVH or functional changes. There were no differences in LV parameters at 40 weeks between mice with transgene inhibition from 20 weeks and mice with continuous transgene expression., Conclusions: These findings highlight the critical role of the early postnatal period in HCM pathogenesis and suggest that mutant sarcomeres manifest irreversible cardiomyocyte defects that induce LVH. In HCM, mutation-silencing therapies are likely to be ineffective for hypertrophy regression and would have to be administered very early in life to prevent hypertrophy development., (Copyright © 2015 American College of Cardiology Foundation. Published by Elsevier Inc. All rights reserved.)

Animal genomes contain a code for construction of the body plan from a fertilized egg. Understanding how genome information is deciphered to create the complex multilayered regulatory systems that drive organismal development, and which become altered in disease, is one of the greatest challenges in the biological sciences. The development of methods that effectively represent and communicate the complexity inherent in gene regulatory networks remains a major barrier. This review introduces the philosophy of systems biology and discusses recent progress in understanding the development of the heart at a systems biology level., (Copyright © 2014 Cold Spring Harbor Laboratory Press; all rights reserved.)

Background: We have recently described a method for the construction of an informative gene expression correlation landscape for a single tissue, longissimus muscle (LM) of cattle, using a small number (less than a hundred) of diverse samples. Does this approach facilitate interspecies comparison of networks?, Findings: Using gene expression datasets from LM samples from a single postnatal time point for high and low muscling sheep, and from a developmental time course (prenatal to postnatal) for normal sheep and sheep exhibiting the Callipyge muscling phenotype gene expression correlations were calculated across subsets of the data comparable to the bovine analysis. An "Always Correlated" gene expression landscape was constructed by integrating the correlations from the subsets of data and was compared to the equivalent landscape for bovine LM muscle. Whilst at the high level apparently equivalent modules were identified in the two species, at the detailed level overlap between genes in the equivalent modules was limited and generally not significant. Indeed, only 395 genes and 18 edges were in common between the two landscapes., Conclusions: Since it is unlikely that the equivalent muscles of two closely related species are as different as this analysis suggests, within tissue gene expression correlations appear to be very sensitive to the samples chosen for their construction, compounded by the different platforms used. Thus users need to be very cautious in interpretation of the differences. In future experiments, attention will be required to ensure equivalent experimental designs and use cross-species gene expression platform to enable the identification of true differences between different species.

Maintenance of cardiac structure and Z-disc signaling are key factors responsible for protecting the heart in a setting of stress, but how these processes are regulated is not well defined. We recently demonstrated that PI3K(p110α) protects the heart against myocardial infarction. The aim of this study was to determine whether PI3K(p110α) directly regulates components of the Z-disc and cardiac structure. To address this question, a unique three-dimensional virtual muscle model was applied to gene expression data from transgenic mice with increased or decreased PI3K(p110α) activity under basal conditions (sham) and in a setting of myocardial infarction to display the location of structural proteins. Key findings from this analysis were then validated experimentally. The three-dimensional virtual muscle model visually highlighted reciprocally regulated transcripts associated with PI3K activation that encoded key components of the Z-disc and costamere, including melusin. Studies were performed to assess whether PI3K and melusin interact in the heart. Here, we identify a novel melusin-PI3K interaction that generates lipid kinase activity. The direct impact of PI3K(p110α) on myocyte structure was assessed by treating neonatal rat ventricular myocytes with PI3K(p110α) inhibitors and examining the myofiber morphology of hearts from PI3K transgenic mice. Results demonstrate that PI3K is critical for myofiber maturation and Z-disc alignment. In summary, PI3K regulates the expression of genes essential for cardiac structure and Z-disc signaling, interacts with melusin, and is critical for Z-disc alignment.

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

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)., Results: Sets of generic striated muscle costamere, Z-disk and filament proteins were constructed from the literature and protein-interaction databases. Expression profiles of the genes encoding these proteins were obtained from mouse C2C12 cells undergoing myogenesis in vitro, as well as a mouse tissue survey dataset. Visualisation of the expression data in VMus3D yielded novel observations with significant relationships between the spatial location and the temporal expression profiles of the structural protein products of these genes. A muscle specificity index was calculated based on muscle expression relative to the median expression in all tissues and, as expected, genes with the highest muscle specificity were also expressed most dynamically during differentiation. Interestingly, most genes encoding costamere as well as some Z-disk proteins appeared to be broadly expressed across most tissues and showed little change in expression during muscle differentiation, in line with the broader cellular role described for some of these proteins., Conclusion: By studying gene expression patterns from a structural perspective we have demonstrated that not all genes encoding proteins that are part of muscle specific structures are simply up-regulated during muscle cell differentiation. Indeed, a group of genes whose expression program appears to be minimally affected by the differentiation process, code for proteins participating in vital skeletal muscle structures. Expression alone is a poor metric of gene behaviour. Instead, the "connectivity model of muscle development" is proposed as a mechanism for muscle development: whereby the closer to the myofibril core of muscle cells, the greater the gene expression changes during muscle differentiation and the greater the muscle specificity.