Gene expression patterns during intramuscular fat development in cattle.

Deposition of intramuscular fat, or "marbling," in beef cattle contributes significantly to meat quality variables, including juiciness, flavor, and tenderness. The accumulation of intramuscular fat is largely influenced by the genetic background of cattle, as well as their age and nutrition. To identify genes that can be used as early biomarkers for the prediction of marbling capacity, we studied the muscle transcriptome of 2 cattle crossbreeds with contrasting intramuscular fat content. The transcriptomes of marbling LM tissue of heifers from Wagyu x Hereford (WxH; n = 6) and Piedmontese x Hereford (PxH; n = 7) crosses were profiled by using a combination of complementary DNA microarray and quantitative reverse transcription-PCR. Five biopsies of LM were taken from each animal at approximately 3, 7, 12, 20, and 25 mo from birth. Tissue was also collected from the LM of each animal at slaughter (approximately 30 mo). Microarray experiments, conducted on the first 3 biopsies of 2 animals from each crossbreed, identified 97 differentially expressed genes. The gene expression results indicated that the LM transcriptome of animals with high marbling potential (WxH) could be reliably distinguished from less marbled animals (PxH) when the animals were as young as 7 mo of age. At this early age, one cannot reliably determine meaningful differences in intramuscular fat deposition. We observed greater expression of a set of adipogenesis- and lipogenesis-related genes in the LM of young WxH animals compared with their PxH contemporaries. In contrast, genes highly expressed in PxH animals were associated with mitochondrial oxidative activity. Further quantitative reverse transcription-PCR experiments revealed that the messenger RNA of 6 of the lipogenesis-related genes also peaked at the age of 20 to 25 mo in WxH animals. The messenger RNA expression of ADIPOQ, SCD, and THRSP was highly correlated with intramuscular fat content of an individual in WxH animals. Our study provides clear evidence of early molecular changes associated with marbling and also identifies specific time frames when intramuscular fat development in cattle muscle can be detected by using gene expression. This information could be used by animal scientists to design optimal nutrition for high marbling potential. In addition, the genes found to be highly expressed during development of marbling could be used to develop genetic markers or biomarkers to assist with beef production strategies.

A rapid method for computationally inferring transcriptome coverage and microarray sensitivity.

MOTIVATION: There are many different gene expression technologies, including cDNA and oligo-based microarrays, SAGE and MPSS. For each organism of interest, coverage of the transcriptome and the genome will be different. We address the question of what level of coverage is required to exploit the sensitivity of the different technologies, and what is the sensitivity of the different approaches in the experimental study. RESULTS: We estimate the transcriptome coverage by randomly sampling transcripts from a pre-defined tag-to-gene mapping function. For a given microarray experiment, we locate the thresholds in intensities that define the distribution of transcript abundance. These values are compared against the distribution obtained by applying the same thresholds to the intensities from differentially expressed genes. The ratio of these two distributions meets at the equilibrium defining sensitivity. We conclude that a collection of approximately 340,000 sequences is adequate for microarrays, but not large enough for maximum utilization of tag-based technologies. In the absence of large-scale sequencing, the majority of the tags detected by the latter approaches will remain unidentified until the genome sequence is available.

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 animal's 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.

Mapping of 12 bovine ribosomal protein genes using a bovine radiation hybrid panel.

Twelve bovine ribosomal protein genes, for which sequence data had been acquired from complementary deoxyribonucleic acid (cDNA) clones isolated from a cattle skin cDNA library, were mapped. As ribosomal protein genes are a group of highly conserved house keeping genes, specific primers were designed to span the intron-exon splice sites and to amplify intronic sequences, in order to obtain bovine-specific polymerase chain reaction (PCR) products. Two of 12 ribosomal protein genes were genotyped in this way and the remaining 10 were mapped using additional primers designed from within the intron. Eleven previously unmapped ribosomal protein genes were localized and one previously reported ribosomal protein gene localization was confirmed. The 12 ribosomal protein genes mapped in this study are spread over 10 chromosomes, including the X chromosome. The locations show conservation of comparative map position in cattle and human.

Genome organization and transcription strategy in the complex GNS-L intergenic region of bovine ephemeral fever rhabdovirus.

A 1622 nucleotide region of the bovine ephemeral fever virus (BEFV) genome, located between the second glycoprotein (GNS) gene and the polymerase (L) gene, has been cloned and sequenced in Australian (BB7721) and Chinese (Beijing-1) isolates of the virus. In the Australian isolate, the region contains five long open reading frames (ORFs) organized into three coding regions (alpha, beta and gamma), each of which are bound by a consensus transcription initiation and transcription termination-polyadenylation-like sequences. The alpha coding region contains three long ORFs (alpha 1, alpha 2 and alpha 3). The alpha 1 ORF encodes a 10.6 kDa polypeptide which contains hydrophobic and highly basic regions characteristic of a viroporin. The alpha 2 ORF encodes a 13.7 kDa polypeptide and overlaps the alpha 3 ORF which encodes a 5.7 kDa polypeptide. The beta coding region contains a single long ORF encoding a polypeptide of 12.2 kDa. The gamma coding region, which does not occur in Adelaide River virus (ARV), contains a single long ORF encoding a polypeptide of 13.4 kDa. The Chinese isolate shares 91% nucleotide sequence identity with the Australian isolate. The organization of the alpha, beta and gamma coding regions is preserved and the sequences of the encoded polypeptides are similar to those of BB7721. The major transcription products of the region were identified in BB7721 as polycistronic alpha (alpha 1-alpha 2-alpha 3) and beta-gamma mRNAs. Sequence similarities in the BEFV alpha-beta and beta-gamma gene junctions, and the gamma-L and beta-L gene junctions of BEFV and ARV, suggest that the gamma gene may have evolved from the beta-gene by sequence duplication.

Vaccinia virus-expressed bovine ephemeral fever virus G but not G(NS) glycoprotein induces neutralizing antibodies and protects against experimental infection.

Two related glycoproteins (G and G(NS)) encoded in the bovine ephemeral fever virus (BEFV) genome were expressed from recombinant vaccinia viruses (rVV). Both proteins were detected in lysates of rVV-infected cells by labelling with D-[6-3H]glucosamine or by immuno-blotting. The recombinant G protein (mol. mass 79 kDa) appeared slightly smaller than the native G protein but reacted with monoclonal antibodies directed against all defined neutralizing antigenic sites (G1, G2, G3a, G3b and G4). The recombinant G(NS) protein (mol. mass 90kDa) was identical in size to the native G(NS) protein and failed to react by immuno-fluorescence with anti-G protein monoclonal or poly-clonal antibodies. Antisera raised in rabbits against rVV-G or rVV-G(NS) both reacted strongly by immuno-fluorescence and immuno-electron microscopy with BEFV-infected cells. The G protein was localized intracellularly in the endoplasmic reticulum/Golgi complex and at the cell surface associated with budding and mature virus particles. The G(NS) protein also localized intracellularly in the endoplasmic reticulum/Golgi complex; however, at the cell surface it was associated with amorphous structures and not with budding or mature virions. Rabbits vaccinated with rVV-G developed high levels of antibodies which neutralized BEFV grown in either mammalian or insect cells. Cattle vaccinated with rVV-G also produced neutralizing antibodies and were protected against experimental BEFV infection. In contrast, rVV-G(NS) vaccinated rabbits and cattle failed to produce neutralizing antibodies and, after challenge, BEFV was isolated from two-thirds of the vaccinated cattle.

Structural and antigenic analysis of the nucleoprotein of bovine ephemeral fever rhabdovirus.

The nucleotide sequence of the bovine ephemeral fever virus (BEFV) genome has been determined from the 3' terminus to the end of the nucleoprotein (N) gene. The 3' leader sequence comprises 50 nucleotides and shares a common terminal three nucleotides (3'-UGC-) and a downstream U-rich domain with vesicular stomatitis virus (VSV) and rabies virus. The N gene comprises 1328 nucleotides from the transcription initiation consensus sequence (AACAGG) to the conserved transcription termination-poly(A) sequence [CATG(A)7] and encodes a polypeptide of 431 amino acids with an estimated M(r) of 49,159 and a pI of 5.4. The deduced amino acid sequence of the BEFV N protein is similar to those of other mammalian rhabdoviruses and is more closely related in sequence to vesiculoviruses (VSV Indiana and New Jersey, Piry, Chandipura) than to lyssaviruses (rabies and Mokola). An almost full-length clone, 1301 bp in length, of the BEFV N gene and clones derived from 5'-terminal (559 bp) and 3'-terminal (742 bp) fragments were expressed in Escherichia coli as glutathione-S-transferase fusion proteins. A panel of 12 BEFV N protein-specific monoclonal antibodies was shown to react in immunoblots with fusion proteins containing the almost full-length N protein and the C-terminal fragment, but not the N-terminal fragment. Two of these antibodies also reacted with baculovirus-expressed rabies virus N protein. Polyclonal mouse ascitic fluids derived from BEFV, rabies virus and several other related viruses were also shown to cross-react in immunoblots with purified preparations of rabies virus and BEFV N proteins.

Complex genome organization in the GNS-L intergenic region of Adelaide River rhabdovirus.

A 2341-nucleotide region of the Adelaide River virus (ARV) genome, located immediately downstream of the second glycoprotein (GNS) gene, has been cloned and sequenced. The region contains four long open reading frames (ORFs), the last of which represents a 1088-nucleotide fragment at the start of the ARV L gene. Between the GNS and L genes are two coding regions, separated by a single nucleotide (C), and each bounded by recognized transcription initiation (AACAG) and termination/polyadenylation (CATG[A]7) sequences. The first coding region comprises 682 nucleotides and contains two long ORFs (alpha 1 and alpha 2) which are in the same reading frame but separated by two consecutive stop codons. The alpha 1 ORF encodes a 12,545-Da polypeptide which contains highly hydrophobic and highly basic domains. The alpha 2 ORF includes a potential initiation codon 18 nucleotides downstream of the tandem stop codons and encodes a polypeptide of 11,951 Da. In ARV-infected cells, the alpha region is transcribed primarily as a long 4.7-kb polycistronic mRNA containing the G, GNS, alpha 1, and alpha ORFs. Direct sequence analysis of the mRNA indicated that the tandem stop codons between the alpha 1 and alpha 2 ORFs are retained in the transcript. The second coding region contains a single long ORF (beta) comprising 493 nucleotides which encodes a polypeptide with a calculated pl of 6.614 and molecular weight of 17,102 Da. The putative beta protein is similar in size to a protein which has been reported as a minor component of virions. The beta gene is transcribed as a 0.65-kb monocistronic mRNA for which the putative transcription termination/polyadenylation signal overlaps the L gene by 22 nucleotides.

Comparison of a dengue-2 virus and its candidate vaccine derivative: sequence relationships with the flaviviruses and other viruses.

A comparison of the sequence of the dengue-2 16681 virus with that of the candidate vaccine strain (16681-PDK53) derived from it identified 53 of the 10,723 nucleotides which differed between the strains. Nucleotide changes occurred in genes coding for all virion and nonvirion proteins, and in the 5' and 3' untranslated regions. Twenty-seven of the nucleotide changes resulted in amino acid alterations. The greatest amino acid sequence differences in the virion proteins occurred in prM (2.20%; 2/91 amino acids) followed by the M protein (1.33%; 1/75 amino acids), the C protein (0.88%; 1/114 amino acid), and the E protein (0.61%; 3/495 amino acids). Differences in the amino acid sequence of nonvirion proteins ranged from 1.51% (6/398 amino acids) in NS4 to 0.33% (3/900 amino acids) in NS5. The encoded protein sequences of 16681-PDK53 were also compared with the published sequences of other flaviviruses to obtain a detailed classification of 17 flaviviruses using the neighbor-joining tree method. The analyses of the sequence data produced dendrograms which supported the traditional groupings based on serological evidence, and they suggested that the flaviviruses have evolved by divergent mutational change and there was no evidence of genetic recombination between members of the group. Comparisons of the sequences of the flavivirus polymerase and helicase-like proteins (NS5 and NS3, respectively) with those from other viruses yielded a classification of the flaviviruses indicating that the primary division of the flaviviruses was between those transmitted by mosquitoes and those transmitted by ticks.

NS 1 gene sequences from eight dengue-2 viruses and their evolutionary relationships with other dengue-2 viruses.

The nucleotide sequences of the NS 1 genes from five Thai and three Sri Lankan dengue-2 viruses were determined by sequencing the viral RNA using synthetic oligonucleotide primers. The results were shown to be similar to four published dengue-2 NS 1 sequences and the classification of these genes was compared with the one obtained for the envelope genes of the same viruses. The classification was similar and showed that the Thai isolates could be divided into two separate groups and that the Sri Lankan isolates were distinct. We found no correlation between disease severity, serological response (1 degree or 2 degrees), or year of isolation and various aspects of NS 1 protein sequence variation; and no particular amino acid changes were correlated with virulence. The sequences were combined with those published and classified elsewhere to provide a comprehensive E/NS 1 gene taxonomy of dengue-2 virus isolates.