The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Discovery of the genes in response to white spot syndrome virus (WSSV) infection in Fenneropenaeus chinensis through cDNA microarray

We used microarray technology to study differentially expressed genes in white spot syndrome virus (WSSV)-infected shrimp. A total of 3136 cDNA targets, including 1578 unique genes from a cephalothorax cDNA library and 1536 cDNA clones from reverse and forward suppression subtractive hybridization (SSH) libraries of Fenneropenaeus chinensis, plus 14 negative and 8 blank control clones, were spotted onto a 18 x 18 mm area of NH2-modified glass slides. Gene expression patterns in the cephalothorax of shrimp at 6 h after WSSV injection and moribund shrimp naturally infected by WSSV were analyzed. A total of 105 elements on the arrays showed a similar regulation pattern in artificially infected shrimp and naturally infected moribund shrimp; parts of the results were confirmed by semiquantitative reverse transcriptase-polymerase chain reaction (RT-PCR). The up-regulated expression of immune-related genes, including heat shock proteins (HSP70 and HSP90), trehalose-phosphate synthase (TPS), ubiquitin C, and so forth, were observed when shrimp were challenged with WSSV. Genes including myosin LC2, ATP synthase A chain, and arginine kinase were found to be down-regulated after WSSV infection. The expression of housekeeping genes such as actin, elongation factor, and tubulin is not stable, and so these genes are not suitable as internal standards for semiquantitative RT-PCR when shrimp are challenged by WSSV. As a substitute, we found that triosephosphate isomerase (TPI) was an ideal candidate of interstandards in this situation.

Identification of 5-Fluorouracil-inducible Target Genes Using cDNA Microarray Profiling

The fluoropyrimidine 5-Fluorouracil (5-FU) is widely used in the treatment of cancer. To identify novel downstream mediators of tumor cell response to 5-FU, we used DNA microarray technology to identify genes that are transcriptionally activated by 5-FU treatment in the MCF-7 breast cancer cell line. Of 2400 genes analyzed, 619 were up-regulated by >3-fold. Highly up-regulated genes (>6-fold) with signal intensities of >3000 were analyzed by Northern blot. Genes that were consistently found to be up-regulated were spermine/spermidine acetyl transferase (SSAT), annexin II, thymosin-beta-10, chaperonin-10, and MAT-8. Treatment of MCF-7 cells with the antifolate tomudex and DNA-damaging agent oxaliplatin also resulted in up-regulation of each of these targets. The 5-FU-induced activation of MAT-8, thymosin-beta-10, and chaperonin-10 was abrogated by inactivation of p53 in MCF-7 cells, whereas induction of SSAT and annexin II was significantly reduced in the absence of p53. Moreover, each of these genes contained more than one potential p53-binding site, suggesting that p53 may play an important regulatory role in 5-FU-induced expression of these genes. In addition, we found that basal expression levels of SSAT, annexin II, thymosin beta-10, and chaperonin-10 were increased (by approximately 2-3-fold), and MAT-8 expression dramatically increased (by approximately 10-fold) in a 5-FU-resistant colorectal cancer cell line (H630-R10) compared with the parental H630 cell line, suggesting these genes may be useful biomarkers of resistance. These results demonstrate the potential of DNA microarrays to identify novel genes involved in mediating the response of tumor cells to chemotherapy.

Identification of eukaryotic mRNAs that are translated at reduced cap binding complex eIF4F concentrations using a cDNA microarray

Although most eukaryotic mRNAs need a functional cap binding complex eIF4F for efficient 5′ end- dependent scanning to initiate translation, picornaviral, hepatitis C viral, and a few cellular RNAs have been shown to be translated by internal ribosome entry, a mechanism that can operate in the presence of low levels of functional eIF4F. To identify cellular mRNAs that can be translated when eIF4F is depleted or in low abundance and that, therefore, may contain internal ribosome entry sites, mRNAs that remained associated with polysomes were isolated from human cells after infection with poliovirus and were identified by using a cDNA microarray. Approximately 200 of the 7000 mRNAs analyzed remained associated with polysomes under these conditions. Among the gene products encoded by these polysome-associated mRNAs were immediate-early transcription factors, kinases, and phosphatases of the mitogen-activated protein kinase pathways and several protooncogenes, including c-myc and Pim-1. In addition, the mRNA encoding Cyr61, a secreted factor that can promote angiogenesis and tumor growth, was selectively mobilized into polysomes when eIF4F concentrations were reduced, although its overall abundance changed only slightly. Subsequent tests confirmed the presence of internal ribosome entry sites in the 5′ noncoding regions of both Cyr61 and Pim-1 mRNAs. Overall, this study suggests that diverse mRNAs whose gene products have been implicated in a variety of stress responses, including inflammation, angiogenesis, and the response to serum, can use translational initiation mechanisms that require little or no intact cap binding protein complex eIF4F.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Possible involvement of plant ABC transporters in cadmium detoxification: a cDNA sub-microarray approach

As a nontolerant plant to a large number of toxic compounds, Arabidopsis thaliana is a suitable model to study regulation of genes involved in response to heavy metals. Using a cDNA-microarray approach, we identified some ABC transporters that are differentially regulated after cadmium treatments, making them putative candidates for being involved in Cd sequestration and redistribution in plants. Regarding yeast and fission yeast, in which Cd is able to form complexes either with glutathione (GSH) or phytochelatins (PC) subsequently transported into vacuoles via ABC transporters, it is also very likely that some plant ABC transporters are able to transport GS2–Cd or PC–Cd complexes into subcellular compartments or outside of the cell. The characterization of such transporters is of great interest for developing molecular biology approaches in phytoremediation.

Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large

Conifers are resistant to attack from a large number of potential herbivores or pathogens. Previous molecular and biochemical characterization of selected conifer defence systems support a model of multigenic, constitutive and induced defences that act on invading insects via physical, chemical, biochemical or ecological (multitrophic) mechanisms. However, the genomic foundation of the complex defence and resistance mechanisms of conifers is largely unknown. As part of a genomics strategy to characterize inducible defences and possible resistance mechanisms of conifers against insect herbivory, we developed a cDNA microarray building upon a new spruce (Picea spp.) expressed sequence tag resource. This first-generation spruce cDNA microarray contains 9720 cDNA elements representing c. 5500 unique genes. We used this array to monitor gene expression in Sitka spruce (Picea sitchensis) bark in response to herbivory by white pine weevils (Pissodes strobi, Curculionidae) or wounding, and in young shoot tips in response to western spruce budworm (Choristoneura occidentalis, Lepidopterae) feeding. Weevils are stem-boring insects that feed on phloem, while budworms are foliage feeding larvae that consume needles and young shoot tips. Both insect species and wounding treatment caused substantial changes of the host plant transcriptome detected in each case by differential gene expression of several thousand array elements at 1 or 2 d after the onset of treatment. Overall, there was considerable overlap among differentially expressed gene sets from these three stress treatments. Functional classification of the induced transcripts revealed genes with roles in general plant defence, octadecanoid and ethylene signalling, transport, secondary metabolism, and transcriptional regulation. Several genes involved in primary metabolic processes such as photosynthesis were down-regulated upon insect feeding or wounding, fitting with the concept of dynamic resource allocation in plant defence. Refined expression analysis using gene-specific primers and real-time PCR for selected transcripts was in agreement with microarray results for most genes tested. This study provides the first large-scale survey of insect-induced defence transcripts in a gymnosperm and provides a platform for functional investigation of plant-insect interactions in spruce. Induction of spruce genes of octadecanoid and ethylene signalling, terpenoid biosynthesis, and phenolic secondary metabolism are discussed in more detail.

Identification of novel genes expressed during rhabdomyosarcoma differentiation using cDNA microarrays

Rhabdomyosarcomas (RMS) are highly aggressive tumors that are thought to arise as a consequence of the regulatory disruption of the growth and differentiation of skeletal muscle progenitor cells. Normal myogenesis is characterized by the expression of the myogenic regulatory factor gene family but, despite their expression in RMS, these tumor cells fail to complete the latter stages of myogenesis. The RMS cell line RD-A was treated with 12-O-tetradecanoylphorbol-13-acetate to induce differentiation and cultured for 10 days. RNA was extracted on days 1, 3, 6, 8 and 10. A human skeletal muscle cDNA microarray was developed and used to analyze the global gene expression of RMS tumors over the time-course of differentiation. As a comparison, the genes identified were subsequently examined during the differentiated primary human skeletal muscle cultures. Prothymosin alpha (PTMA), and translocase of inner mitochondrial membrane 10 (Tim10), two genes not previously implicated in RMS, showed reduced expression during differentiation. Marked differences in the expression of PTMA and Tim10 were observed during the differentiation of human primary skeletal muscle cells. These results identify several new genes with potential roles in the myogenic arrest present in rhabdomyosarcoma. PTMA expression in RMS biopsy samples might prove to be an effective diagnostic marker for this disease.