The association of single nucleotide polymorphisms with endometrial cancer risk and prognosis

Endometrial cancer is one of the most common female diseases in developed nations and is the most commonly diagnosed gynaecological cancer in Australia. The disease is commonly classified by histology: endometrioid or non-endometrioid endometrial cancer. While non-endometrioid endometrial cancers are accepted to be high-grade, aggressive cancers, endometrioid cancers (comprising 80% of all endometrial cancers diagnosed) generally carry a favourable patient prognosis. However, endometrioid endometrial cancer patients endure significant morbidity due to surgery and radiotherapy used for disease treatment, and patients with recurrent disease have a 5-year survival rate of less than 50%. Genetic analysis of women with endometrial cancer could uncover novel markers associated with disease risk and/or prognosis, which could then be used to identify women at high risk and for the use of specialised treatments. Proteases are widely accepted to play an important role in the development and progression of cancer. This PhD project hypothesised that SNPs from two protease gene families, the matrix metalloproteases (MMPs, including their tissue inhibitors, TIMPs) and the tissue kallikrein-related peptidases (KLKs) would be associated with endometrial cancer susceptibility and/or prognosis. In the first part of this study, optimisation of the genotyping techniques was performed. Results from previously published endometrial cancer genetic association studies were attempted to be validated in a large, multicentre replication set (maximum cases n = 2,888, controls n = 4,483, 3 studies). The rs11224561 progesterone receptor SNP (PGR, A/G) was observed to be associated with increased endometrial cancer risk (per A allele OR 1.31, 95% CI 1.12-1.53; p-trend = 0.001), a result which was initially reported among a Chinese sample set. Previously reported associations for the remaining 8 SNPs investigated for this section of the PhD study were not confirmed, thereby reinforcing the importance of validation of genetic association studies. To examine the effect of SNPs from the MMP and KLK families on endometrial cancer risk, we selected the most significantly associated MMP and KLK SNPs from genome-wide association study analysis (GWAS) to be genotyped in the GWAS replication set (cases n = 4,725, controls n = 9,803, 13 studies). The significance of the MMP24 rs932562 SNP was unchanged after incorporation of the stage 2 samples (Stage 1 per allele OR 1.18, p = 0.002; Combined Stage 1 and 2 OR 1.09, p = 0.002). The rs10426 SNP, located 3' to KLK10 was predicted by bioinformatic analysis to effect miRNA binding. This SNP was observed in the GWAS stage 1 result to exhibit a recessive effect on endometrial cancer risk, a result which was not validated in the stage 2 sample set (Stage 1 OR 1.44, p = 0.007; Combined Stage 1 and 2 OR 1.14, p = 0.08). Investigation of the regions imputed surrounding the MMP, TIMP and KLK genes did not reveal any significant targets for further analysis. Analysis of the case data from the endometrial cancer GWAS to identify genetic variation associated with cancer grade did not reveal SNPs from the MMP, TIMP or KLK genes to be statistically significant. However, the representation of SNPs from the MMP, TIMP and KLK families by the GWAS genotyping platform used in this PhD project was examined and observed to be very low, with the genetic variation of four genes (MMP23A, MMP23B, MMP28 and TIMP1) not captured at all by this technique. This suggests that comprehensive candidate gene association studies will be required to assess the role of SNPs from these genes with endometrial cancer risk and prognosis. Meta-analysis of gene expression microarray datasets curated as part of this PhD study identified a number of MMP, TIMP and KLK genes to display differential expression by endometrial cancer status (MMP2, MMP10, MMP11, MMP13, MMP19, MMP25 and KLK1) and histology (MMP2, MMP11, MMP12, MMP26, MMP28, TIMP2, TIMP3, KLK6, KLK7, KLK11 and KLK12). In light of these findings these genes should be prioritised for future targeted genetic association studies. Two SNPs located 43.5 Mb apart on chromosome 15 were observed from the GWAS analysis to be associated with increased endometrial cancer grade, results that were validated in silico in two independent datasets. One of these SNPs, rs8035725 is located in the 5' untranslated region of a MYC promoter binding protein DENND4A (Stage 1 OR 1.15, p = 9.85 x 10P -5 P, combined Stage 1 and in silico validation OR 1.13, p = 5.24 x 10P -6 P). This SNP has previously been reported to alter the expression of PTPLAD1, a gene involved in the synthesis of very long fatty acid chains and in the Rac1 signaling pathway. Meta-analysis of gene expression microarray data found PTPLAD1 to display increased expression in the aggressive non-endometrioid histology compared with endometrioid endometrial cancer, suggesting that the causal SNP underlying the observed genetic association may influence expression of this gene. Neither rs8035725 nor significant SNPs identified by imputation were predicted bioinformatically to affect transcription factor binding sites, indicating that further studies are required to assess their potential effect on other regulatory elements. The other grade- associated SNP, rs6606792, is located upstream of an inferred pseudogene, ELMO2P1 (Stage 1 OR 1.12, p = 5 x 10P -5 P; combined Stage 1 and in silico validation OR 1.09, p = 3.56 x 10P -5 P). Imputation of the ±1 Mb region surrounding this SNP revealed a cluster of significantly associated variants which are predicted to abolish various transcription factor binding sites, and would be expected to decrease gene expression. ELMO2P1 was not included on the microarray platforms collected for this PhD, and so its expression could not be investigated. However, the high sequence homology of ELMO2P1 with ELMO2, a gene important to cell motility, indicates that ELMO2 could be the parent gene for ELMO2P1 and as such, ELMO2P1 could function to regulate the expression of ELMO2. Increased expression of ELMO2 was seen to be associated with increasing endometrial cancer grade, as well as with aggressive endometrial cancer histological subtypes by microarray meta-analysis. Thus, it is hypothesised that SNPs in linkage disequilibrium with rs6606792 decrease the transcription of ELMO2P1, reducing the regulatory effect of ELMO2P1 on ELMO2 expression. Consequently, ELMO2 expression is increased, cell motility is enhanced leading to an aggressive endometrial cancer phenotype. In summary, these findings have identified several areas of research for further study. The results presented in this thesis provide evidence that a SNP in PGR is associated with risk of developing endometrial cancer. This PhD study also reports two independent loci on chromosome 15 to be associated with increased endometrial cancer grade, and furthermore, genes associated with these SNPs to be differentially expressed according in aggressive subtypes and/or by grade. The studies reported in this thesis support the need for comprehensive SNP association studies on prioritised MMP, TIMP and KLK genes in large sample sets. Until these studies are performed, the role of MMP, TIMP and KLK genetic variation remains unclear. Overall, this PhD study has contributed to the understanding of genetic variation involvement in endometrial cancer susceptibility and prognosis. Importantly, the genetic regions highlighted in this study could lead to the identification of novel gene targets to better understand the biology of endometrial cancer and also aid in the development of therapeutics directed at treating this disease.

A novel saliva-based microRNA biomarker panel to detect head and neck cancers

Background MicroRNAs (miRNAs) are known to play an important role in cancer development by post-transcriptionally affecting the expression of critical genes. The aims of this study were two-fold: (i) to develop a robust method to isolate miRNAs from small volumes of saliva and (ii) to develop a panel of saliva-based diagnostic biomarkers for the detection of head and neck squamous cell carcinoma (HNSCC). Methods Five differentially expressed miRNAs were selected from miScript™ miRNA microarray data generated using saliva from five HNSCC patients and five healthy controls. Their differential expression was subsequently confirmed by RT-qPCR using saliva samples from healthy controls (n = 56) and HNSCC patients (n = 56). These samples were divided into two different cohorts, i.e., a first confirmatory cohort (n = 21) and a second independent validation cohort (n = 35), to narrow down the miRNA diagnostic panel to three miRNAs: miR-9, miR-134 and miR-191. This diagnostic panel was independently validated using HNSCC miRNA expression data from The Cancer Genome Atlas (TCGA), encompassing 334 tumours and 39 adjacent normal tissues. Receiver operating characteristic (ROC) curve analysis was performed to assess the diagnostic capacity of the panel. Results On average 60 ng/μL miRNA was isolated from 200 μL of saliva. Overall a good correlation was observed between the microarray data and the RT-qPCR data. We found that miR-9 (P <0.0001), miR-134 (P <0.0001) and miR-191 (P <0.001) were differentially expressed between saliva from HNSCC patients and healthy controls, and that these miRNAs provided a good discriminative capacity with area under the curve (AUC) values of 0.85 (P <0.0001), 0.74 (P < 0.001) and 0.98 (P < 0.0001), respectively. In addition, we found that the salivary miRNA data showed a good correlation with the TCGA miRNA data, thereby providing an independent validation. Conclusions We show that we have developed a reliable method to isolate miRNAs from small volumes of saliva, and that the saliva-derived miRNAs miR-9, miR-134 and miR-191 may serve as novel biomarkers to reliably detect HNSCC. © 2014 International Society for Cellular Oncology.

The role of DNA repair pathways in cisplatin resistant lung cancer

Platinum chemotherapeutic agents such as cisplatin are currently used in the treatment of various malignancies such as lung cancer. However, their efficacy is significantly hindered by the development of resistance during treatment. While a number of factors have been reported that contribute to the onset of this resistance phenotype, alterations in the DNA repair capacity of damaged cells is now recognised as an important factor in mediating this phenomenon. The mode of action of cisplatin has been linked to its ability to crosslink purine bases on the DNA, thereby interfering with DNA repair mechanisms and inducing DNA damage. Following DNA damage, cells respond by activating a DNA-damage response that either leads to repair of the lesion by the cell thereby promoting resistance to the drug, or cell death via activation of the apoptotic response. Therefore, DNA repair is a vital target to improving cancer therapy and reduce the resistance of tumour cells to DNA damaging agents currently used in the treatment of cancer patients. To date, despite the numerous findings that differential expression of components of the various DNA repair pathways correlate with response to cisplatin, translation of such findings in the clinical setting are still warranted. The identification of alterations in specific proteins and pathways that contribute to these unique DNA repair pathways in cisplatin resistant cancer cells may potentially lead to a renewed interest in the development of rational novel therapies for cisplatin resistant cancers, in particular, lung cancer.

RNASeqBrowser: A genome browser for simultaneous visualization of raw strand specific RNAseq reads and UCSC genome browser custom tracks

Background Strand specific RNAseq data is now more common in RNAseq projects. Visualizing RNAseq data has become an important matter in Analysis of sequencing data. The most widely used visualization tool is the UCSC genome browser that introduced the custom track concept that enabled researchers to simultaneously visualize gene expression at a particular locus from multiple experiments. Our objective of the software tool is to provide friendly interface for visualization of RNAseq datasets. Results This paper introduces a visualization tool (RNASeqBrowser) that incorporates and extends the functionality of the UCSC genome browser. For example, RNASeqBrowser simultaneously displays read coverage, SNPs, InDels and raw read tracks with other BED and wiggle tracks -- all being dynamically built from the BAM file. Paired reads are also connected in the browser to enable easier identification of novel exon/intron borders and chimaeric transcripts. Strand specific RNAseq data is also supported by RNASeqBrowser that displays reads above (positive strand transcript) or below (negative strand transcripts) a central line. Finally, RNASeqBrowser was designed for ease of use for users with few bioinformatic skills, and incorporates the features of many genome browsers into one platform. Conclusions The features of RNASeqBrowser: (1) RNASeqBrowser integrates UCSC genome browser and NGS visualization tools such as IGV. It extends the functionality of the UCSC genome browser by adding several new types of tracks to show NGS data such as individual raw reads, SNPs and InDels. (2) RNASeqBrowser can dynamically generate RNA secondary structure. It is useful for identifying non-coding RNA such as miRNA. (3) Overlaying NGS wiggle data is helpful in displaying differential expression and is simple to implement in RNASeqBrowser. (4) NGS data accumulates a lot of raw reads. Thus, RNASeqBrowser collapses exact duplicate reads to reduce visualization space. Normal PC’s can show many windows of NGS individual raw reads without much delay. (5) Multiple popup windows of individual raw reads provide users with more viewing space. This avoids existing approaches (such as IGV) which squeeze all raw reads into one window. This will be helpful for visualizing multiple datasets simultaneously. RNASeqBrowser and its manual are freely available at http://www.australianprostatecentre.org/research/software/rnaseqbrowser webcite or http://sourceforge.net/projects/rnaseqbrowser/ webcite

Proteogenomics of selective susceptibility to endotoxin using circulating acute phase biomarkers and bioassay development in sheep: a review

Scientists have injected endotoxin into animals to investigate and understand various pathologies and novel therapies for several decades. Recent observations have shown that there is selective susceptibility to Escherichia coli lipopolysaccharide (LPS) endotoxin in sheep, despite having similar breed characteristics. The reason behind this difference is unknown, and has prompted studies aiming to explain the variation by proteogenomic characterisation of circulating acute phase biomarkers. It is hypothesised that genetic trait, biochemical, immunological and inflammation marker patterns contribute in defining and predicting mammalian response to LPS. This review discusses the effects of endotoxin and host responses, genetic basis of innate defences, activation of the acute phase response (APR) following experimental LPS challenge, and the current approaches employed in detecting novel biomarkers including acute phase proteins (APP) and micro-ribonucleic acids (miRNAs) in serum or plasma. miRNAs are novel targets for elucidating molecular mechanisms of disease because of their differential expression during pathological, and in healthy states. Changes in miRNA profiles during a disease challenge may be reflected in plasma. Studies show that gel-based two-dimensional electrophoresis (2-DE) coupled with either matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF MS) or liquid chromatography-mass spectrometry (LC-MS/MS) are currently the most used methods for proteome characterisation. Further evidence suggests that proteomic investigations are preferentially shifting from 2-DE to non-gel based LC-MS/MS coupled with data extraction by sequential window acquisition of all theoretical fragment-ion spectra (SWATH) approaches that are able to identify a wider range of proteins. Enzyme-linked immunosorbent assay (ELISA), quantitative real-time polymerase chain reaction (qRT-PCR), and most recently proteomic methods have been used to quantify low abundance proteins such as cytokines. qRT-PCR and next generation sequencing (NGS) are used for the characterisation of miRNA. Proteogenomic approaches for detecting APP and novel miRNA profiling are essential in understanding the selective resistance to endotoxin in sheep. The results of these methods could help in understanding similar pathology in humans. It might also be helpful in the development of physiological and diagnostic screening assays for determining experimental inclusion and endpoints, and in clinical trials in future

Genome-wide association study for sight-threatening diabetic retinopathy reveals association with genetic variation near the GRB2 gene

Aims/hypothesis Diabetic retinopathy is a serious complication of diabetes mellitus and can lead to blindness. A genetic component, in addition to traditional risk factors, has been well described although strong genetic factors have not yet been identified. Here, we aimed to identify novel genetic risk factors for sight-threatening diabetic retinopathy using a genome-wide association study. Methods Retinopathy was assessed in white Australians with type 2 diabetes mellitus. Genome-wide association analysis was conducted for comparison of cases of sight-threatening diabetic retinopathy (n = 336) with diabetic controls with no retinopathy (n = 508). Top ranking single nucleotide polymorphisms were typed in a type 2 diabetes replication cohort, a type 1 diabetes cohort and an Indian type 2 cohort. A mouse model of proliferative retinopathy was used to assess differential expression of the nearby candidate gene GRB2 by immunohistochemistry and quantitative western blot. Results The top ranked variant was rs3805931 with p = 2.66 × 10−7, but no association was found in the replication cohort. Only rs9896052 (p = 6.55 × 10−5) was associated with sight-threatening diabetic retinopathy in both the type 2 (p = 0.035) and the type 1 (p = 0.041) replication cohorts, as well as in the Indian cohort (p = 0.016). The study-wide meta-analysis reached genome-wide significance (p = 4.15 × 10−8). The GRB2 gene is located downstream of this variant and a mouse model of retinopathy showed increased GRB2 expression in the retina. Conclusions/interpretation Genetic variation near GRB2 on chromosome 17q25.1 is associated with sight-threatening diabetic retinopathy. Several genes in this region are promising candidates and in particular GRB2 is upregulated during retinal stress and neovascularisation.

Representations and properties of para-Bose oscillator operators. I. Energy position and momentum eigenstates

Para-Bose commutation relations are related to the SL(2,R) Lie algebra. The irreducible representation [script D]alpha of the para-Bose system is obtained as the direct sum Dbeta[direct-sum]Dbeta+1/2 of the representations of the SL(2,R) Lie algebra. The position and momentum eigenstates are then obtained in this representation [script D]alpha, using the matrix mechanical method. The orthogonality, completeness, and the overlap of these eigenstates are derived. The momentum eigenstates are also derived using the wave mechanical method by specifying the domain of the definition of the momentum operator in addition to giving it a formal differential expression. By a careful consideration in this manner we find that the two apparently different solutions obtained by Ohnuki and Kamefuchi in this context are actually unitarily equivalent. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Orphan nuclear receptor subfamily NR4A their interplay with other nuclear receptors and functions in osteoblasts

Nurr1, NGFI-B and Nor1 (NR4A2, NR4A1 and NR4A3, respectively) belong to the NR4A subfamily of nuclear receptors. The NR4A receptors are orphan nuclear receptors which means that activating or repressing ligands for these receptors have not been found. NR4A expression is rapidly induced in response to various stimuli including growth factors and the parathyroid hormone (PTH). The studies concerning the NR4A receptors in the central nervous system have demonstrated that they have a major role in the development and function of the dopaminergic neurons of the midbrain and in regulating hypothalamus-pituitary-adrenal-axis. However, the peripheral functions of the NR4A family are largely unknown. Cultured mouse primary osteoblasts, a preosteoblastic cell line and several osteoblastic cell lines were used to investigate the role of NR4A receptors in osteoblasts. NR4A receptors were shown to directly bind to and activate the promoter of the osteopontin gene (OPN) in osteoblastic cells, thus regulating its expression. OPN is a major bone matrix protein expressed throughout the differentiation of preosteoblastic cells into osteoblasts. The activation of the OPN promoter was shown to be dependent on the activation function-1 located in the N-terminal part of Nurr1 and to occur in both monomeric and RXR heterodimeric forms of NR4A receptors. Furthermore, PTH was shown to upregulate OPN expression through the NR4A family. It was also demonstrated that the fibroblast growth factor-8b (FGF-8b) induces the expression of NR4A receptors in osteoblasts as immediate early genes. This induction involved phosphatidylinositol-3 kinase, protein kinase C, and mitogen activated protein kinase, which are all major pathways of FGF signalling. Nurr1 and NGFI-B were shown to induce the proliferation of preosteoblastic cells and to reduce their apoptosis. FGF-8b was shown to stimulate the proliferation of osteoblastic cells through the NR4A receptors. These results suggest that NR4A receptors have a role both in the differentiation of osteoblasts and in the proliferation and apoptosis of preosteoblast. The NR4A receptors were found to bind to the same response element on OPN as the members of the NR3B family of orphan receptors do. Mutual repression was observed between the NR4A receptors and the NR3B receptors. This repression was shown to be dependent on the DNA-binding domains of both receptor families, but to result neither from the competition of DNA binding nor from the competition for coactivators. As the repression was dependent on the relative expression levels of the NR4As and NR3Bs, it seems likely that the ratio of the receptors mediates their activity on their response elements. Rapid induction of the NR4As in response to various stimuli and differential expression of the NR3Bs can effectively control the gene activation by the NR4A receptors. NR4A receptors can bind DNA as monomers, and Nurr1 and NGFI-B can form permissive heterodimers with the retinoid X receptor (RXR). Permissive heterodimers can be activated with RXR agonists, unlike non-permissive heterodimers, which are formed by RXR and retinoic acid receptor or thyroid hormone receptor (RAR and TR, respectively). Non-permissive heterodimers can only be activated by the agonists of the heterodimerizing partner. The mechanisms behind differential response to RXR agonists have remained unresolved. As there are no activating or repressing ligands for the NR4A receptors, it would be important to find out, how they are regulated. Permissiviness of Nurr1/RXR heterodimers was linked to the N-terminal part of Nurr1 ligand-binding domain. This region has previously been shown to mediate the interaction between NRs and corepressors. Non-permissive RAR and TR, permissive Nurr1 and NGFI-B, and RXR were overexpressed with corepressors silencing mediator for retinoic acid and thyroid hormone receptors (SMRT), and with nuclear receptor corepressor in several cell lines. Nurr1 and NGFI-B were found to be repressed by SMRT. The interaction of RXR heterodimers with corepressors was weak in permissive heterodimers and much stronger in non-permissive heterodimers. Non-permissive heterodimers also released corepressors only in response to the agonist of the heterodimeric partner of RXR. In the permissive Nurr1/RXR heterodimer, however, SMRT was released following the treatment with RXR agonists. Corepressor release in response to ligands was found to differentiate permissive heterodimers from non-permissive ones. Corepressors were thus connected to the regulation of NR4A functions. In summary, the studies presented here linked the NR4A family of orphan nuclear receptors to the regulation of osteoblasts. Nurr1 and NGFI-B were found to control the proliferation and apoptosis of preosteoblasts. The studies also demonstrated that cross-talk with the NR3B receptors controls the activity of these orphan receptors. The results clarified the mechanism of permissiviness of RXR-heterodimers. New information was obtained on the regulation and functions of NR4A receptors, for which the ligands are unknown.

Functional evaluation of genetic variants associated with endometriosis near GREB1

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.

Latent class and genetic analysis does not support migraine with aura and migraine without aura as separate entities

Latent class and genetic analyses were used to identify subgroups of migraine sufferers in a community sample of 6,265 Australian twins (55% female) aged 25-36 who had completed an interview based on International Headache Society (IHS) criteria. Consistent with prevalence rates from other population-based studies, 703 (20%) female and 250 (9%) male twins satisfied the IHS criteria for migraine without aura (MO), and of these, 432 (13%) female and 166 (6%) male twins satisfied the criteria for migraine with aura (MA) as indicated by visual symptoms. Latent class analysis (LCA) of IHS symptoms identified three major symptomatic classes, representing 1) a mild form of recurrent nonmigrainous headache, 2) a moderately severe form of migraine, typically without visual aura symptoms (although 40% of individuals in this class were positive for aura), and 3) a severe form of migraine typically with visual aura symptoms (although 24% of individuals were negative for aura). Using the LCA classification, many more individuals were considered affected to some degree than when using IHS criteria (35% vs. 13%). Furthermore, genetic model fitting indicated a greater genetic contribution to migraine using the LCA classification (heritability, h(2)=0.40; 95% CI, 0.29-0.46) compared with the IHS classification (h(2)=0.36; 95% CI, 0.22-0.42). Exploratory latent class modeling, fitting up to 10 classes, did not identify classes corresponding to either the IHS MO or MA classification. Our data indicate the existence of a continuum of severity, with MA more severe but not etiologically distinct from MO. In searching for predisposing genes, we should therefore expect to find some genes that may underlie all major recurrent headache subtypes, with modifying genetic or environmental factors that may lead to differential expression of the liability for migraine.

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.

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.

Pathogenic Mechanisms of Polycystic Lipomembranous Osteodysplasia with Sclerosing Leukoencephalopathy (PLOSL)

PATHOGENIC MECHANISMS OF PLOSL Polycystic lipomembranous osteodysplasia with sclerosing leukoencephalopathy (PLOSL), also known as Nasu-Hakola disease, is a recessively inherited disease of brain and bone. PLOSL manifests as early-onset progressive dementia and bone fractures. Mutations in the TYROBP (DAP12) and TREM2 genes have been identified as the primary cause of PLOSL. DAP12 and TREM2 encode important signalling molecules in cells of the innate immune system. The mechanism by which loss-of-function of the DAP12/TREM2 signalling complex leads to PLOSL is currently unknown. The aim of this thesis work was to gain insight into the pathogenic mechanisms behind PLOSL. To first identify the central nervous system (CNS) cell types that express both Dap12 and Trem2, the expression patterns of Dap12 and Trem2 in mouse CNS were analyzed. Dap12 and Trem2 expression was seen from embryonic stage to adulthood and microglial cells and oligodendrocytes were identified as the major Dap12/Trem2 producing cells of the CNS. To subsequently identify the pathways and biological processes associated with DAP12/TREM2 mediated signalling in human cells, genome wide transcript analysis of in vitro differentiated dendritic cells (DCs) of PLOSL patients representing functional knockouts of either DAP12 or TREM2 was performed. Both DAP12 and TREM2 deficient cells differentiated into DCs and responded to pathogenic stimuli. However, the DCs showed morphological differences compared to control cells due to defects in the actin filaments. Transcript profiles of the patient DCs showed differential expression of genes involved in immune response and for genes earlier associated with other disorders of the CNS as well as genes involved in the remodeling of bone, linking the findings with the tissue phenotype of PLOSL patients. To analyze the effect of Dap12 deficiency in the CNS, genome wide expression analysis of Dap12 deficient mouse brain and Dap12 deficient microglia as well as functional analysis of Dap12 deficient microglia was performed. Regulation of several pathways involved in synaptic function and transcripts coding for the myelin components was seen in Dap12 knockout mice. Decreased migration, morphological changes and shortened lifespan of the Dap12 knockout microglia was further observed. Taken together, this thesis work showed that both Dap12 and Trem2 are expressed by CNS microglia and that Dap12 deficiency results in functional defects of these cells. Lack of Dap12 in the CNS also leads to synaptic abnormalities even before pathological changes are seen in the tissue level.This work further showed that loss-of-function of DAP12 or TREM2 leads to changes in morphology and gene expression in human dendritic cells. These data underline the functional diversity of the molecules of the innate immune system and implies their significant contribution also in demyelinating CNS disorders, including those resulting in dementia.

The role of MADS and TCP transcription factors in Gerbera hybrida flower development

Angiosperms represent a huge diversity in floral structures. Thus, they provide an attractive target for comparative developmental genetics studies. Research on flower development has focused on few main model plants, and studies on these species have revealed the importance of transcription factors, such as MADS-box and TCP genes, for regulating the floral form. The MADS-box genes determine floral organ identities, whereas the TCP genes are known to regulate flower shape and the number of floral organs. In this study, I have concentrated on these two gene families and their role in regulating flower development in Gerbera hybrida, a species belonging to the large sunflower family (Asteraceae). The Gerbera inflorescence is comprised of hundreds of tightly clustered flowers that differ in their size, shape and function according to their position in the inflorescence. The presence of distinct flower types tells Gerbera apart from the common model species that bear only single kinds of flowers in their inflorescences. The marginally located ray flowers have large bilaterally symmetrical petals and non-functional stamens. The centrally located disc flowers are smaller, have less pronounced bilateral symmetry and carry functional stamens. Early stages of flower development were studied in Gerbera to understand the differentiation of flower types better. After morphological analysis, we compared gene expression between ray and disc flowers to reveal transcriptional differences in flower types. Interestingly, MADS-box genes showed differential expression, suggesting that they might take part in defining flower types by forming flower-type-specific regulatory complexes. Functional analysis of a CYCLOIDEA-like TCP gene GhCYC2 provided evidence that TCP transcription factors are involved in flower type differentiation in Gerbera. The expression of GhCYC2 is ray-flower-specific at early stages of development and activated only later in disc flowers. Overexpression of GhCYC2 in transgenic Gerbera-lines causes disc flowers to obtain ray-flower-like characters, such as elongated petals and disrupted stamen development. The expression pattern and transgenic phenotypes further suggest that GhCYC2 may shape ray flowers by promoting organ fusion. Cooperation of GhCYC2 with other Gerbera CYC-like TCP genes is most likely needed for proper flower type specification, and by this means for shaping the elaborate inflorescence structure. Gerbera flower development was also approached by characterizing B class MADS-box genes, which in the main model plants are known regulators of petal and stamen identity. The four Gerbera B class genes were phylogenetically grouped into three clades; GGLO1 into the PI/GLO clade, GDEF2 and GDEF3 into the euAP3 clade and GDEF1 into the TM6 clade. Putative orthologs for GDEF2 and GDEF3 were identified in other Asteraceae species, which suggests that they appeared through an Asteraceae-specific duplication. Functional analyses indicated that GGLO1 and GDEF2 perform conventional B-function as they determine petal and stamen identities. Our studies on GDEF1 represent the first functional analysis of a TM6-like gene outside the Solanaceae lineage and provide further evidence for the role of TM6 clade members in specifying stamen development. Overall, the Gerbera B class genes showed both commonalities and diversifications with the conventional B-function described in the main model plants.

Transcription of Human Resistin Gene Involves an Interaction of Sp1 with Peroxisome Proliferator-Activating Receptor Gamma (PPAR gamma)

Background: Resistin is a cysteine rich protein, mainly expressed and secreted by circulating human mononuclear cells. While several factors responsible for transcription of mouse resistin gene have been identified, not much is known about the factors responsible for the differential expression of human resistin.Methodology/Principal Finding: We show that the minimal promoter of human resistin lies within similar to 80 bp sequence upstream of the transcriptional start site (-240) whereas binding sites for cRel, CCAAT enhancer binding protein alpha (C/EBP-alpha), activating transcription factor 2 (ATF-2) and activator protein 1 (AP-1) transcription factors, important for induced expression, are present within sequences up to -619. Specificity Protein 1(Sp1) binding site (-276 to -295) is also present and an interaction of Sp1 with peroxisome proliferator activating receptor gamma (PPAR gamma) is necessary for constitutive expression in U937 cells. Indeed co-immunoprecipitation assay demonstrated a direct physical interaction of Sp1 with PPAR gamma in whole cell extracts of U937 cells. Phorbol myristate acetate (PMA) upregulated the expression of resistin mRNA in U937 cells by increasing the recruitment of Sp1, ATF-2 and PPAR gamma on the resistin gene promoter. Furthermore, PMA stimulation of U937 cells resulted in the disruption of Sp1 and PPAR gamma interaction. Chromatin immunoprecipitation (ChIP) assay confirmed the recruitment of transcription factors phospho ATF-2, Sp1, Sp3, PPAR gamma, chromatin modifier histone deacetylase 1 (HDAC1) and the acetylated form of histone H3 but not cRel, C/EBP-alpha and phospho c-Jun during resistingene transcription.Conclusion: Our findings suggest a complex interplay of Sp1 and PPAR gamma along with other transcription factors that drives the expression of resistin in human monocytic U937 cells.