On the simultaneous use of clinical and microarray expression data in the cluster analysis of tissue samples

This paper considers a model-based approach to the clustering of tissue samples of a very large number of genes from microarray experiments. It is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. Frequently in practice, there are also clinical data available on those cases on which the tissue samples have been obtained. Here we investigate how to use the clinical data in conjunction with the microarray gene expression data to cluster the tissue samples. We propose two mixture model-based approaches in which the number of components in the mixture model corresponds to the number of clusters to be imposed on the tissue samples. One approach specifies the components of the mixture model to be the conditional distributions of the microarray data given the clinical data with the mixing proportions also conditioned on the latter data. Another takes the components of the mixture model to represent the joint distributions of the clinical and microarray data. The approaches are demonstrated on some breast cancer data, as studied recently in van't Veer et al. (2002).

Statistical analysis and meta-analysis of microarray data

The microarray technology provides a high-throughput technique to study gene expression. Microarrays can help us diagnose different types of cancers, understand biological processes, assess host responses to drugs and pathogens, find markers for specific diseases, and much more. Microarray experiments generate large amounts of data. Thus, effective data processing and analysis are critical for making reliable inferences from the data. ^ The first part of dissertation addresses the problem of finding an optimal set of genes (biomarkers) to classify a set of samples as diseased or normal. Three statistical gene selection methods (GS, GS-NR, and GS-PCA) were developed to identify a set of genes that best differentiate between samples. A comparative study on different classification tools was performed and the best combinations of gene selection and classifiers for multi-class cancer classification were identified. For most of the benchmarking cancer data sets, the gene selection method proposed in this dissertation, GS, outperformed other gene selection methods. The classifiers based on Random Forests, neural network ensembles, and K-nearest neighbor (KNN) showed consistently god performance. A striking commonality among these classifiers is that they all use a committee-based approach, suggesting that ensemble classification methods are superior. ^ The same biological problem may be studied at different research labs and/or performed using different lab protocols or samples. In such situations, it is important to combine results from these efforts. The second part of the dissertation addresses the problem of pooling the results from different independent experiments to obtain improved results. Four statistical pooling techniques (Fisher inverse chi-square method, Logit method. Stouffer's Z transform method, and Liptak-Stouffer weighted Z-method) were investigated in this dissertation. These pooling techniques were applied to the problem of identifying cell cycle-regulated genes in two different yeast species. As a result, improved sets of cell cycle-regulated genes were identified. The last part of dissertation explores the effectiveness of wavelet data transforms for the task of clustering. Discrete wavelet transforms, with an appropriate choice of wavelet bases, were shown to be effective in producing clusters that were biologically more meaningful. ^

Alternative fungicides for controlling husk spot caused by Pseudocercospora macadamiae in macadamia

Husk spot, caused by Pseudocercospora macadamiae is a major fungal disease of macadamia in Australia. Chemicals to control the disease are limited and frequent failure to control the disease is a major concern to growers. The overall goal of this research was to improve the chemical control strategy of P. macadamiae through the provision of fungicides with different modes of action to carbendazim, which is the current industry standard. Husk spot incidence, premature fruit abscission, kernel quality and yield were evaluated following application of different fungicide products in replicated field experiments at three different sites. Results showed significant differences in disease incidence and premature fruit abscission between fungicide treatments, field sites and years. Generally, disease incidence and premature fruit abscission on trees treated with fungicide were significantly (P < 0.05) lower than the untreated control. Pyraclostrobin conferred significantly better protection than trifloxystrobin, reducing disease severity by 70% compared with a 50% reduction by trifloxystrobin. The pyraclostrobin treatment had a similar efficacy to the current industry standard (70% reduction cf. 73% reduction by tank-mixed carbendazim and copper). Higher amounts of immature kernels occurred in the untreated control, followed by difenoconazole and trifloxystrobin. Diseased fruit accounted for 78% of premature fruit abscission, which indicates that husk spot enhances fruit abscission in macadamia. Our results suggest that pyraclostrobin provided similar efficacy to the industry standard and could, therefore, play a key role in the management of husk spot.

Rice responses to soil management in a rice-based cropping system in the semi-arid tropics of southern Lombok, Eastern Indonesia

This paper is the first of a series that investigates whether new cropping systems with permanent raised beds (PRBs) or Flat land could be successfully used to increase farmers' incomes from rainfed crops in Lombok in Eastern Indonesia. This paper discusses the rice phase of the cropping system. Low grain yields of dry-seeded rice (Oryza sativa) grown on Flat land on Vertisols in the rainfed region of southern Lombok, Eastern Indonesia, are probably mainly due to (a) erratic rainfall (870-1220 mm/yr), with water often limiting at sensitive growth stages, (b) consistently high temperatures (average maximum - 31 C), and (c) low solar radiation. Farmers are therefore poor, and labour is hard and costly, as all operations are manual. Two replicated field experiments were run at Wakan (annual rainfall = 868 mm) and Kawo (1215 mm) for 3 years (2001/2002 to 2003/2004) on Vertisols in southern Lombok. Dry-seeded rice was grown in 4 treatments with or without manual tillage on (a) PRBs, 1.2 m wide, 200 mm high, separated by furrows 300 mm wide, 200 mill deep, with no rice sown in the well-graded furrows, and (b) well-graded Flat land. Excess surface water was harvested from each treatment and used for irrigation after the vegetative stage of the rice. All operations were manual. There were no differences between treatments in grain yield of rice (mean grain yield = 681 g/m(2)) which could be partly explained by total number of tillers/hill and mean panicle length, but not number of productive tillers/hill, plant height or weight of 1000 grains. When the data from both treatments on PRBs and from both treatments on Flat land, each year at each site were analysed, there were also no differences in grain yield of rice (g/m(2)). When rainfall in the wet season up to harvest was over 1000 mm (Year 2; Wakan, Kawo), or plants were water-stressed during crop establishment (Year 1; Wakan) or during grain-fill (Year 3: Kawo), there were significant differences in grain yield (g/1.5 m(2)) between treatments; generally the grain yield (g/1.5 m(2)) on PRBs with or without tillage was less than that on Flat land with or without tillage. However, when the data from both treatments on PRBs and from both treatments on Flat land, each year at each site, were analysed, the greater grain yield of dry-seeded rice on Flat land (mean yield 1 092 g/1.5 m(2)) than that on PRBs (mean 815 g/1.5 m(2)) was mainly because there were 25% more plants on Flat land. Overall when the data in the 2 outer rows and the 2 inner rows on PRBs were each combined, there was a higher number of productive tillers in the combined outer rows (mean 20.7 tillers/hill) compared with that in the combined inner rows on each PRB (mean 18.2 tillers/hill). However, there were no differences in grain yield between combined rows (mean 142 g/m row). Hence with a gap of 500 mm (the distance between the outer rows of plants on adjacent raised beds), plants did not compensate in grain yield for missing plants in furrows. This suggests that rice (a) also sown in furrows, or (b) sown in 7 rows with narrower row-spacing, or (c) sown in 6 rows with slightly wider row-spacing, and narrower gap between outer rows on adjacent beds, may further increase grain yield (g/1.5 m(2)) in this system of PRBs. The growth and the grain yield (y in g/m(2)) of rainfed rice (with rainfall on-site the only source of water for irrigation) depended mainly on the rainfall (x in mm) in the wet season up to harvest (due either to site or year) with y = 1. 1x -308; r(2) = 0.54; p < 0.005. However, 280 mm (i.e. 32%) of the rainfall was not directly used to produce grain (i.e. when y = 0 g/m(2)). Manual tillage did not affect growth and grain yield of rice (g/m(2); g/1.5 m(2)), either on PRB or on Flat land.

Resistance to the root-lesion nematode Pratylenchus thornei in wheat landraces and cultivars from the West Asia and North Africa (WANA) region.

Resistance to the root-lesion nematode Pratylenchus thornei was sought in wheat from the West Asia and North Africa (WANA) region in the Watkins Collection (148 bread and 139 durum wheat accessions) and the McIntosh Collection (59 bread and 43 durum wheat accessions). It was considered that landraces from this region, encompassing the centres of origin of wheat and where P. thornei also occurs, could be valuable sources of resistance for use in wheat breeding. Resistance was determined by number of P. thornei/kg soil after the growth of the plants in replicated glasshouse experiments. On average, durum accessions produced significantly lower numbers of P. thornei than bread wheat accessions in both the Watkins and McIntosh Collections. Selected accessions with low P. thornei numbers were re-tested and 13 bread wheat and 10 durum accessions were identified with nematode numbers not significantly different from GS50a, a partially resistant bread wheat line used as a reference standard. These resistant accessions, which originated in Iran, Iraq, Syria, Egypt, Sudan, Morocco, and Tunisia, represent a resource of resistance genes in the primary wheat gene pool, which could be used in Australian wheat breeding programs to reduce the economic loss from P. thornei.

The impacts and management of foxes Vulpes vulpes in Australia

1. The successful introduction of the red fox Vulpes vulpes into Australia in the 1870s has had dramatic and deleterious impacts on both native fauna and agricultural production. Historical accounts detail how the arrival of foxes in many areas coincided with the local demise of native fauna. Recent analyses suggest that native fauna can be successfully reintroduced to their former ranges only if foxes have been controlled, and several replicated removal experiments have confirmed that foxes are the major agents of extirpation of native fauna. Predation is the primary cause of losses, but competition and transmission of disease may be important for some species. 2. In agricultural landscapes, fox predation on lambs can cause losses of 1–30%; variation is due to flock size, health and management, as well as differences in the timing and duration of lambing and the density of foxes. 3. Fox control measures include trapping, shooting, den fumigation and exclusion fencing; baiting using the toxin 1080 is the most commonly employed method. Depending on the baiting strategy, habitat and area covered, baiting can reduce fox activity by 50–97%. We review patterns of baiting in a large sheep-grazing region in central New South Wales, and propose guidelines to increase landholder awareness of baiting strategies, to concentrate and coordinate bait use, and to maximize the cost-effectiveness of baiting programs. 4. The variable reduction in fox density within the baited area, together with the ability of the fox to recolonize rapidly, suggest that current baiting practices in eastern Australia are often ineffective, and that reforms are required. These might include increasing landholder awareness and involvement in group control programs, and the use of more efficient broadscale techniques, such as aerial baiting.

QTLs for water absorption and flour yield identified in the doubled haploid wheat population Lang/QT8766

Two key quality traits in milling wheat are flour yield (FY) and water absorption (WA). Ideally, breeders would prefer to use markers to select promising lines rather than time consuming rheology tests. In this study, we measured FY and WA on a wheat mapping population (Lang/QT8766) of 162 individuals grown in two replicated field experiments at three locations over 2 years. We also carried out near infrared reflectance spectroscopy (NIRS) predictions on the grain for these traits to see if NIRS phenotypic data could provide useful mapping results when compared to the reference phenotypic data. Several common QTLs were identified for FY and WA by both sets of data. The QTL on chromosome 4D was a consistently recurring QTL region for both traits. The QTL on chromosome 2A was positively linked to protein content which was supported by genetic correlation data. The results also indicated it was possible to obtain useful phenotypic data for mapping FY and WA using NIRS data. This would save time and costs as NIRS is quicker and cheaper than current rheology methods.

Methods to improve gene signal: Application to cDNA microarrays

Microarrays are high throughput biological assays that allow the screening of thousands of genes for their expression. The main idea behind microarrays is to compute for each gene a unique signal that is directly proportional to the quantity of mRNA that was hybridized on the chip. A large number of steps and errors associated with each step make the generated expression signal noisy. As a result, microarray data need to be carefully pre-processed before their analysis can be assumed to lead to reliable and biologically relevant conclusions. This thesis focuses on developing methods for improving gene signal and further utilizing this improved signal for higher level analysis. To achieve this, first, approaches for designing microarray experiments using various optimality criteria, considering both biological and technical replicates, are described. A carefully designed experiment leads to signal with low noise, as the effect of unwanted variations is minimized and the precision of the estimates of the parameters of interest are maximized. Second, a system for improving the gene signal by using three scans at varying scanner sensitivities is developed. A novel Bayesian latent intensity model is then applied on these three sets of expression values, corresponding to the three scans, to estimate the suitably calibrated true signal of genes. Third, a novel image segmentation approach that segregates the fluorescent signal from the undesired noise is developed using an additional dye, SYBR green RNA II. This technique helped in identifying signal only with respect to the hybridized DNA, and signal corresponding to dust, scratch, spilling of dye, and other noises, are avoided. Fourth, an integrated statistical model is developed, where signal correction, systematic array effects, dye effects, and differential expression, are modelled jointly as opposed to a sequential application of several methods of analysis. The methods described in here have been tested only for cDNA microarrays, but can also, with some modifications, be applied to other high-throughput technologies. Keywords: High-throughput technology, microarray, cDNA, multiple scans, Bayesian hierarchical models, image analysis, experimental design, MCMC, WinBUGS.

Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

Gene Expression : From Microarrays to Functional Genomics

The time of the large sequencing projects has enabled unprecedented possibilities of investigating more complex aspects of living organisms. Among the high-throughput technologies based on the genomic sequences, the DNA microarrays are widely used for many purposes, including the measurement of the relative quantity of the messenger RNAs. However, the reliability of microarrays has been strongly doubted as robust analysis of the complex microarray output data has been developed only after the technology had already been spread in the community. An objective of this study consisted of increasing the performance of microarrays, and was measured by the successful validation of the results by independent techniques. To this end, emphasis has been given to the possibility of selecting candidate genes with remarkable biological significance within specific experimental design. Along with literature evidence, the re-annotation of the probes and model-based normalization algorithms were found to be beneficial when analyzing Affymetrix GeneChip data. Typically, the analysis of microarrays aims at selecting genes whose expression is significantly different in different conditions followed by grouping them in functional categories, enabling a biological interpretation of the results. Another approach investigates the global differences in the expression of functionally related groups of genes. Here, this technique has been effective in discovering patterns related to temporal changes during infection of human cells. Another aspect explored in this thesis is related to the possibility of combining independent gene expression data for creating a catalog of genes that are selectively expressed in healthy human tissues. Not all the genes present in human cells are active; some involved in basic activities (named housekeeping genes) are expressed ubiquitously. Other genes (named tissue-selective genes) provide more specific functions and they are expressed preferably in certain cell types or tissues. Defining the tissue-selective genes is also important as these genes can cause disease with phenotype in the tissues where they are expressed. The hypothesis that gene expression could be used as a measure of the relatedness of the tissues has been also proved. Microarray experiments provide long lists of candidate genes that are often difficult to interpret and prioritize. Extending the power of microarray results is possible by inferring the relationships of genes under certain conditions. Gene transcription is constantly regulated by the coordinated binding of proteins, named transcription factors, to specific portions of the its promoter sequence. In this study, the analysis of promoters from groups of candidate genes has been utilized for predicting gene networks and highlighting modules of transcription factors playing a central role in the regulation of their transcription. Specific modules have been found regulating the expression of genes selectively expressed in the hippocampus, an area of the brain having a central role in the Major Depression Disorder. Similarly, gene networks derived from microarray results have elucidated aspects of the development of the mesencephalon, another region of the brain involved in Parkinson Disease.

Gene expression profile of epithelial cells and mesenchymal cells derived from limbal explant culture

Purpose: Limbal stem cell deficiency is a challenging clinical problem and the current treatment involves replenishing the depleted limbal stem cell (LSC) pool by either limbal tissue transplantation or use of cultivated limbal epithelial cells (LEC). Our experience of cultivating the LEC on denuded human amniotic membrane using a feeder cell free method, led to identification of mesenchymal cells of limbus (MC-L), which showed phenotypic resemblance to bone marrow derived mesenchymal stem cells (MSC-BM). To understand the transcriptional profile of these cells, microarray experiments were carried out.Methods: RNA was isolated from cultured LEC, MC-L and MSC-BM and microarray experiments were carried out by using Agilent chip (4x44 k). The microarray data was validated by using Realtime and semiquntitative reverse transcription polymerase chain reaction. Results: The microarray analysis revealed specific gene signature of LEC and MC-L, and also their complementary role related to cytokine and growth factor profile, thus supporting the nurturing roles of the MC-L. We have also observed similar and differential gene expression between MC-L and MSC-BM.Conclusions: This study represents the first extensive gene expression analysis of limbal explant culture derived epithelial and mesenchymal cells and as such reveals new insight into the biology, ontogeny, and in vivo function of these cells.

Role of GDNF and its Cross-Talk with Other Growth Factors in the Dopaminergic System

Parkinson´s Disease (PD) is a neurodegenerative movement disorder resulting from loss of dopaminergic (DA) neurons in substantia nigra (SN). Possible causative treatment strategies for PD include neurotrophic factors, which protect and in some cases restore the function of dopaminergic neurons. Glial cell line-derived neurotrophic factor (GDNF) family of neurotrophic factors have been to date the most promising candidates for treatment of PD, demonstrating both neuroprotective and neurorestorative properties. We have investigated the role of GDNF in the rodent dopaminergic system and its possible crosstalk with other growth factors. We characterized the GDNF-induced gene expression changes by DNA microarray analysis in different neuronal systems, including in vitro cultured Neuro2A cells treated with GDNF, as well as midbrains from GDNF heterozygous (Hz) knockout mice. These microarray experiments, resulted in the identification of GDNF-induced genes, which were also confirmed by other methods. Further analysis of the dopaminergic system of GDNF Hz mice demonstrated about 40% reduction in GDNF levels, revealed increased intracellular dopamine concentrations and FosB/DeltaFosB expression in striatal areas. These animals did not show any significant changes in behavioural analysis of acute and repeated cocaine administration on locomotor activity, nor did they exhibit any changes in dopamine output following treatment with acute cocaine. We further analysed the significance of GDNF receptor RET signalling in dopaminergic system of MEN2B knock-in animals with constitutively active Ret. The MEN2B animals showed a robust increase in extracellular dopamine and its metabolite levels in striatum, increased tyrosine hydroxylase (TH) and dopamine transporter (DAT) protein levels by immunohistochemical staining and Western blotting, as well as increased Th mRNA levels in SN. MEN2B mice had increased number of DA neurons in SN by about 25% and they also exhibited increased sensitivity to the stimulatory effects of cocaine. We also developed a semi-throughput in vitro micro-island assay for the quantification of neuronal survival and TH levels by computer-assisted methodology from limited amounts of tissue. This assay can be applied for the initial screening for dopaminotrophic molecules, as well as chemical drug library screening. It is applicable to any neuronal system for the screening of neurotrophic molecules. Since our microarray experiments revealed possible GDNF-VEGF-C crosstalk we further concentrated on studying the neurotrophic effects of VEGF-C. We showed that VEGF-C acts as a neurotrophic molecule for the DA neurons both in vitro and in vivo, however without additive effect when used together with GDNF. The neuroprotective effect for VEGF-C in vivo in rat 6-OHDA model of PD was demonstrated. The possible signalling mechanisms of VEGF-C in the nervous system were investigated - infusion of VEGF-C to rat brain induced ERK activation, however no direct activation of RET signalling in vitro was found. VEGF-C treatment of rat striatum lead to up-regulation of VEGFR-1-3, indicating that VEGF-C can regulate the expression level of its own receptor. VEGF-C dopaminotrophic activity in vivo was further supported by increased vascular tissue in the neuroprotection experiments.

Computational interrogation of cis-regulatory elements of genes that are common targets of luteotropin and luteolysin in the primate corpus luteum

The rapid recent increase in microarray-based gene expression studies in the corpus luteum (CL) utilizing macaque models gathered increasing volume of data in publically accessible microarray expression databases. Examining gene pathways in different functional states of CL may help to understand the factors that control luteal function and hence human fertility. Co-regulation of genes in microarray experiments may imply common transcriptional regulation by sequence-specific DNA-binding transcriptional factors. We have computationally analyzed the transcription factor binding sites (TFBS) in a previously reported macaque luteal microarray gene set (n = 15) that are common targets of luteotropin (luteinizing hormone (LH) and human chorionic gonadotropin (hCG)) and luteolysin (prostaglandin (PG) F-2 alpha). This in silico approach can reveal transcriptional networks that control these important genes which are representative of the interplay between luteotropic and luteolytic factors in the control of luteal function. Our computational analyses revealed 6 matrix families whose binding sites are significantly over-represented in promoters of these genes. The roles of these factors are discussed, which might help to understand the transcriptional regulatory network in the control of luteal function. These factors might be promising experimental targets for investigation of human luteal insufficiency. (C) 2012 Elsevier B.V. All rights reserved.

Ferrochelatase is a conserved downstream target of the blue light-sensing White collar complex in fungi.

Light is a universal signal perceived by organisms, including fungi, in which light regulates common and unique biological processes depending on the species. Previous research has established that conserved proteins, originally called White collar 1 and 2 from the ascomycete Neurospora crassa, regulate UV/blue light sensing. Homologous proteins function in distant relatives of N. crassa, including the basidiomycetes and zygomycetes, which diverged as long as a billion years ago. Here we conducted microarray experiments on the basidiomycete fungus Cryptococcus neoformans to identify light-regulated genes. Surprisingly, only a single gene was induced by light above the commonly used twofold threshold. This gene, HEM15, is predicted to encode a ferrochelatase that catalyses the final step in haem biosynthesis from highly photoreactive porphyrins. The C. neoformans gene complements a Saccharomyces cerevisiae hem15Delta strain and is essential for viability, and the Hem15 protein localizes to mitochondria, three lines of evidence that the gene encodes ferrochelatase. Regulation of HEM15 by light suggests a mechanism by which bwc1/bwc2 mutants are photosensitive and exhibit reduced virulence. We show that ferrochelatase is also light-regulated in a white collar-dependent fashion in N. crassa and the zygomycete Phycomyces blakesleeanus, indicating that ferrochelatase is an ancient target of photoregulation in the fungal kingdom.

Effect of pooling samples on the efficiency of comparative studies using microarrays

Motivation: Many biomedical experiments are carried out by pooling individual biological samples. However, pooling samples can potentially hide biological variance and give false confidence concerning the data significance. In the context of microarray experiments for detecting differentially expressed genes, recent publications have addressed the problem of the efficiency of sample pooling, and some approximate formulas were provided for the power and sample size calculations. It is desirable to have exact formulas for these calculations and have the approximate results checked against the exact ones. We show that the difference between the approximate and the exact results can be large.