Transcriptomic and proteomic profiling of two porcine tissues using high-throughput technologies

Comparative transcriptome and proteome analyses of the longissimus dorsi muscle for explaining the difference between donkey meat and other meats

Domestic donkeys (Equus asinus) have been maintained mainly for service purposes in the past. Nowadays, there is an increasing interest in donkey milk and meat production in several countries, including China. Donkey meat is highly consumed because of its nutritional value and unique flavor. However, genomic studies on donkey meat are limited. Therefore, in this study, we aimed to examine the molecular difference of longissimus dorsi muscles of donkey, cow, and goat. RNA sequencing and Proteome sequencing technology were used to analyze the transcriptome and proteome of the longissimus dorsi muscle of donkey, cow, and goat. A total of 1338 and 1780 differentially expressed genes (DEGs) were identified in donkey meat compared with that in cow and goat meat, respectively. Most of the DEGs were involved in biological processes, including small GTPase-mediated signal transduction, protein ubiquitination, protein glycosylation, and MAP kinase tyrosine/serine/threonine phosphatase activity. Additionally, 764 and 1024 differentially expressed proteins (DEPs) were identified in cow vs. donkey, and goat vs. donkey, respectively; these DEPs were mainly involved in metabolism. Genetic variation and regulatory factors can combine as a database to provide more valuable molecular information for further analysis.

Salivary proteomics of healthy dogs: An in depth catalog

Objective

To provide an in-depth catalog of the salivary proteome and endogenous peptidome of healthy dogs, evaluate proteins and peptides with antimicrobial properties, and compare the most common salivary proteins and peptides between different breed phylogeny groups.

Methods

36 healthy dogs without evidence of periodontal disease representing four breed phylogeny groups, based upon single nucleotide polymorphism haplotypes (ancient, herding/sighthound, and two miscellaneous groups). Saliva collected from dogs was pooled by phylogeny group and analyzed using nanoscale liquid chromatography-tandem mass spectrometry. Resulting tandem mass spectra were compared to databases for identification of endogenous peptides and inferred proteins.

Results

2,491 proteins and endogenous peptides were found in the saliva of healthy dogs with no periodontal disease. All dog phylogeny groups’ saliva was rich in proteins and peptides with antimicrobial functions. The ancient breeds group was distinct in that it contained unique proteins and was missing many proteins and peptides present in the other groups.

Conclusions and clinical relevance

Using a sophisticated nanoscale liquid chromatography-tandem mass spectrometry, we were able to identify 10-fold more salivary proteins than previously reported in dogs. Seven of the top 10 most abundant proteins or peptides serve immune functions and many more with various antimicrobial mechanisms were found. This is the most comprehensive analysis of healthy canine saliva to date, and will provide the groundwork for future studies analyzing salivary proteins and endogenous peptides in disease states.

Tissue-Based Mapping of the Fathead Minnow (Pimephales promelas) Transcriptome and Proteome

Omics approaches are broadly used to explore endocrine and toxicity-related pathways and functions. Nevertheless, there is still a significant gap in knowledge in terms of understanding the endocrine system and its numerous connections and intricate feedback loops, especially in non-model organisms. The fathead minnow (Pimephales promelas) is a widely used small fish model for aquatic toxicology and regulatory testing, particularly in North America. A draft genome has been published, but the amount of available genomic or transcriptomic information is still far behind that of other more broadly studied species, such as the zebrafish. Here, we used a proteogenomics approach to survey the tissue-specific proteome and transcriptome profiles in adult male fathead minnow. To do so, we generated a draft transcriptome using short and long sequencing reads from liver, testis, brain, heart, gill, head kidney, trunk kidney, and gastrointestinal tract. We identified 30,378 different putative transcripts overall, with the assembled contigs ranging in size from 264 to over 9,720 nts. Over 17,000 transcripts were >1,000 nts, suggesting a robust transcriptome that can be used to interpret RNA sequencing data in the future. We also performed RNA sequencing and proteomics analysis on four tissues, including the telencephalon, hypothalamus, liver, and gastrointestinal tract of male fish. Transcripts ranged from 0 to 600,000 copies per gene and a large portion were expressed in a tissue-specific manner. Specifically, the telencephalon and hypothalamus shared the most expressed genes, while the gastrointestinal tract and the liver were quite distinct. Using protein profiling techniques, we identified a total of 4,045 proteins in the four tissues investigated, and their tissue-specific expression pattern correlated with the transcripts at the pathway level. Similarly to the findings with the transcriptomic data, the hypothalamus and telencephalon had the highest degree of similarity in the proteins detected. The main purpose of this analysis was to generate tissue-specific omics data in order to support future aquatic ecotoxicogenomic and endocrine-related studies as well as to improve our understanding of the fathead minnow as an ecological model.

The Pig PeptideAtlas: A resource for systems biology in animal production and biomedicine

Biological research of Sus scrofa, the domestic pig, is of immediate relevance for food production sciences, and for developing pig as a model organism for human biomedical research. Publicly available data repositories play a fundamental role for all biological sciences, and protein data repositories are in particular essential for the successful development of new proteomic methods. Cumulative proteome data repositories, including the PeptideAtlas, provide the means for targeted proteomics, system-wide observations, and cross-species observational studies, but pigs have so far been underrepresented in existing repositories. We here present a significantly improved build of the Pig PeptideAtlas, which includes pig proteome data from 25 tissues and three body fluid types mapped to 7139 canonical proteins. The content of the Pig PeptideAtlas reflects actively ongoing research within the veterinary proteomics domain, and this article demonstrates how the expression of isoform-unique peptides can be observed across distinct tissues and body fluids. The Pig PeptideAtlas is a unique resource for use in animal proteome research, particularly biomarker discovery and for preliminary design of SRM assays, which are equally important for progress in research that supports farm animal production and veterinary health, as for developing pig models with relevance to human health research.

Integration of genomic and proteomic data to identify candidate genes in HT-29 cells after incubation with Bifidobacterium bifidum ATCC 29521

As the predominant group inhabiting the human gastrointestinal tract, bifidobacteria play a vital role in human nutrition, therapeutics, and health by shaping and maintaining the gut ecosystem, reducing blood cholesterol, and promoting the supply of nutrients. The interaction between bacterial cells and human intestinal epithelial cell lines has been studied for decades in an attempt to understand the mechanisms of action. These studies, however, have been limited by lack of genomic and proteomic database to aid in achieving comprehensive understanding of these mechanisms at molecular levels. Microarray data (GSE: 74119) coupled with isobaric tags for relative and absolute quantitation (iTRAQ) were performed to detect differentially expressed genes and proteins in HT-29 cells after incubation with Bifidobacterium bifidum. Real-time quantitative PCR, gene ontology, and Kyoto Encyclopedia of Genes and Genomes analyses were further conducted for mRNA validation, functional annotation, and pathway identification, respectively. According to the results of microarray, 1,717 differentially expressed genes, including 1,693 upregulated and 24 downregulated genes, were selected and classified by the gene ontology database. The iTRAQ analysis identified 43 differentially expressed proteins, where 29 proteins were upregulated and 14 proteins were downregulated. Eighty-two candidate genes showing consistent differences with microarray and iTRAQ were further validated in HT-29 and Caco-2 cells by real-time quantitative PCR. Nine of the top genes showing interesting results with high confidence were further investigated in vivo in mice intestine samples. Integration of genomic and proteomic data provides an approach to identify candidate genes that are more likely to function in ubiquitin-mediated proteolysis, positive regulation of apoptosis, membrane proteins, and transferase catalysis. These findings might contribute to our understanding of molecular mechanisms regulating the interaction between probiotics and intestinal epithelial cell lines.

Integrative analysis of transcriptomics and proteomics of skeletal muscles of the Chinese indigenous Shaziling pig compared with the Yorkshire breed

Background

The Shaziling pig (Sus scrofa) is a well-known indigenous breed in China. One of its main advantages over European breeds is its high meat quality. However, little genetic information is available for the Shaziling pig. To screen for differentially expressed genes and proteins that might be responsible for the meat quality, the longissimus dorsi muscles from Shaziling and Yorkshire pig breeds were investigated using an integrative analysis of transcriptomics and proteomics, involving high-throughput sequencing, the two-dimensional gel electrophoresis, and mass spectrometry.

Results

Sequencing produced 79,320 unigenes by de novo assembly, and 488 differentially expressed genes in the longissimus dorsi muscle of Shaziling pig compared with the Yorkshire breed were identified. Gene Ontology term enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes analysis showed that the gene products were mainly involved in metabolism, protein binding, and regulation of skeletal muscle development. At the protein level, 23 differentially expressed proteins were identified, which were potentially associated with fatty acid metabolism, the glycolytic pathway, and skeletal muscle growth. Eight differentially expressed genes were confirmed by real-time PCR. These results give an insight into the mechanisms underlying the formation of skeletal muscle in the Shaziling pig.

Conclusions

Certain differentially expressed genes and proteins are involved in fatty acid metabolism, intramuscular fat deposition, and skeletal muscle growth in the Shaziling pig. These results provide candidate genes for improving meat quality and will promote further transcriptomic research in Shaziling pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0389-y) contains supplementary material, which is available to authorized users.

A quantitative proteomics-based signature of platinum sensitivity in ovarian cancer cell lines

Although DNA encodes the molecular instructions that underlie the control of cell function, it is the proteins that are primarily responsible for implementing those instructions. Therefore quantitative analyses of the proteome would be expected to yield insights into important candidates for the detection and treatment of disease. We present an iTRAQ (isobaric tag for relative and absolute quantification)-based proteomic analysis of ten ovarian cancer cell lines and two normal ovarian surface epithelial cell lines. We profiled the abundance of 2659 cellular proteins of which 1273 were common to all 12 cell lines. Of the 1273, 75 proteins exhibited elevated expression and 164 proteins had diminished expression in the cancerous cells compared with the normal cell lines. The iTRAQ expression profiles allowed us to segregate cell lines based upon sensitivity and resistance to carboplatin. Importantly, we observed no substantial correlation between protein abundance and RNA expression or epigenetic DNA methylation data. Furthermore, we could not discriminate between sensitivity and resistance to carboplatin on the basis of RNA expression and DNA methylation data alone. The present study illustrates the importance of proteomics-based discovery for defining the basis for the carboplatin response in ovarian cancer and highlights candidate proteins, particularly involved in cellular redox regulation, homologous recombination and DNA damage repair, which otherwise could not have been predicted from whole genome and expression data sources alone.

Chernobyl seed project. Advances in the identification of differentially abundant proteins in a radio-contaminated environment

Plants have the ability to grow and successfully reproduce in radio-contaminated environments, which has been highlighted by nuclear accidents at Chernobyl (1986) and Fukushima (2011). The main aim of this article is to summarize the advances of the Chernobyl seed project which has the purpose to provide proteomic characterization of plants grown in the Chernobyl area. We present a summary of comparative proteomic studies on soybean and flax seeds harvested from radio-contaminated Chernobyl areas during two successive generations. Using experimental design developed for radio-contaminated areas, altered abundances of glycine betaine, seed storage proteins, and proteins associated with carbon assimilation into fatty acids were detected. Similar studies in Fukushima radio-contaminated areas might complement these data. The results from these Chernobyl experiments can be viewed in a user-friendly format at a dedicated web-based database freely available at http://www.chernobylproteomics.sav.sk.

Metagenomics using next-generation sequencing

Traditionally, microbial genome sequencing has been restricted to the small number of species that can be grown in pure culture. The progressive development of culture-independent methods over the last 15 years now allows researchers to sequence microbial communities directly from environmental samples. This approach is commonly referred to as "metagenomics" or "community genomics". However, the term metagenomics is applied liberally in the literature to describe any culture-independent analysis of microbial communities. Here, we define metagenomics as shotgun ("random") sequencing of the genomic DNA of a sample taken directly from the environment. The metagenome can be thought of as a sampling of the collective genome of the microbial community. We outline the considerations and analyses that should be undertaken to ensure the success of a metagenomic sequencing project, including the choice of sequencing platform and methods for assembly, binning, annotation, and comparative analysis.

LC-MS/MS analysis of visceral and subcutaneous adipose tissue proteomes in young goats with focus on innate immunity and inflammation related proteins

The endocrine role of adipose tissue and its involvement in several physiological and pathological processes are well recognized. Studies on human, mouse and rat adipose tissues have made clear that subcutaneous and visceral deposits play different roles, which is also reflected by different protein and gene expression patterns. In ruminants, fat tissues play important biological roles not only for animal health, but also for quality and gain in meat and milk production. Yet very few studies have explored the ruminant adipose tissue proteomes. The aim of our study was to compare subcutaneous and visceral adipose tissues of goat, focusing on proteins involved in immune and inflammatory response. A 2-D LC-MS/MS approach followed by cluster analysis shows a clear distinction between subcutaneous and visceral fat tissue proteomes, and qualitative RT-PCR based analysis of 30 potential adipokines further confirmed the individual expression patterns of 26 of these, including 7 whose mRNA expression was observed for the first time in adipose tissues. This study provides a first description of adipose tissue proteomes in goat, and presents observations on novel proteins related to metabolic and inflammatory pathways. The mass spectrometry data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD000564.

BIOLOGICAL SIGNIFICANCE

The proteomic analysis of different subcutaneous and visceral adipose tissue deposits showed tissue specific differences in protein expressions of well known as well as novel adipokines. This highlights the importance of sampling site when studying adipose tissue's metabolic roles. The protein expression characteristics of adipose tissues was evaluated by quantitative RT-PCR, and confirmed that adipose tissues play a central role in controlling inflammation, detoxification and coagulation pathways, as well as regulation of body fat mobilization in dairy animals. These findings are of particular interest in farm animals where health and production traits are important for animal welfare and for economic gains.

Expression of proteins associated with adipocyte lipolysis was significantly changed in the adipose tissues of the obese spontaneously hypertensive/NDmcr-cp rat

BACKGROUND

The etiology of the metabolic syndrome is complex, and is determined by the interplay of both genetic and environmental factors. The present study was designed to identify genes and proteins in the adipose tissues with altered expression in the spontaneously hypertensive/NIH -corpulent rat, SHR/NDmcr-cp (CP) and to find possible molecular targets associated with the pathogenesis or progression of obesity related to the metabolic syndrome.

METHODS

We extracted RNAs and proteins from the epididymal adipose tissues in CP, SHR/Lean (Lean), and Wistar Kyoto (WKY) rats and performed microarray analysis and two-dimensional difference in gel electrophoresis (2D-DIGE) linked to a matrix-assisted laser desorption ionization time-of-flight tandem mass spectrometry (MALDI-TOF/TOF MS).

RESULTS

The results showed different mRNA and protein expression levels in the adipose tissue: oligo DNA microarray identified 33 genes that were significantly (P < 0.01) up-regulated and 17 genes significantly down-regulated in CP compared with WKY and Lean rats at both 6 and 25 weeks of age. The affected genes-proteins were associated with lipolytic enzymes stimulated by peroxisome proliferator-activated receptor (PPAR) signaling. Further analysis using the 2D-DIGE connected with MALDI-TOF/TOF analysis, the expression of monoglyceride lipase (MGLL) was significantly up-regulated and that of carboxylesterase 3 (CES3) was significantly down-regulated in 6- and 25-week-old CP compared with age-matched control (WKY and Lean rats).

CONCLUSIONS

Our results suggest the possible involvement of proteins associated with adipocyte lipolysis in obesity related to the metabolic syndrome.

Characterization of Liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis

Background

Liaoning cashmere goat is a famous goat breed for cashmere wool. In order to increase the transcriptome data and accelerate genetic improvement for this breed, we performed denovo transcriptome sequencing to generate the first expressed sequence tag dataset for the Liaoning cashmere goat, using next-generation sequencing technology.

Results

Transcriptome sequencing of Liaoning cashmere goat on a Roche 454 platform yielded 804,601 high-quality reads. Clustering and assembly of these reads produced a non-redundant set of 117,854 unigenes, comprising 13,194 isotigs and 104,660 singletons. Based on similarity searches with known proteins, 17,356 unigenes were assigned to 6,700 GO categories, and the terms were summarized into three main GO categories and 59 sub-categories. 3,548 and 46,778 unigenes had significant similarity to existing sequences in the KEGG and COG databases, respectively. Comparative analysis revealed that 42,254 unigenes were aligned to 17,532 different sequences in NCBI non-redundant nucleotide databases. 97,236 (82.51%) unigenes were mapped to the 30 goat chromosomes. 35,551 (30.17%) unigenes were matched to 11,438 reported goat protein-coding genes. The remaining non-matched unigenes were further compared with cattle and human reference genes, 67 putative new goat genes were discovered. Additionally, 2,781 potential simple sequence repeats were initially identified from all unigenes.

Conclusion

The transcriptome of Liaoning cashmere goat was deep sequenced, denovo assembled, and annotated, providing abundant data to better understand the Liaoning cashmere goat transcriptome. The potential simple sequence repeats provide a material basis for future genetic linkage and quantitative trait loci analyses.

A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus)

Background

Genomic resources for plant and animal species that are under exploitation primarily for human consumption are increasingly important, among other things, for understanding physiological processes and for establishing adequate genetic selection programs. Current available techniques for high-throughput sequencing have been implemented in a number of species, including fish, to obtain a proper description of the transcriptome. The objective of this study was to generate a comprehensive transcriptomic database in turbot, a highly priced farmed fish species in Europe, with potential expansion to other areas of the world, for which there are unsolved production bottlenecks, to understand better reproductive- and immune-related functions. This information is essential to implement marker assisted selection programs useful for the turbot industry.

Results

Expressed sequence tags were generated by Sanger sequencing of cDNA libraries from different immune-related tissues after several parasitic challenges. The resulting database (“Turbot 2 database”) was enlarged with sequences generated from a 454 sequencing run of brain-hypophysis-gonadal axis-derived RNA obtained from turbot at different development stages. The assembly of Sanger and 454 sequences generated 52,427 consensus sequences (“Turbot 3 database”), of which 23,661 were successfully annotated. A total of 1,410 sequences were confirmed to be related to reproduction and key genes involved in sex differentiation and maturation were identified for the first time in turbot (AR, AMH, SRY-related genes, CYP19A, ZPGs, STAR FSHR, etc.). Similarly, 2,241 sequences were related to the immune system and several novel key immune genes were identified (BCL, TRAF, NCK, CD28 and TOLLIP, among others). The number of genes of many relevant reproduction- and immune-related pathways present in the database was 50–90% of the total gene count of each pathway. In addition, 1,237 microsatellites and 7,362 single nucleotide polymorphisms (SNPs) were also compiled. Further, 2,976 putative natural antisense transcripts (NATs) including microRNAs were also identified.

Conclusions

The combined sequencing strategies employed here significantly increased the turbot genomic resources available, including 34,400 novel sequences. The generated database contains a larger number of genes relevant for reproduction- and immune-associated studies, with an excellent coverage of most genes present in many relevant physiological pathways. This database also allowed the identification of many microsatellites and SNP markers that will be very useful for population and genome screening and a valuable aid in marker assisted selection programs.

Pressure-overload hypertrophy of the developing heart reveals activation of divergent gene and protein pathways in the left and right ventricular myocardium

Right ventricular (RV) and left ventricular (LV) myocardium differ in their pathophysiological response to pressure-overload hypertrophy. In this report we use microarray and proteomic analyses to identify pathways modulated by LV-aortic banding (AOB) and RV-pulmonary artery banding (PAB) in the immature heart. Newborn New Zealand White rabbits underwent banding of the descending thoracic aorta [LV-AOB; n = 6]. RV-PAB was achieved by banding the pulmonary artery (n = 6). Controls (n = 6 each) were sham-manipulated. After 4 (LV-AOB) and 6 (RV-PAB) wk recovery, the hearts were removed and matched RNA and proteins samples were isolated for microarray and proteomic analysis. Microarray and proteomic data demonstrate that in LV-AOB there is increased transcript expression levels for oxidative phosphorylation, mitochondria energy pathways, actin, ILK, hypoxia, calcium, and protein kinase-A signaling and increased protein expression levels of proteins for cellular macromolecular complex assembly and oxidative phosphorylation. In RV-PAB there is also an increased transcript expression levels for cardiac oxidative phosphorylation but increased protein expression levels for structural constituents of muscle, cardiac muscle tissue development, and calcium handling. These results identify divergent transcript and protein expression profiles in LV-AOB and RV-PAB and provide new insight into the biological basis of ventricular specific hypertrophy. The identification of these pathways should allow for the development of specific therapeutic interventions for targeted treatment and amelioration of LV-AOB and RV-PAB to ameliorate morbidity and mortality.

Altered gene and protein expression in liver of the obese spontaneously hypertensive/NDmcr-cp rat

BACKGROUND

It is difficult to study the mechanisms of the metabolic syndrome in humans due to the heterogeneous genetic background and lifestyle. The present study investigated changes in the gene and protein profiles in an animal model of the metabolic syndrome to identify the molecular targets associated with the pathogenesis and progression of obesity related to the metabolic syndrome.

METHODS

We extracted mRNAs and proteins from the liver tissues of 6- and 25-week-old spontaneously hypertensive/NIH -corpulent rat SHR/NDmcr-cp (CP), SHR/Lean (Lean) and Wistar Kyoto rats (WKY) and performed microarray analysis and two-dimensional difference in gel electrophoresis (2D-DIGE) linked to a matrix-assisted laser desorption ionization time-of-flight tandem mass spectrometry (MALDI-TOF/TOF MS).

RESULTS

The microarray analysis identified 25 significantly up-regulated genes (P < 0.01; log10 > 1) and 31 significantly down-regulated genes (P < 0.01; log10 < -1) in 6- and 25-week-old CP compared with WKY and Lean. Several of these genes are known to be involved in important biological processes such as electron transporter activity, electron transport, lipid metabolism, ion transport, transferase, and ion channel activity. MALDI-TOF/TOF MS identified 31 proteins with ±1.2 fold change (P < 0.05) in 6- and 25-week-old CP, compared with age-matched WKY and Lean. The up-regulated proteins are involved in metabolic processes, biological regulation, catalytic activity, and binding, while the down-regulated proteins are involved in endoplasmic reticulum stress-related unfolded protein response.

CONCLUSION

Genes with significant changes in their expression in transcriptomic analysis matched very few of the proteins identified in proteomics analysis. However, annotated functional classifications might provide an important reference resource to understand the pathogenesis of obesity associated with the metabolic syndrome.

Nuclear Factor kappa B is central to Marek's disease herpesvirus induced neoplastic transformation of CD30 expressing lymphocytes in-vivo

Background

Marek’s Disease (MD) is a hyperproliferative, lymphomatous, neoplastic disease of chickens caused by the oncogenic Gallid herpesvirus type 2 (GaHV-2; MDV). Like several human lymphomas the neoplastic MD lymphoma cells overexpress the CD30 antigen (CD30^hi) and are in minority, while the non-neoplastic cells (CD30^lo) form the majority of population. MD is a unique natural in-vivo model of human CD30^hi lymphomas with both natural CD30^hi lymphomagenesis and spontaneous regression. The exact mechanism of neoplastic transformation from CD30^lo expressing phenotype to CD30^hi expressing neoplastic phenotype is unknown. Here, using microarray, proteomics and Systems Biology modeling; we compare the global gene expression of CD30^lo and CD30^hi cells to identify key pathways of neoplastic transformation. We propose and test a specific mechanism of neoplastic transformation, and genetic resistance, involving the MDV oncogene Meq, host gene products of the Nuclear Factor Kappa B (NF-κB) family and CD30; we also identify a novel Meq protein interactome.

Results

Our results show that a) CD30^lo lymphocytes are pre-neoplastic precursors and not merely reactive lymphocytes; b) multiple transformation mechanisms exist and are potentially controlled by Meq; c) Meq can drive a feed-forward cycle that induces CD30 transcription, increases CD30 signaling which activates NF-κB, and, in turn, increases Meq transcription; d) Meq transcriptional repression or activation of the CD30 promoter generally correlates with polymorphisms in the CD30 promoter distinguishing MD-lymphoma resistant and susceptible chicken genotypes e) MDV oncoprotein Meq interacts with proteins involved in physiological processes central to lymphomagenesis.

Conclusions

In the context of the MD lymphoma microenvironment (and potentially in other CD30^hi lymphomas as well), our results show that the neoplastic transformation is a continuum and the non-neoplastic cells are actually pre-neoplastic precursor cells and not merely immune bystanders. We also show that NF-κB is a central player in MDV induced neoplastic transformation of CD30-expressing lymphocytes in vivo. Our results provide insights into molecular mechanisms of neoplastic transformation in MD specifically and also herpesvirus induced lymphoma in general.

Microarray and proteomic analysis of the cardioprotective effects of cold blood cardioplegia in the mature and aged male and female

Recently we have shown that the cardioprotection afforded by cardioplegia is modulated by age and gender and is significantly decreased in the aged female. In this report we use microarray and proteomic analyses to identify transcriptomic and proteomic alterations affecting cardioprotection using cold blood cardioplegia in the mature and aged male and female heart. Mature and aged male and female New Zealand White rabbits were used for in situ blood perfused cardiopulmonary bypass. Control hearts received 30 min sham ischemia and 120 min sham reperfusion. Global ischemia (GI) hearts received 30 min of GI achieved by cross-clamping of the aorta. Cardioplegia (CP) hearts received cold blood cardioplegia prior to GI. Following 30 min of GI the hearts were reperfused for 120 min and then used for RNA and protein isolation. Microarray and proteomic analyses were performed. Functional enrichment analysis showed that mitochondrial dysfunction, oxidative phosphorylation and calcium signaling pathways were significantly enriched in all experimental groups. Glycolysis/gluconeogenesis and the pentose phosphate pathway were significantly changed in the aged male only (P < 0.05), while glyoxylate/dicarboxylate metabolism was significant in the aged female only (P < 0.05). Our data show that specific pathways associated with the mitochondrion modulate cardioprotection with CP in the aged and specifically in the aged female. The alteration of these pathways significantly contributes to decreased myocardial functional recovery and myonecrosis following ischemia and may be modulated to allow for enhanced cardioprotection in the aged and specifically in the aged female.

Gene expression and behaviour in mouse models of HD

Huntington's disease (HD) is an autosomal dominant neurodegenerative disease, resulting in expansion of the CAG repeat in exon 1 of the HTT gene. The resulting mutant huntingtin protein has been implicated in the disruption of a variety of cellular functions, including transcription. Mouse models of HD have been central to the development of our understanding of gene expression changes in this disease, and are now beginning to elucidate the relationship between gene expression and behaviour. Here, we review current mouse models of HD and their characterisation in terms of gene expression. In addition, we look at how this can inform behaviours observed in mouse models of disease. The relationship between gene expression and behaviour in mouse models of HD is important, as this will further our knowledge of disease progression and its underlying molecular events, highlight new treatment targets, and potentially provide new biomarkers for therapeutic trials.

Saliva proteomics as an emerging, non-invasive tool to study livestock physiology, nutrition and diseases

Saliva is an extraordinary fluid in terms of research and diagnostic possibilities. Its composition in electrolytes, hormones and especially its proteome contains information about feeding status, nutritional requirements and adaptations to diet and environment, and also about health status of animals. It is easy to collect on a non-invasive and routine basis without any need for special training. Therefore, the analysis of salivary proteomes is going to emerge into a field of high interest with the future goal to maintain and improve livestock productivity and welfare. Moreover, the comprehensive analysis and identification of salivary proteins and peptides in whole and glandular saliva is a necessary pre-requisite to identify animal disease biomarkers and a powerful tool to better understand animal physiology. This review focuses on the different approaches used to study the salivary proteomes of farm animals, in respect to the physiology of nutrition and food perception in relation to food choices. The potential of animal saliva as a source of disease biomarkers will also be pointed out. Special emphasis is laid on the 'ruminating triad' - cattle, goat and sheep - as well as swine as major species of animal production in Western and Southern Europe.

Insight into molecular pathways of retinal metabolism, associated with vitellogenesis in zebrafish

Retinal is the main retinoid stored in oviparous eggs of fish, amphibians, and reptiles, reaching the oocytes in association with vitellogenins, the yolk precursor proteins. During early presegmentation stages of zebrafish embryos, retinal is metabolized to retinoic acid (RA), which regulates genes involved in cell proliferation, differentiation, and tissue function and is therefore essential for normal embryonic development. While synthesis of vitellogenin and its regulation by 17β-estradiol (E(2)) were extensively investigated, pathways for retinal synthesis remain obscure. We determined the expression pattern of 46 candidate genes, aiming at identifying enzymes associated with retinal synthesis, ascertaining whether they were regulated by E(2), and finding pathways that could fulfill the demand for retinoids during vitellogenesis. Genes associated with retinal synthesis were upregulated in liver (rdh10, rdh13, sdr) and surprisingly also in intestine (rdh13) and ovary (rdh1, sdr), concomitantly with higher gene expression and synthesis of vitellogenins in liver but also in extrahepatic tissues, shown here for the first time. Vitellogenin synthesis in the ovary was regulated by E(2). Gene expression studies suggest that elevated retinal synthesis in liver, intestine, and ovary also depends on cleavage of carotenoids (by Bcdo2 or Bmco1), but in the ovary it may also be contingent on higher uptake of retinol from the circulatory system (via Stra6) and retinol synthesis from retinyl esters (by Lpl). Decrease in oxidation (by Raldh2 or Raldh3) of retinal to RA and/or degradation of RA (by Cyp26a1) may also facilitate higher hepatic retinal levels. Together, these processes enable meeting the putative demands of retinal for binding to vitellogenins. Bioinformatic tools reveal multiple hormone response elements in the studied genes, suggesting complex and intricate regulation of these processes.

Transcriptome changes in age-related macular degeneration

Age-related macular degeneration (AMD) is a debilitating, common cause of visual impairment. While the last decade has seen great progress in understanding the pathophysiology of AMD, the molecular changes that occur in eyes with AMD are still poorly understood. In the current issue of Genome Medicine, Newman and colleagues present the first systematic transcriptional profile analysis of AMD-affected tissues, providing a comprehensive set of expression data for different regions (macula versus periphery), tissues (retina versus retinal pigment epithelium (RPE)/choroid), and disease state (control versus early or advanced AMD). Their findings will serve as a foundation for additional systems-level research into the pathogenesis of this blinding disease.

Biotechnologies for the management of genetic resources for food and agriculture

In recent years, the land area under agriculture has declined as also has the rate of growth in agricultural productivity while the demand for food continues to escalate. The world population now stands at 7 billion and is expected to reach 9 billion in 2045. A broad range of agricultural genetic diversity needs to be available and utilized in order to feed this growing population. Climate change is an added threat to biodiversity that will significantly impact genetic resources for food and agriculture (GRFA) and food production. There is no simple, all-encompassing solution to the challenges of increasing productivity while conserving genetic diversity. Sustainable management of GRFA requires a multipronged approach, and as outlined in the paper, biotechnologies can provide powerful tools for the management of GRFA. These tools vary in complexity from those that are relatively simple to those that are more sophisticated. Further, advances in biotechnologies are occurring at a rapid pace and provide novel opportunities for more effective and efficient management of GRFA. Biotechnology applications must be integrated with ongoing conventional breeding and development programs in order to succeed. Additionally, the generation, adaptation, and adoption of biotechnologies require a consistent level of financial and human resources and appropriate policies need to be in place. These issues were also recognized by Member States at the FAO international technical conference on Agricultural Biotechnologies for Developing Countries (ABDC-10), which took place in March 2010 in Mexico. At the end of the conference, the Member States reached a number of key conclusions, agreeing, inter alia, that developing countries should significantly increase sustained investments in capacity building and the development and use of biotechnologies to maintain the natural resource base; that effective and enabling national biotechnology policies and science-based regulatory frameworks can facilitate the development and appropriate use of biotechnologies in developing countries; and that FAO and other relevant international organizations and donors should significantly increase their efforts to support the strengthening of national capacities in the development and appropriate use of pro-poor agricultural biotechnologies.

Comparative transcriptomics and gene expression in larval tiger salamander (Ambystoma tigrinum) gill and lung tissues as revealed by pyrosequencing

Biologists are beginning to unravel the complexities of gene expression in model organisms by studying the transcriptome, the complement of genes that are transcribed in a given tissue. It is unclear, however, if findings from model systems apply to non-model organisms because of environmental effects on gene expression. Furthermore, there have been few efforts to quantify how transcriptome or gene expression varies across individuals and across tissues in natural environments. Herein, we describe transcriptomic profiling of gene expression in lung and gill tissue of three larval tiger salamanders. We do so with a hierarchical experimental design that captures variation in expression among genes, among tissues, and among individuals. Using 454 pyrosequencing, we produced high-quality sequence data of 59 megabases and assembled ~200,000 reads into 19,501 contigs. These contigs BLASTed to 3,599 transcripts, of which 721 were expressed in both tissues, 1,668 were unique to gill, and 1,210 unique to lung. Our data showed tissue-specific patterns in gene expression level with variation among transcripts and individuals. We identified genes and gene ontology terms related to respiration and compared their relative expression levels between gill and lung tissues. We also found evidence of exogenous genes associated with larval salamanders, and we identified ~1400 potential molecular markers (microsatellites and single nucleotide polymorphisms) that are associated with expressed genes. Given the tissue-specific differences we observed in transcriptomes, these data reinforce the idea that changes in gene expression serve as a primary mechanism underlying phenotypic plasticity.

Analysis of muscle and ovary transcriptome of Sus scrofa: assembly, annotation and marker discovery

Pig (Sus scrofa) is an important organism for both agricultural and medical purpose. This study aims to investigate the S. scrofa transcriptome by the use of Roche 454 pyrosequencing. We obtained a total of 558 743 and 528 260 reads for the back-leg muscle and ovary tissue each. The overall 1 087 003 reads give rise to 421 767 341 bp total residues averaging 388 bp per read. The de novo assemblies yielded 11 057 contigs and 60 270 singletons for the back-leg muscle, 12 204 contigs and 70 192 singletons for the ovary and 18 938 contigs and 102 361 singletons for combined tissues. The overall GC content of S. scrofa transcriptome is 42.3% for assembled contigs. Alternative splicing was found within 4394 contigs, giving rise to 1267 isogroups or genes. A total of 56 589 transcripts are involved in molecular function (40 916), biological process (38 563), cellular component (35 787) by further gene ontology analyses. Comparison analyses showed that 336 and 553 genes had significant higher expression in the back-leg muscle and ovary each. In addition, we obtained a total of 24 214 single-nucleotide polymorphisms and 11 928 simple sequence repeats. These results contribute to the understanding of the genetic makeup of S. scrofa transcriptome and provide useful information for functional genomic research in future.

Identification of tissue-specific, abiotic stress-responsive gene expression patterns in wine grape (Vitis vinifera L.) based on curation and mining of large-scale EST data sets

BACKGROUND

Abiotic stresses, such as water deficit and soil salinity, result in changes in physiology, nutrient use, and vegetative growth in vines, and ultimately, yield and flavor in berries of wine grape, Vitis vinifera L. Large-scale expressed sequence tags (ESTs) were generated, curated, and analyzed to identify major genetic determinants responsible for stress-adaptive responses. Although roots serve as the first site of perception and/or injury for many types of abiotic stress, EST sequencing in root tissues of wine grape exposed to abiotic stresses has been extremely limited to date. To overcome this limitation, large-scale EST sequencing was conducted from root tissues exposed to multiple abiotic stresses.

RESULTS

A total of 62,236 expressed sequence tags (ESTs) were generated from leaf, berry, and root tissues from vines subjected to abiotic stresses and compared with 32,286 ESTs sequenced from 20 public cDNA libraries. Curation to correct annotation errors, clustering and assembly of the berry and leaf ESTs with currently available V. vinifera full-length transcripts and ESTs yielded a total of 13,278 unique sequences, with 2302 singletons and 10,976 mapped to V. vinifera gene models. Of these, 739 transcripts were found to have significant differential expression in stressed leaves and berries including 250 genes not described previously as being abiotic stress responsive. In a second analysis of 16,452 ESTs from a normalized root cDNA library derived from roots exposed to multiple, short-term, abiotic stresses, 135 genes with root-enriched expression patterns were identified on the basis of their relative EST abundance in roots relative to other tissues.

CONCLUSIONS

The large-scale analysis of relative EST frequency counts among a diverse collection of 23 different cDNA libraries from leaf, berry, and root tissues of wine grape exposed to a variety of abiotic stress conditions revealed distinct, tissue-specific expression patterns, previously unrecognized stress-induced genes, and many novel genes with root-enriched mRNA expression for improving our understanding of root biology and manipulation of rootstock traits in wine grape. mRNA abundance estimates based on EST library-enriched expression patterns showed only modest correlations between microarray and quantitative, real-time reverse transcription-polymerase chain reaction (qRT-PCR) methods highlighting the need for deep-sequencing expression profiling methods.

Comparative proteomics of seed maturation in oilseeds reveals differences in intermediary metabolism

Proteomics is increasingly being used to understand enzyme expression and regulatory mechanisms involved in the accumulation of storage reserves in crops with sequenced genomes. During the past six years, considerable progress has been made to characterize proteomes of both mature and developing seeds, particularly oilseeds - plants which accumulate principally oil and protein as storage reserves. This review summarizes the emerging proteomics data, with emphasis on seed filling in soy, rapeseed, castor and Arabidopsis as each of these oilseeds were analyzed using very similar proteomic strategies. These parallel studies provide a comprehensive view of source-sink relationships, specifically sucrose assimilation into organic acid intermediates for de novo amino acid and fatty acid synthesis. We map these biochemical processes for seed maturation and illustrate the differences and similarities among the four oilseeds. For example, while the four oilseeds appear capable of producing cytosolic phosphoenolpyruvate as the principal carbon intermediate, soybean and castor also express malic enzymes and malate dehydrogenase, together capable of producing malate that has been previously proposed to be the major intermediate for fatty acid synthesis in castor. We discuss these and other differences in the context of intermediary metabolism for developing oilseeds.

Toward a complete in silico, multi-layered embryonic stem cell regulatory network

Recent efforts in systematically profiling embryonic stem (ES) cells have yielded a wealth of high-throughput data. Complementarily, emerging databases and computational tools facilitate ES cell studies and further pave the way toward the in silico reconstruction of regulatory networks encompassing multiple molecular layers. Here, we briefly survey databases, algorithms, and software tools used to organize and analyze high-throughput experimental data collected to study mammalian cellular systems with a focus on ES cells. The vision of using heterogeneous data to reconstruct a complete multi-layered ES cell regulatory network is discussed. This review also provides an accompanying manually extracted dataset of different types of regulatory interactions from low-throughput experimental ES cell studies available at http://amp.pharm.mssm.edu/iscmid/literature.

Integration of deep transcriptome and proteome analyses reveals the components of alkaloid metabolism in opium poppy cell cultures

BACKGROUND

Papaver somniferum (opium poppy) is the source for several pharmaceutical benzylisoquinoline alkaloids including morphine, the codeine and sanguinarine. In response to treatment with a fungal elicitor, the biosynthesis and accumulation of sanguinarine is induced along with other plant defense responses in opium poppy cell cultures. The transcriptional induction of alkaloid metabolism in cultured cells provides an opportunity to identify components of this process via the integration of deep transcriptome and proteome databases generated using next-generation technologies.

RESULTS

A cDNA library was prepared for opium poppy cell cultures treated with a fungal elicitor for 10 h. Using 454 GS-FLX Titanium pyrosequencing, 427,369 expressed sequence tags (ESTs) with an average length of 462 bp were generated. Assembly of these sequences yielded 93,723 unigenes, of which 23,753 were assigned Gene Ontology annotations. Transcripts encoding all known sanguinarine biosynthetic enzymes were identified in the EST database, 5 of which were represented among the 50 most abundant transcripts. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) of total protein extracts from cell cultures treated with a fungal elicitor for 50 h facilitated the identification of 1,004 proteins. Proteins were fractionated by one-dimensional SDS-PAGE and digested with trypsin prior to LC-MS/MS analysis. Query of an opium poppy-specific EST database substantially enhanced peptide identification. Eight out of 10 known sanguinarine biosynthetic enzymes and many relevant primary metabolic enzymes were represented in the peptide database.

CONCLUSIONS

The integration of deep transcriptome and proteome analyses provides an effective platform to catalogue the components of secondary metabolism, and to identify genes encoding uncharacterized enzymes. The establishment of corresponding transcript and protein databases generated by next-generation technologies in a system with a well-defined metabolite profile facilitates an improved linkage between genes, enzymes, and pathway components. The proteome database represents the most relevant alkaloid-producing enzymes, compared with the much deeper and more complete transcriptome library. The transcript database contained full-length mRNAs encoding most alkaloid biosynthetic enzymes, which is a key requirement for the functional characterization of novel gene candidates.

Computational identification of root hair-specific genes in Arabidopsis

Activated cortical domains (ACDs) are regions of the plant cell cortex performing localized membrane turnover, delimited by concerted action of the cortical cytoskeleton and endomembrane compartments. Arabidopsis thaliana rhizodermis consists of two cell types differing by a single ACD (trichoblasts, carrying tip-growing root hairs, and hairless atrichoblasts), providing a model for the study of ACD determination. We compiled a set of genes specifically upregulated in root hairs from published transcriptome data, and compared it with a "virtual Arabidopsis root hair proteome", i.e. a list of computationally identified homologs of proteins from the published soybean root hair proteome. Both data sets were enriched in genes and proteins associated with root hairs in functional studies, but there was little overlap between the transcriptome and the proteome: the former captured gene products specific to root hairs, while the latter selected those abundant in root hairs but not necessarily specific to them. Decisive steps in ACD specification may be performed by signaling proteins of high expression specifity and low abundance. Nevertheless, 73 genes specifically transcribed in Arabidopsis trichoblasts or root hairs encode homologs of abundant root hair proteins from soybean. Most of them encode "housekeeping" proteins required for rapid tip growth. However, among the "candidates" is also a generative actin isoform, ACT11. Preliminary characterization of an act11 mutant allele indeed suggests a hitherto unexpected role for this gene in root and root hair development.

Breast cancer genomics: normal tissue and cancer markers

Mammography is a powerful screening tool for early detection of breast cancer, but it has limitations in terms of both specificity and sensitivity. Imaging tools such as MRI that complement mammography are too costly to serve as first-line screens. Recently, progress has been made on blood markers, particularly microRNAs and proteins. There are new methods for protein marker discovery directly in blood, but they are limited in the number of patients that can be examined. An alternative is to discover markers as transcripts in tissues, followed by development of blood protein tests for those that perform best. To identify genes that are overexpressed in malignancy it is paramount to include normal control tissues from healthy individuals. Here we report the identification of potential breast cancer markers, including some that are overexpressed in aggressive disease.

Proteomic survey of the cestode Mesocestoides corti during the first 24 hours of strobilar development

Despite the fact that cestodes represent major etiological agents of both human and domestic animal diseases, little is known about the molecular aspects of cestode development. In this work, Mesocestoides corti, a model cestode species, was studied from the early development of its larval form (tetrathyridium) into adult worms (strobilation) using different proteomic approaches. The protein profiles of M. corti tetrathyridia induced or not induced to undergo strobilation were compared. Proteomic mapping by two-dimensional gel electrophoresis showed the resolution of 248 and 154 spots from tetrathyridia that were subjected or not subjected to strobilation induction, respectively, allowing for the detection of at least nine spots exclusive to each group. Spot analysis by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) or MALDI-TOF MS/MS identified four reference proteins (six spots). LC-MS/MS analyses of protein extracts identified 66 proteins, eight of which were found exclusively in non-induced tetrathyridia, while 13 were found exclusively in strobilation-induced tetrathyridia. Among the proteins exclusively identified in strobilation-induced worms, there was a predominance of proteins with functions relating to chaperone activity and protein synthesis and turnover. Quantitative differential expression analysis between M. corti tetrathyridia prior to and after strobilation induction revealed six proteins upregulated in strobilation-induced worms; these proteins were involved in metabolic pathways, cell proliferation, and cytoskeletal rearrangement. Overall, despite the absence of a sequenced M. corti genome, using sequences from other platyhelminthes, we were able to establish comprehensive protein profiles for tetrathyridia prior to and after strobilation induction and identify several proteins potentially involved in the early events leading to strobilation.

Response of porcine intestinal in vitro organ culture tissues following exposure to Lactobacillus plantarum JC1 and Salmonella enterica serovar Typhimurium SL1344

The development of novel intervention strategies for the control of zoonoses caused by bacteria such as Salmonella spp. in livestock requires appropriate experimental models to assess their suitability. Here, a novel porcine intestinal in vitro organ culture (IVOC) model utilizing cell crown (CC) technology (CCIVOC) (Scaffdex) was developed. The CCIVOC model was employed to investigate the characteristics of association of S. enterica serovar Typhimurium strain SL1344 with porcine intestinal tissue following exposure to a Lactobacillus plantarum strain. The association of bacteria to host cells was examined by light microscopy and electron microscopy (EM) after appropriate treatments and staining, while changes in the proteome of porcine jejunal tissues were investigated using quantitative label-free proteomics. Exposure of porcine intestinal mucosal tissues to L. plantarum JC1 did not reduce the numbers of S. Typhimurium bacteria associating to the tissues but was associated with significant (P < 0.005) reductions in the percentages of areas of intestinal IVOC tissues giving positive staining results for acidic mucins. Conversely, the quantity of neutrally charged mucins present within the goblet cells of the IVOC tissues increased significantly (P < 0.05). In addition, tubulin-α was expressed at high levels following inoculation of jejunal IVOC tissues with L. plantarum. Although L. plantarum JC1 did not reduce the association of S. Typhimurium strain SL1344 to the jejunal IVOC tissues, detection of increased acidic mucin secretion, host cytoskeletal rearrangements, and proteins involved in the porcine immune response demonstrated that this strain of L. plantarum may contribute to protecting the pig from infections by S. Typhimurium or other pathogens.

Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon

BACKGROUND

Several recent studies have demonstrated the use of Roche 454 sequencing technology for de novo transcriptome analysis. Low error rates and high coverage also allow for effective SNP discovery and genetic diversity estimates. However, genetically diverse datasets, such as those sourced from natural populations, pose challenges for assembly programs and subsequent analysis. Further, estimating the effectiveness of transcript discovery using Roche 454 transcriptome data is still a difficult task.

RESULTS

Using the Roche 454 FLX Titanium platform, we sequenced and assembled larval transcriptomes for two butterfly species: the Propertius duskywing, Erynnis propertius (Lepidoptera: Hesperiidae) and the Anise swallowtail, Papilio zelicaon (Lepidoptera: Papilionidae). The Expressed Sequence Tags (ESTs) generated represent a diverse sample drawn from multiple populations, developmental stages, and stress treatments. Despite this diversity, > 95% of the ESTs assembled into long (> 714 bp on average) and highly covered (> 9.6x on average) contigs. To estimate the effectiveness of transcript discovery, we compared the number of bases in the hit region of unigenes (contigs and singletons) to the length of the best match silkworm (Bombyx mori) protein--this "ortholog hit ratio" gives a close estimate on the amount of the transcript discovered relative to a model lepidopteran genome. For each species, we tested two assembly programs and two parameter sets; although CAP3 is commonly used for such data, the assemblies produced by Celera Assembler with modified parameters were chosen over those produced by CAP3 based on contig and singleton counts as well as ortholog hit ratio analysis. In the final assemblies, 1,413 E. propertius and 1,940 P. zelicaon unigenes had a ratio > 0.8; 2,866 E. propertius and 4,015 P. zelicaon unigenes had a ratio > 0.5.

CONCLUSIONS

Ultimately, these assemblies and SNP data will be used to generate microarrays for ecoinformatics examining climate change tolerance of different natural populations. These studies will benefit from high quality assemblies with few singletons (less than 26% of bases for each assembled transcriptome are present in unassembled singleton ESTs) and effective transcript discovery (over 6,500 of our putative orthologs cover at least 50% of the corresponding model silkworm gene).

Systems analysis of seed filling in Arabidopsis: using general linear modeling to assess concordance of transcript and protein expression

Previous systems analyses in plants have focused on a single developmental stage or time point, although it is often important to additionally consider time-index changes. During seed development a cascade of events occurs within a relatively brief time scale. We have collected protein and transcript expression data from five sequential stages of Arabidopsis (Arabidopsis thaliana) seed development encompassing the period of reserve polymer accumulation. Protein expression profiling employed two-dimensional gel electrophoresis coupled with tandem mass spectrometry, while transcript profiling used oligonucleotide microarrays. Analyses in biological triplicate yielded robust expression information for 523 proteins and 22,746 genes across the five developmental stages, and established 319 protein/transcript pairs for subsequent pattern analysis. General linear modeling was used to evaluate the protein/transcript expression patterns. Overall, application of this statistical assessment technique showed concurrence for a slight majority (56%) of expression pairs. Many specific examples of discordant protein/transcript expression patterns were detected, suggesting that this approach will be useful in revealing examples of posttranscriptional regulation.

Microarrays and high-throughput transcriptomic analysis in species with incomplete availability of genomic sequences

Microarrays produce a measurement of gene expression based on the relative measures of dye intensities that correspond to the amount of target RNA. This technology is fast developing and its application is expanding from Homo sapiens to a wide number of species, where enough information on sequences and annotations exist. Anyway, the number of species for which a dedicated platform exists is not high. The use of heterologous array hybridization, screening for gene expression in one species using an array developed for another one, is still quite frequent, even though cross-species microarray hybridization has raised many arguments. Some methods which are high throughput and do not rely on knowledge of the DNA/RNA sequence exist, namely serial analysis of gene expression (SAGE), Massively Parallel Signature Sequencing (MPSS) and deep sequencing of full transcriptome. Although very powerful, particularly the latter, they are still quite costly and cumbersome methods. In some species where genome sequences are largely unknown, several anonymous sequences are deposited in gene banks as a result of Expressed Sequence Tags (ESTs) sequencing projects. The ESTs databases represent a valuable knowledge that can be exploited with some bioinformatic effort to build species-specific microarrays. We present here a method of high-density in situ synthesized microarrays starting from available EST sequences in, Ovis aries. Our data indicate that the method is very efficient and can be easily extended to other species of which genetic sequences are present in public databases, but neglected so far with advanced devices like microarrays. As a perspective, the approach can be applied also to species of which no sequences are available to date, thanks to high-throughput deep sequencing methods.