Gene expression profiling: insights into skeletal muscle growth and development

Differential isoform expression and selective muscle involvement in muscular dystrophies

Despite the expression of the mutated gene in all muscles, selective muscles are involved in genetic muscular dystrophies. Different muscular dystrophies show characteristic patterns of fatty degenerative changes by muscle imaging, even to the extent that the patterns have been used for diagnostic purposes. However, the underlying molecular mechanisms explaining the selective involvement of muscles are not known. To test the hypothesis that different muscles may express variable amounts of different isoforms of muscle genes, we applied a custom-designed exon microarray containing probes for 57 muscle-specific genes to assay the transcriptional profiles in sets of human adult lower limb skeletal muscles. Quantitative real-time PCR and whole transcriptome sequencing were used to further analyze the results. Our results demonstrate significant variations in isoform and gene expression levels in anatomically different muscles. Comparison of the known patterns of selective involvement of certain muscles in two autosomal dominant titinopathies and one autosomal dominant myosinopathy, with the isoform and gene expression results, shows a correlation between the specific muscles involved and significant differences in the level of expression of the affected gene and exons in these same muscles compared with some other selected muscles. Our results suggest that differential expression levels of muscle genes and isoforms are one determinant in the selectivity of muscle involvement in muscular dystrophies.

Sequencing and Characterization of Divergent Marbling Levels in the Beef Cattle (Longissimus dorsi Muscle) Transcriptome

Marbling is an important trait regarding the quality of beef. Analysis of beef cattle transcriptome and its expression profile data are essential to extend the genetic information resources and would support further studies on beef cattle. RNA sequencing was performed in beef cattle using the Illumina High-Seq2000 platform. Approximately 251.58 million clean reads were generated from a high marbling (H) group and low marbling (L) group. Approximately 80.12% of the 19,994 bovine genes (protein coding) were detected in all samples, and 749 genes exhibited differential expression between the H and L groups based on fold change (>1.5-fold, p<0.05). Multiple gene ontology terms and biological pathways were found significantly enriched among the differentially expressed genes. The transcriptome data will facilitate future functional studies on marbling formation in beef cattle and may be applied to improve breeding programs for cattle and closely related mammals.

Effect of sexual maturation on muscle gene expression of rainbow trout: RNA-Seq approach

Muscle degradation occurs as a response to various physiological states that are regulated by specific molecular mechanisms. Previously, we characterized the metabolic changes of muscle deterioration of the female rainbow trout at full sexual maturity and spawning (Salem et al., Physiol. Genomics 2006;28:33–45; J. Proteomics 2010;73:778–789). Muscle deterioration in this model represents nutrient mobilization as a response to the energetic overdemands of the egg/ovarian growth phase. Our recent studies showed that most of the changes in muscle growth and quality start 2–3 months before spawning. Gravid fish exhibited reduced intramuscular fat that is lower in saturated and monounsaturated fatty acids and higher in polyunsaturated fatty acids compared to sterile fish. In this study, RNA-Seq was used to explain the mechanisms underlying changes during this phase of sexual maturity. Furthermore, to minimize changes due to nutrient deficits, fish were fed on a high-plane of nutrition. The RNA-Seq technique identified a gene expression signature that is consistent with metabolic changes of gravid fish. Gravid fish exhibited increased abundance of transcripts in metabolic pathways of fatty acid degradation and up-regulated expression of genes involved in biosynthesis of unsaturated fatty acids. In addition, increased expression of genes involved in the citric acid cycle and oxidative phosphorylation was observed for gravid fish. This muscle transcriptomic signature of fish fed on a high nutritional plane is quite distinct from that previously described for fish at terminal stages of maturity and suggest that female rainbow trout approaching spawning, on high nutritional planes, likely mobilize intramuscular fat rather than protein to support gonadal maturation.

Underlying functional genomics of fat deposition in adipose tissue

The objective of this study was to gain insight into the underlying mechanisms of fat deposition. Two sheep breeds with large fat-tail (Lori-Bakhtiari) and with thin-tail (Zel) were used as models. To determine important and key candidate lipid metabolism related genes, comparative genomic approaches were employed. Gene expression profiles of adipose tissues were analyzed in human, pig, and cattle by express sequence tag (EST) analysis. EST analysis determined 65, 102 and 125 transcripts in human, pig and cattle respectively with at least 10 fold over-expression in the adipose tissue. Based on our comparative functional genomic analysis, seven genes were more abundant and common in investigated mammalian adipose tissues promising a conserved novel gene network in mammalian lipid metabolism. The candidate genes including fatty acid binding protein 4 (FABP4), fatty acid synthase (FASN), Stearoyl-CoA desaturase (SCD) and Lipoprotein lipase (LPL) were selected for further gene expression investigation within two sheep breeds. The real time PCR results showed that among the genes tested, FABP4 was expressed at higher levels than the others. The expression of FABP4 was significantly higher in the fat-tail of Lori-Bakhtiari than in the fat-tail and visceral adipose tissues of Zel (P<0.05). The findings suggest that the FABP4 gene expression in the fat-tail is an important index of fat deposition.

Variability-specific differential gene expression across reproductive stages in sows

Differential gene expression analyses typically focus on departures across mathematical expectations (i.e. mean) from two or more groups of microarrays, without considering alternative patterns of departure. Nevertheless, recent studies in humans and great apes have suggested that differential gene expression could also be characterized in terms of heterogeneous dispersion patterns. This must be viewed as a very interesting genetic phenomenon clearly linked to the regulation mechanisms of gene transcription. Unfortunately, we completely lack information about the incidence and relevance of dispersion-specific differential gene expression in livestock species, although a specific Bayes factor (BF) for testing this kind of differential gene expression (i.e. within-probe heteroskedasticity) has been recently developed. Within this context, our main objective was to characterize the incidence of dispersion-specific differential gene expression in pigs and, if possible, providing the first evidence of this phenomenon in a livestock species. We evaluated dispersion-specific differential gene expression on ovary, uterus and hypophysis samples from 22 F2 Iberian × Meishan sows, where a total of 15,252 probes were interrogated. For each tissue, heteroskedasticity of probe-specific residual variances was evaluated by three pairwise comparisons involving three physiological stages, that is, heat, 15 days of pregnancy and 45 days of pregnancy. Between 2.9% and 37.4% of the analyzed probes provided statistical evidence of within-tissue across-physiological stages dispersion-specific differential gene expression (BF >1), and between 0.1% and 3.0% of them reported decisive evidence (BF >100). It is important to highlight that <8% of the heteroskedastic probes were also linked to differential gene expression in terms of departures among the probe-specific mathematical expectation of each physiological stage. This discarded the disturbance of scale effects in a high percentage of probes and suggested that probe-specific heteroskedasticity must be viewed as an independent phenomenon within the context of differential gene expression. As a whole, our results report a remarkable incidence of dispersion-specific differential gene expression across the whole genome of the pig, establishing a very interesting starting point for further studies focused on deciphering the genetic mechanisms underlying heteroskedasticity.

Characterization of the longissimus lumborum transcriptome response to adding propionate to the diet of growing Angus beef steers

Development of management paradigms that enhance the rate of gain and qualitative characteristics of beef carcass development has the potential to impact production and nutrient use efficiency but also mitigate losses to the environment. We used eight Black Angus beef steers (272.5 ± 17.6 kg initial body wt) fed a forage-based pelleted diet alone (n = 4) or supplemented with sodium propionate included (n = 4) for 42 days. High-quality RNA was extracted from the longissimus lumborum and subjected to transcriptome sequencing using RNA-seq technology. Trimmed reads were aligned to the bovine reference genome (Btau4.0, release 63) and uniquely mapped reads from control and propionate treatment groups were subject to further analysis using edgeR. Candidates were filtered to account for multiple testing and differentially expressed genes (153 at a false discovery rate of <5%) were analyzed using Gene Ontology (GO) analysis (GOseq) to select terms where enrichment had occurred. Significant GO terms included regulation of cholesterol transport, regulation of sterol transport, and cellular modified amino acid metabolic process. Furthermore, the top four identified gene networks included lipid metabolism, small molecule biochemistry, carbohydrate metabolism, and molecular transport-related categories. Notably, changes in lipid metabolism specific genes reflect both increased oxidative and lipid synthetic capacities. Metabolism-related gene changes are reflective of expected enhancements in lean tissue accretion patterns exhibited in steers where high ruminal propionate relative to other short chain fatty acids is observed. Propionate feeding induced increased N retention in rapidly growing Angus cattle, and the observed alterations in LL tissue lipid metabolism-related gene networks are consistent with enhanced cell formation and function (protein synthesis, and lipogenic vs. lipolytic activities).

Episcopic 3D Imaging Methods: Tools for Researching Gene Function

This work aims at describing episcopic 3D imaging methods and at discussing how these methods can contribute to researching the genetic mechanisms driving embryogenesis and tissue remodelling, and the genesis of pathologies. Several episcopic 3D imaging methods exist. The most advanced are capable of generating high-resolution volume data (voxel sizes from 0.5x0.5x1 µm upwards) of small to large embryos of model organisms and tissue samples. Beside anatomy and tissue architecture, gene expression and gene product patterns can be three dimensionally analyzed in their precise anatomical and histological context with the aid of whole mount in situ hybridization or whole mount immunohistochemical staining techniques. Episcopic 3D imaging techniques were and are employed for analyzing the precise morphological phenotype of experimentally malformed, randomly produced, or genetically engineered embryos of biomedical model organisms. It has been shown that episcopic 3D imaging also fits for describing the spatial distribution of genes and gene products during embryogenesis, and that it can be used for analyzing tissue samples of adult model animals and humans. The latter offers the possibility to use episcopic 3D imaging techniques for researching the causality and treatment of pathologies or for staging cancer. Such applications, however, are not yet routine and currently only preliminary results are available. We conclude that, although episcopic 3D imaging is in its very beginnings, it represents an upcoming methodology, which in short terms will become an indispensable tool for researching the genetic regulation of embryo development as well as the genesis of malformations and diseases.

Molecular genetic studies of gene identification for sarcopenia

Sarcopenia, which is characterized by a progressive decrease of skeletal muscle mass and function with aging, is closely related to several common diseases (such as cardiovascular and airway diseases) and functional impairment/disability. Strong genetic determination has been reported for muscle mass and muscle strength, two most commonly recognized and studied risk phenotypes for sarcopenia, with heritability ranging from 30 to 85% for muscle strength and 45-90% for muscle mass. Sarcopenia has been the subject of increasing genetic research over the past decade. This review is designed to comprehensively summarize the most important and representative molecular genetic studies designed to identify genetic factors associated with sarcopenia. We have methodically reviewed whole-genome linkage studies in humans, quantitative trait loci mapping in animal models, candidate gene association studies, newly reported genome-wide association studies, DNA microarrays and microRNA studies of sarcopenia or related skeletal muscle phenotypes. The major results of each study are tabulated for easy comparison and reference. The findings of representative studies are discussed with respect to their influence on our present understanding of the genetics of sarcopenia. This is a comprehensive review of molecular genetic studies of gene identification for sarcopenia, and an overarching theme for this review is that the currently accumulating results are tentative and occasionally inconsistent and should be interpreted with caution pending further investigation. Consequently, this overview should enhance recognition of the need to validate/replicate the genetic variants underlying sarcopenia in large human cohorts and animal. We believe that further progress in understanding the genetic etiology of sarcopenia will provide valuable insights into important fundamental biological mechanisms underlying muscle physiology that will ultimately lead to improved ability to recognize individuals at risk for developing sarcopenia and our ability to treat this debilitating condition.

Muscle transcriptome profiling in divergent phenotype swine breeds during growth using microarray and RT-PCR tools

Using an array consisting of 10 665 70-mer oligonucleotide probes, the longissimus dorsi muscle tissue expression during growth in nine pigs belonging to Casertana (CT), an autochthonous breed characterized by slow growth and a massive accumulation of backfat, was compared with that of two cosmopolitan breeds, Large White (LW) and a crossbreed (CB; Duroc × Landrace × Large White). The results were validated by real-time PCR. All animals were of the same age and were raised under the same environmental conditions. Muscle tissues were collected at 3, 6, 9 and 11 months of age, and a total of 173 genes showed significant differential expression between CT and the cosmopolitan genetic types at 3 months of age. Time series cluster analysis indicated that the CT breed had a different pattern of gene expression compared with that of the LW and the CB. Four of the eight clusters highlighted the gene differences between CT and the other two breeds, which were further supported by statistical analyses: clusters 4 and 5 contained a total of 71 genes that were underexpressed at 3 months of age, and cluster 3 and cluster 7 included 28 and 42 genes respectively that were overexpressed at 3 months of age. As expected, differentially expressed genes belonged to the category of genes coding for contractile fibres and transcription factors involved in muscle development and differentiation. These findings highlight muscle expression genes during pig growth and are useful to understand the genetic meaning of the different developmental rates.

Evaluation of post mortem stability of porcine skeletal muscle RNA

The objective of this study was to evaluate the effect of postmortem times on the quality of porcine skeletal muscle total RNA in order to consider the possibility to use postmortem material for gene expression analysis. Samples of Musculus semimembranosus were collected at 20min, 2h, 6h, 24h and 48h postmortem from the left legs of four commercial heavy pigs. Total RNA was analysed by agarose gel electrophoresis stained with ethidium bromide and by microfluidic capillary electrophoresis on an Agilent 2100 Bioanalyzer instrument obtaining 28S:18S rRNA peak ratios and RIN values. The average RIN values of the analysed samples were 7.45±0.13, 7.43±0.15, 7.45±0.10, 7.33±0.15 and 3.95±0.58 for the same postmortem times, respectively, indicating that RNA degradation was present at 48h postmortem. In a similar experiment, carried out by other authors on beef cattle muscle total RNA extracted at different postmortem times, RNA was stable up to 8days after death as indicated by intact 28S and 18S rRNA bands. Thus, differences among species or other environmental factors might affect the level of RNA degradation. In the porcine postmortem samples, qualitative assessment of GAPDH transcripts by PCR amplification of different cDNA fragments indicated that postmortem stages did not affect the possibility of analysing this housekeeping gene. Thus, postmortem porcine skeletal muscle can be an useful tissue to obtain gene expression based information.

The Peter Pan paradigm

Genetic and environmental agents that disrupt organogenesis are numerous and well described. Less well established, however, is the role of delay in the developmental processes that yield functionally immature tissues at birth. Evidence is mounting that organs do not continue to develop postnatally in the context of these organogenesis insults, condemning the patient to utilize under-developed tissues for adult processes. These poorly differentiated organs may appear histologically normal at birth but with age may deteriorate revealing progressive or adult-onset pathology. The genetic and molecular underpinning of the proposed paradigm reveals the need for a comprehensive systems biology approach to evaluate the role of maternal-fetal environment on organogenesis."You may delay, but time will not" Benjamin Franklin, USA Founding Father.

Transcriptomic approach to the study of osmoregulation in the European eel Anguilla anguilla

In euryhaline teleosts, osmoregulation is a fundamental and dynamic process that is essential for the maintenance of ion and water balance, especially when fish migrate between fresh water (FW) and sea water (SW) environments. The European eel has proved to be an excellent model species to study the molecular and physiological adaptations associated with this osmoregulatory plasticity. The life cycle of the European eel includes two migratory periods, the second being the migration of FW eels back to the Sargasso Sea for reproduction. Various anatomical and physiological changes allow the successful transition to SW. The aim of this study was to use a microarray approach to screen the osmoregulatory tissues of the eel for changes in gene expression following acclimation to SW. Tissues were sampled from fish at selected intervals over a 5-mo period following FW/SW transfer, and RNA was isolated. Suppressive subtractive hybridization was used for enrichment of differentially expressed genes. Microarrays comprising 6,144 cDNAs from brain, gill, intestine, and kidney libraries were hybridized with appropriate targets and analyzed; 229 differentially expressed clones with unique sequences were identified. These clones represented the sequences for 95 known genes, with the remaining sequences (59%) being unknown. The results of the microarray analysis were validated by quantification of 28 differentially expressed genes by Northern blotting. A number of the differentially expressed genes were already known to be involved in osmoregulation, but the functional roles of many others, not normally associated with ion or water transport, remain to be characterized.

Salinity adaptation and gene profiling analysis in the European eel (Anguilla anguilla) using microarray technology

The life cycle of the European eel (Anguilla anguilla) includes two long migratory periods, when the newly hatched leptocephali larvae drift on ocean currents from the Sargasso Sea to the shores of Western Europe and then again up to 30 years later when adult eels swim back to their place of birth for reproductive purposes. Prior to the migration from fresh water (FW) to sea water (SW) adult yellow eels undergo various anatomical and physiological adaptations (silvering) which promote sexual development and aid the transition to increased environmental salinities. The aim of this study was to identify and characterise changes in gene expression within the major osmoregulatory tissues of the eel which enable these fish to make the physiological adaptations required for transfer to SW environments. In particular, changes in the expression of the FW-adapting hormone prolactin were correlated with differential expression of known osmoregulatory important genes within the gill, intestine and kidney following the acclimation of eels to SW. Various tissues were sampled from individual fish at selected intervals over a 5-month period following FW/SW transfer and RNA was isolated. Suppressive subtractive hybridization (SSH) was used for enrichment of differentially expressed genes. Microarrays comprising 6144 cDNAs spotted in triplicate, from brain, gill, intestine and kidney libraries (1536 randomly selected clones per tissue library), were hybridized with appropriate targets and analysed. Microarray results were validated using known genes implicated in osmoregulation, such as prolactin, growth hormone, Na, K-ATPase and some unknown genes, the role of which in osmoregulation needs to be elucidated.

Advances in swine transcriptomics

The past five years have seen a tremendous rise in porcine transcriptomic data. Available porcine Expressed Sequence Tags (ESTs) have expanded greatly, with over 623,000 ESTs deposited in Genbank. ESTs have been used to expand the pig-human comparative maps, but such data has also been used in many ways to understand pig gene expression. Several methods have been used to identify genes differentially expressed (DE) in specific tissues or cell types under different treatments. These include open screening methods such as suppression subtractive hybridization, differential display, serial analysis of gene expression, and EST sequence frequency, as well as closed methods that measure expression of a defined set of sequences such as hybridization to membrane arrays and microarrays. The use of microarrays to begin large-scale transcriptome analysis has been recently reported, using either specialized or broad-coverage arrays. This review covers published results using the above techniques in the pig, as well as unpublished data provided by the research community, and reports on unpublished Affymetrix data from our group. Published and unpublished bioinformatics efforts are discussed, including recent work by our group to integrate two broad-coverage microarray platforms. We conclude by predicting experiments that will become possible with new anticipated tools and data, including the porcine genome sequence. We emphasize that the need for bioinformatics infrastructure to efficiently store and analyze the expanding amounts of gene expression data is critical, and that this deficit has emerged as a limiting factor for acceleration of genomic understanding in the pig.