Yonghong Wang

Gene expression studies of developing bovine longissimus muscle from two different beef cattle breeds

BackgroundThe muscle fiber number and fiber composition of muscle is largely determined during prenatal development. In order to discover genes that are involved in determining adult muscle phenotypes, we studied the gene expression profile of developing fetal bovine longissimus muscle from animals with two different genetic backgrounds using a bovine cDNA microarray. Fetal longissimus muscle was sampled at 4 stages of myogenesis and muscle maturation: primary myogenesis (d 60), secondary myogenesis (d 135), as well as beginning (d 195) and final stages (birth) of functional differentiation of muscle fibers. All fetuses and newborns (total n = 24) were from Hereford dams and crossed with either Wagyu (high intramuscular fat) or Piedmontese (GDF8 mutant) sires, genotypes that vary markedly in muscle and compositional characteristics later in postnatal life.ResultsWe obtained expression profiles of three individuals for each time point and genotype to allow comparisons across time and between sire breeds. Quantitative reverse transcription-PCR analysis of RNA from developing longissimus muscle was able to validate the differential expression patterns observed for a selection of differentially expressed genes, with one exception. We detected large-scale changes in temporal gene expression between the four developmental stages in genes coding for extracellular matrix and for muscle fiber structural and metabolic proteins. FSTL1 and IGFBP5 were two genes implicated in growth and differentiation that showed developmentally regulated expression levels in fetal muscle. An abundantly expressed gene with no functional annotation was found to be developmentally regulated in the same manner as muscle structural proteins. We also observed differences in gene expression profiles between the two different sire breeds. Wagyu-sired calves showed higher expression of fatty acid binding protein 5 (FABP5) RNA at birth. The developing longissimus muscle of fetuses carrying the Piedmontese mutation shows an emphasis on glycolytic muscle biochemistry and a large-scale up-regulation of the translational machinery at birth. We also document evidence for timing differences in differentiation events between the two breeds.ConclusionTaken together, these findings provide a detailed description of molecular events accompanying skeletal muscle differentiation in the bovine, as well as gene expression differences that may underpin the phenotype differences between the two breeds. In addition, this study has highlighted a non-coding RNA, which is abundantly expressed and developmentally regulated in bovine fetal muscle.

Gene expression profiling of muscle tissue in Brahman steers during nutritional restriction

Expression profiling using microarrays allows for the detailed characterization of the gene networks that regulate an animals response to environmental stresses. During nutritional restriction, processes such as protein turnover, connective tissue remodeling, and muscle atrophy take place in the skeletal muscle of the animal. These processes and their regulation are of interest in the context of managing livestock for optimal production efficiency and product quality. Here we expand on recent research applying complementary DNA (cDNA) microarray technology to the study of the effect of nutritional restriction on bovine skeletal muscle. Using a custom cDNA microarray of 9,274 probes from cattle muscle and s.c. fat libraries, we examined the differential gene expression profile of the LM from 10 Brahman steers under three different dietary treatments. The statistical approach was based on mixed-model ANOVA and model-based clustering of the BLUP solutions for the gene x diet interaction effect. From the results, we defined a transcript profile of 156 differentially expressed array elements between the weight loss and weight gain diet substrates. After sequence and annotation analyses, the 57 upregulated elements represented 29 unique genes, and the 99 downregulated elements represented 28 unique genes. Most of these co-regulated genes cluster into groups with distinct biological function related to protein turnover and cytoskeletal metabolism and contribute to our mechanistic understanding of the processes associated with remodeling of muscle tissue in response to nutritional stress.

The genome of bovine ephemeral fever rhabdovirus contains two related glycoprotein genes.

A 3789 nucleotide region of the bovine ephemeral fever virus (BEFV) genome, located 1.65 kb downstream of the N gene, has been cloned and sequenced. The region contains two long open reading frames (ORFs) which are bounded by putative consensus (AACAGG) and polyadenylation (CATG[A]7) sequences and are separated by an intergenic region of 53 nucleotides. Discrete mRNAs corresponding to each ORF have been identified. The first ORF encodes a polypeptide comprising 623 residues which was identified by peptide sequencing as the virion G protein. The deduced amino acid sequence of the G protein includes putative signal and transmembrane domains and five potential glycosylation sites. The second ORF encodes a polypeptide of 586 amino acids which also has characteristics of a rhabdovirus glycoprotein, including putative signal and transmembrane domains and eight potential glycosylation sites, and appears to correspond to a 90-kDa nonstructural glycoprotein (GNS) identified in BEFV-infected cells (Walker et al. [1991] J. Gen. Virol. 72, 67-74). A database search indicated that both the G and GNS proteins share significant amino acid sequence homology with other rhabdovirus G proteins and with each other. Highest homology scores for each protein were with sigma virus and vesicular stomatitis virus serotypes.

Simultaneous identification of differential gene expression and connectivity in inflammation, adipogenesis and cancer

MOTIVATION Biological differences between classes are reflected in transcriptional changes which in turn affect the levels by which essential genes are individually expressed and collectively connected. The purpose of this communication is to introduce an analytical procedure to simultaneously identify genes that are differentially expressed (DE) as well as differentially connected (DC) in two or more classes of interest. RESULTS Our procedure is based on a two-step approach: First, mixed-model equations are applied to obtain the normalized expression levels of each gene in each class treatment. These normalized expressions form the basis to compute a measure of (possible) DE as well as the correlation structure existing among genes. Second, a two-component mixture of bi-variate distributions is fitted to identify the component that encapsulates those genes that are DE and/or DC. We demonstrate our approach using three distinct datasets including a human systemic inflammation oligonucleotide data; a spotted cDNA data dealing with bovine in vitro adipogenesis and SAGE database on cancerous and normal tissue samples.

Proteins of bovine ephemeral fever virus.

The proteins of bovine ephemeral fever virus (BEFV) were examined in purified virions and in infected BHK-21 cells. Five structural proteins were named L (180K), G (81K), N (52K), M1 (43K) and M2 (29K). The 81K G protein incorporated [3H]glucosamine, was removed from virions by treatment with Triton X-100 and bound monoclonal antibodies which were both neutralizing and protective. Treatment of virions with Triton X-100 and 0.2 to 1.0 M-NaCl progressively released L, M1 and M2. The N protein remained associated with nucleocapsids in up to 2.5 M-NaCl. The glycoprotein (G), nucleoprotein (N) and matrix protein (M2) were phosphorylated. In BEFV-infected BHK-21 cells, five virus-induced proteins were detected from 12 h post-infection. The L, N, M1 and M2 proteins corresponded to those detected in virions whereas the G protein existed in two forms. In tunicamycin-treated cells these occurred as 67K and 71K non-glycosylated precursors. In the absence of tunicamycin, 77K and 79K glycosylated forms were further modified to produce the 81K virion G protein and a 90K cell-associated form. Five viral proteins were also detected in cells infected with the closely related Berrimah virus; the Berrimah virus G protein was also present in two forms.

Inferring the transcriptional landscape of bovine skeletal muscle by integrating co-expression networks

Background Despite modern technologies and novel computational approaches, decoding causal transcriptional regulation remains challenging. This is particularly true for less well studied organisms and when only gene expression data is available. In muscle a small number of well characterised transcription factors are proposed to regulate development. Therefore, muscle appears to be a tractable system for proposing new computational approaches. Methodology/Principal Findings Here we report a simple algorithm that asks “which transcriptional regulator has the highest average absolute co-expression correlation to the genes in a co-expression module?” It correctly infers a number of known causal regulators of fundamental biological processes, including cell cycle activity (E2F1), glycolysis (HLF), mitochondrial transcription (TFB2M), adipogenesis (PIAS1), neuronal development (TLX3), immune function (IRF1) and vasculogenesis (SOX17), within a skeletal muscle context. However, none of the canonical pro-myogenic transcription factors (MYOD1, MYOG, MYF5, MYF6 and MEF2C) were linked to muscle structural gene expression modules. Co-expression values were computed using developing bovine muscle from 60 days post conception (early foetal) to 30 months post natal (adulthood) for two breeds of cattle, in addition to a nutritional comparison with a third breed. A number of transcriptional landscapes were constructed and integrated into an always correlated landscape. One notable feature was a ‘metabolic axis’ formed from glycolysis genes at one end, nuclear-encoded mitochondrial protein genes at the other, and centrally tethered by mitochondrially-encoded mitochondrial protein genes. Conclusions/Significance The new module-to-regulator algorithm complements our recently described Regulatory Impact Factor analysis. Together with a simple examination of a co-expression modules contents, these three gene expression approaches are starting to illuminate the in vivo transcriptional regulation of skeletal muscle development.

Gene expression profiling of bovine in vitro adipogenesis using a cDNA microarray

The gene expression profile of bovine bone marrow stromal cells undergoing adipogenesis was established using a custom cDNA microarray. Cells that were treated with adipogenic stimulants and those that were not were collected at each of the six time points, and gene expression differences between the treated and untreated samples within each time point were compared using a microarray. Statistical analyses revealed that 158 genes showed a minimum fold change of 2 in at least one of the five post-differentiation time points. These genes are involved in various cellular pathways and functions, including lipogenesis, glycolysis, cytoskeleton remodelling, extracellular matrix, transcription as well as various signalling pathways such as insulin, calcium and wingless signalling. The experiment also identified 17 differentially expressed (DE) microarray elements with no assigned function. Quantitative real-time PCR was employed to validate eight DE genes, and the PCR data were found to reproduce the microarray data for these eight genes. Subsequent gene ontology annotation was able to provide a global overview of the molecular function of DE genes during adipogenesis. This analysis was able to indicate the importance of different gene categories at various stages of adipogenic conversion, thereby providing further insights into the molecular changes during bovine adipogenesis.

Chronic exposure to anabolic steroids induces the muscle expression of oxytocin and a more than fiftyfold increase in circulating oxytocin in cattle

Molecular mechanisms in skeletal muscle associated with anabolic steroid treatment of cattle are unclear and we aimed to characterize transcriptional changes. Cattle were chronically exposed (68 ± 20 days) to a steroid hormone implant containing 200 mg trenbolone acetate and 20 mg estradiol (Revalor-H). Biopsy samples from 48 cattle (half treated) from longissimus dorsi (LD) muscle under local anesthesia were collected. Gene expression levels were profiled by microarray, covering 16,944 unique bovine genes: 121 genes were differentially expressed (DE) due to the implant (99.99% posterior probability of not being false positives). Among DE genes, a decrease in expression of a number of fat metabolism-associated genes, likely reflecting the lipid storage activity of intramuscular adipocytes, was observed. The expression of IGF1 and genes related to the extracellular matrix, slow twitch fibers, and cell cycle (including SOX8, a satellite cell marker) was increased in the treated muscle. Unexpectedly, a very large 21- (microarray) to 97 (real time quantitative PCR)-fold higher expression of the mRNA encoding the neuropeptide hormone oxytocin was observed in treated muscle. We also observed an ∼50-fold higher level of circulating oxytocin in the plasma of treated animals at the time of biopsy. Using a coexpression network strategy OXTR was identified as more likely than IGF1R to be a major mediator of the muscle response to Revalor-H. A re-investigation of in vivo cattle LD muscle samples during early to mid-fetal development identified a >128-fold increased expression of OXT, coincident with myofiber differentiation and fusion. We propose that oxytocin may be involved in mediating the anabolic effects of Revalor-H treatment.

Gene expression profiling of Hereford Shorthorn cattle following challenge with Boophilus microplus tick larvae

The ability of cattle to resist tick infestations is partly genetically determined. In order to better define the nature of Bos taurus resistance to the cattle tick Boophilus microplus, skin gene expression was studied using a cattle skin derived cDNA microarray. Expression profiles were determined in skin biopsies sampled from three highly tick resistant animals (HR) and two animals with lower tick resistance (LR) at time 0, immediately before challenge, and again 24 h after challenge. The analysis of the resulting expression data addressed two biological questions: first, for any animal exposed to ticks, which genes are differentially expressed in the 24 h following challenge; and second, which genes are differentially expressed between animals of high and low resistance at 24 h after challenge? In total, 214 genes were found to be differentially expressed in response to larval challenge across all the animals. Seventy-two genes were upregulated and 76 were downregulated at 24 h after challenge. Genes with significantly altered gene expression levels following tick infestation were predominantly keratin genes or mitochondrial genes, as well as odorant binding protein (OBP) and Bos taurus major allergen BDA20. In addition, we identified 66 genes with differential expression between HR and LR animals at 24 h. Of these, genes representing the extracellular matrix and immunoglobulin gene expression pathways were overrepresented. Three differentially expressed genes, OBP, Bos taurus major allergen BDA20 and dendritic cell protein HFL-B5 were further analysed by quantitative reverse transcription PCR (qRT-PCR). The qRT-PCR assay results closely mirrored the expression profiles found in the microarray experiment.