Cloning and characterization of pig muscle cDNAs by an expressed sequence tag approach

Full-length enriched cDNA library construction from tissues related to energy metabolism in pigs

Genome sequencing of the pig is being accelerated because of its importance as an evolutionary and biomedical model animal as well as a major livestock animal. However, information on expressed porcine genes is insufficient to allow annotation and use of the genomic information. A series of expressed sequence tags of 5' ends of five full-length enriched cDNA libraries (SUSFLECKs) were functionally characterized. SUSFLECKs were constructed from porcine abdominal fat, induced fat cells, loin muscle, liver, and pituitary gland, and were composed of non-normalized and normalized libraries. A total of 55,658 ESTs that were sequenced once from the 5' ends of clones were produced and assembled into 17,684 unique sequences with 7,736 contigs and 9,948 singletons. In Gene Ontology analysis, two significant biological process leaf nodes were found: gluconeogenesis and translation elongation. In functional domain analysis based on the Pfam database, the beta transducin repeat domain of WD40 protein was the most frequently occurring domain. Twelve genes, including SLC25A6, EEF1G, EEF1A1, COX1, ACTA1, SLA, and ANXA2, were significantly more abundant in fat tissues than in loin muscle, liver, and pituitary gland in the SUSFLECKs. These characteristics of SUSFLECKs determined by EST analysis can provide important insight to discover the functional pathways in gene networks and to expand our understanding of energy metabolism in the pig.

A pig multi-tissue normalised cDNA library: large-scale sequencing, cluster analysis and 9K micro-array resource generation

Background

Domestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.

The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species.

Results

We constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs.

Conclusion

In the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

Generation and analysis of large-scale expressed sequence tags (ESTs) from a full-length enriched cDNA library of porcine backfat tissue

Background

Genome research in farm animals will expand our basic knowledge of the genetic control of complex traits, and the results will be applied in the livestock industry to improve meat quality and productivity, as well as to reduce the incidence of disease. A combination of quantitative trait locus mapping and microarray analysis is a useful approach to reduce the overall effort needed to identify genes associated with quantitative traits of interest.

Results

We constructed a full-length enriched cDNA library from porcine backfat tissue. The estimated average size of the cDNA inserts was 1.7 kb, and the cDNA fullness ratio was 70%. In total, we deposited 16,110 high-quality sequences in the dbEST division of GenBank (accession numbers: DT319652-DT335761). For all the expressed sequence tags (ESTs), approximately 10.9 Mb of porcine sequence were generated with an average length of 674 bp per EST (range: 200–952 bp). Clustering and assembly of these ESTs resulted in a total of 5,008 unique sequences with 1,776 contigs (35.46%) and 3,232 singleton (65.54%) ESTs. From a total of 5,008 unique sequences, 3,154 (62.98%) were similar to other sequences, and 1,854 (37.02%) were identified as having no hit or low identity (<95%) and 60% coverage in The Institute for Genomic Research (TIGR) gene index of Sus scrofa. Gene ontology (GO) annotation of unique sequences showed that approximately 31.7, 32.3, and 30.8% were assigned molecular function, biological process, and cellular component GO terms, respectively. A total of 1,854 putative novel transcripts resulted after comparison and filtering with the TIGR SsGI; these included a large percentage of singletons (80.64%) and a small proportion of contigs (13.36%).

Conclusion

The sequence data generated in this study will provide valuable information for studying expression profiles using EST-based microarrays and assist in the condensation of current pig TCs into clusters representing longer stretches of cDNA sequences. The isolation of genes expressed in backfat tissue is the first step toward a better understanding of backfat tissue on a genomic basis.

Complementary DNA macroarray analyses of differential gene expression in porcine fetal and postnatal muscle

To study differential gene expression in porcine skeletal muscle, a porcine complementary DNA (cDNA) macroarray was produced that contained 327 expressed sequence tags (EST) derived from whole embryo and adult skeletal muscle, and differential display PCR products from fetal and postnatal muscle. Total RNA from four muscle samples, 75- and 105-d fetal hind limb muscles, and 1- and 7-wk postnatal semitendinosus muscle was used to make radiolabeled targets for duplicate hybridization to the macroarray membranes in an initial screen for expression. All EST that gave clear signals (n = 238) were then re-arrayed, and hybridization was conducted with additional biological replication of samples in the 75-d and 1-wk ages. Signal intensity for each gene was normalized to signal intensity measured at control spots on each membrane, which consisted of total cDNA from liver, lung, spleen, and skeletal muscle. Both normalized ratio levels and a mixed linear model analyses were used to identify genes differentially expressed among the muscle samples. Results showed 28 genes had differences in expression level greater than twofold between the 75-d fetal and 1-wk muscle RNA samples. All 28 genes were also identified as genes with significantly different (P < 0.01) expression using a mixed linear model analysis. Nineteen of these 28 genes had significant matches (basic local alignment search tool [BLAST] score > 100; P < 0.01) to known genes, two matched genes encoding human hypothetical proteins, and seven had no significant matches to Genbank nonredundant and dbEST (database of expressed sequence tags) entries. These results were confirmed for representative genes with RNA blot analysis of seven developmental time points, including RNA from the same muscle samples tested previously in the macroarray. The RNA blot results confirmed the macroarray results for all selected genes, demonstrating that the macroarray technique used in this study is accurate and reproducible. An unknown muscle clone (M218) with a slightly less than twofold increase in expression from the 75-d to the 1-wk age (1 wk/75 d = 1.94; P = 0.0114) was also shown to differ between these two ages using RNA blot analysis, demonstrating the methods used to identify differentially expressed genes may be conservative. The association between expression patterns of vimentin and desmin was also investigated. Results indicate the switch in intermediate filament protein from vimentin to desmin occurs primarily at the level of transcription and/or RNA processing.

Generation of expressed sequence tags from a normalized porcine skeletal muscle cDNA library

Recent developments in microarray technologies permit scientists to analyze expression of thousands of genes simultaneously in diverse biological systems. In an effort to provide integrated resources for application of microarray technologies to studies of skeletal muscle growth and development in swine, we have constructed a normalized cDNA library from porcine skeletal muscle. The effectiveness of normalization was evaluated by DNA sequencing of clones randomly picked from the library before and after normalization, and also by Southern blot hybridization using probes representing abundant transcripts. Our data suggests that the normalization procedure successfully reduced the highly abundant cDNA species in the normalized library. To date, a total of 782 EST (expressed sequence tag) sequences have been generated from this normalized library (687 ESTs) and the original library (95 ESTs). The sequence information of these ESTs plus their BLAST results has been made available through a web accessible database (http://nbfgc.msu.edu). Cluster analysis of the data indicates that a total of 742 unique sequences are present in this collection. BLASTN search of the 742 EST sequences against the public database (dbEST) revealed that 139 had no significant matches (E-value > 10(-15)) to porcine ESTs already entered in the database, suggesting the possibility of their specific expression in porcine skeletal muscle. Generation of non-redundant ESTs from this library will allow us to construct cDNA microarrays for identification of gene expression changes that regulate muscle growth and affect meat quality in swine.

Isolation of expressed sequence tags of skeletal muscle of neonatal healthy and splay leg piglets and mapping by somatic cell hybrid analysis

We have isolated 14 differentially displayed and 10 further expressed sequence tags (ESTs) from Musculus biceps femoris of newborn healthy and splay leg piglets. By comparison with EMBL/GenBank data we could identify nine porcine homologues to human genes (TATA box binding protein associated factor B TAF1B; B-cell CLL/lymphoma 7B BCL7B; pyruvate dehydrogenase kinase, isoenzyme 4 PDK4; ribosomal protein S10 RPS10; SPARC-like 1 SPARCL1; epithelial protein lost in neoplasm beta EPLIN; N-myc downstream-regulated gene 2 NDRG2; pleiomorphic adenoma gene like 2 PLAGL and, BCL-2 associated transcription factor short form BTFS). Eight fragments correspond to uncharacterized ESTs and 7 ESTs had no significant match with database sequences. These data provide the first expression profiles in skeletal muscle of neonatal piglets and are a basis for candidate gene investigations for congenital splay leg in piglets. Eleven ESTs were physically mapped to porcine chromosomes 1, 3, 4, 5, 7, 8, 9 and 10 and contribute to the comparative map of humans and pigs.

Analysis of expressed sequence tags of porcine skeletal muscle

Porcine skeletal muscle genes play a major role in determining muscle growth and meat quality. Therefore, to progress towards a better understanding of the genetic factors influencing these traits, the first step is to characterize the genes expressed in skeletal muscle tissue in pig. To this aim, we constructed a porcine biceps femoris muscle cDNA library and sequenced 111 randomly isolated clones. By FASTA analysis, we identified 72 unique clones: 47 showed homology to previously identified genes in human or other mammals, 20 matched uncharacterized expressed sequence tags (ESTs), two showed no significant matches to sequences already present in DNA databases, and three other clones containing only repetitive elements were excluded from further analysis. Mitochondrial genes (16.2%), myosin heavy chain genes (9%) and the actin alpha skeletal muscle gene (9%) were the most abundant transcripts. Among the 47 identified genes several muscle-specific or predominant sequences expressed in skeletal muscle were found. The sequences of the clones matching uncharacterized human, mouse or porcine ESTs were tested by GRAIL in order to identify putative coding regions. The results of our analysis allowed the establishment of a first list of genes expressed in porcine skeletal muscle.