Skip to main content

ArrayD: A general purpose software for Microarray design

Abstract

Background

Microarray is a high-throughput technology to study expression of thousands of genes in parallel. A critical aspect of microarray production is the design aimed at space optimization while maximizing the number of gene probes and their replicates to be spotted.

Results

We have developed a software called 'ArrayD' that offers various alternative design solutions for an array given a set of user requirements. The user feeds the following inputs: type of source plates to be used, number of gene probes to be printed, number of replicates and number of pins to be used for printing. The solutions are stored in a text file. The choice of a design solution to be used will be governed by the spotting chemistry to be used and the accuracy of the robot.

Conclusions

ArrayD is a software for standard cartesian robots. The software aids users in preparing a judicious and elegant design. ArrayD is universally applicable and is available at http://www.igib.res.in/scientists/arrayd/arrayd.html.

Background

Microarray is a popularly used high-throughput technology to investigate gene expression of thousands of genes simultaneously at the level of mRNA. Ever since the development of this technology [1–3], transcriptional profiling at the genomic level has been employed to address numerous issues in biology and in medicine [4–8]. It is likely that microarrays will continue to be used to explore various biological phenomena. The basic underlying principle involves spotting DNA fragments either derived from polymerase chain reaction or preparation of plasmids or oligonucleotides at high density (~10,000–25,000 spots on a glass slide of 25 mm × 75 mm) representing the probes of the genes under study. The surface on which the DNA fragments or oligonucleotides are spotted is usually glass slides coated with poly-L-lysine or amino alkyl silane that serve to improve adherence of DNA to the surface. Uniform spotting at high density requires robotic operation and a variety of robots are now available for spotting [9]. The robots employed for the preparation of microarrays are of the cartesian type with movement in x-y-z direction.

A critical aspect of microarray production is the design considering space optimization to produce high-density arrays for a given set of samples and replicates. The softwares generally supplied with robotic spotters translate user input parameters into a set of instructions in robot language for printing arrays. These softwares do not offer design capabilities in which spotting parameters and grid configurations can be chosen for a given set of samples and replicates. Presently various solutions have to be derived manually in most academic laboratories. We have developed a user-friendly software 'ArrayD' that can be used by experts and novice alike to fill this gap to simplify and aid in rapid design. ArrayD offers a variety of design solutions given a set of requirements: Number of gene probes, number of replicates, and the source plate (384 well or 96 well). Because the algorithm implemented in ArrayD is inherently simple and uses fundamental principles of robot operation, the design solutions offered by ArrayD are universally applicable to any system. The choice of a design solution would be governed by the spotting chemistry and the humidity used in addition to elegant appearance. The hallmark of ArrayD is its overall simplicity and the variety of alternative designs it offers for users to decide on choosing the appropriate spotting parameters. The multiple design solutions offered by ArrayD provides a wide range of arrays from compact to loosely spaced spots as well as convenient grid patterning, which can be user selected for printing.

Implementation

ArrayD program is developed in C and can be compiled and operated on UNIX V5.1, IRIX 5.1 and Red Hat Linux 7.0 (or higher) operating systems. A companion computer program ArraySolution was developed in Perl (Practical Extraction Report Language) version 5.6.1 and can be implemented on any UNIX or Linux operating system.

Inputs to be defined for ArrayD

  1. (a)

    Type of source plate to be used

The input parameter toggles between either 96 well or 384 well plate.

  1. (b)

    Number of gene probes to be printed

The number of gene probes including positive and negative controls and blanks.

  1. (c)

    Number of replicates

The number of replicates of each gene probe. Although the most common pattern chosen is duplicates, users can choose any number of replicates. Replicates are assumed to be printed in the Y axis.

  1. (d)

    Number of pins to be used for printing

This parameter relates to time taken for printing the slides and the number of spots arrayed per slide. The number of pins in X-axis and Y-axis need to be specified. The type of pins used is assumed to be stealth pins, which are widely used. It is not necessary to specify pin type for ArrayD. Instead, this aspect is considered in the printing software according to the pin type used for implementing a particular design.

Results and discussion

A general microarray design layout is displayed in Figure 1. ArrayD accepts standard slide dimensions (25 mm × 75 mm), conceptualizes the spotting area to be 50 mm × 22 mm to provide space for barcode labeling and for appropriate placement of coverslip over the print area.

Figure 1
figure 1

A general layout of microarray. In this example, a 2 × 2 pin configuration for 192 samples and 384 spots (duplicate) were considered. The origin is marked. Each pin prints one grid. One Super grid comprises of 4 grids, two in each axis. The number of super grids in X-direction is 2 and the number of super grids in Y-direction is 1. The spot pattern of a portion of one grid is zoomed for clarity. Replicates are spotted adjacent to each other in Y-direction. The number of samples in each grid is 24 and the total number of spots is 48. Two spots are further zoomed to show the diameter of the spot and the inter-spot distance.

The reference direction of the robot for picking probes from source plate is left → right followed by top → down; the printing direction is top → down followed by left → right. Replicates are considered to be spotted in Y-axis (Figure 1). After the user has entered the parameters, the software generates a text file called 'solution.txt' that carries possible alternative array design parameters for the given input. The algorithm implemented in ArrayD is displayed in Figure 2.

Figure 2
figure 2

The flowchart of algorithm implemented in ArrayD to compute all possible design solutions for a given input parameters. The program first validates the given input parameters and then calculates the grid configuration. Note that in the case of 384 well plate type the pin-to-pin distance is 4500 microns on the print-head and therefore the max-grid-size is set at 4000 microns, which is 500 microns less than the upper limit (4500 microns). Similarly for 96 well plate type, the max-grid-size is set at 8500 microns (500 microns less than the upper limit 9000 microns). Maximum number of pins in the print head is taken as 48 and conforms to most printing robots. The spot distance database has inter-spot distances of 300, 250, 220, 200, 190, 180 and 170 microns. Users can expand this database. The variables used in the flowchart: pin_x: pin number in X-direction; pin_y: pin number in Y-direction; R: number of replicates; D: pin-to-pin distance; d: inter-spot distance; SGx (max): maximum possible number of super grids in X-direction; SGy (max): maximum possible number of super grids in Y-direction; GT: Total number of grids in the slide; Gx: Total Number of Grids in X-direction; Gy: Total Number of grids in Y-direction; ST: Maximum possible number of samples in each grid; Sx: Maximum possible number of samples in each grid in X-direction; Sy: Maximum possible number of samples in each grid in Y-direction; Xs: number of samples in each grid in X-direction for a solution; Ys: number of samples in each grid in Y-direction for a solution; Sg: Total number of samples per grid for a solution; A: Total area covered by the microarray on the slide.

The program first validates the input given by the user for appropriate number of pins in each direction and the plate type to be used. For a valid input, ArrayD calculates maximum possible number of super grids in X (or Y) direction based on the coverslip dimensions, pin number in X (or Y) direction and pin-to-pin distance (Figure 2). The coverslip dimensions have been set in the program as 50 mm × 22 mm for the longest size coverslip that can be effectively used during hybridization. The pin-to-pin distance is fixed at 4500 microns in the print head for 384 well plates and at 9000 microns for 96 well plates.

ArrayD uses a predefined inter-spot distance database. Design solutions of ArrayD encompass various inter-spot distances that would be compatible with different spotting chemistries and conditions of humidity. We have used inter-spot distances of 170 μm, 180 μm, 190 μm, 200 μm, 220 μm, 250 μm and 300 μm. This database can be expanded to incorporate even lower inter-spot distances for use with other spotting chemistries by simple modification. We chose 170 μm as least distance based on several trial experiments in our laboratory using 50% DMSO as spotting solution and SMP3 pin type. In our experience, a minimum inter-spot distance of 200 microns works best with 50% DMSO at 40% – 50% humidity at 25°C.

The options for the inter-spot distance currently offered by the software can work successfully for SMP2, SMP2B, SMP2XB, SMP2.5, SMP2.5B, SMP2.5XB, SMP3, SMP3B, SMP3XB, SMP4, SMP4B, SMP4XB, SMP5, SMP5B and SMP5XB stealth pins (See Table 1). For each possible super grid configuration, the number of grids in each direction is optimized based on the number of gene probes (samples) input by the user as shown in figure 2.

Table 1 Stealth pins for implementing various design solutions offered by ArrayD

Design solutions offered by the program

Alternative array designs for a given set of input parameters are ranked on the basis of 'Distance area ratio' that describes the area covered by the array for each design. The array design spanning least area is ranked highest. This strategy allows for applying the labeled target sparingly. Subsequently, an easy report in tabular form can be generated by feeding the output data file from ArrayD into the companion Perl program 'ArraySolution.pl', which classifies array solutions into 'Square', 'Rectangle (Horizontal bar)', or 90° rotated 'Rectangle (Vertical column)' based on the geometry of a given design solution. If the number of grids are equal in both the direction we have a 'Square' design. In all other cases we obtain a 'Rectangle' design, which can be either of two types: the long side of the array is parallel to the length (Horizontal) or the width (Vertical) of the slide. The output of ArraySolution is a tab-delimited text file called 'filename.solution' where filename corresponds to the input name of the file carrying design solutions. The tabular report consist of Number of super grids in X – direction, Number of super grids in Y – direction, Number of spots per grid in X – direction, Number of spots per grid in Y – direction, Distance between two spots (in microns), Distance Area ratio and geometry of design (Square or Rectangle). This can aid users to decide on a particular design solution based on space optimization and elegant appearance.

An example of a sample run is provided in Figure 3. The number of gene probes (including controls and blanks) to be spotted using a 2 by 2 pin configuration in X-Y axis is fed as 2304 (Figure 3). The gene probes have to be spotted in duplicates so the total number of spots on the slide would be 4608. The program provides 67 different array designs for various inter-spot distances, number of grids and number of spots in each grid. Two examples of different solutions are presented in Figure 4. In this example, the first solution is ranked highest with inter-spot distance of 170 microns and a 24 × 24 grid pattern with 4 grids in Y axis and 2 grids in X axis. An alternative solution provides a design with a higher inter-spot distance of 200 microns and 18 × 16 grid pattern with 4 grids each in X axis and Y axis. The first solution can be used in conditions when humidity is low and the spotting solution does not absorb moisture and spread after printing. The second solution is more appropriate for printing samples in 50% DMSO. The classification of all design solutions based on the geometries obtained from ArraySolution is displayed in Table 2 [see Additional file 1].

Figure 3
figure 3

Screenshot of input parameters to ArrayD. Five inputs (Plate type (384 or 96), total number of samples, total number of replicates, pin number in X-axis and pin number in Y-axis) are fed to the program sequentially.

Figure 4
figure 4

A sample run of ArrayD. The inputs given to the program were: type of plate-384; number of samples-2304; number of replicates-2; number of pins in X-axis-2; number of pins in Y-axis-2. A total of 67 solutions were offered. Two solutions are shown: top ranking solution with inter-spot distance 170 μm and an alternative design solution with inter-spot distance 200 μm. The ArrayD output was subsequently fed to the companion program 'ArraySolution' to classify each design based on its geometry and the report generated is displayed in Table 2 [see Additional file 1].

Conclusion

We have developed a simple and rapid software ArrayD that offers various design solutions of designing microarrays for a specific set of user-defined requirements.

Availability and requirements

The source code and the executable file for ArrayD and ArraySolution programs are freely available and can be downloaded from our website [11]. The source code can be compiled and executed on Unix v 5.1, or IRIX v 5.1 or Red Hat Linux v 7.0 (or higher). The executable files can be downloaded for Windows platform (Windows 98/NT/XP/2000). Further information can be requested by sending e-mail to ramu@igib.res.in or ramucbt@yahoo.com.

References

  1. Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 1995, 270: 467–70.

    Article  CAS  PubMed  Google Scholar 

  2. Schena M, Shalon D, Heller R, Chai A, Brown PO, Davis RW: Parallel human genome analysis: microarray-based expression monitoring of 1000 genes. Proc Natl Acad Sci U S A 1996, 93: 10614–9. 10.1073/pnas.93.20.10614

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Schena M: Genome analysis with gene expression microarrays. Bioessays 1996, 18: 427–31.

    Article  CAS  PubMed  Google Scholar 

  4. Gerhold DL, Jensen RV, Gullans SR: Better therapeutics through microarrays. Nature Genetics 2002, 32(suppl):547–552. 10.1038/ng1042

    Article  CAS  PubMed  Google Scholar 

  5. Petricoin III EF, Hackett JL, Lesko LJ, Puri RK, Gutman SI, Chumakov K, Woodcock J, Feigal DW, Zoon KC Jr, Sistare FD: Medical applications of microarray technologies: a regulatory science perspective. Nature Genetics 2002, 32(suppl):474–479. 10.1038/ng1029

    Article  Google Scholar 

  6. Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, Bloomfield CD, Lander ES: Molecular classification of cancer class discovery and class prediction by gene expression monitoring. Science 1999, 286: 531–7. 10.1126/science.286.5439.531

    Article  CAS  PubMed  Google Scholar 

  7. Cummings CA, Relman DA.: Using DNA microarrays to study host-microbe interactions. Emerg Infect Dis 2000, 6: 513–25.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Quackenbush J: Computational analysis of microarray data. Nat Rev Genet 2001, 2: 418–27. 10.1038/35076576

    Article  CAS  PubMed  Google Scholar 

  9. Schena M: Microarray analysis. Hoboken, New Jersey: Wiley and Sons, Inc; 2002.

    Google Scholar 

  10. Arrayit Web site[http://arrayit.com/Products/Printing/Stealth/stealth.html]

  11. ArrayD Web site[http://www.igib.res.in/scientists/arrayd/arrayd.html]

Download references

Acknowledgements

AS, GPS and VKS thank Council of Scientific and Industrial Research (CSIR) for providing financial assistance. We thank Dipayan Dasgupta for his help in compilation of the software and Mamta Khandelwal for her help in uploading the software on the web site. We thank the assistance of Technosol.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Srinivasan Ramachandran.

Additional information

Authors' contributions

AS explained the operational details to software expert, did the experimentation, maintenance and testing of the software. GPS did the basic software writing and implementation of 'ArrayD'. VKS prepared the code for the companion program 'ArraySolution' to classify the solutions on the basis of design geometry. SR is the Group Leader generating demand, sourcing and linking people, explaining the concepts, testing, critical examination, presentation of data, providing salary through grants.

Electronic supplementary material

12859_2004_258_MOESM1_ESM.pdf

Additional File 1: Report generated by ArraySolution. A total of 67 different design solutions were classified into three categories of Square, Rectangle (Horizontal) and Rectangle (Vertical) as mentioned in text. The detailed design parameters including number of supergrids in X and Y direction, number of spots per grid in X and Y direction, distance between the spots, distance area ratio and geometry of design are provided for preparing elegant microarrays. (PDF 5 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Sharma, A., Srivastava, G.P., Sharma, V.K. et al. ArrayD: A general purpose software for Microarray design. BMC Bioinformatics 5, 142 (2004). https://doi.org/10.1186/1471-2105-5-142

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-5-142

Keywords