REPRODUCTION REVIEW A guide to issues in microarray analysis: application to endometrial biology Christine A White1,2 and Lois A Salamonsen1 1 Prince Henry’s Institute of Medical Research, PO Box 5152, Clayton, Victoria, 3168, Australia and 2Dept of Obstetrics & Gynaecology, Monash University, Clayton, Victoria, 3168, Australia Correspondence should be addressed to C A White; Email: [email protected] Abstract Within the last decade, the development of DNA microarray technology has enabled the simultaneous measurement of thousands of gene transcripts in a biological sample. Conducting a microarray study is a multi-step process; starting with a welldefined biological question, moving through experimental design, target RNA preparation, microarray hybridisation, image acquisition and data analysis – finishing with a biological interpretation requiring further study. Advances continue to be made in microarray quality and methods of statistical analysis, improving the reliability and therefore appeal of microarray analysis for a wide range of biological questions. The purpose of this review is to provide both an introduction to microarray methodology, as well as a practical guide to the use of microarrays for gene expression analysis, using endometrial biology as an example of the applications of this technology. While recommendations are based on previous experience in our laboratory, this review also summarises the methods currently considered to be best practice in the field. Reproduction (2005) 130 1–13 Principles of microarray analysis As for many other techniques used in molecular biology, microarrays rely on the complementarity of the DNA duplex, i.e. that the two strands will always reassemble with base pairing A to T and C to G. In addition, singlestranded DNA will bind strongly to a solid support, where it is available for hybridisation with complementary DNA (cDNA). In a DNA microarray experiment, many genespecific ‘probes’ are immobilised on a solid support (usually nylon membrane or glass) and the array is exposed to labelled cDNA ‘targets’ derived from one or more biological samples (Schena et al. 1995). Nylon membrane arrays are normally hybridised with a single cDNA population, labelled with a radioactive or chemiluminescent tag, so that the intensity of the signal generated by each bound probe indicates the abundance of that transcript in the sample. Transcript abundance in different samples (e.g. treatment and control) can then be compared across serial or parallel array hybridisations. The advantage of the glass DNA microarray is that cDNA from two or more biological samples can be labelled with different fluorescent dyes and competitively hybridised, so that the relative abundance of gene transcripts can be determined by the fluorescent signal obtained (Figure 1). This review will focus on the techniques and data analysis associated with the two most common cDNA microarray q 2005 Society for Reproduction and Fertility ISSN 1470–1626 (paper) 1741–7899 (online) platforms; dual colour fluorescence on glass and radioactively labelled nylon membranes. A summary of some of the key issues involved in microarray analysis is provided in Table 1. Choice of microarray As most laboratories are not equipped with robotic printers, the investigator must usually obtain microarrays from commercial or other sources. There are a range of mouse, rat and human microarrays available, and they fall into two broad categories; cDNA and high-density synthetic oligonucleotide (reviewed in Barrett & Kawasaki 2003). As described above, cDNA microarrays can be based on either nylon membrane or dual colour fluorescence on glass. Both formats of cDNA microarray involve the deposition of purified and/or PCR-amplified DNA in solution onto the solid support at defined locations using a robotic printer (Figure 1). The quality of the microarray is therefore dependent on the performance of the print tips, which must deliver reproducible volumes and uniform spot sizes to enable effective data analysis. The main advantages of cDNA microarrays are their relatively low cost (, US $150 per slide) and greater flexibility in terms of producing and spotting custom-made clone sets. Manufacturers of high-density synthetic oligonucleotide microarrays, such as Affymetrix (Santa Clara, CA, USA; DOI: 10.1530/rep.1.00685 Online version via www.reproduction-online.org 2 C A White and L A Salamonsen Table 1 Some key issues involved in microarray analysis Parameter Issue † Consider the biological question(s) and the ability to achieve statistical significance † Seek expert statistical advice during the early planning stages † Microarray experiments have multiple sources of variation and must be carefully controlled † Biological and technical replication are essential † Sample pooling should be avoided if accurate sample synchronisation is not possible † Microarray analysis of purified cells will only reveal genes expressed by these cells, but removal from the in vivo microenvironment may alter gene expression † There are limitations in the use of both whole tissue and purified cells, which may necessitate the use of microdissection and RNA amplification techniques † When using clinical samples, detailed patient history and tissue histopathology are critical to the interpretation of gene expression profiles Target RNA preparation † The quality of the target RNA is one of the most important factors in the success or failure of a microarray experiment Data analysis † While critical to the outcome of a microarray experiment, statistical analysis of microarray data is not well understood by many biologists and expert advice should be sought Data validation † The biomedical research community does not yet accept that microarray data can stand alone without independent validation † The investigator must decide which genes to examine further, and those with larger fold changes and statistical significance are often the best candidates † To describe a biological event or system, gene expression data obtained by microarray analysis must be extended to the study of protein products Experimental design http://www.affymetrix.com), use photolithography and solid-phase DNA synthesis to generate synthetic 25 base polynucleotides (25mers) directly on the glass surface (Lipschutz et al. 1999). Each gene is represented by 11 to 20 different 25 mers, in ‘perfect match’ or ‘mismatch’ sequence pairs. Probes can be generated representing a unique part of a gene transcript, enabling discrimination between closely related genes or splice variants and the ‘mismatch’ sequences provide an internal control for every gene. Longer oligonucleotide (50 to 100 mers) microarrays are also available, which provide even greater hybridisation specificity (Barrett & Kawasaki 2003). Oligonucleotide microarrays are hybridised with a single fluorescently labelled sample and gene expression in different samples compared across multiple microarrays. The main disadvantage of these microarrays is their high cost (up to Reproduction (2005) 130 1–13 US$800 per slide), so their appeal is likely to increase as they become more economical. Commercially-available microarrays are printed or synthesised with a particular clone set (reviewed in Bowtell 1999) and these differ in the proportion of known genes and expressed sequence tags (ESTs). Some ESTs correspond to a segment of a known gene, but most represent partially sequenced novel genes. A number of groups including those at Merck, Washington University, the IMAGE (Integrated Analysis of Genomes and their Expression) Consortium and the Cancer Genome Anatomy Project (CGAP) have been responsible for sequencing over one million human ESTs (Bowtell 1999). The online database dbEST (a division of GenBank; http://www.ncbi.nlm.nih.gov/dbEST/ index.html) houses all these EST sequences and the automated process known as UniGene assigns overlapping sequences to a single cluster, which may or may not have a known identity (http://www.ncbi.nlm.nih.gov/UniGene/ index.html). As the human genome sequence is effectively complete (Venter et al. 2001) it is expected that all ESTs will progressively be assigned an identity. Until then, careful consideration should be given to whether identifying differentially expressed ESTs is a priority in any particular microarray experiment. If not, then using a more tailored array specific to a cellular process or pathway may be more appropriate and cost-effective. Experimental design The most important considerations in microarray experimental design are the biological question under study and the ability to achieve statistical significance (reviewed in Churchill 2002, Kerr & Churchill 2001, Smyth et al. 2003, Yang & Speed 2002). Answering the biological question may require identification of downstream target genes, and may also involve a time course. Designing a microarray experiment with sufficient statistical power requires input from a statistician or bioinformatician with experience in microarray technology. The right statistical advice during early planning can save vast amounts of time and money. The reason for their complexity is that microarray experiments have multiple sources of variation, each of which must be considered in the experimental design (reviewed in Churchill 2002, Chen et al. 2004). Firstly, there is biological variation between animals or patients. Secondly, technical variation arises from the RNA extraction, reverse transcription, label incorporation and hybridisation steps. Thirdly, measurement errors occur due to differences in hybridisation efficiencies between spots, between different print-tip groups across the array and between slides in the same and different print runs. Technical variation and measurement error can also interact. For example, the scanning properties of the fluorescent dyes can vary with the spot intensity and spatial position on the slide (Smyth & Speed 2003). A conventional power analysis requires prior knowledge of the variance of individual measurements, the www.reproduction-online.org Guide to issues in microarray analysis 3 Figure 1 Dual colour fluorescence cDNA microarray analysis. Clones of interest (probes) are amplified by PCR and printed onto treated glass slides using a robotic printer. Total RNA samples extracted from treated and control cells/tissues (targets) are reverse transcribed and labelled with either Cy3 (green) or Cy5 (red). The samples are combined and competitively hybridised to the microarray under stringent conditions. Following washing to remove non-hybridised target, laser excitation is applied and the emissions measured in each colour channel. Specialised software is used to attach gene names, fluorescence intensity values and intensity ratios to each spot, which are then exported for advanced statistical analysis to identify differentially expressed genes. magnitude of the effect to be detected, the acceptable false-positive rate and the desired ‘power’ of the calculation; that is, the probability of detecting an effect of the specified or greater magnitude (Yang & Speed 2002, Yang et al. 2003, Chen et al. 2004). In a microarray experiment, two of these components are unknown; both the variance of the expression ratio measurements and the magnitude of the effects of interest will be different for every gene on the microarray. To overcome this, power calculations can be performed using the median variance across all of the genes in a previous microarray hybridisation (Yang & Speed 2002). It is critical that microarray experiments are carefully controlled, particularly when using dual colour fluorescence microarrays in which the endpoint is a ratio of expression between two or more samples. As in any experiment, treatment controls must be carefully incorporated into the study design. To ensure that there is only one source of experimental variation, consistency must also be applied to tissue collection, processing and RNA extraction, as well as the microarray hybridisations. Even with a single variable, such as a differentiation stimulus, it is possible to end up comparing cells or tissues in completely different physiological states. In this situation, differentially expressed genes will likely be the consequence, rather than the cause, of the differences in phenotype. This problem can be minimised by using carefully controlled inducible systems and examining early rather than later time points. When an experiment involves comparisons across multiple dual colour fluorescence microarrays, there are www.reproduction-online.org a number of possible design matrices (Figure 2). Hybridising an appropriate reference sample to each microarray (Figure 2A) can provide a consistent control across multiple slides. The ideal reference contains all possible mRNA transcripts present in the experimental samples, so that fluorescence ratio measurements less than zero cannot occur. This is usually achieved by generating a pool of RNA from multiple samples of the tissue or cell type under study. Another approach is to create a reference mixture of all the PCR products spotted on the microarray (Sterrenburg et al. 2002). Although there are some advantages to using a reference sample, it is more precise and economical to make the critical comparisons directly on the same microarray (Figure 2B; Kerr & Churchill 2001, Churchill 2002). As a reference design requires more complicated statistical analysis, it should only be used for a well-defined purpose. It may also be useful to use a common reference if a large number of experimental samples are to be collected and analysed over a long period of time. A saturated design (Figure 2C) may be used when an experiment has more than two treatments, and all comparisons are of interest in answering the biological question. The efficiency of time course experiments can be maximised by a loop design (Figure 2D). Regardless of the design matrix, care should be taken that all pairings are biologically relevant and controlled. For example, pairs could be wild type and knockout littermates, or isolated cells from the same endometrial biopsy with and without treatment. Replication is critical in a microarray experiment, as it enables the data to be effectively analysed using formal Reproduction (2005) 130 1–13 4 C A White and L A Salamonsen A B REF A1 A1 B1 REF A2 A2 B2 REF A3 A3 B3 REF B1 REF B2 REF B3 C B1 A2 A1 B3 B2 A3 D Time 0 Time 1 Time 2 Time 3 Figure 2 Microarray experimental designs (adapted from Churchill 2002, Smyth et al. 2003). Letters refer to different treatments/genotypes and subscripts indicate biological replicates. Each arrow represents one microarray, with the arrow pointing away from the Cy3 (green) labelled sample and towards the Cy5 (red) labelled sample. Double arrows indicate dye-swap pairs. A, Indirect comparison with a common reference; B, Direct comparison; C, Saturated design and D, Time course experiment loop design. statistical methods, and the results of that analysis to be broadly applicable to the sampled population. Given the sources of variation described above, replication should be both biological and technical. Biological replication is the use of RNA samples from multiple animals or patients. Technical replication (i.e. repeated measures) includes the presence of duplicate DNA probes on the microarray and hybridisation of the same RNA sample on multiple microarrays. As biological variability is usually greater than technical variability, the hybridisation of independent RNA samples should be prioritised. Replication allows the investigator to identify and remove false-positives and false-negatives, so only reproducible data are considered for further analysis. A single hybridisation may be justified if the aim is to generate hypotheses for further testing, but advanced statistical analysis will be more productive on 3 or more replicates per treatment (reviewed in Sasik et al. 2004). In the case of dual colour fluorescence microarrays, dye swap replicates should also be performed to account for unequal dye incorporation and quenching (Figure 2). To avoid or minimise bias, the assignment of dye labels should be randomised. Reproduction (2005) 130 1–13 Pooling samples is often considered when RNA is in limited supply, or to minimise the effects of biological variation. In addition, it has been demonstrated that pooling RNA from an increased number of subjects can reduce the number of microarrays required, without any loss of precision (Kendziorski et al. 2003, Peng et al. 2003). Sample pooling should be avoided when it is not possible to accurately synchronise the samples (Sasik et al. 2004). For example, when using pseudopregnant mice on the same day after vaginal plug detection, there is likely to be variation in the time at which mating occurred and therefore in the physiological state of the uterus. In this situation, it is better to maintain sample independence and find common gene expression features at the data analysis stage. Careful consideration should be given to the use of whole tissue or purified homogeneous cell populations in a microarray study. The advantage of using whole tissue is that there is a greater amount of RNA available for technical replicates and subsequent validation studies (see ‘Data validation’ below). It is important to consider, however, that tissues contain a range of different cell types. Whole endometrial biopsies, for example, contain luminal www.reproduction-online.org Guide to issues in microarray analysis and glandular epithelium, non-decidualised and decidualised stroma, endothelial cells, smooth muscle cells and leukocytes. If the cell type of interest makes up a small proportion of the total tissue, whole tissue gene expression data may not be particularly informative. Microarray analysis of purified cells will only reveal genes expressed by those cells, but they will also have been removed from their in vivo microenvironment and cultured under conditions which are likely to alter gene expression. The limitations of both approaches need to be balanced against the aims of the experiment and perhaps additional technologies such as laser capture microdissection (EmmertBuck et al. 1996) considered. New RNA amplification methods (see ‘Target RNA preparation’ below) have improved the feasibility of using microdissected tissue for gene expression studies and this approach has the advantage of maintaining close to in vivo cellular context. The parallel measurement of other biological parameters can be used to assist the interpretation of microarray data. For example, variation in prolactin secretion levels from different preparations of decidualising human endometrial stromal cells may correlate with variations in the expression levels of other genes. When using clinical samples, patient history and tissue histopathology are also critical in the final interpretation of gene expression profiles. Target RNA preparation The target in microarray analysis is a labelled population of cDNAs representing the mRNA repertoire of a cell type or tissue. The purity of the extracted total RNA is one of the most critical factors in the success or failure of a microarray experiment. Residual chloroform, phenol or ethanol from the extraction process can interfere in the efficiency of reverse transcription to cDNA, and other contaminants such as cellular protein, lipids and carbohydrates can cause non-specific binding of fluorescent cDNAs to the glass surface (Duggan et al. 1999). For most tissues, including human endometrium, optimal quality total RNA can be extracted using TRIzol Reagent (Invitrogen, Carlsbad, CA, USA), followed by column purification using the RNeasy Minikit (Qiagen, Hilden, Germany). Highly pure total RNA can be extracted from isolated cells using the RNeasy kit alone. Both of these methods should be followed by removal of genomic DNA. It is recommended that a number of different RNA extraction methods be compared for the tissue or cell type of interest, to maximise RNA quality and yield. Particularly in microarray experiments, it is critical to accurately measure RNA quality and quantity, to minimise variation and therefore improve labelling and hybridisation consistency. The standard UV spectrophotometer is useful for an initial estimate of RNA quality and quantity. Optical density (OD) can be measured at 230, 260 and 280 nm and RNA purity considered acceptable at values of OD260/OD280 1.8 –2.0 and OD230/OD280 , 1.0. Agarose gel electrophoresis may then be used to further www.reproduction-online.org 5 confirm RNA integrity. However, other systems such as the RiboGreen Assay (Molecular Probes, Eugene, OR, USA) or the Agilent 2100 Bioanalyser (Agilent Technologies UK Ltd, Cheadle, Cheshire, UK) are more sensitive and accurate. The Agilent Bioanalyser requires only 50 – 500 ng of total RNA and produces a detailed electrophoretogram which will reveal any RNA degradation or genomic DNA contamination. The amount of RNA required per hybridisation is the greatest limitation to the use of this technology, particularly when the tissue of interest is in limited supply, or when using isolated cell populations. Until very recently, it was recommended that 50 –200 mg total RNA per sample be used for each hybridisation to generate a sufficient signal (Duggan et al. 1999). However, improved RNA purification, fluorescent labelling methods and hybridisation conditions have reduced this requirement to 5 – 10 mg for both glass and nylon membrane arrays. As this amount of RNA may still be difficult to obtain in some systems, a number of different RNA amplification approaches have been developed. The most commonly used method is T7 polymerase in vitro transcription (IVT; van Gelder et al. 1990). While this method can significantly reduce the RNA requirement for each hybridisation, it is also expensive, time-consuming and labour-intensive. In addition, multiple rounds of amplification may be required, which decreases the linearity of the amplification and may result in a cDNA target which is no longer representative of the original sample (Petalidis et al. 2003). Newly developed PCR-based cDNA amplification techniques can decrease the amount of starting total RNA required to 200 ng, while maintaining amplification linearity (Petalidis et al. 2003). This method vastly improves the feasibility of glass microarray studies on clinical samples. Microarray hybridisation Effective hybridisation of the target to the microarray is essential in obtaining high quality data. Coverslips with raised Teflon edging (Lifter Slips; Erie Scientific, Portsmouth, NH, USA) are a useful way of ensuring the full hybridisation volume maintains contact with the arrayed probes. Hybridisation chambers, such as those available from Corning (Corning, NY, USA), allow the microarrays to be submerged in water of a set temperature, greatly improving hybridisation consistency across the slide. Commercially available buffers such as ExpressHyb Hybridisation Solution (Clontech, Palo Alto, CA, USA) are very effective for nylon membrane microarray hybridisation (Evans et al. 2003). The use of non-specific blocking agents including Cot-1 DNA, polyadenylic acid and salmon sperm DNA ensures that the signal detected from each spot is specific to the particular probe sequence and background is minimised. Most microarrays include positive and negative (printing buffer alone) hybridisation controls, as well as spiked controls and RT efficacy controls. Internal control spots can be useful for assessing Reproduction (2005) 130 1–13 6 C A White and L A Salamonsen data quality, but it is not recommended that they be used to standardise the data, as a housekeeping gene would be used in quantitative RT-PCR (Churchill 2002). Optimal target quality and hybridisation conditions will ensure that the maximum microarray sensitivity is achieved, allowing even very low abundance genes to be detected and differential expression determined. Image acquisition Nylon membrane cDNA microarrays hybridised with a radioactive 33P-labelled probe are scanned with a phosphorimager screen. Commercially available software such as Imagene (BioDiscovery, CA, USA) is then used to align the specific grid of arrayed DNA spots and quantify the signal intensity at each location. Scanners and software used for fluorescent image acquisition from glass cDNA microarrays are more complex, and there are many different systems available (see http://www.biocompare.com). The Axon GenePix 4000B scanner has two lasers which simultaneously excite a small region of the glass surface (, 100 mm2), at a focal plane pre-set by the user (2 50 mm to þ200 mm relative to the slide). The entire image is obtained by moving the laser lens across the glass slide. Light emitted at the wavelengths of the fluorescent labels (532 nm for Cy3 and 635 nm for Cy5) is converted to an electrical signal with a photomultiplier tube (PMT). These signals are then displayed as a 16-bit tagged image file format (TIFF) image (Cy5 coloured red and Cy3 coloured green) and given numerical values. Scanning and processing images from glass microarrays requires the investigator to perform a number of manual tasks, each of which demands a high level of technical knowledge. The PMT of the GenePix scanner needs to be adjusted to maximise its dynamic range (0 to 65 536 pixels), to prevent signal saturation and balance the intensities of the two excitation wavelengths (Forster et al. 2003). Both red and green foreground and background intensities are measured for each spot. The foreground intensity for the spot is given as the mean intensity of all foreground pixels and this is assumed to be proportional to the number of complementary mRNA molecules present in the sample. The method of background calculation differs between software packages, but local background correction is preferred over global methods (Kim et al. 2002). The default method used by GenePix software is local background subtraction, in which a different background value is computed for each individual spot and the median background pixel intensity is used for correction purposes. Following image acquisition, the user must align an appropriate grid containing spot identities to the image, as well as identify artefacts of the hybridisation process so that they can be removed from subsequent analyses. As the settings used for background calculation, background thresholds and ratio calculation can greatly influence data quality, the investigator should be aware of the implications of using each of the different methods. Reproduction (2005) 130 1–13 Importantly, the methods used for image acquisition can be optimised from slide to slide, but those used for image quantification should be identical for all slides in the experiment (Forster et al. 2003). Data analysis Following image acquisition and conversion of the image into spot intensity and intensity ratio measures, this large body of data must be stored in spreadsheet form for further analysis. Depending on the number of microarrays processed, a high capacity system such as an SQL server may be required. Data analysis is perhaps the most critical aspect of a microarray experiment and the least understood by the majority of biologists. Many different analytical approaches have been developed to achieve sensitivity in detecting gene changes while also providing a measure of statistical significance and likelihood of error. Measures of statistical significance can be made for a single microarray, or across multiple biological and technical replicate microarrays. The complexity of the data and statistical analysis requires the use of sophisticated visualisation and analysis software. Graphical displays are useful in determining the overall success of a microarray experiment (Smyth et al. 2003). The red-green image produced during scanning can detect any problems with colour balance, hybridisation, spatial effects, spot quality or artefacts such as scratches and dust. The original red-green image can also be used to check differential expression of a particular gene. Before they are plotted or analysed, the raw intensity data are always log-transformed (log2) to spread the values more evenly across the scale from 0 to 65 535 pixels. If any negative values for red (R) or green (G) foreground intensity have arisen due to high spot background, these will be removed from the analysis on a log scale. Using these log-transformed values, an informative visualisation tool is the MA-plot (Dudoit et al. 2002b). This scatterplot has M-values (R/G ratio log-transformed to M ¼ log2R/G) on the vertical axis and A-values (spot intensity expressed as p A ¼ log2 ðR £ GÞ) on the horizontal axis. Particularly when using large microarrays with thousands of probes, the majority of genes should not be differentially expressed. An MA-plot of good quality microarray data should therefore have an elongated comet shape centred around M ¼ 0 (i.e. equal red and green intensities over a wide intensity range). As well as helping to identify spot artefacts and intensity-dependent patterns, MA-plots can be used to display the effects of normalisation on the data (Smyth et al. 2003). Normalisation Normalisation is essential in microarray experiments to adjust the data for systematic non-biological effects arising from technical variation and measurement error (see ‘Experimental design’ above). The aim of normalisation is to remove the effect of this ‘noise’ from the data, while www.reproduction-online.org Guide to issues in microarray analysis still maintaining the ability to detect significantly differentially expressed genes. When using multiple nylon membrane or other single sample microarrays, each of the arrays must be paired with another and normalised or ‘scaled’ to its pair (reviewed in Evans et al. 2003). Dual colour fluorescence microarrays require normalisation to account for differences between microarrays, print-tips groups and fluorescent dye channels (reviewed in Smyth & Speed 2003). There is no universally accepted method of microarray data normalisation, and a description and comparison of all available methods is beyond the scope of this review. Overall, the literature supports the use of intensity-dependent normalisation methods, such as printtip loess (local weighted regression) normalisation (Dudoit et al. 2002b, Yang et al. 2002, Park et al. 2003, Smyth & Speed 2003). This method is capable of removing biases without altering the structure of the data. Essentially, printtip loess normalisation corrects the M-values (log2R/G ratios) for non-biological spatial and intensity effects. Statistical analysis Clustering was one of the first methods used to impose order on microarray data (Eisen et al. 1998). This method involves grouping genes on the basis of similar expression patterns, with the assumption that each cluster of genes is co-ordinately regulated, perhaps as part of the same signaling pathway. Clustering can be useful in assigning potential functions to unidentified genes and ESTs, which can then be tested in further studies. Related methods such as supervised clustering, principle component analysis, self-organising maps and linear discriminant analysis are also widely used to discover patterns of gene expression common to a particular physiological state. The aim of a microarray experiment is usually to identify differentially expressed genes, with a measure of statistical significance (reviewed in Dudoit et al. 2002b, Cui & Churchill 2003). Most microarray experiments are designed with only one categorical factor (eg. treatment or genotype), so the statistical analysis is based on the paired t-test. Experiments with multiple categorical factors (eg. genotype and time) require methods based on the analysis of variance (ANOVA). Once the data are appropriately normalised, it is common practice to consider a univariate testing problem for each gene and calculate t-statistics (Dudoit et al. 2002b). The t-statistic tests the null hypothesis of equal mean expression levels in the two samples (e.g. treatment and control). Another useful indicator of differential expression is the B-statistic (Lonnstedt & Speed 2002), which is an estimate of the odds that the gene is differentially expressed. The challenge in assigning statistical significance to a differentially expressed gene is that the often thousands of genes on a microarray result in a high level of multiple testing. Determining the false discovery rate is the most powerful method of controlling for multiple testing (Tusher et al. 2001), but this can also be achieved using adjusted P values (Dudoit et al. 2002b). www.reproduction-online.org 7 Time course experiments require even more specialised statistical analysis (Cui & Churchill 2003) and should only be conducted if the primary biological question is one of time dependence. Just as diagnostic MA-plots can be invaluable for visualising trends in raw and normalised data, plots of values obtained during statistical analysis are also useful. Both fold change (difference in gene abundance between two samples) and significance measures can be represented graphically in a ‘volcano plot’ (Cui & Churchill 2003), with the log odds of differential expression on the vertical axis and the mean M-value (log2R/G ratio) on the horizontal axis. Genes with statistically significant differential expression will appear above a horizontal threshold line and those with large fold changes (up- or downregulated) will lie to the far left or right. Differentially expressed genes identified by the B-statistic will appear in the upper left or right quadrants. There are many different software packages available for performing normalisation, statistical analysis and visualisation with single and dual sample microarrays. Some of the more widely used packages include Cyber-T (Baldi & Long 2001), SAM (Tusher et al. 2001), BRB-ArrayTools (http://linus.nci.nih.gov/BRB-ArrayTools), QVALUE (Storey & Tibshirani 2003) and Focus (Cole et al. 2003). The statistical language R (Ihaka & Gentleman 1996, http://www. r-project.org) has also been used successfully for the analysis of microarray data (Dudoit et al. 2002a) and indeed many of the other packages are based on R commands. Bioconductor (http://www.bioconductor.org) provides a more user-friendly interface for the R statistical language. Although they are often easier for biologists to use, care should be exercised in the choice of commercially available software packages. Some are excellent for data visualisation and normalisation, but cannot assign measures of statistical significance within and across multiple microarrays, or do not handle time course data. Data validation The biomedical research community does not yet accept that microarray data can stand alone, without independent validation (reviewed in Rockett & Hellmann 2004). There are a number of reasons for this caution, including the relatively recent development of the technology, the lack of standard operating procedures and the potential for errors to exist in the data (Knight 2001, Kothapalli et al. 2002). Hybridisation errors may occur due to crosshybridisation between transcripts of high homology and data may be misleading if mis-annotation of probe sequences has occurred. Even in well-maintained clone sets, it is estimated that 1 –5% of clones do not contain the correct sequence (Knight 2001). In addition, as the statistical tests used for microarray data analysis are yet to be standardised, often several methods are used and the resulting data requires further validation. Best practice in microarray analysis can achieve less than 5 –10% Reproduction (2005) 130 1–13 8 C A White and L A Salamonsen variation in signal intensity from replicate probes on the same microarray, and around 10 –30% variation between corresponding probe signal intensities on different microarrays (reviewed in Stears et al. 2003). Despite this, there is an expectation that an additional mRNA quantification method will be used to confirm the differential expression of the genes of interest (Firestein & Pisetsky 2002). The first important task for the investigator is to decide which genes to investigate further. From experience, genes displaying a large fold change (. 2) and statistical significance are the best candidates for validation. Before embarking on additional studies, it is good practice to review the primary red-green image data to confirm differential expression of these genes, and the spotted DNA sequence may also be checked for correct annotation. Comparing data with that obtained from other microarray studies on the same system can also provide ‘in silico’ validation and increase confidence in the data set as a whole (Chuaqui et al. 2002). Quantitative real-time RT-PCR is commonly used to confirm mRNA levels, as it has higher sensitivity and lower RNA requirements than Northern blot. Previous studies have demonstrated that genes with relatively high expression and at least 2-fold regulation are likely to be validated using real-time RT-PCR (Rajeevan et al. 2001). The advantage of Northern blot and RNase protection assay is that they provide a quantitative measure as well as reveal the number and size of transcripts detected by the particular spotted DNA sequence. Quantitative data obtained with microarray and Northern blot are comparable, with Northern blot slightly more sensitive in detecting differential expression compared with microarray (Taniguchi et al. 2001). In complex tissues such as the endometrium, defining the cellular localisation of mRNA expression using in situ hybridisation can provide important functional information. As it is almost impossible to differentiate between primary and secondary gene expression effects in microarray data, further testing may be required to define the molecular interactions occurring. While mRNA reflects the functional state of the cell, it is the proteins which ultimately carry out the instructions of the genome. Translation of mRNA into protein may be controlled independently of transcription and proteins may undergo post-translational modifications that alter their function. To describe a biological event or system, therefore, gene expression data obtained by microarray analysis must be extended to the study of protein products. Particularly if target RNA has been prepared from whole tissue, characterising the cellular distribution of the corresponding protein by immunostaining or tissue array is critical to understanding the function of a gene. Protein quantification by Western blot or ELISA will indicate whether transcription and translation are co-ordinately regulated. Defining the functions of differentially expressed genes may be considered the ultimate validation of microarray Reproduction (2005) 130 1–13 data. Functional studies may include in vitro experiments using dominant-negative mutants or RNA interference, or in vivo experiments using antisense morpholino oligonucleotides, knockout or conditional knockout technologies. Though the experiments may be carried out some time later, each level of data validation (mRNA, protein and function) should be considered at the microarray experimental design stage, to allow additional controlled samples to be obtained. Endometrial gene expression analysis In the last few years, both cDNA and oligonucleotide microarray technology have been successfully applied to the study of endometrial gene expression (reviewed in Giudice 2003). The endometrium is a uniquely dynamic tissue, with the capacity to undergo dramatic remodelling in response to cyclic variations in steroid hormones and local autocrine and paracrine factors. The mRNA complement of each of its different cell types is altered during different phases of the menstrual cycle, with the onset of decidualisation and in response to an implanting embryo, as well as in pathological conditions such as endometriosis, abnormal bleeding, infection or cancer. Gene expression profiling has the capacity to identify new targets for the manipulation of fertility and the diagnosis and treatment of endometrial abnormalities. A number of endometrial gene expression studies have been discussed in recent reviews (Giudice 2003, 2004, Horcajadas et al. 2004), so rather than providing a detailed description of their findings, the experimental design features of these studies have been summarised in Table 2. With only three exceptions (Popovici et al. 2000, Martin et al. 2002, Okada et al. 2003), all of these studies included validation of a small number of differentially expressed genes (usually less than 10) by an independent mRNA quantification method (Northern blot, semi-quantitative or quantitative RT-PCR). Less than half also included cellular localisation studies (in situ hybridisation and/or immunohistochemistry). The molecular events of endometrial stromal cell decidualisation are still not well defined, so microarray analysis is ideally suited to identify genes with important regulatory roles in this process. A number of factors including combined oestrogen and progesterone or cAMP can be used to induce decidualisation in vitro and the comparison of decidualised with non-decidualised cells has revealed many of the genes which are likely to contribute to this process (Popovici et al. 2000, Brar et al. 2001, Tierney et al. 2003). The two studies investigating the window of implantation (Kao et al. 2002, Borthwick et al. 2003) were identical in experimental design except that the microarrays used in the first study were hybridised with RNA samples from individual endometrial biopsies, whereas the second study used pooled samples. It is important to note that there is no consensus on the use of pooled or individual samples www.reproduction-online.org Guide to issues in microarray analysis 9 Table 2 Summary of studies using DNA microarray analysis to investigate gene expression in the endometrium during normal processes or in response to stimuli Species Cell/tissue type Process/stimulus Reference Human Isolated endometrial stromal cells Decidualisation Endometrial biopsies Progesterone Pre-receptive vs receptive Popovici et al. 2000 Brar et al. 2001 Tierney et al. 2003 Okada et al. 2003 Carson et al. 2002 Martin et al. 2002 Dominguez et al. 2003 Riesewijk et al. 2003 Horcajadas et al. 2004 Kao et al. 2002 Borthwick et al. 2003 Ponnampalam et al. 2004 Mirkin et al. 2004 Horcajadas et al. 2005 Mutter et al. 2001 Risinger et al. 2003 Moreno-Bueno et al. 2003 Cao et al. 2004 Saidi et al. 2004 Ferguson et al. 2004 Ferguson et al. 2005 Eyster et al. 2002 Lebovic et al. 2002 Kao et al. 2003 Arimoto et al. 2003 Matsuzaki et al. 2004 Yanaihara et al. 2005 Catalano et al. 2003 Punyadeera et al. 2005 Tan et al. 2003 Yoshioka et al. 2000 Reese et al. 2001 Cheon et al. 2002 Curtis-Hewitt et al. 2003 Ho Hong et al. 2004 Watanabe et al. 2003 Yao et al. 2003 Poggi et al. 2003 Naciff et al. 2002 Wu et al. 2003 Ace & Okulicz, 2004 Tynan et al. 2005 Ishiwata et al. 2003 Window of implantation vs late proliferative phase Across normal menstrual cycle Normal vs IVF cycles Endometrial cancer Endometriosis Laser capture microdissection Endometrial explants Epithelial vs stromal cells Progestin antagonist (RU486) Menstrual vs late proliferative phase; oestrogen Across normal oestrous cycle Pre-vs post-implantation Implantation vs interimplantation sites Progestin antagonist (RU486) Oestrogen Mouse Whole uterus Rat Dissected decidua Whole uterus Hoxa-10 deficiency Alcohol Ovariectomy and oestrogen treatment Rhesus monkey Cynomolgus monkey Cow Endometrial biopsies Endometrial biopsies Endometrial biopsies Pre-receptive vs receptive Progestin antagonist (RU486) Pregnant vs non-pregnant (see ‘Experimental design’ above). Many of the differentially expressed genes identified by Borthwick et al. (2003) were the same as those reported by Kao et al. (2002), suggesting that different experimental designs may be equally valid. Similarly, there was some consensus between the study by Kao et al. (2002) and that of Carson et al. (2002) examining the transition of the endometrium into receptivity. The study by Riesewijk et al. (2003) was slightly different from that of Carson et al. (2002) in that consecutive endometrial biopsies were taken from the same patient (rather than from different patients) and the timing of the biopsies was more precise. There was good agreement in the data obtained by Riesewijk et al. (2003) and Kao et al. (2002), but less between the two more similar studies (Carson et al. 2002, Riesewijk et al. 2003). Discrepancies in data sets are not unexpected, considering the www.reproduction-online.org many sources of variation in microarray experiments, and the differences in methodology used by different laboratories. Any genes which show consistent differential expression under similar experimental conditions can therefore be considered almost certain to play an important biological role. Genes involved in endometrial receptivity and implantation have also been examined using the progesterone receptor antagonist RU486 (Cheon et al. 2002, Catalano et al. 2003, Tynan et al. 2005). As RU486 is known to inhibit implantation in mice and humans, its downstream target genes are likely to be involved in normal implantation, and have been identified in the whole mouse uterus (Cheon et al. 2002), in human endometrial explants (Catalano et al. 2003) and in cynomolgus monkey endometrial biopsies (Tynan et al. 2005). While no genes were found to be regulated by RU486 in all Reproduction (2005) 130 1–13 10 C A White and L A Salamonsen three species, some were identified as downregulated with both RU486 treatment (Cheon et al. 2002) and in post-implantation compared with pre-implantation mouse uterus (Yoshioka et al. 2000). Reese et al. (2001) examined genes involved in mouse implantation using a combined approach of implantation versus interimplantation sites and activated versus delayed implantation. Interestingly, many of the genes regulated in both models were associated with the maternal immune response. Genes found to be regulated in both mouse and human endometrium with the onset of decidualisation and/or receptivity are attractive targets for the manipulation of implantation mechanisms conserved across species. Microarrays have also been utilised to identify potential markers of endometrial pathologies. Studies exploring differential gene expression in endometriotic lesions versus eutopic endometrium (Eyster et al. 2002, Lebovic et al. 2002, Arimoto et al. 2003) have revealed dysregulation of a number of genes in endometriotic tissue (reviewed in Giudice 2003), which may prove to have functional roles in this disease. Interestingly, a comparison of gene expression in eutopic endometrium from women with and without endometriosis (Kao et al. 2003) has shown that the endometrium of women with endometriosis has an altered transcriptional profile to that of women without the disease. As endometriosis is often associated with infertility, genes with altered expression in endometriosis patients may be involved in endometrial receptivity and embryo implantation (reviewed in Giudice et al. 2002). Similarly, genes found by microarray to be differentially expressed in endometrial tumours compared with normal endometrium (Mutter et al. 2001, Saidi et al. 2004) are likely to provide diagnostic markers and treatment targets in the future (reviewed in Giudice 2003). Another powerful application of microarray technology is the classification of tumour types by their gene expression profiles and a number of studies have successfully utilised this approach in endometrial cancer (Moreno-Bueno et al. 2003, Risinger et al. 2003, Cao et al. 2004, Ferguson et al. 2004, 2005). Molecular classification of tumours using microarray technology has the potential to greatly enhance patient management and improve treatment and prognosis. Conclusions Conducting a microarray experiment requires extensive planning and a high level of technical ability. Each of the steps involved, from experimental design through to data analysis, requires careful consideration and greatly benefits from collaboration with a bioinformatician. As the volume of microarray data continues to expand, there is an increasing need for accepted standard operating procedures and coordinated data deposition. The Microarray Gene Expression Data (MGED) Society exists to encourage MIAME (Minimum Information About a Reproduction (2005) 130 1–13 Microarray Experiment) compliance, so that microarray data from different laboratories and platforms can be unambiguously interpreted and independently verified. Indeed publication in a number of leading journals, including Nature, is now dependent upon meeting MIAME standards (Rockett & Hellmann 2004). Standardisation of the field could also be achieved by including a Standard Gene Set (SGS) on all microarray platforms (Fryer et al. 2002). Publicly-available databases are also required to enable large-scale data mining. For example, mRNA levels measured in many human tissues using the Affymetrix system are now available online (HuGE Index; http://www.hugeindex. org) and an endometrial database containing gene expression data from many of the studies mentioned above can be accessed at http://endometrium.bcm.tmc. edu/edr. Despite many recent advances, microarray analysis should not be considered the end-point of an investigation, but rather as a tool to assist in the formulation of hypotheses. With improved microarray quality, standardised data analysis methods and integration with proteomic approaches, gene expression profiling will be an extremely effective tool towards understanding the biology of the reproductive system and in developing diagnostic tests and therapeutic strategies for reproductive abnormalities. Acknowledgements The authors would like to thank Dr Garry Myers, Dr Gordon Smyth, Prof Terry Speed and Dr Andrew Sharkey for training and helpful discussions and Sue Panckridge for assistance in the preparation of Figure 1. LAS is supported by the National Health and Medical Research Council of Australia (grants #241000 and #143798) and CAW by an Australian Postgraduate Award. The authors declare that there is no conflict of interest that would prejudice the impartiality of this scientific work. References Ace CI & Okulicz WC 2004 Microarray profiling of progesteroneregulated endometrial genes during the rhesus monkey secretory phase. Reproductive Biology & Endocrinology 2 54. Arimoto T, Katagiri T, Oda K, Tsunoda T, Yasugi T, Osuga Y, Yoshikawa H, Nishii O, Yano T, Taketani Y & Nakamura Y 2003 Genome-wide cDNA microarray analysis of gene-expression profiles involved in ovarian endometriosis. International Journal of Oncology 22 551 –560. Baldi P & Long A 2001 A Bayesian framework for the analysis of microarray expression data: regularized t-test and statistical inferences of gene changes. Bioinformatics 17 509–519. Barrett JC & Kawasaki ES 2003 Microarrays: the use of oligonucleotides and cDNA for the analysis of gene expression. Drug Discovery Today 8 134– 141. Borthwick JM, Charnock-Jones DS, Tom BD, Hull ML, Teirney R, Phillips SC & Smith SK 2003 Determination of the transcript profile of human endometrium. Molecular Human Reproduction 9 19–33. www.reproduction-online.org Guide to issues in microarray analysis Bowtell DDL 1999 Options available - from start to finish - for obtaining expression data by microarray. Nature Genetics 21 Suppl 25– 32. Brar AK, Handwerger S, Kessler CA & Aronow BJ 2001 Gene induction and categorical reprogramming during in vitro human endometrial fibroblast decidualization. Physiological Genomics 7 135– 148. Cao QJ, Belbin T, Socci N, Balan R, Prystowsky MB, Childs G & Jones JG 2004 Distinctive gene expression profiles by cDNA microarrays in endometrioid and serous carcinomas of the endometrium. International Journal of Gynecological Pathology 23 321–329. Carson DD, Lagow E, Thathiah A, Al-Shami R, Farach-Carson MC, Vernon M, Yuan L, Fritz MA & Lessey B 2002 Changes in gene expression during the early to mid-luteal (receptive phase) transition in human endometrium detected by high-density microarray screening. Molecular Human Reproduction 8 871–879. Catalano RD, Yanaihara A, Evans AL, Rocha D, Prentice A, Saidi S, Print CG, Charnock-Jones DS, Sharkey AM & Smith SK 2003 The effect of RU486 on the gene expression profile in an endometrial explant model. Molecular Human Reproduction 9 465–473. Chen JJ, Delongchamp RR, Tsai C-A, Hsueh H-M, Sistare F, Thompson KL, Desai VG & Fuscoe JC 2004 Analysis of variance components in gene expression data. Bioinformatics 20 1436–1446. Cheon Y-P, Li Q, Xu X, DeMayo FJ, Bagchi IC & Bagchi MK 2002 A genomic approach to identify novel progesterone receptor regulated pathways in the uterus during implantation. Molecular Endocrinology 16 2853– 2871. Chuaqui RF, Bonner RF, Best CJ, Gillespie JW, Flaig MJ, Hewitt SM, Phillips JL, Krizman DB, Tangrea MA & Ahram M et al. 2002 Postanalysis follow-up and validation of microarray experiments. Nature Genetics 32 Suppl 509–514. Churchill GA 2002 Fundamentals of experimental design for cDNA microarrays. Nature Genetics 32 Suppl 490–495. Cole SW, Galic Z & Zack JA 2003 Controlling false negative errors in microarray differential expression analysis: a PRIM approach. Bioinformatics 19 1808–1816. Cui X & Churchill GA 2003 Statistical tests for differential expression in cDNA microarray experiments. Genome Biology 4 210. Curtis-Hewitt S, Deroo BJ, Hansen K, Collins J, Grissom S, Afshari CA & Korach KS 2003 Estrogen receptor-dependent genomic responses in the uterus mirror the biphasic physiological response to estrogen. Molecular Endocrinology 17 2070– 2083. Dominguez F, Avila S, Cervero A, Martin J, Pellicer A, Castrillo JL & Simon C 2003 A combined approach for gene discovery identifies insulin-like growth factor-binding protein-related protein 1 as a new gene implicated in human endometrial receptivity. Journal of Clinical Endocrinology & Metabolism 88 1849–1857. Dudoit S, Yang YH & Bolstad B 2002a Using R for the analysis of DNA microarray data. R News 2 24– 32. Dudoit S, Yang YH, Callow MJ & Speed TP 2002b Statistical methods for identifying differentially expressed genes in replicated cDNA microarray experiments. Statistica Sinica 12 111–139. Duggan DJ, Bittner M, Chen Y, Meltzer P & Trent JM 1999 Expression profiling using cDNA microarrays. Nature Genetics 21 Suppl 10– 14. Eisen MB, Spellman PT, Brown PO & Botstein D 1998 Cluster analysis and display of genome-wide expression patterns. PNAS 95 14863–14868. Emmert-Buck MR, Bonner RF, Smith PD, Chuaqui RF, Zhuang Z, Goldstein SR, Weiss RA & Liotta LA 1996 Laser capture microdissection. Science 274 998–1001. Evans AL, Sharkey AS, Saidi SA, Print CG, Catalano RD, Smith SK & Charnock-Jones DS 2003 Generation and use of a tailored gene array to investigate vascular biology. Angiogenesis 6 93–104. www.reproduction-online.org 11 Eyster KM, Boles AL, Brannian JD & Hansen KA 2002 DNA microarray analysis of gene expression markers of endometriosis. Fertility & Sterility 77 38–42. Ferguson SE, Olshen AB, Viale A, Awtrey CS, Barakat RR & Boyd J 2004 Gene expression profiling of tamoxifen-associated uterine cancers: evidence for two molecular classes of endometrial carcinoma. Gynecological Oncology 92 719– 725. Ferguson SE, Olshen AB, Viale A, Barakat RR & Boyd J 2005 Stratification of intermediate-risk endometrial cancer patients into groups at high risk or low risk for recurrence based on tumor gene expression profiles. Clinical Cancer Research 11 2252–2257. Firestein GS & Pisetsky DS 2002 DNA microarrays: boundless technology or bound by technology? Guidelines for studies using microarray technology. Arthritis & Rheumatism 46 859–861. Forster T, Roy D & Ghazal P 2003 Experiments using microarray technology: limitations and standard operating procedures. Journal of Endocrinology 178 195 –204. Fryer RM, Randall J, Yoshida T, Hsaio L-L, Blumenstock J, Jensen KE, Dimofte T, Jensen RV & Gullans SR 2002 Global analysis of gene expression: methods, interpretation, and pitfalls. Experimental Nephrology 10 64–74. Giudice LC, Telles TL, Lobo S & Kao LC 2002 The molecular basis for implantation failure in endometriosis: on the road to discovery. Annals of the New York Academy of Sciences 955 252–264. Giudice LC 2003 Elucidating endometrial function in the post-genomic era. Human Reproduction Update 9 223 –235. Giudice LC 2004 Microarray expression profiling reveals candidate genes for human uterine receptivity. American Journal of Pharmacogenomics 4 299–312. Ho Hong S, Young Nah H, Yoon Lee J, Chan Gye M, Hoon Kim C & Kyoo Kim M 2004 Analysis of estrogen-regulated genes in mouse uterus using cDNA microarray and laser capture microdissection. Journal of Endocrinology 181 157–167. Horcajadas JA, Riesewijk A, Martin J, Cervero A, Mosselman S, Pellicer A & Simon C 2004 Global gene expression profiling of human endometrial receptivity. Journal of Reproductive Immunology 63 41–49. Horcajadas JA, Riesewijk A, Polman J, van Os R, Pellicer A, Mosselman S & Simon C 2005 Effect of controlled ovarian hyperstimulation in IVF on endometrial gene expression profiles. Molecular Human Reproduction 11 195– 205. Ihaka R & Gentleman R 1996 R: a language for data analysis and graphics. Journal of Computational & Graphical Statistics 5 299–314. Ishiwata H, Katsuma S, Kizaki K, Patel OV, Nakano H, Takahashi T, Imai K, Hirasawa A, Shiojima S, Ikawa H, Suzuki Y, Tsujimoto G, Izaike Y, Todoroki J & Hashizume K 2003 Characterization of gene expression profiles in early bovine pregnancy using a custom cDNA microarray. Molecular Reproduction & Development 65 9 –18. Kao LC, Germeyer A, Tulac S, Lobo S, Yang JP, Taylor RN, Osteen K, Lessey BA & Giudice LC 2003 Expression profiling of endometrium from women with endometriosis reveals candidate genes for disease-based implantation failure and infertility. Endocrinology 144 2870– 2881. Kao LC, Tulac S, Lobo S, Imani B, Yang JP, Germeyer A, Osteen K, Taylor RN, Lessey BA & Giudice LC 2002 Global gene profiling in human endometrium during the window of implantation. Endocrinology 143 2119–2138. Kendziorski CM, Zhang Y, Lan H & Attie AD 2003 The efficiency of pooling mRNA in microarray experiments. Biostatistics 4 465–477. Kerr MK & Churchill GA 2001 Statistical design and the analysis of gene expression microarray data. Genetics Research 77 123–128. Kim JH, Shin DM & Lee YS 2002 Effect of local background intensities in the normalization of cDNA microarray data with a skewed expression profiles. Experimental Molecular Medicine 34 224–232. Reproduction (2005) 130 1–13 12 C A White and L A Salamonsen Knight J 2001 When the chips are down. Nature 410 860–861. Kothapalli R, Yoder SJ, Mane S & Loughran TP Jr 2002 Microarray results: how accurate are they? BMC Bioinformatics 3 22. Lebovic DI, Baldocchi RA, Mueller MD & Taylor RN 2002 Altered expression of a cell-cycle suppressor gene, Tob-1, in endometriotic cells by cDNA array analyses. Fertility & Sterility 78 849–854. Lipschutz RJ, Fodor SPA, Gingeras TR & Lockhart DJ 1999 High density synthetic oligonucleotide arrays. Nature Genetics 21 Suppl 20–24. Lonnstedt I & Speed TP 2002 Replicated microarray data. Statistica Sinica 12 31–46. Matsuzaki S, Canis M, Vaurs-Barriere C, Pouly JL, Boespflug-Tanguy O, Penault-Llorca F, Dechelotte P, Dastugue B, Okamura K & Mage G 2004 DNA microarray analysis of gene expression profiles in deep endometriosis using laser capture microdissection. Molecular Human Reproduction 10 719–728. Martin J, Dominguez F, Avila S, Castrillo JL, Remohi J, Pellicer A & Simon C 2002 Human endometrial receptivity: gene regulation. Journal of Reproductive Immunology 55 131 –139. Mirkin S, Nikas G, Hsiu J-G, Diaz J & Oehninger S 2004 Gene expression profiles and structural/functional features of the peri-implantation endometrium in natural and gonadotropin-stimulated cycles. Journal of Clinical Endocrinology & Metabolism 89 5742–5752. Moreno-Bueno G, Sanchez-Estevez C, Cassia R, Rodriguez-Perales S, Diaz-Uriarte R, Dominguez O, Hardisson D, Andujar M, Prat J, Matias-Guiu X, Cigudosa JC & Palacios J 2003 Differential gene expression profile in endometrioid and nonendometrioid endometrial carcinoma: STK15 is frequently overexpressed and amplified in nonendometrioid carcinomas. Cancer Research 63 5697–5702. Mutter GL, Baak JPA, Fitzgerald JT, Gray R, Neuberg D, Kust GA, Gentleman R, Gullans SR, Wei LJ & Wilcox G 2001 Global expression changes of constitutive and hormonally regulated genes during endometrial neoplastic transformation. Gynecological Oncology 83 177 –185. Naciff JM, Jump ML, Torontali SM, Carr GJ, Tiesman JP, Overmann GJ & Daston GP 2002 Gene expression profile induced by 17alpha-ethynyl estradiol, bisphenol A, and genistein in the developing female reproductive system of the rat. Toxicological Sciences 68 184–199. Okada H, Nakajima T, Yoshimura T, Yasuda K & Kanzaki H 2003 Microarray analysis of genes controlled by progesterone in human endometrial stromal cells in vitro. Gynecological Endocrinology 17 271–280. Park T, Yi S-G, Kang S-H, Lee S-Y, Lee Y-S & Simon R 2003 Evaluation of normalization methods for microarray data. BMC Bioinformatics 4 33. Peng X, Wood CL, Blalock EM, Chen KC, Landfield PW & Stromberg AJ 2003 Statistical implications of pooling RNA samples for microarray experiments. BMC Bioinformatics 4 26. Petalidis L, Bhattacharyya S, Morris GA, Collins VP, Freeman TC & Lyons PA 2003 Global amplification of mRNA by template-switching PCR: linearity and application to microarray analysis. Nucleic Acids Research 31 e142. Poggi SH, Goodwin KM, Hill JM, Brenneman DE, Tendi E, Schninelli S & Spong CY 2003 Differential expression of c-fos in a mouse model of fetal alcohol syndrome. American Journal of Obstetrics & Gynecology 189 786–789. Ponnampalam AP, Weston GC, Trajstman AC, Susil B & Rogers PA 2004 Molecular classification of human endometrial cycle stages by transcriptional profiling. Molecular Human Reproduction 10 879–893. Popovici RM, Kao LC & Giudice LC 2000 Discovery of new inducible genes in in vitro decidualized human endometrial stromal cells using microarray technology. Endocrinology 141 3510–3513. Punyadeera C, Dassen H, Klomp J, Dunselman G, Kamps R, Dijcks F, Ederveen A, de Goeij A & Groothuis P 2005 Reproduction (2005) 130 1–13 Oestrogen-modulated gene expression in the human endometrium. Cellular & Molecular Life Sciences 62 239–250. Rajeevan MS, Vernon SD, Taysavang N & Unger ER 2001 Validation of array-based gene expression profiles by real-time (kinetic) RTPCR. Journal of Molecular Diagnostics 3 26–31. Reese J, Das SK, Paria BC, Lim H, Song H, Matsumoto H, Knudtson KL, DuBois RN & Dey SK 2001 Global gene expression analysis to identify molecular markers of uterine receptivity and embryo implantation. Journal of Biological Chemistry 276 44137–44145. Riesewijk A, Martin J, van Os R, Horcajadas JA, Polman J, Pellicer A, Mosselman S & Simon C 2003 Gene expression profiling of human endometrial receptivity on days LH þ 2 versus LH þ 7 by microarray technology. Molecular Human Reproduction 9 253 –264. Risinger JI, Maxwell GL, Chandramouli GV, Jazaeri A, Aprelikova O, Patterson T, Berchuck A & Barrett JC 2003 Microarray analysis reveals distinct gene expression profiles among different histologic types of endometrial cancer. Cancer Research 63 6–11. Rockett JC & Hellmann GM 2004 Confirming microarray data - is it really necessary? Genomics 83 541– 549. Saidi SA, Holland CM, Kreil DP, MacKay DJ, Charnock-Jones DS, Print CG & Smith SK 2004 Independent component analysis of microarray data in the study of endometrial cancer. Oncogene 23 6677–6683. Sasik R, Woelk CH & Corbeil J 2004 Microarray truths and consequences. Journal of Molecular Endocrinology 33 1–9. Schena M, Shalon D, Davis RW & Brown PO 1995 Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 270 467– 470. Smyth GK & Speed T 2003 Normalization of cDNA microarray data. Methods 31 265 –273. Smyth GK, Yang YH & Speed T 2003 Statistical issues in cDNA microarray data analysis. Methods in Molecular Biology 224 111– 136. Stears RL, Martinsky T & Schena M 2003 Trends in microarray analysis. Nature Medicine 9 140 –145. Sterrenburg E, Turk R, Boer JM, van Ommen GB & den Dunnen JT 2002 A common reference for cDNA microarray hybridizations. Nucleic Acids Research 30 e116. Storey JD & Tibshirani R 2003 Statistical significance for genomewide studies. PNAS 100 9440–9445. Tan YF, Li FX, Piao YS, Sun XY & Wang YL 2003 Global gene profiling analysis of mouse uterus during the oestrous cycle. Reproduction 126 171–182. Taniguchi M, Miura K, Iwao H & Yamanaka S 2001 Quantitative assessment of DNA microarrays - comparison with Northern blot analyses. Genomics 71 34–39. Tierney EP, Tulac S, Huang S-TJ & Giudice LC 2003 Activation of the protein kinase A pathway in human endometrial stromal cells reveals sequential categorical gene regulation. Physiological Genomics 16 47–66. Tusher VG, Tibshirani R & Chu G 2001 Significance analysis of microarrays applied to the ionizing radiation response. PNAS 98 5116–5121. Tynan S, Pacia E, Haynes-Johnson D, Lawrence D, D’Andrea MR, Guo JZ, Lundeen S & Allan G 2005 The putative tumor suppressor deleted in malignant brain tumors 1 (DMBT1) is an estrogen-regulated gene in rodent and primate endometrial epithelium. Endocrinology 146 1066–1073. van Gelder RN, von Zastrow ME, Yool A, Dement WC, Barchas JD & Eberwine JH 1990 Amplified RNA synthesized from limited quantities of heterogeneous cDNA. PNAS 87 1663–1667. Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA & Holt RA et al. 2001 The sequence of the human genome. Science 291 1304– 1351. Watanabe H, Suzuki A, Kobayashi M, Takahashi E, Itamoto M, Lubahn DB, Handa H & Iguchi T 2003 Analysis of temporal changes in the expression of estrogen-regulated genes in the uterus. Journal of Molecular Endocrinology 30 347–358. www.reproduction-online.org Guide to issues in microarray analysis Wu X, Pang ST, Sahlin L, Blanck A, Norstedt G & Flores-Morales A 2003 Gene expression profiling of the effects of castration and estrogen treatment in the rat uterus. Biology of Reproduction 69 1308–1317. Yanaihara A, Otsuka Y, Iwasaki S, Aida T, Tachikawa T, Irie T & Okai T 2005 Differences in gene expression in the proliferative human endometrium. Fertility & Sterility 83 1206–1215. Yang MCK, Yang JJ, McIndoe RA & She JX 2003 Microarray experimental design: power and sample size considerations. Physiological Genomics 16 24–28. Yang YH & Speed T 2002 Design issues for cDNA microarray experiments. Nature Reviews Genetics 3 579–588. Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J & Speed TP 2002 Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Research 30 e15. www.reproduction-online.org 13 Yao MW, Lim H, Schust DJ, Choe SE, Farago A, Ding Y, Michaud S, Church GM & Maas RL 2003 Gene expression profiling reveals progesterone-mediated cell cycle and immunoregulatory roles of Hoxa-10 in the preimplantation uterus. Molecular Endocrinology 17 610 –627. Yoshioka K, Matsuda F, Takakura K, Noda Y, Imakawa K & Sakai S 2000 Determination of genes involved in the process of implantation: application of GeneChip to scan 6500 genes. Biochemical & Biophysical Research Communications 272 531 –538. Received 10 February 2005 First decision 12 April 2005 Revised manuscript received 27 April 2005 Accepted 3 May 2005 Reproduction (2005) 130 1–13
© Copyright 2024