Bayesian survival analysis using gene expression

This thesis developed and applied Bayesian models for the analysis of survival data. The gene expression was considered as explanatory variables within the Bayesian survival model which can be considered the new contribution in the analysis of such data. The censoring factor that is inherent of survival data has also been addressed in terms of its impact on the fitting of a finite mixture of Weibull distribution with and without covariates. To investigate this, simulation study were carried out under several censoring percentages. Censoring percentage as high as 80% is acceptable here as the work involved high dimensional data. Lastly the Bayesian model averaging approach was developed to incorporate model uncertainty in the prediction of survival.

A flexible Bayesian model for describing temporal variability of N2O emissions from an Australian pasture

Soil-based emissions of nitrous oxide (N2O), a well-known greenhouse gas, have been associated with changes in soil water-filled pore space (WFPS) and soil temperature in many previous studies. However, it is acknowledged that the environment-N2O relationship is complex and still relatively poorly unknown. In this article, we employed a Bayesian model selection approach (Reversible jump Markov chain Monte Carlo) to develop a data-informed model of the relationship between daily N2O emissions and daily WFPS and soil temperature measurements between March 2007 and February 2009 from a soil under pasture in Queensland, Australia, taking seasonal factors and time-lagged effects into account. The model indicates a very strong relationship between a hybrid seasonal structure and daily N2O emission, with the latter substantially increased in summer. Given the other variables in the model, daily soil WFPS, lagged by a week, had a negative influence on daily N2O; there was evidence of a nonlinear positive relationship between daily soil WFPS and daily N2O emission; and daily soil temperature tended to have a linear positive relationship with daily N2O emission when daily soil temperature was above a threshold of approximately 19°C. We suggest that this flexible Bayesian modeling approach could facilitate greater understanding of the shape of the covariate-N2O flux relation and detection of effect thresholds in the natural temporal variation of environmental variables on N2O emission.

Gene expression studies in multiple sclerosis

Multiple sclerosis (MS) is a serious neurological disorder affecting young Caucasian individuals, usually with an age of onset at 18 to 40 years old. Females account for approximately 60x of MS cases and the manifestation and course of the disease is highly variable from patient to patient. The disorder is characterised by the development of plaques within the central nervous system (CNS). Many gene expression studies have been undertaken to look at the specific patterns of gene transcript levels in MS. Human tissues and experimental mice were used in these gene-profiling studies and a very valuable and interesting set of data has resulted from these various expression studies. In general, genes showing variable expression include mainly immunological and inflammatory genes, stress and antioxidant genes, as well as metabolic and central nervous system markers. Of particular interest are a number of genes localised to susceptible loci previously shown to be in linkage with MS. However due to the clinical complexity of the disease, the heterogeneity of the tissues used in expression studies, as well as the variable DNA chips/membranes used for the gene profiling, it is difficult to interpret the available information. Although this information is essential for the understanding of the pathogenesis of MS, it is difficult to decipher and define the gene pathways involved in the disorder. Experiments in gene expression profiling in MS have been numerous and lists of candidates are now available for analysis. Researchers have investigated gene expression in peripheral mononuclear white blood cells (PBMCs), in MS animal models Experimental Allergic Encephalomyelitis (EAE) and post mortem MS brain tissues. This review will focus on the results of these studies.

In vitro and in vivo MMP gene expression localisation by In Situ-RT-PCR in cell culture and paraffin embedded human breast cancer cell line xenografts

Background Members of the matrix metalloproteinase (MMP) family of proteases are required for the degradation of the basement membrane and extracellular matrix in both normal and pathological conditions. In vitro, MT1-MMP (MMP-14, membrane type-1-MMP) expression is higher in more invasive human breast cancer (HBC) cell lines, whilst in vivo its expression has been associated with the stroma surrounding breast tumours. MMP-1 (interstitial collagenase) has been associated with MDA-MB-231 invasion in vitro, while MMP-3 (stromelysin-1) has been localised around invasive cells of breast tumours in vivo. As MMPs are not stored intracellularly, the ability to localise their expression to their cells of origin is difficult. Methods We utilised the unique in situ-reverse transcription-polymerase chain reaction (IS-RT-PCR) methodology to localise the in vitro and in vivo gene expression of MT1-MMP, MMP-1 and MMP-3 in human breast cancer. In vitro, MMP induction was examined in the MDA-MB-231 and MCF-7 HBC cell lines following exposure to Concanavalin A (Con A). In vivo, we examined their expression in archival paraffin embedded xenografts derived from a range of HBC cell lines of varied invasive and metastatic potential. Mouse xenografts are heterogenous, containing neoplastic human parenchyma with mouse stroma and vasculature and provide a reproducible in vivo model system correlated to the human disease state. Results In vitro, exposure to Con A increased MT1-MMP gene expression in MDA-MB-231 cells and decreased MT1-MMP gene expression in MCF-7 cells. MMP-1 and MMP-3 gene expression remained unchanged in both cell lines. In vivo, stromal cells recruited into each xenograft demonstrated differences in localised levels of MMP gene expression. Specifically, MDA-MB-231, MDA-MB-435 and Hs578T HBC cell lines are able to influence MMP gene expression in the surrounding stroma. Conclusion We have demonstrated the applicability and sensitivity of IS-RT-PCR for the examination of MMP gene expression both in vitro and in vivo. Induction of MMP gene expression in both the epithelial tumour cells and surrounding stromal cells is associated with increased metastatic potential. Our data demonstrate the contribution of the stroma to epithelial MMP gene expression, and highlight the complexity of the role of MMPs in the stromal-epithelial interactions within breast carcinoma.

Age-related changes in cardiac adenosine receptor expression

Adenosine is an important cardioprotective agent that works via several adenosine receptor (ADOR) subtypes to regulate cardiovascular activity. It is well established that functional responses to adenosine decline with age. What is unclear, though, is whether these changes occur at the receptor, second messenger or translational level. In this study we determined the effect of age on cardiac adenosine receptor expression using the housekeeping gene 18S rRNA versus the adenosine A2B receptor gene as internal controls. Absolute quantification showed that no age-related changes occurred in the expression of 18S rRNA or adenosine A2B receptor internal control genes. Subsequently, relative analysis of the adenosine receptor subtypes using 18S rRNA found a significant age-related reduction in the expression of the adenosine A1 receptor (5.5-fold), with no changes in the expression of the adenosine A2A, A2B and A3 receptors. When using the expression of the adenosine A2B receptor as the internal control gene, a significant down regulation of both the adenosine A1 (5.4-fold) and A2A (2.2-fold) receptors with no change in the expression of adenosine A3 receptor was found. Therefore, the high level of expression of the 18S rRNA housekeeping gene was found to mask a significant change in expression of the adenosine A2A receptor with age. Ultimately, these findings show an age-related reduction in adenosine A1 and A2A receptor expression in rat heart.

Stimulation of MMP-11 (stromelysin-3) expression in mouse fibroblasts by cytokines, collagen and co-culture with human breast cancer cell lines

Background Matrix metalloproteinases (MMPs) are central to degradation of the extracellular matrix and basement membrane during both normal and carcinogenic tissue remodeling. MT1-MMP (MMP-14) and stromelysin-3 (MMP-11) are two members of the MMP family of proteolytic enzymes that have been specifically implicated in breast cancer progression. Expressed in stromal fibroblasts adjacent to epithelial tumour cells, the mechanism of MT1-MMP and MMP-11 induction remains unknown. Methods To investigate possible mechanisms of induction, we examined the effects of a number of plausible regulatory agents and treatments that may physiologically influence MMP expression during tumour progression. Thus NIH3T3 and primary mouse embryonic fibroblasts (MEFs) were: a) treated with the cytokines IL-1β, IL-2, IL-6, IL-8 and TGF-β for 3, 6, 12, 24, and 48 hours; b) grown on collagens I, IV and V; c) treated with fibronectin, con-A and matrigel; and d) co-cultured with a range of HBC (human breast cancer) cell lines of varied invasive and metastatic potential. Results Competitive quantitative RT-PCR indicated that MMP-11 expression was stimulated to a level greater than 100%, by 48 hour treatments of IL-1β, IL-2, TGF-β, fibronectin and collagen V. No other substantial changes in expression of MMP-11 or MT1-MMP in either tested fibroblast culture, under any treatment conditions, were observed. Conclusion We have demonstrated significant MMP-11 stimulation in mouse fibroblasts using cytokines, matrix constituents and HBC cell lines, and also some inhibition of MT1-MMP. Our data suggest that the regulation of these genes in the complex stromal-epithelial interactions that occur in human breast carcinoma, is influenced by several mechanisms.

An assessment of MMP and TIMP gene expression in cell lines and stroma : tumour differences in microdissected breast cancer biopsies

To examine matrix metalloproteinase (MMP) and tissue inhibitor of metalloproteinases (TIMP) mRNA levels in archival breast cancer biopsies, we employed microdissection to separate tumour tissue from the surrounding breast tissue, or stroma and RT-PCR to determine gross qualitative and small quantitative differences in the patterns of expression. In this study, a significant correlation (p < 0.05, by Mann-Whitney U analysis) between TIMP-2 expression and lymph node involvement was identified, while MMP-11 and TIMP-1 expression patterning also significantly (p < 0.05) differed between those tumours showing calcification and those that did not. When compared by Spearmans’ ρ correlation analysis, a significant association (p < 0.05, ρ = 0.404) was identified in the pattern of MMP-2 and MMP-9 gene expression. In this study, the use of microdissection and a systematic strategy of RT-PCR analysis have allowed us to investigate localized MMP and MMP inhibitor expression within breast tumours. We have identified patterns of gene expression that may further reveal aspects of breast carcinogenesis, and a robust method for examining changes in clinically important genes using archival biopsies and across stroma-tumour boundaries.

Quantitative and qualitative changes in gene expression patterns characterize the activity of plaques in multiple sclerosis

Multiple sclerosis (MS) is a complex autoimmune disorder of the CNS with both genetic and environmental contributing factors. Clinical symptoms are broadly characterized by initial onset, and progressive debilitating neurological impairment. In this study, RNA from MS chronic active and MS acute lesions was extracted, and compared with patient matched normal white matter by fluorescent cDNA microarray hybridization analysis. This resulted in the identification of 139 genes that were differentially regulated in MS plaque tissue compared to normal tissue. Of these, 69 genes showed a common pattern of expression in the chronic active and acute plaque tissues investigated (Pvalue<0.0001, ρ=0.73, by Spearman's ρ analysis); while 70 transcripts were uniquely differentially expressed (≥1.5-fold) in either acute or chronic active tissues. These results included known markers of MS such as the myelin basic protein (MBP) and glutathione S-transferase (GST) M1, nerve growth factors, such as nerve injury-induced protein 1 (NINJ1), X-ray and excision DNA repair factors (XRCC9 and ERCC5) and X-linked genes such as the ribosomal protein, RPS4X. Primers were then designed for seven array-selected genes, including transferrin (TF), superoxide dismutase 1 (SOD1), glutathione peroxidase 1 (GPX1), GSTP1, crystallin, alpha-B (CRYAB), phosphomannomutase 1 (PMM1) and tubulin β-5 (TBB5), and real time quantitative (Q)-PCR analysis was performed. The results of comparative Q-PCR analysis correlated significantly with those obtained by array analysis (r=0.75, Pvalue<0.01, by Pearson's bivariate correlation). Both chronic active and acute plaques shared the majority of factors identified suggesting that quantitative, rather than gross qualitative differences in gene expression pattern may define the progression from acute to chronic active plaques in MS.

The measurement of adenosine and estrogen receptor expression in rat brains following ovariectomy using quantitative PCR analysis

In our laboratory we have developed a quantitative-polymerase chain reaction (Q-PCR) strategy to examine the differential expression of adenosine receptor (ADOR), A(1), A(2A), A(2B) and A(3), and estrogen receptors (ER) alpha and beta. Brain and uterine mRNA were first used to optimise specific amplification conditions prior to SYBR Green I real time analysis of receptor subtype expression. SYBR Green I provided a convenient and sensitive means of examining specific PCR amplification product in real time, and allowed the generation of standard curves from which relative receptor abundance could be determined. Real time Q-PCR analysis was then performed, to examine changes in receptor expression levels in brains of adult female Wistar rats 3-month post ovariectomy. Comparison with sham-operated age-matched control rats demonstrated both comparative and absolute-copy number changes in receptor levels. Evaluation of both analytical methods investigated 18S rRNA as an internal reference for comparative gene expression analysis in the brain. The results of this study revealed preferential repression of ADORA(2A) (>4-fold down) and consistent (>2-fold) down-regulation of ADORA(1), ADORA(3), and ER-beta, following ovariectomy. No change was found in ADORA(2B) or ER-alpha. Analysis of absolute copy number in this study revealed a correlation between receptor expression in response to ovariectomy, and relative receptor subtype abundance in the brain.

Differential gene expression in breast cancer cell lines and stroma-tumor differences in microdissected breast cancer biopsies revealed by display array analysis

To examine gene-expression patterning in late-stage breast cancer biopsies, we used a microdissection technique to separate tumor from the surrounding breast tissue or stroma. A DD-PCR protocol was then used to amplify expressed products, which were resolved using PAGE and used as probe to hybridize with representative human arrays and cDNA libraries. The probe derived from the tumor–stroma comparison was hybridized with a gene array and an arrayed cDNA library derived from a GCT of bone; 21 known genes or expressed sequence tags were detected, of which 17 showed differential expression. These included factors associated with epithelial to mesenchymal transition (vimentin), the cargo selection protein (TIP47) and the signal transducer and activator of transcription (STAT3). Northern blot analysis was used to confirm those genes also expressed by representative breast cancer cell lines. Notably, 6 genes of unknown function were restricted to tumor while the majority of stroma-associated genes were known. When applied to transformed breast cancer cell lines (MDA-MB-435 and T47D) that are known to have different metastatic potential, DD array analysis revealed a further 20 genes; 17 of these genes showed differential expression. Use of microdissection and the DD-PCR array protocol allowed us to identify factors whose localized expression within the breast may play a role in abnormal breast development or breast carcinogenesis.

The expression of HAT and HDAC enzymes in asthma airways

Asthma is chronic inflammatory disease of the lower airways that is both, genetically inherited and environmentally influenced. This project investigated how molecular mechanisms known to be influenced both genetically and environmentally, contribute to the onset of asthma.

Transcription and expression of Plasmodium falciparum histidine-rich proteins in different stages and strains : implications for rapid diagnostic tests

Background: Although rapid diagnostic tests (RDTs) for Plasmodium falciparum infection that target histidine rich protein 2 (PfHRP2) are generally sensitive, their performance has been reported to be variable. One possible explanation for variable test performance is differences in expression level of PfHRP in different parasite isolates. Methods: Total RNA and protein were extracted from synchronised cultures of 7 P. falciparum lines over 5 time points of the life cycle, and from synchronised ring stages of 10 falciparum lines. Using quantitative real-time polymerase chain reaction, Western blot analysis and ELISA we investigated variations in the transcription and protein levels of pfhrp2, pfhrp3 and PfHRP respectively in the different parasite lines, over the parasite intraerythrocytic life cycle. Results: Transcription of pfhrp2 and pfhrp3 in different parasite lines over the parasite life cycle was observed to vary relative to the control parasite K1. In some parasite lines very low transcription of these genes was observed. The peak transcription was observed in ring-stage parasites. Pfhrp2 transcription was observed to be consistently higher than pfhrp3 transcription within parasite lines. The intraerythrocytic lifecycle stage at which the peak level of protein was present varied across strains. Total protein levels were more constant relative to total mRNA transcription, however a maximum 24 fold difference in expression at ring-stage parasites relative to the K1 strain was observed. Conclusions: The levels of transcription of pfhrp2 and pfhrp3, and protein expression of PfHRP varied between different P. falciparum strains. This variation may impact on the detection sensitivity of PfHRP2-detecting RDTs.

Mining temporal and spatial travel regularities for transit planning

Smart Card data from Automated Fare Collection system has been considered as a promising source of information for transit planning. However, literature has been limited to mining travel patterns from transit users and suggesting the potential of using this information. This paper proposes a method for mining spatial regular origins-destinations and temporal habitual travelling time from transit users. These travel regularity are discussed as being useful for transit planning. After reconstructing the travel itineraries, three levels of Density-Based Spatial Clustering of Application with Noise (DBSCAN) have been utilised to retrieve travel regularity of each of each frequent transit users. Analyses of passenger classifications and personal travel time variability estimation are performed as the examples of using travel regularity in transit planning. The methodology introduced in this paper is of interest for transit authorities in planning and managements

Influence of nitrogen fertiliser application and timing on greenhouse gas emissions from a lychee (Litchi chinensis) orchard in humid subtropical Australia

Nitrous oxide emissions from intensive, fertilised agricultural systems have been identified as significant contributors to both Australia's and the global greenhouse gas (GHG) budget. This is expected to increase as rates of agriculture intensification and land use change accelerate to support population growth and food production. Limited data exists on N2O trace gas fluxes from subtropical or tropical tree cropping soils critical for the development of effective mitigation strategies.This study aimed to quantify GHG emissions over two consecutive years (March 2007 to March 2009) from a 30 year (lychee) orchard in the humid subtropical region of Australia. GHG fluxes were measured using a combination of high temporal resolution automated sampling and manually sampled chambers. No fertiliser was added to the plots during the 2007 measurement season. A split application of nitrogen fertiliser (urea) was added at the rate of 265kgNha-1 during the autumn and spring of 2008. Emissions of N2O were influenced by rainfall events and seasonal temperatures during 2007 and the fertilisation events in 2008. Annual N2O emissions from the lychee canopy increased from 1.7kgN2O-Nha-1yr-1 for 2007, to 7.6kgN2O-Nha-1yr-1 following fertiliser application in 2008. This represented an emission factor of 1.56%, corrected for background emissions. The timing of the split application was found to be critical to N2O emissions, with over twice as much lost following an application in spring (2.44%) compared to autumn (EF: 1.10%). This research suggests that avoiding fertiliser application during the hot and moist spring/summer period can reduce N2O losses without compromising yields.

High resolution melt analysis : a novel method for studying the genetic relatedness of Pseudomonas aeruginosa isolates from clinical and environmental sources

A novel method was developed for studying the genetic relatedness of Pseudomonas aeruginosa isolates from clinical and environmental sources. This bacterium is ubiquitous in the natural environment and is an important pathogen known to infect Cystic Fibrosis (CF) patients. The transmission route of strains has not yet been defined; current theories include acquisition from an environmental source or through patient-to-patient spread. A highly discriminatory, bioinformatics based, DNA typing method was developed to investigate the relatedness of clinical and environmental isolates. This study found a similarity between the environmental and several CF clonal strains and also highlighted occurrence of environmental P. aeruginosa strains in CF infections.