Yeast microarrays for genome wide parallel genetic and gene expression analysis

We have developed high-density DNA microarrays of yeast ORFs. These microarrays can monitor hybridization to ORFs for applications such as quantitative differential gene expression analysis and screening for sequence polymorphisms. Automated scripts retrieved sequence information from public databases to locate predicted ORFs and select appropriate primers for amplification. The primers were used to amplify yeast ORFs in 96-well plates, and the resulting products were arrayed using an automated micro arraying device. Arrays containing up to 2,479 yeast ORFs were printed on a single slide. The hybridization of fluorescently labeled samples to the array were detected and quantitated with a laser confocal scanning microscope. Applications of the microarrays are shown for genetic and gene expression analysis at the whole genome level.

Identification of diverse nerve growth factor-regulated genes by serial analysis of gene expression (SAGE) profiling

Neurotrophic factors such as nerve growth factor (NGF) promote a wide variety of responses in neurons, including differentiation, survival, plasticity, and repair. Such actions often require changes in gene expression. To identify the regulated genes and thereby to more fully understand the NGF mechanism, we carried out serial analysis of gene expression (SAGE) profiling of transcripts derived from rat PC12 cells before and after NGF-promoted neuronal differentiation. Multiple criteria supported the reliability of the profile. Approximately 157,000 SAGE tags were analyzed, representing at least 21,000 unique transcripts. Of these, nearly 800 were regulated by 6-fold or more in response to NGF. Approximately 150 of the regulated transcripts have been matched to named genes, the majority of which were not previously known to be NGF-responsive. Functional categorization of the regulated genes provides insight into the complex, integrated mechanism by which NGF promotes its multiple actions. It is anticipated that as genomic sequence information accrues the data derived here will continue to provide information about neurotrophic factor mechanisms.

BodyMap incorporated PCR-based expression profiling data and a gene ranking system

BodyMap is a human and mouse gene expression database that is based on site-directed 3′-expressed sequence tags generated at Osaka University. To date, it contains more than 300 000 tag sequences from 64 human and 39 mouse tissues. For the recent release, the precise anatomical expression patterns for more than half of the human gene entries were generated by introduced amplified fragment length polymorphism (iAFLP), which is a PCR-based high-throughput expression profiling method. The iAFLP data incorporated into BodyMap describe the relative contents of more than 12 000 transcripts across 30 tissue RNAs. In addition, a newly developed gene ranking system helps users obtain lists of genes that have desired expression patterns according to their significance. BodyMap supports complete transfer of unique data sets and provides analysis that is accessible through the WWW at http://bodymap.ims.u-tokyo.ac.jp.

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.

Systematic expression profiling of the mouse transcriptome using RIKEN cDNA mircoarrays

The number of known mRNA transcripts in the mouse has been greatly expanded by the RIKEN Mouse Gene Encyclopedia project. Validation of their reproducible expression in a tissue is an important contribution to the study of functional genomics. In this report, we determine the expression profile of 57,931 clones on 20 mouse tissues using cDNA microarrays. Of these 57,931 clones, 22,928 clones correspond to the FANTOM2 clone set. The set represents 20,234 transcriptional units (TUs) out of 33,409 TUs in the FANTOM2 set. We identified 7206 separate clones that satisfied stringent criteria for tissue-specific expression. Gene Ontology terms were assigned for these 7206 clones, and the proportion of 'molecular function' ontology for each tissue-specific clone was examined. These data will provide insights into the function of each tissue. Tissue-specific gene expression profiles obtained using our cDNA microarrays were also compared with the data extracted from the GNF Expression Atlas based on Affymetrix microarrays. One major outcome of the RIKEN transcriptome analysis is the identification of numerous nonprotein-coding mRNAs. The expression profile was also used to obtain evidence of expression for putative noncoding RNAs. In addition, 1926 clones (70%) of 2768 clones that were categorized as unknown EST, and 1969 (58%) clones of 3388 clones that were categorized as unclassifiable were also shown to be reproducibly expressed.

A simple implementation of a normal mixture approach to differential gene expression in multiclass microarrays

Motivation: An important problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. We provide a straightforward and easily implemented method for estimating the posterior probability that an individual gene is null. The problem can be expressed in a two-component mixture framework, using an empirical Bayes approach. Current methods of implementing this approach either have some limitations due to the minimal assumptions made or with more specific assumptions are computationally intensive. Results: By converting to a z-score the value of the test statistic used to test the significance of each gene, we propose a simple two-component normal mixture that models adequately the distribution of this score. The usefulness of our approach is demonstrated on three real datasets.

MediPlEx - a tool to combine in silico & experimental gene expression profiles of the model legume Medicago truncatula

Background: Expressed Sequence Tags (ESTs) are in general used to gain a first insight into gene activities from a species of interest. Subsequently, and typically based on a combination of EST and genome sequences, microarray-based expression analyses are performed for a variety of conditions. In some cases, a multitude of EST and microarray experiments are conducted for one species, covering different tissues, cell states, and cell types. Under these circumstances, the challenge arises to combine results derived from the different expression profiling strategies, with the goal to uncover novel information on the basis of the integrated datasets. Findings: Using our new analysis tool, MediPlEx (MEDIcago truncatula multiPLe EXpression analysis), expression data from EST experiments, oligonucleotide microarrays and Affymetrix GeneChips® can be combined and analyzed, leading to a novel approach to integrated transcriptome analysis. We have validated our tool via the identification of a set of well-characterized AM-specific and AM-induced marker genes, identified by MediPlEx on the basis of in silico and experimental gene expression profiles from roots colonized with AM fungi. Conclusions: MediPlEx offers an integrated analysis pipeline for different sets of expression data generated for the model legume Medicago truncatula. As expected, in silico and experimental gene expression data that cover the same biological condition correlate well. The collection of differentially expressed genes identified via MediPlEx provides a starting point for functional studies in plant mutants.

Wolbachia uses host microRNAs to manipulate host gene expression and facilitate colonization of the dengue vector Aedes aegypti

The obligate endosymbiont Wolbachia pipientis is found in a wide range of invertebrates where they are best known for manipulating host reproduction. Recent studies have shown that Wolbachia also can modulate the lifespan of host insects and interfere with the development of human pathogens in mosquito vectors. Despite considerable study, very little is known about the molecular interactions between Wolbachia and its hosts that might mediate these effects. Using microarrays, we show that the microRNA (miRNA) profile of the mosquito, Aedes aegypti, is significantly altered by the wMelPop-CLA strain of W. pipientis. We found that a host miRNA (aae-miR-2940) is induced after Wolbachia infection in both mosquitoes and cell lines. One target of aae-miR-2940 is the Ae. aegypti metalloprotease gene. Interestingly, expression of the target gene was induced after Wolbachia infection, ectopic expression of the miRNA independent of Wolbachia, or transfection of an artificial mimic of the miRNA into mosquito cells. We also confirmed the interaction of aae-miR-2940 with the target sequences using GFP as a reporter gene. Silencing of the metalloprotease gene in both Wolbachia-infected cells and adult mosquitoes led to a significant reduction in Wolbachia density, as did inhibition of the miRNA in cells. These results indicate that manipulation of the mosquito metalloprotease gene via aae-miR-2940 is crucial for efficient maintenance of the endosymbiont. This report shows how Wolbachia alters the host miRNA profile and provides insight into the mechanisms of host manipulation used by this widespread endosymbiont.

Straminipilan protists : a study of their taxonomic position, morphology and abiotic stress-mediated gene expression

Thraustochytrids have become of considerable industrial and scientific interest in the past decade due to their health benefits. They have been proven to be the principle source in marine and estuarine fish diets with high percentage of long chain (LC) or polyunsaturated fatty acids (PUFA). Therefore, the oil extracted from fish for human document.forms[0].elements[13].select();consumption is rich in PUFA with high omega-3 fatty acid content. Docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) of all of the omega-3 fatty acids, are considered beneficial essential oils for humans with a wide range of health benefits. These include brain and neural development in infants, general wellbeing of adults and drug delivery through precursor molecules. They have become one of the most extensively studied organisms for industrial oil preparations as PUFA extraction from fish becomes less profitable. Many forms of these Thraustochytrid oils are being trialled for human consumption all over the world. In Australia, there has been little research performed on these organisms in the past ten years. A few Australian studies have been conducted in the form of comparative studies related to PUFA production within the related genera, but not focussed on their identification or cellular and genomic characterisation. Therefore, the main aim of this study was to investigate the morphological and genetic characteristics of Australian Thraustochytrids in order to aid in their identification and characterisation, as well as to better understand the effect of environmental conditions in the regulation of PUFA production. It was also noted that there was a knowledge gap in the preservation and total genomic DNA extraction of these organisms for the purposes of scientific research. The cryopreservation of these organisms for studies around the world follows existing generic methods. However, it is well understood that many of these generic methods attract not only high costs for chemicals, but also uses considerable storage space and other resources, all of which can be improved with new or modified approaches. In this context, a simple and inexpensive bead preservation method is described, without compromising the storage shelf life. We also describe, for the first time, the effects of culture age on the successful cryopreservation of Thraustochytrids. It was evident in the literature that DNA and RNA extractions for molecular and genetic studies of Thraustochytrids follow the classical phenol-chloroform extraction methods. It was also observed that modern protocols failed to avoid the use of phenol-chloroform rather than improving preparation and cell disruption. In order to provide a high quantity and quality DNA extraction, a modified protocol has been introduced that employs the use of modern commercial extraction kits and standard laboratory equipment. Thraustochytrids have been shown to be highly conserved in their 18S rDNA gene sequences, which is used as the current standard for identification. It was demonstrated that the 18S rDNA gene sequence limits the recognition of closely related genera or within the genera from each member. Therefore, it was proposed that another profile, such as a randomly amplified polymorphic DNA (RAPD) based profiling system, be tested for use in the characterisation of Thraustochytrids. The RAPD profiles were shown to provide a unique DNA fingerprint for each isolate and small variations in their genome were able to be detected. This method involved the use of a minimum number of standard arbitrary primers and with an increase in the number of different primers used, a very high discrimination between organisms could be achieved. However, the method was not suitable for taxonomic purposes because the results did not correlate with other taxonomic features such as morphology. Another knowledge gap was found with respect to Australian Thraustochytrid growth characteristics, in that these had not been recorded and published. In order to rectify this, a record of colony and microscopic features of 12 selected isolates was performed. The results of preliminary studies indicated that further microbiological and biochemical studies are needed for full characterisation of these organisms. This information is of great importance to bio-prospecting of new Thraustochytrids from Australian ecosystems and would allow for their accurate identification, and so permit the prediction of their PUFA capability by comparison with related genera/species. It was well recognized that environmental stress plays a role in the PUFA production and is mainly due to the reactive oxygen species as abiotic stress (Chiou et al., 2001; Okuyama et al., 2008; Shabala et al., 2009; Shabala et al., 2001). In this aspect, this study makes the first attempt towards better understanding of this phenomenon by way of the use of real-time PCR for the detection of environmental effects on the regulation of PUFA production. Three main environmental conditions including temperature, pH and oxygen availability were monitored as stress inducers. In summary, this study provides novel approaches for the preservation and handling of Thraustochytrids, their molecular biological features, taxonomy, characterisation and responses to environmental factors with respect to their oil production enzymes. The information produced from this study will prove to be vital for both industrial and scientific investigations in the future.

Differential expression profiling of components associated with exoskeletal hardening in crustaceans

Background: Exoskeletal hardening in crustaceans can be attributed to mineralization and sclerotization of the organic matrix. Glycoproteins have been implicated in the calcification process of many matrices. Sclerotization, on the other hand, is catalysed by phenoloxidases, which also play a role in melanization and the immunological response in arthropods. Custom cDNA microarrays from Portunus pelagicus were used to identify genes possibly associated with the activation pathways involved in these processes. Results: Two genes potentially involved in the recognition of glycosylation, the C-type lectin receptor and the mannose-binding protein, were found to display molt cycle-related differential expression profiles. C-type lectin receptor up-regulation was found to coincide with periods associated with new uncalcified cuticle formation, while the up-regulation of mannose-binding protein occurred only in the post-molt stage, during which calcification takes place, implicating both in the regulation of calcification. Genes presumed to be involved in the phenoloxidase activation pathway that facilitates sclerotization also displayed molt cycle-related differential expression profiles. Members of the serine protease superfamily, trypsin-like and chymotrypsin-like, were up-regulated in the intermolt stage when compared to post-molt, while trypsin-like was also up-regulated in pre-molt compared to ecdysis. Additionally, up-regulation in pre- and intermolt stages was observed by transcripts encoding other phenoloxidase activators including the putative antibacterial protein carcinin-like, and clotting protein precursor-like. Furthermore, hemocyanin, itself with phenoloxidase activity, displayed an identical expression pattern to that of the phenoloxidase activators, i.e. up-regulation in pre- and intermolt. Conclusion: Cuticle hardening in crustaceans is a complex process that is precisely timed to occur in the post-molt stage of the molt cycle. We have identified differential expression patterns of several genes that are believed to be involved in biomineralization and sclerotization and propose possible regulatory mechanisms for these processes based on their expression profiles, such as the potential involvement of C-type lectin receptors and mannose binding protein in the regulation of calcification.

Gene expression evidence for off-target effects caused by RNA interference-mediated gene silencing of Ubiquitin-63E in the cattle tick Rhipicephalus microplus.

Knowledge of cattle tick (Rhipicephalus (Boophilus) microplus; Acari: Ixodidae) molecular and cellular pathways has been hampered by the lack of an annotated genome. In addition, most of the tick expressed sequence tags (ESTs) available to date consist of similar to 50% unassigned sequences without predicted functions. The most common approach to address this has been the application of RNA interference (RNAi) methods to investigate genes and their pathways. This approach has been widely adopted in tick research despite minimal knowledge of the tick RNAi pathway and double-stranded RNA (dsRNA) uptake mechanisms. A strong knockdown phenotype of adult female ticks had previously been observed using a 594 bp dsRNA targeting the cattle tick homologue for the Drosophila Ubiquitin-63E gene leading to nil or deformed eggs. A NimbleGen cattle tick custom microarray based on the BmiGI.V2 database of R. microplus ESTs was used to evaluate the expression of mRNAs harvested from ticks treated with the tick Ubiquitin-63E 594 bp dsRNA compared with controls. A total of 144 ESTs including TC6372 (Ubiquitin-63E) were down-regulated with 136 ESTs up-regulated following treatment. The results obtained substantiated the knockdown phenotype with ESTs identified as being associated with ubiquitin proteolysis as well as oogenesis, embryogenesis, fatty acid synthesis and stress responses. A bioinformatics analysis was undertaken to predict off-target effects (OTE) resulting from the in silico dicing of the 594 bp Ubiquitin-63E dsRNA which identified 10 down-regulated ESTs (including TC6372) within the list of differentially expressed probes on the microarrays. Subsequent knockdown experiments utilising 196 and 109 bp dsRNAs, and a cocktail of short hairpin RNAs (shRNA) targeting Ubiquitin-63E, demonstrated similar phenotypes for the dsRNAs but nil effect following shRNA treatment. Quantitative reverse transcriptase PCR analysis confirmed differential expression of TC6372 and selected ESTs. Our study demonstrated the minimisation of predicted OTEs in the shorter dsRNA treatments (similar to 100-200 bp) and the usefulness of microarrays to study knockdown phenotypes.

Genomic and functional profiling of gastric cancer

Helicobacter pylori infection is a risk factor for gastric cancer, which is a major health issue worldwide. Gastric cancer has a poor prognosis due to the unnoticeable progression of the disease and surgery is the only available treatment in gastric cancer. Therefore, gastric cancer patients would greatly benefit from identifying biomarker genes that would improve diagnostic and prognostic prediction and provide targets for molecular therapies. DNA copy number amplifications are the hallmarks of cancers in various anatomical locations. Mechanisms of amplification predict that DNA double-strand breaks occur at the margins of the amplified region. The first objective of this thesis was to identify the genes that were differentially expressed in H. pylori infection as well as the transcription factors and signal transduction pathways that were associated with the gene expression changes. The second objective was to identify putative biomarker genes in gastric cancer with correlated expression and copy number, and the last objective was to characterize cancers based on DNA copy number amplifications. DNA microarrays, an in vitro model and real-time polymerase chain reaction were used to measure gene expression changes in H. pylori infected AGS cells. In order to identify the transcription factors and signal transduction pathways that were activated after H. pylori infection, gene expression profiling data from the H. pylori experiments and a bioinformatics approach accompanied by experimental validation were used. Genome-wide expression and copy number microarray analysis of clinical gastric cancer samples and immunohistochemistry on tissue microarray were used to identify putative gastric cancer genes. Data mining and machine learning techniques were applied to study amplifications in a cross-section of cancers. FOS and various stress response genes were regulated by H. pylori infection. H. pylori regulated genes were enriched in the chromosomal regions that are frequently changed in gastric cancer, suggesting that molecular pathways of gastric cancer and premalignant H. pylori infection that induces gastritis are interconnected. 16 transcription factors were identified as being associated with H. pylori infection induced changes in gene expression. NF-κB transcription factor and p50 and p65 subunits were verified using elecrophoretic mobility shift assays. ERBB2 and other genes located in 17q12- q21 were found to be up-regulated in association with copy number amplification in gastric cancer. Cancers with similar cell type and origin clustered together based on the genomic localization of the amplifications. Cancer genes and large genes were co-localized with amplified regions and fragile sites, telomeres, centromeres and light chromosome bands were enriched at the amplification boundaries. H. pylori activated transcription factors and signal transduction pathways function in cellular mechanisms that might be capable of promoting carcinogenesis of the stomach. Intestinal and diffuse type gastric cancers showed distinct molecular genetic profiles. Integration of gene expression and copy number microarray data allowed the identification of genes that might be involved in gastric carcinogenesis and have clinical relevance. Gene amplifications were demonstrated to be non-random genomic instabilities. Cell lineage, properties of precursor stem cells, tissue microenvironment and genomic map localization of specific oncogenes define the site specificity of DNA amplifications, whereas labile genomic features define the structures of amplicons. These conclusions suggest that the definition of genomic changes in cancer is based on the interplay between the cancer cell and the tumor microenvironment.

Construction of a global map of human gene expression : the process, tools and analysis

This thesis studies human gene expression space using high throughput gene expression data from DNA microarrays. In molecular biology, high throughput techniques allow numerical measurements of expression of tens of thousands of genes simultaneously. In a single study, this data is traditionally obtained from a limited number of sample types with a small number of replicates. For organism-wide analysis, this data has been largely unavailable and the global structure of human transcriptome has remained unknown. This thesis introduces a human transcriptome map of different biological entities and analysis of its general structure. The map is constructed from gene expression data from the two largest public microarray data repositories, GEO and ArrayExpress. The creation of this map contributed to the development of ArrayExpress by identifying and retrofitting the previously unusable and missing data and by improving the access to its data. It also contributed to creation of several new tools for microarray data manipulation and establishment of data exchange between GEO and ArrayExpress. The data integration for the global map required creation of a new large ontology of human cell types, disease states, organism parts and cell lines. The ontology was used in a new text mining and decision tree based method for automatic conversion of human readable free text microarray data annotations into categorised format. The data comparability and minimisation of the systematic measurement errors that are characteristic to each lab- oratory in this large cross-laboratories integrated dataset, was ensured by computation of a range of microarray data quality metrics and exclusion of incomparable data. The structure of a global map of human gene expression was then explored by principal component analysis and hierarchical clustering using heuristics and help from another purpose built sample ontology. A preface and motivation to the construction and analysis of a global map of human gene expression is given by analysis of two microarray datasets of human malignant melanoma. The analysis of these sets incorporate indirect comparison of statistical methods for finding differentially expressed genes and point to the need to study gene expression on a global level.

Expression profiling of lymph nodes in tuberculosis patients reveal inflammatory milieu at site of infection

Extrapulmonary manifestations constitute 15 to 20% of tuberculosis cases, with lymph node tuberculosis (LNTB) as the most common form of infection. However, diagnosis and treatment advances are hindered by lack of understanding of LNTB biology. To identify host response, Mycobacterium tuberculosis infected lymph nodes from LNTB patients were studied by means of transcriptomics and quantitative proteomics analyses. The selected targets obtained by comparative analyses were validated by quantitative PCR and immunohistochemistry. This approach provided expression data for 8,728 transcripts and 102 proteins, differentially regulated in the infected human lymph node. Enhanced inflammation with upregulation of T-helper1-related genes, combined with marked dysregulation of matrix metalloproteinases, indicates tissue damage due to high immunoactivity at infected niche. This expression signature was accompanied by significant upregulation of an immunoregulatory gene, leukotriene A4 hydrolase, at both transcript and protein levels. Comparative transcriptional analyses revealed LNTB-specific perturbations. In contrast to pulmonary TB-associated increase in lipid metabolism, genes involved in fatty-acid metabolism were found to be downregulated in LNTB suggesting differential lipid metabolic signature. This study investigates the tissue molecular signature of LNTB patients for the first time and presents findings that indicate the possible mechanism of disease pathology through dysregulation of inflammatory and tissue-repair processes.