The EADGENE Microarray Data Analysis Workshop (open access publication).

Microarray analyses have become an important tool in animal genomics. While their use is becoming widespread, there is still a lot of ongoing research regarding the analysis of microarray data. In the context of a European Network of Excellence, 31 researchers representing 14 research groups from 10 countries performed and discussed the statistical analyses of real and simulated 2-colour microarray data that were distributed among participants. The real data consisted of 48 microarrays from a disease challenge experiment in dairy cattle, while the simulated data consisted of 10 microarrays from a direct comparison of two treatments (dye-balanced). While there was broader agreement with regards to methods of microarray normalisation and significance testing, there were major differences with regards to quality control. The quality control approaches varied from none, through using statistical weights, to omitting a large number of spots or omitting entire slides. Surprisingly, these very different approaches gave quite similar results when applied to the simulated data, although not all participating groups analysed both real and simulated data. The workshop was very successful in facilitating interaction between scientists with a diverse background but a common interest in microarray analyses.

Analysis of the real EADGENE data set: multivariate approaches and post analysis (open access publication).

The aim of this paper was to describe, and when possible compare, the multivariate methods used by the participants in the EADGENE WP1.4 workshop. The first approach was for class discovery and class prediction using evidence from the data at hand. Several teams used hierarchical clustering (HC) or principal component analysis (PCA) to identify groups of differentially expressed genes with a similar expression pattern over time points and infective agent (E. coli or S. aureus). The main result from these analyses was that HC and PCA were able to separate tissue samples taken at 24 h following E. coli infection from the other samples. The second approach identified groups of differentially co-expressed genes, by identifying clusters of genes highly correlated when animals were infected with E. coli but not correlated more than expected by chance when the infective pathogen was S. aureus. The third approach looked at differential expression of predefined gene sets. Gene sets were defined based on information retrieved from biological databases such as Gene Ontology. Based on these annotation sources the teams used either the GlobalTest or the Fisher exact test to identify differentially expressed gene sets. The main result from these analyses was that gene sets involved in immune defence responses were differentially expressed.

Analysis of a simulated microarray dataset: comparison of methods for data normalisation and detection of differential expression (open access publication).

Microarrays allow researchers to measure the expression of thousands of genes in a single experiment. Before statistical comparisons can be made, the data must be assessed for quality and normalisation procedures must be applied, of which many have been proposed. Methods of comparing the normalised data are also abundant, and no clear consensus has yet been reached. The purpose of this paper was to compare those methods used by the EADGENE network on a very noisy simulated data set. With the a priori knowledge of which genes are differentially expressed, it is possible to compare the success of each approach quantitatively. Use of an intensity-dependent normalisation procedure was common, as was correction for multiple testing. Most variety in performance resulted from differing approaches to data quality and the use of different statistical tests. Very few of the methods used any kind of background correction. A number of approaches achieved a success rate of 95% or above, with relatively small numbers of false positives and negatives. Applying stringent spot selection criteria and elimination of data did not improve the false positive rate and greatly increased the false negative rate. However, most approaches performed well, and it is encouraging that widely available techniques can achieve such good results on a very noisy data set.

Analysis of the real EADGENE data set: comparison of methods and guidelines for data normalisation and selection of differentially expressed genes (open access publication).

A large variety of methods has been proposed in the literature for microarray data analysis. The aim of this paper was to present techniques used by the EADGENE (European Animal Disease Genomics Network of Excellence) WP1.4 participants for data quality control, normalisation and statistical methods for the detection of differentially expressed genes in order to provide some more general data analysis guidelines. All the workshop participants were given a real data set obtained in an EADGENE funded microarray study looking at the gene expression changes following artificial infection with two different mastitis causing bacteria: Escherichia coli and Staphylococcus aureus. It was reassuring to see that most of the teams found the same main biological results. In fact, most of the differentially expressed genes were found for infection by E. coli between uninfected and 24 h challenged udder quarters. Very little transcriptional variation was observed for the bacteria S. aureus. Lists of differentially expressed genes found by the different research teams were, however, quite dependent on the method used, especially concerning the data quality control step. These analyses also emphasised a biological problem of cross-talk between infected and uninfected quarters which will have to be dealt with for further microarray studies.

Biochemical pathways analysis of microarray results: regulation of myogenesis in pigs.

BACKGROUND: Combining microarray results and biological pathway information will add insight into biological processes. Pathway information is widely available in databases through the internet. Mammalian muscle formation has been previously studied using microarray technology in pigs because these animals are an interesting animal model for muscle formation due to selection for increased muscle mass. Results indicated regulation of the expression of genes involved in proliferation and differentiation of myoblasts, and energy metabolism. The aim of the present study was to analyse microarrays studying myogenesis in pigs. It was necessary to develop methods to search biochemical pathways databases. RESULTS: PERL scripts were developed that used the names of the genes on the microarray to search databases. Synonyms of gene names were added to the list by searching the Gene Ontology database. The KEGG database was searched for pathway information using this updated gene list. The KEGG database returned 88 pathways. Most genes were found in a single pathway, but others were found in up to seven pathways. Combining the pathways and the microarray information 21 pathways showed sufficient information content for further analysis. These pathways were related to regulation of several steps in myogenesis and energy metabolism. Pathways regulating myoblast proliferation and muscle fibre formation were described. Furthermore, two networks of pathways describing the formation of the myoblast cytoskeleton and regulation of the energy metabolism during myogenesis were presented. CONCLUSION: Combining microarray results and pathways information available through the internet provide biological insight in how the process of porcine myogenesis is regulated.

Transcriptome expression profiles in prenatal pigs in relation to myogenesis.

Myogenesis, the formation of muscle fibers, is a complex process. Pigs have been selected for efficient muscle growth for the past decades making them interesting to study myogenesis. We studied expression profiles of genes known to affect myogenesis, muscle structural proteins, and energy metabolism in prenatal pigs from 14 to 91 days of gestation. Primary and secondary muscle fiber formation takes place during days 30-60 and 54-90 of gestation, respectively. Differential expression and expression levels of the genes were studied using microarray technology. Gene activation and repression profiles were studied counting the number of spots with detectable signal. The number of spots for muscle tissue structural protein genes showing upregulated expression increased constantly from day 14 until day 91 of gestation indicating continued activation of genes during this period. The mRNA expression level of the genes showed a peak around day 35 of gestation. The expression levels of genes affecting myogenic differentiation (stimulating and inhibiting) showed a peak at day 35 of gestation. The number of spots for differentiation-stimulating genes showing differential expression reaches a first peak around day 35 of gestation and a nadir at day 49 of gestation while the number of spots for differentiation-inhibiting genes reaches a nadir at day 35 of gestation. Myogenic differentiation seems less a matter of the expression level of genes affecting differentiation, but depends on the balance between the number of significantly activated genes for stimulating and inhibiting differentiation. Genes stimulating myoblast proliferation showed a small peak expression prior to day 35 of gestation indicating myoblast proliferation before differentiation. The number of spots and the expression levels of genes for glycolysis and ATP-metabolism are at a nadir around days 35 and 49-63 of gestation suggesting that the energy metabolism is low during fusion of myoblasts into multinucleated muscle fibers.