Proteostasis as a fundamental principle of Tau immunotherapy.

The microtubule-associated protein Tau is a driver of neuronal dysfunction in Alzheimer's disease and other tauopathies. In this process, Tau initially undergoes subtle changes to its abundance, subcellular localisation and a vast array of post-translational modifications including phosphorylation, that progressively result in the protein's somatodendritic accumulation and dysregulation of multiple Tau-dependent cellular processes. Given the various loss- and gain-of-functions of Tau in disease and the brain-wide changes in the proteome that characterise tauopathies, we asked whether targeting Tau would restore the alterations in proteostasis observed in disease. Therefore, by phage display, we generated a novel pan-Tau antibody, RNJ1, that preferentially binds human Tau and neutralises proteopathic seeding activity in multiple cell lines, and benchmarked it against a clinically tested pan-Tau antibody, HJ8.5 (murine version of tilavonemab). We then evaluated both antibodies, alone and in combination, in the K3 tauopathy mouse model, showing reduced Tau pathology and improvements in neuronal function following 14 weekly treatments, without obtaining synergy for the combination. These effects were more pronounced in female mice. To investigate the molecular mechanisms contributing to improvements in neuronal function, we employed quantitative proteomics, phosphoproteomics and kinase prediction analysis to first establish alterations in K3 mice relative to WT controls at the proteome level. In female K3 mice, we found 342 differentially abundant proteins, which are predominantly involved in metabolic and microtubule-associated processes, strengthening previously reported findings of defects in several functional domains in multiple tauopathy models. We next asked whether antibody-mediated Tau target engagement indirectly affects levels of deregulated proteins in the K3 model. Importantly, both immunotherapies, in particular RNJ1, induced abundance shifts towards a restoration to wild-type levels (proteostasis). A total of 257 of 342 (∼75%) proteins altered in K3 were closer in abundance to WT levels after RNJ1 treatment, and 73% after HJ8.5 treatment. However, the magnitude of these changes was less pronounced than that observed with RNJ1, as reflected by a far smaller number of differentially abundant proteins. Furthermore, analysis of the phosphoproteome showed an even stronger restoration effect with RNJ1, with ∼82% of altered phosphopeptides in K3 showing a shift to WT levels, and 75% with HJ8.5. Gene set over-representation analysis (ORA) further confirmed that proteins undergoing restoration are involved in biological pathways affected in K3 mice. Together, our study suggests that a Tau immunotherapy-induced restoration of proteostasis links target engagement and treatment efficacy.

Global quantitative proteomic analysis of aged mouse hippocampus.

Understanding the molecular changes associated with the aged brain forms the basis for developing potential strategies for slowing cognitive decline associated with normal aging. Focusing on the hippocampus, a critical brain region involved in learning and memory, we employed tandem mass tag methodology to investigate global proteomic changes that occur in advanced-aged (20-month) versus young (3-month) C57BL/6 male mice. Our analysis revealed the upregulation of 236 proteins in the old hippocampal proteome, including those enriched within several age-related processes, such as the adaptive immune response and molecular metabolic pathways, whereas downregulated proteins (88 in total) are mainly involved in axonogenesis and growth cone-related processes. Categorizing proteins by cell-type enrichment in the brain identified a general upregulation of proteins preferentially expressed in microglia, astrocytes, and oligodendrocytes. In contrast, proteins with neuron-specific expression displayed an overall age-related downregulation. By integrating our proteomic with our previously published transcriptomic data, we discovered a mild but significant positive correlation between mRNA and protein expression changes in the aged hippocampus. Therefore, this proteomic data is a valuable additional resource for further understanding age-related molecular mechanisms.

Early cortical oligodendrocyte precursor cells are transcriptionally distinct and lack synaptic connections.

Oligodendrocyte precursor cells (OPCs) generate oligodendrocytes, a process that may be tuned by neuronal activity, possibly via synaptic connections to OPCs. However, a developmental role of synaptic signaling to OPCs has so far not been shown unequivocally. To address this question, we comparatively analyzed functional and molecular characteristics of highly proliferative and migratory OPCs in the embryonic brain. Embryonic OPCs in mice (E18.5) shared the expression of voltage-gated ion channels and their dendritic morphology with postnatal OPCs, but almost completely lacked functional synaptic currents. Transcriptomic profiling of PDGFRα+ OPCs revealed a limited abundance of genes coding for postsynaptic signaling and synaptogenic cell adhesion molecules in the embryonic versus the postnatal period. RNA sequencing of single OPCs showed that embryonic synapse-lacking OPCs are found in clusters distinct from postnatal OPCs and with similarities to early progenitors. Furthermore, single-cell transcriptomics demonstrated that synaptic genes are transiently expressed only by postnatal OPCs until they start to differentiate. Taken together, our results indicate that embryonic OPCs represent a unique developmental stage biologically resembling postnatal OPCs but without synaptic input and a transcriptional signature in the continuum between OPCs and neural precursors.

The temporal profile of activity-dependent presynaptic phospho-signalling reveals long-lasting patterns of poststimulus regulation.

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.

Reactive Oxygen Species (ROS)-Activated ATM-Dependent Phosphorylation of Cytoplasmic Substrates Identified by Large-Scale Phosphoproteomics Screen.

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.

Analysis of steric effects in DamID profiling of transcription factor target genes.

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.

CompGO: an R package for comparing and visualizing Gene Ontology enrichment differences between DNA binding experiments.

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.

Multiomics of early epileptogenesis in mice reveals phosphorylation and dephosphorylation-directed growth and synaptic weakening.

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.

Highly dynamic inflammatory and excitability transcriptional profiles in hippocampal CA1 following status epilepticus.

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.

Tau forms synaptic nano-biomolecular condensates controlling the dynamic clustering of recycling synaptic vesicles.

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.

Phosphoinositide 3-kinase (PI3K(p110alpha)) directly regulates key components of the Z-disc and cardiac structure.

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.

A presynaptic phosphosignaling hub for lasting homeostatic plasticity.

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 RIMS1045E. 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.

An Analytically and Diagnostically Sensitive RNA Extraction and RT-qPCR Protocol for Peripheral Blood Mononuclear Cells.

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 × 104 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.

A systematic approach to simultaneously evaluate safety, immunogenicity, and efficacy of novel tuberculosis vaccination strategies.

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.

Uptake of Schistosoma mansoni extracellular vesicles by human endothelial and monocytic cell lines and impact on vascular endothelial cell gene expression.

The ability of the parasitic blood fluke Schistosoma mansoni and other parasitic helminths to manipulate host biology is well recognised, but the mechanisms that underpin these phenomena are not well understood. An emerging paradigm is that helminths transfer their biological cargo to host cells by secretion of extracellular vesicles (EVs). Herein, we show that two populations of S. mansoni secreted EVs - exosome-like vesicles (ELVs) and microvesicles (MVs) - are actively internalised in two distinct human cell lines that reflect the resident cell types encountered by the parasite in vivo: human umbilical vein endothelial cells (HUVECs) and THP-1 monocytes. RNA-sequencing of HUVECs co-cultured with S. mansoni ELVs compared with untreated HUVECs revealed differential expression of genes associated with intravascular parasitism, including vascular endothelial contraction, coagulation, arachidonic acid metabolism and immune cell trafficking and signalling. Finally, we show that antibodies raised against recombinant tetraspanin (TSP) proteins from the surface of S. mansoni EVs significantly blocked EV uptake by both HUVECs and THP-1 monocytes whereas pre-immunisation antibodies did not. To our knowledge, this is the first evidence demonstrating the internalisation of secreted EVs from any helminth into vascular endothelial cells, providing novel insight into the potential mechanisms underlying host-schistosome interactions. The ability of anti-TSP antibodies to block vesicle uptake by host target cells further supports the potential of TSPs as promising antigens for an anti-fluke vaccine. It also suggests a potential mechanism whereby the current candidate human schistosomiasis vaccine, Sm-TSP-2, exerts its protective effect in animal models.

TWIST1 Homodimers and Heterodimers Orchestrate Lineage-Specific Differentiation.

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.

consensusDE: an R package for assessing consensus of multiple RNA-seq algorithms with RUV correction.

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/.

Statistical Analysis of ATM-Dependent Signaling in Quantitative Mass Spectrometry Phosphoproteomics.

Ataxia-telangiectasia mutated (ATM) is a serine/threonine protein kinase, which when perturbed is associated with modified protein signaling that ultimately leads to a range of neurological and DNA repair defects. Recent advances in phospho-proteomics coupled with high-resolution mass-spectrometry provide new opportunities to dissect signaling pathways that ATM utilize under a number of conditions. This chapter begins by providing a brief overview of ATM function, its various regulatory roles and then leads into a workflow focused on the use of the statistical programming language R, together with code, for the identification of ATM-dependent substrates in the cytoplasm. This chapter cannot cover statistical properties in depth nor the range of possible methods in great detail, but instead aims to equip researchers with a set of tools to perform analysis between two conditions through examples with R functions.

Progesterone signalling in broiler skeletal muscle is associated with divergent feed efficiency.

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.