Inverse regulation of the cytosolic Ca²⁺ buffer parvalbumin and mitochondrial volume in muscle cells via SIRT1/PGC-1α axis

What Is Parvalbumin for?

Parvalbumin (PA) is a small, acidic, mostly cytosolic Ca2+-binding protein of the EF-hand superfamily. Structural and physical properties of PA are well studied but recently two highly conserved structural motifs consisting of three amino acids each (clusters I and II), which contribute to the hydrophobic core of the EF-hand domains, have been revealed. Despite several decades of studies, physiological functions of PA are still poorly known. Since no target proteins have been revealed for PA so far, it is believed that PA acts as a slow calcium buffer. Numerous experiments on various muscle systems have shown that PA accelerates the relaxation of fast skeletal muscles. It has been found that oxidation of PA by reactive oxygen species (ROS) is conformation-dependent and one more physiological function of PA in fast muscles could be a protection of these cells from ROS. PA is thought to regulate calcium-dependent metabolic and electric processes within the population of gamma-aminobutyric acid (GABA) neurons. Genetic elimination of PA results in changes in GABAergic synaptic transmission. Mammalian oncomodulin (OM), the β isoform of PA, is expressed mostly in cochlear outer hair cells and in vestibular hair cells. OM knockout mice lose their hearing after 3–4 months. It was suggested that, in sensory cells, OM maintains auditory function, most likely affecting outer hair cells’ motility mechanisms.

Ryanodine receptor leak triggers fiber Ca2+ redistribution to preserve force and elevate basal metabolism in skeletal muscle

Muscle contraction depends on tightly regulated Ca²⁺ release. Aberrant Ca²⁺ leak through ryanodine receptor 1 (RyR1) on the sarcoplasmic reticulum (SR) membrane can lead to heatstroke and malignant hyperthermia (MH) susceptibility, as well as severe myopathy. However, the mechanism by which Ca²⁺ leak drives these pathologies is unknown. Here, we investigate the effects of four mouse genotypes with increasingly severe RyR1 leak in skeletal muscle fibers. We find that RyR1 Ca²⁺ leak initiates a cascade of events that cause precise redistribution of Ca²⁺ among the SR, cytoplasm, and mitochondria through altering the Ca²⁺ permeability of the transverse tubular system membrane. This redistribution of Ca²⁺ allows mice with moderate RyR1 leak to maintain normal function; however, severe RyR1 leak with RYR1 mutations reduces the capacity to generate force. Our results reveal the mechanism underlying force preservation, increased ATP metabolism, and susceptibility to MH in individuals with gain-of-function RYR1 mutations.

Parvalbumin affects skeletal muscle trophism through modulation of mitochondrial calcium uptake

Parvalbumin (PV) is a cytosolic Ca²⁺-binding protein highly expressed in fast skeletal muscle, contributing to an increased relaxation rate. Moreover, PV is an “atrogene” downregulated in most muscle atrophy conditions. Here, we exploit mice lacking PV to explore the link between the two PV functions. Surprisingly, PV ablation partially counteracts muscle loss after denervation. Furthermore, acute PV downregulation is accompanied by hypertrophy and upregulation by atrophy. PV ablation has a minor impact on sarcoplasmic reticulum but is associated with increased mitochondrial Ca²⁺ uptake, mitochondrial size and number, and contacts with Ca²⁺ release sites. Mitochondrial calcium uniporter (MCU) silencing abolishes the hypertrophic effect of PV ablation, suggesting that mitochondrial Ca²⁺ uptake is required for hypertrophy. In turn, an increase of mitochondrial Ca²⁺ is required to enhance expression of the pro-hypertrophy gene PGC-1α4, whose silencing blocks hypertrophy due to PV ablation. These results reveal how PV links cytosolic Ca²⁺ control to mitochondrial adaptations, leading to muscle mass regulation.

•

PV is downregulated during skeletal muscle atrophy, and its levels affect trophism

•

Skeletal muscle mitochondria undergo remodeling in PV knockout mice

•

Mitochondria increase cytosolic Ca²⁺ buffer capacity in PV knockout skeletal muscles

•

Increased mitochondrial Ca²⁺ triggers the PGC-1α4 pathway, inducing muscle growth

NeuroHeal Improves Muscle Regeneration after Injury

Musculoskeletal injuries represent a challenging medical problem. Although the skeletal muscle is able to regenerate and recover after injury, the process engaged with conservative therapy can be inefficient, leading to a high re-injury rate. In addition, the formation of scar tissue implies an alteration of mechanical properties in muscle. There is still a need for new treatments of the injured muscle. NeuroHeal may be one option. Published studies demonstrated that it reduces muscle atrophy due to denervation and disuse. The main objective of the present work was to assess the potential of NeuroHeal to improve muscle regeneration after traumatic injury. Secondary objectives included characterizing the effect of NeuroHeal treatment on satellite cell biology. We used a rat model of sport-induced injury in the gastrocnemius and analyzed the effects of NeuroHeal on functional recovery by means of electrophysiology and tetanic force analysis. These studies were accompanied by immunohistochemistry of the injured muscle to analyze fibrosis, satellite cell state, and fiber type. In addition, we used an in vitro model to determine the effect of NeuroHeal on myoblast biology and partially decipher its mechanism of action. The results showed that NeuroHeal treatment advanced muscle fiber recovery after injury in a preclinical model of muscle injury, and significantly reduced the formation of scar tissue. In vitro, we observed that NeuroHeal accelerated the formation of myotubes. The results pave the way for novel therapeutic avenues for muscle/tendinous disorders.

The Parvalbumin Hypothesis of Autism Spectrum Disorder

The prevalence of autism spectrum disorder (ASD)-a type of neurodevelopmental disorder-is increasing and is around 2% in North America, Asia, and Europe. Besides the known genetic link, environmental, epigenetic, and metabolic factors have been implicated in ASD etiology. Although highly heterogeneous at the behavioral level, ASD comprises a set of core symptoms including impaired communication and social interaction skills as well as stereotyped and repetitive behaviors. This has led to the suggestion that a large part of the ASD phenotype is caused by changes in a few and common set of signaling pathways, the identification of which is a fundamental aim of autism research. Using advanced bioinformatics tools and the abundantly available genetic data, it is possible to classify the large number of ASD-associated genes according to cellular function and pathways. Cellular processes known to be impaired in ASD include gene regulation, synaptic transmission affecting the excitation/inhibition balance, neuronal Ca2+ signaling, development of short-/long-range connectivity (circuits and networks), and mitochondrial function. Such alterations often occur during early postnatal neurodevelopment. Among the neurons most affected in ASD as well as in schizophrenia are those expressing the Ca2+-binding protein parvalbumin (PV). These mainly inhibitory interneurons present in many different brain regions in humans and rodents are characterized by rapid, non-adaptive firing and have a high energy requirement. PV expression is often reduced at both messenger RNA (mRNA) and protein levels in human ASD brain samples and mouse ASD (and schizophrenia) models. Although the human PVALB gene is not a high-ranking susceptibility/risk gene for either disorder and is currently only listed in the SFARI Gene Archive, we propose and present supporting evidence for the Parvalbumin Hypothesis, which posits that decreased PV level is causally related to the etiology of ASD (and possibly schizophrenia).

Parvalbumin-Deficiency Accelerates the Age-Dependent ROS Production in Pvalb Neurons in vivo: Link to Neurodevelopmental Disorders

In neurodevelopmental disorders (NDDs) including autism spectrum disorder (ASD) and schizophrenia, impairment/malfunctioning of a subpopulation of interneurons expressing the calcium-binding protein parvalbumin (PV) -here termed Pvalb neurons- has gradually emerged as a possible cause. These neurons may represent a hub or point-of-convergence in the etiology of NDD. Increased oxidative stress associated with mitochondria impairment in Pvalb neurons is discussed as an essential step in schizophrenia etiology. Since PV downregulation is a common finding in ASD and schizophrenia individuals and PV-deficient (PV-/-) mice show a strong ASD-like behavior phenotype, we investigated the putative link between PV expression, alterations in mitochondria and oxidative stress. In a longitudinal study with 1, 3, and 6-months old PV-/- and wild type mice, oxidative stress was investigated in 9 Pvalb neuron subpopulations in the hippocampus, striatum, somatosensory cortex, medial prefrontal cortex, thalamic reticular nucleus (TRN) and cerebellum. In Pvalb neuron somata in the striatum and TRN, we additionally determined mitochondria volume and distribution at these three time points. In all Pvalb neuron subpopulations, we observed an age-dependent increase in oxidative stress and the increase strongly correlated with PV expression levels, but not with mitochondria density in these Pvalb neurons. Moreover, oxidative stress was elevated in Pvalb neurons of PV-/- mice and the magnitude of the effect was again correlated with PV expression levels in the corresponding wild type Pvalb neuron subpopulations. The PV-dependent effect was insignificant at 1 month and relative differences between WT and PV-/- Pvalb neurons were largest at 3 months. Besides the increase in mitochondria volume in PV's absence in TRN and striatal PV-/- Pvalb neurons fully present already at 1 month, we observed a redistribution of mitochondria from the perinuclear region toward the plasma membrane at all time points. We suggest that in absence of PV, slow Ca2+ buffering normally exerted by PV is compensated by a (mal)adaptive, mostly sub-plasmalemmal increase in mitochondria resulting in increased oxidative stress observed in 3- and 6-months old mice. Since PV-/- mice display core ASD-like symptoms already at 1 month, oxidative stress in Pvalb neurons is not a likely cause for their ASD-related behavior observed at this age.

Parvalbumin immunohistochemical expression in the spectrum of perivascular epithelioid cell (PEC) lesions of the kidney

Parvalbumin is a cytosolic calcium-binding protein expressed in the distal convoluted tubule of the renal nephron. Among epithelial renal tumors, the reactivity for parvalbumin is observed in chromophobe renal cell carcinomas and frequently in oncocytomas. On the other hand, there are no data available on parvalbumin expression in the mesenchymal tumors of the kidney. Therefore, the purpose of this study was to evaluate the expression of parvalbumin in the spectrum of PEC (perivascular epithelioid cells) lesions of the kidney. Sixty-six PEC lesions (37 classic angiomyolipomas, 10 microscopic angiomyolipomas, 7 epithelioid angiomyolipomas/pure epithelioid PEComas, 5 leiomyoma-like angiomyolipomas, 3 lipoma-like angiomyolipomas, 2 intraglomerular lesions, 1 angiomyolipoma with epithelial cysts (AMLEC), and 1 sclerosing angiomyolipoma) were immunohistochemically stained with parvalbumin. Overall, parvalbumin immunostain was found in fifty-six PEC lesions (85%) and absent in the remaining ten cases (15%). Classic angiomyolipomas were positive in almost all cases (97%). Intraglomerular lesions and AMLEC showed parvalbumin immunolabeling as well. None of the 7 epithelioid angiomyolipomas/pure epithelioid PEComas or the only sclerosing angiomyolipoma expressed parvalbumin. In conclusion, we demonstrated the immunolabeling of parvalbumin in almost all PEC lesions of the kidney, but not in the epithelioid angiomyolipoma/pure epithelioid PEComa. This finding could shed light on some biological characteristics observed in the PEC lesions such as the plasticity of their cellular component. Moreover, parvalbumin may be another useful tool in the differential diagnosis among epithelioid angiomyolipoma/pure epithelioid PEComa with other renal eosinophilic tumors, such as oncocytoma and chromophobe renal cell carcinoma.

Cytosolic Ca2+ Buffers Are Inherently Ca2+ Signal Modulators

For precisely regulating intracellular Ca²⁺ signals in a time- and space-dependent manner, cells make use of various components of the "Ca²⁺ signaling toolkit," including Ca²⁺ entry and Ca²⁺ extrusion systems. A class of cytosolic Ca²⁺-binding proteins termed Ca²⁺ buffers serves as modulators of such, mostly short-lived Ca²⁺ signals. Prototypical Ca²⁺ buffers include parvalbumins (α and β isoforms), calbindin-D9k, calbindin-D28k, and calretinin. Although initially considered to function as pure Ca²⁺ buffers, that is, as intracellular Ca²⁺ signal modulators controlling the shape (amplitude, decay, spread) of Ca²⁺ signals, evidence has accumulated that calbindin-D28k and calretinin have additional Ca²⁺ sensor functions. These other functions are brought about by direct interactions with target proteins, thereby modulating their targets' function/activity. Dysregulation of Ca²⁺ buffer expression is associated with several neurologic/neurodevelopmental disorders including autism spectrum disorder (ASD) and schizophrenia. In some cases, the presence of these proteins is presumed to confer a neuroprotective effect, as evidenced in animal models of Parkinson's or Alzheimer's disease.

Parvalbumin expression in oligodendrocyte-like CG4 cells causes a reduction in mitochondrial volume, attenuation in reactive oxygen species production and a decrease in cell processes' length and branching

Forebrain glial cells - ependymal cells and astrocytes -acquire upon injury- a "reactive" phenotype associated with parvalbumin (PV) upregulation. Since free radicals, e.g. reactive oxygen species (ROS) play a role in the pathogenesis of multiple sclerosis, and that PV-upregulation in glial cells is inversely correlated with the level of oxidative stress, we hypothesized that PV-upregulation might also protect oligodendrocytes by decreasing ROS production. Lentiviral transduction techniques allowed for PV overexpression in CG4 oligodendrocyte progenitor cells (OPCs). Depending on the growth medium CG4 cells can be maintained in an OPC-like state, or induced to differentiate into an oligodendrocyte (OLG)-like phenotype. While increased levels of PV had no effect on cell proliferation and invasiveness in vitro, PV decreased the mitochondria volume in CG4 cell bodies, as well as the mitochondrial density in CG4 processes in both OPC-like and OLG-like states. In line with the PV-induced global decrease in mitochondrial volume, elevated PV levels reduced transcript levels of mitochondrial transcription factors involved in mitochondria biogenesis. In differentiated PV-overexpressing CG4 cells with a decreased mitochondrial volume, UV-induced ROS production was lower than in control CG4 cells hinting towards a possible role of PV in counteracting oxidative stress. Unexpectedly, PV also decreased the length of processes in undifferentiated CG4 cells and moreover diminished branching of differentiated CG4 cell processes, strongly correlated with the decreased density of mitochondria in CG4 cell processes. Thus besides conferring a protective role against oxidative stress, PV in a cell autonomous fashion additionally affects process' growth and branching in CG4 cells.

Parvalbumin alters mitochondrial dynamics and affects cell morphology

The Ca2+-binding protein parvalbumin (PV) and mitochondria play important roles in Ca2+ signaling, buffering and sequestration. Antagonistic regulation of PV and mitochondrial volume is observed in in vitro and in vivo model systems. Changes in mitochondrial morphology, mitochondrial volume and dynamics (fusion, fission, mitophagy) resulting from modulation of PV were investigated in MDCK epithelial cells with stable overexpression/downregulation of PV. Increased PV levels resulted in smaller, roundish cells and shorter mitochondria, the latter phenomenon related to reduced fusion rates and decreased expression of genes involved in mitochondrial fusion. PV-overexpressing cells displayed increased mitophagy, a likely cause for the decreased mitochondrial volumes and the smaller overall cell size. Cells showed lower mobility in vitro, paralleled by reduced protrusions. Constitutive PV down-regulation in PV-overexpressing cells reverted mitochondrial morphology and fractional volume to the state present in control MDCK cells, resulting from increased mitochondrial movement and augmented fusion rates. PV-modulated, bi-directional and reversible mitochondrial dynamics are key to regulation of mitochondrial volume.

Six1 homeoprotein drives myofiber type IIA specialization in soleus muscle

BACKGROUND

Adult skeletal muscles are composed of slow and fast myofiber subtypes which each express selective genes required for their specific contractile and metabolic activity. Six homeoproteins are transcription factors regulating muscle cell fate through activation of myogenic regulatory factors and driving fast-type gene expression during embryogenesis.

RESULTS

We show here that Six1 protein accumulates more robustly in the nuclei of adult fast-type muscles than in adult slow-type muscles, this specific enrichment takes place during perinatal growth. Deletion of Six1 in soleus impaired fast-type myofiber specialization during perinatal development, resulting in a slow phenotype and a complete lack of Myosin heavy chain 2A (MyHCIIA) expression. Global transcriptomic analysis of wild-type and Six1 mutant myofibers identified the gene networks controlled by Six1 in adult soleus muscle. This analysis showed that Six1 is required for the expression of numerous genes encoding fast-type sarcomeric proteins, glycolytic enzymes and controlling intracellular calcium homeostasis. Parvalbumin, a key player of calcium buffering, in particular, is a direct target of Six1 in the adult myofiber.

CONCLUSIONS

This analysis revealed that Six1 controls distinct aspects of adult muscle physiology in vivo, and acts as a main determinant of fast-fiber type acquisition and maintenance.

Antagonistic Regulation of Parvalbumin Expression and Mitochondrial Calcium Handling Capacity in Renal Epithelial Cells

Parvalbumin (PV) is a cytosolic Ca²⁺-binding protein acting as a slow-onset Ca²⁺ buffer modulating the shape of Ca²⁺ transients in fast-twitch muscles and a subpopulation of neurons. PV is also expressed in non-excitable cells including distal convoluted tubule (DCT) cells of the kidney, where it might act as an intracellular Ca²⁺ shuttle facilitating transcellular Ca²⁺ resorption. In excitable cells, upregulation of mitochondria in “PV-ergic” cells in PV-/- mice appears to be a general hallmark, evidenced in fast-twitch muscles and cerebellar Purkinje cells. Using Gene Chip Arrays and qRT-PCR, we identified differentially expressed genes in the DCT of PV-/- mice. With a focus on genes implicated in mitochondrial Ca²⁺ transport and membrane potential, uncoupling protein 2 (Ucp2), mitocalcin (Efhd1), mitochondrial calcium uptake 1 (Micu1), mitochondrial calcium uniporter (Mcu), mitochondrial calcium uniporter regulator 1 (Mcur1), cytochrome c oxidase subunit 1 (COX1), and ATP synthase subunit β (Atp5b) were found to be up-upregulated. At the protein level, COX1 was increased by 31 ± 7%, while ATP-synthase subunit β was unchanged. This suggested that these mitochondria were better suited to uphold the electrochemical potential across the mitochondrial membrane, necessary for mitochondrial Ca²⁺ uptake. Ectopic expression of PV in PV-negative Madin-Darby canine kidney (MDCK) cells decreased COX1 and concomitantly mitochondrial volume, while ATP synthase subunit β levels remained unaffected. Suppression of PV by shRNA in PV-expressing MDCK cells led subsequently to an increase in COX1 expression. The collapsing of the mitochondrial membrane potential by the uncoupler CCCP occurred at lower concentrations in PV-expressing MDCK cells than in control cells. In support, a reduction of the relative mitochondrial mass was observed in PV-expressing MDCK cells. Deregulation of the cytoplasmic Ca²⁺ buffer PV in kidney cells was counterbalanced in vivo and in vitro by adjusting the relative mitochondrial volume and modifying the mitochondrial protein composition conceivably to increase their Ca²⁺-buffering/sequestration capacity.

Routes of Ca2+ Shuttling during Ca2+ Oscillations: FOCUS ON THE ROLE OF MITOCHONDRIAL Ca2+ HANDLING AND CYTOSOLIC Ca2+ BUFFERS

In some cell types, Ca²⁺ oscillations are strictly dependent on Ca²⁺ influx across the plasma membrane, whereas in others, oscillations also persist in the absence of Ca²⁺ influx. We observed that, in primary mesothelial cells, the plasmalemmal Ca²⁺ influx played a pivotal role. However, when the Ca²⁺ transport across the plasma membrane by the “lanthanum insulation method” was blocked prior to the induction of the serum-induced Ca²⁺ oscillations, mitochondrial Ca²⁺ transport was found to be able to substitute for the plasmalemmal Ca²⁺ exchange function, thus rendering the oscillations independent of extracellular Ca²⁺. However, in a physiological situation, the Ca²⁺-buffering capacity of mitochondria was found not to be essential for Ca²⁺ oscillations. Moreover, brief spontaneous Ca²⁺ changes were observed in the mitochondrial Ca²⁺ concentration without apparent changes in the cytosolic Ca²⁺ concentration, indicating the presence of a mitochondrial autonomous Ca²⁺ signaling mechanism. In the presence of calretinin, a Ca²⁺-buffering protein, the amplitude of cytosolic spikes during oscillations was decreased, and the amount of Ca²⁺ ions taken up by mitochondria was reduced. Thus, the increased calretinin expression observed in mesothelioma cells and in certain colon cancer might be correlated to the increased resistance of these tumor cells to proapoptotic/pronecrotic signals. We identified and characterized (experimentally and by modeling) three Ca²⁺ shuttling pathways in primary mesothelial cells during Ca²⁺ oscillations: Ca²⁺ shuttled between (i) the endoplasmic reticulum (ER) and mitochondria, (ii) the ER and the extracellular space, and (iii) the ER and cytoplasmic Ca²⁺ buffers.

Contractile and metabolic properties of engineered skeletal muscle derived from slow and fast phenotype mouse muscle

Satellite cells derived from fast and slow muscles have been shown to adopt contractile and metabolic properties of their parent muscle. Mouse muscle shows less distinctive fiber-type profiles than rat or rabbit muscle. Therefore, in this study we sought to determine whether three-dimensional muscle constructs engineered from slow soleus (SOL) and fast tibialis anterior (TA) from mice would adopt the contractile and metabolic properties of their parent muscle. Time-to-peak tension (TPT) and half-relaxation time (1/2RT) was significantly slower in SOL constructs. In agreement with TPT, TA constructs contained significantly higher levels of fast myosin heavy chain (MHC) and fast troponin C, I, and T isoforms. Fast SERCA protein, both slow and fast calsequestrin isoforms and parvalbumin were found at higher levels in TA constructs. SOL constructs were more fatigue resistant and contained higher levels of the mitochondrial proteins SDH and ATP synthase and the fatty acid transporter CPT-1. SOL constructs contained lower levels of the glycolytic enzyme phosphofructokinase but higher levels of the β-oxidation enzymes LCAD and VLCAD suggesting greater fat oxidation. Despite no changes in PGC-1α protein, SOL constructs contained higher levels of SIRT1 and PRC. TA constructs contained higher levels of the slow-fiber program repressor SOX6 and the six transcriptional complex (STC) proteins Eya1 and Six4 which may underlie the higher in fast-fiber and lower slow-fiber program proteins. Overall, we have found that muscles engineered from predominantly slow and fast mouse muscle retain contractile and metabolic properties of their native muscle.

Comparison of growth-related traits and gene expression profiles between the offspring of neomale (XX) and normal male (XY) rainbow trout

All-female lines of fish are created by crossing sex reversed (XX genotype) males with normal females. All-female lines avoid the deleterious phenotypic effects that are typical of precocious maturation in males. To determine whether all-female and mixed sex populations of rainbow trout (Oncorhynchus mykiss) differ in performance, we compared the growth and gene expression profiles in progeny groups produced by crossing a XX male and a XY male to the same five females. Body weight and length were measured in the resulting all-female (XX) and mixed sex (XX/XY) offspring groups. Microarray experiments with liver and white muscle were used to determine if the gene expression profiles of large and small XX offspring differ from those in large and small XX/XY offspring. We detected no significant differences in body length and weight between offspring groups but XX offspring were significantly less variable in the value of these traits. A large number of upregulated genes were shared between the large XX and large XX/XY offspring; the small XX and small XX/XY offspring also shared similar expression profiles. No GO category differences were seen in the liver or between the large XX and large XX/XY offspring in the muscle. The greatest differences between the small XX and small XX/XY offspring were in the genes assigned to the "small molecule metabolic process" and "cellular metabolic process" GO level 3 categories. Similarly, genes within these categories as well as the category "macromolecule metabolic process" were more highly expressed in small compared to large XX fish.

Differential gene expression in small and large rainbow trout derived from two seasonal spawning groups

Background

Growth in fishes is regulated via many environmental and physiological factors and is shaped by the genetic background of each individual. Previous microarray studies of salmonid growth have examined fish experiencing either muscle wastage or accelerated growth patterns following refeeding, or the influence of growth hormone and transgenesis. This study determines the gene expression profiles of genetically unmanipulated large and small fish from a domesticated salmonid strain reared on a typical feeding regime. Gene expression profiles of white muscle and liver from rainbow trout (Oncorhynchus mykiss) from two seasonal spawning groups (September and December lots) within a single strain were examined when the fish were 15 months of age to assess the influence of season (late fall vs. onset of spring) and body size (large vs. small).

Results

Although IGFBP1 gene expression was up-regulated in the livers of small fish in both seasonal lots, few expression differences were detected in the liver overall. Faster growing Dec. fish showed a greater number of differences in white muscle expression compared to Sept. fish. Significant differences in the GO Generic Level 3 categories ‘response to external stimulus’, ‘establishment of localization’, and ‘response to stress’ were detected in white muscle tissue between large and small fish. Larger fish showed up-regulation of cytoskeletal component genes while many genes related to myofibril components of muscle tissue were up-regulated in small fish. Most of the genes up-regulated in large fish within the ‘response to stress’ category are involved in immunity while in small fish most of these gene functions are related to apoptosis.

Conclusions

A higher proportion of genes in white muscle compared to liver showed similar patterns of up- or down-regulation within the same size class across seasons supporting their utility as biomarkers for growth in rainbow trout. Differences between large and small Sept. fish in the ‘response to stress’ and ‘response to external stimulus’ categories for white muscle tissue, suggests that smaller fish have a greater inability to handle stress compared to the large fish. Sampling season had a significant impact on the expression of genes related to the growth process in rainbow trout.