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.

Characterization of genomic diversity in extraintestinal pathogenic Escherichia coli (ExPEC) and development of a diagnostic DNA microarray for the differentiation of ExPEC isolates causing urinary tract infections

Extraintestinal pathogenic Escherichia coli (ExPEC) represent a diverse group of strains of E. coli, which infect extraintestinal sites, such as the urinary tract, the bloodstream, the meninges, the peritoneal cavity, and the lungs. Urinary tract infections (UTIs) caused by uropathogenic E. coli (UPEC), the major subgroup of ExPEC, are among the most prevalent microbial diseases world wide and a substantial burden for public health care systems. UTIs are responsible for serious morbidity and mortality in the elderly, in young children, and in immune-compromised and hospitalized patients. ExPEC strains are different, both from genetic and clinical perspectives, from commensal E. coli strains belonging to the normal intestinal flora and from intestinal pathogenic E. coli strains causing diarrhea. ExPEC strains are characterized by a broad range of alternate virulence factors, such as adhesins, toxins, and iron accumulation systems. Unlike diarrheagenic E. coli, whose distinctive virulence determinants evoke characteristic diarrheagenic symptoms and signs, ExPEC strains are exceedingly heterogeneous and are known to possess no specific virulence factors or a set of factors, which are obligatory for the infection of a certain extraintestinal site (e. g. the urinary tract). The ExPEC genomes are highly diverse mosaic structures in permanent flux. These strains have obtained a significant amount of DNA (predictably up to 25% of the genomes) through acquisition of foreign DNA from diverse related or non-related donor species by lateral transfer of mobile genetic elements, including pathogenicity islands (PAIs), plasmids, phages, transposons, and insertion elements. The ability of ExPEC strains to cause disease is mainly derived from this horizontally acquired gene pool; the extragenous DNA facilitates rapid adaptation of the pathogen to changing conditions and hence the extent of the spectrum of sites that can be infected. However, neither the amount of unique DNA in different ExPEC strains (or UPEC strains) nor the mechanisms lying behind the observed genomic mobility are known. Due to this extreme heterogeneity of the UPEC and ExPEC populations in general, the routine surveillance of ExPEC is exceedingly difficult. In this project, we presented a novel virulence gene algorithm (VGA) for the estimation of the extraintestinal virulence potential (VP, pathogenicity risk) of clinically relevant ExPECs and fecal E. coli isolates. The VGA was based on a DNA microarray specific for the ExPEC phenotype (ExPEC pathoarray). This array contained 77 DNA probes homologous with known (e.g. adhesion factors, iron accumulation systems, and toxins) and putative (e.g. genes predictably involved in adhesion, iron uptake, or in metabolic functions) ExPEC virulence determinants. In total, 25 of DNA probes homologous with known virulence factors and 36 of DNA probes representing putative extraintestinal virulence determinants were found at significantly higher frequency in virulent ExPEC isolates than in commensal E. coli strains. We showed that the ExPEC pathoarray and the VGA could be readily used for the differentiation of highly virulent ExPECs both from less virulent ExPEC clones and from commensal E. coli strains as well. Implementing the VGA in a group of unknown ExPECs (n=53) and fecal E. coli isolates (n=37), 83% of strains were correctly identified as extraintestinal virulent or commensal E. coli. Conversely, 15% of clinical ExPECs and 19% of fecal E. coli strains failed to raster into their respective pathogenic and non-pathogenic groups. Clinical data and virulence gene profiles of these strains warranted the estimated VPs; UPEC strains with atypically low risk-ratios were largely isolated from patients with certain medical history, including diabetes mellitus or catheterization, or from elderly patients. In addition, fecal E. coli strains with VPs characteristic for ExPEC were shown to represent the diagnostically important fraction of resident strains of the gut flora with a high potential of causing extraintestinal infections. Interestingly, a large fraction of DNA probes associated with the ExPEC phenotype corresponded to novel DNA sequences without any known function in UTIs and thus represented new genetic markers for the extraintestinal virulence. These DNA probes included unknown DNA sequences originating from the genomic subtractions of four clinical ExPEC isolates as well as from five novel cosmid sequences identified in the UPEC strains HE300 and JS299. The characterized cosmid sequences (pJS332, pJS448, pJS666, pJS700, and pJS706) revealed complex modular DNA structures with known and unknown DNA fragments arranged in a puzzle-like manner and integrated into the common E. coli genomic backbone. Furthermore, cosmid pJS332 of the UPEC strain HE300, which carried a chromosomal virulence gene cluster (iroBCDEN) encoding the salmochelin siderophore system, was shown to be part of a transmissible plasmid of Salmonella enterica. Taken together, the results of this project pointed towards the assumptions that first, (i) homologous recombination, even within coding genes, contributes to the observed mosaicism of ExPEC genomes and secondly, (ii) besides en block transfer of large DNA regions (e.g. chromosomal PAIs) also rearrangements of small DNA modules provide a means of genomic plasticity. The data presented in this project supplemented previous whole genome sequencing projects of E. coli and indicated that each E. coli genome displays a unique assemblage of individual mosaic structures, which enable these strains to successfully colonize and infect different anatomical sites.

Identification of genes associated with tumorigenesis of meibomian cell carcinoma by microarray analysis

Meibomian cell carcinoma (MCC) is a malignant tumor of the meibomian glands located in the eyelids. No information exists on the cytogenctic and genetic aspects of MCC. There is no report on the gene expression profile of MCC. Thus there is a need, for both scientific and clinical reasons, to identify genes and pathways that are involved in the development and progression of MCC. We analyzed the gene expression profile of MCC by the microarray technique. Forty-four genes were upregulated and 149 genes were downregulated in MCC. Differential expression data were confirmed for 5 genes by semiquantitative RT-PCR in MCC tumors: GTF2H4, RBM12, UBE2D3, DDX17, and LZTS1. We found dysregulation of two major pathways in MCC: MAPK and JAK/STAT. Clusters of genes on chromosomes 1, 12, and 19 were dysregUlated in MCC. The data presented here will facilitate the identification of specific markers and therapeutic targets for the treatment of MCC patients. (c) 2007 Elsevier Inc. All rights reserved.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Microbial identification by detection of ligation probes on DNA microarray

Microbes in natural and artificial environments as well as in the human body are a key part of the functional properties of these complex systems. The presence or absence of certain microbial taxa is a correlate of functional status like risk of disease or course of metabolic processes of a microbial community. As microbes are highly diverse and mostly notcultivable, molecular markers like gene sequences are a potential basis for detection and identification of key types. The goal of this thesis was to study molecular methods for identification of microbial DNA in order to develop a tool for analysis of environmental and clinical DNA samples. Particular emphasis was placed on specificity of detection which is a major challenge when analyzing complex microbial communities. The approach taken in this study was the application and optimization of enzymatic ligation of DNA probes coupled with microarray read-out for high-throughput microbial profiling. The results show that fungal phylotypes and human papillomavirus genotypes could be accurately identified from pools of PCR amplicons generated from purified sample DNA. Approximately 1 ng/μl of sample DNA was needed for representative PCR amplification as measured by comparisons between clone sequencing and microarray. A minimum of 0,25 amol/μl of PCR amplicons was detectable from amongst 5 ng/μl of background DNA, suggesting that the detection limit of the test comprising of ligation reaction followed by microarray read-out was approximately 0,04%. Detection from sample DNA directly was shown to be feasible with probes forming a circular molecule upon ligation followed by PCR amplification of the probe. In this approach, the minimum detectable relative amount of target genome was found to be 1% of all genomes in the sample as estimated from 454 deep sequencing results. Signal-to-noise of contact printed microarrays could be improved by using an internal microarray hybridization control oligonucleotide probe together with a computational algorithm. The algorithm was based on identification of a bias in the microarray data and correction of the bias as shown by simulated and real data. The results further suggest semiquantitative detection to be possible by ligation detection, allowing estimation of target abundance in a sample. However, in practise, comprehensive sequence information of full length rRNA genes is needed to support probe design with complex samples. This study shows that DNA microarray has the potential for an accurate microbial diagnostic platform to take advantage of increasing sequence data and to replace traditional, less efficient methods that still dominate routine testing in laboratories. The data suggests that ligation reaction based microarray assay can be optimized to a degree that allows good signal-tonoise and semiquantitative detection.

Novel Functions of (p)ppGpp and Cyclic di-GMP in Mycobacterial Physiology Revealed by Phenotype Microarray Analysis of Wild-Type and Isogenic Strains of Mycobacterium smegmatis

The bacterial second messengers (p)ppGpp and bis-(3'-5')-cyclic dimeric GMP (c-di-GMP) regulate important functions, such as transcription, virulence, biofilm formation, and quorum sensing. In mycobacteria, they regulate long-term survival during starvation, pathogenicity, and dormancy. Recently, a Pseudomonas aeruginosa strain lacking (p) ppGpp was shown to be sensitive to multiple classes of antibiotics and defective in biofilm formation. We were interested to find out whether Mycobacterium smegmatis strains lacking the gene for either (p)ppGpp synthesis (Delta rel(Msm)) or c-di-GMP synthesis (Delta dcpA) would display similar phenotypes. We used phenotype microarray technology to compare the growth of the wild-type and the knockout strains in the presence of several antibiotics. Surprisingly, the Delta rel(Msm) and Delta dcpA strains showed enhanced survival in the presence of many antibiotics, but they were defective in biofilm formation. These strains also displayed altered surface properties, like impaired sliding motility, rough colony morphology, and increased aggregation in liquid cultures. Biofilm formation and surface properties are associated with the presence of glycopeptidolipids (GPLs) in the cell walls of M. smegmatis. Thin-layer chromatography analysis of various cell wall fractions revealed that the levels of GPLs and polar lipids were reduced in the knockout strains. As a result, the cell walls of the knockout strains were significantly more hydrophobic than those of the wild type and the complemented strains. We hypothesize that reduced levels of GPLs and polar lipids may contribute to the antibiotic resistance shown by the knockout strains. Altogether, our data suggest that (p)ppGpp and c-di-GMP may be involved in the metabolism of glycopeptidolipids and polar lipids in M. smegmatis.

Independent component analysis of microarray data in the study of endometrial cancer.

Gene microarray technology is highly effective in screening for differential gene expression and has hence become a popular tool in the molecular investigation of cancer. When applied to tumours, molecular characteristics may be correlated with clinical features such as response to chemotherapy. Exploitation of the huge amount of data generated by microarrays is difficult, however, and constitutes a major challenge in the advancement of this methodology. Independent component analysis (ICA), a modern statistical method, allows us to better understand data in such complex and noisy measurement environments. The technique has the potential to significantly increase the quality of the resulting data and improve the biological validity of subsequent analysis. We performed microarray experiments on 31 postmenopausal endometrial biopsies, comprising 11 benign and 20 malignant samples. We compared ICA to the established methods of principal component analysis (PCA), Cyber-T, and SAM. We show that ICA generated patterns that clearly characterized the malignant samples studied, in contrast to PCA. Moreover, ICA improved the biological validity of the genes identified as differentially expressed in endometrial carcinoma, compared to those found by Cyber-T and SAM. In particular, several genes involved in lipid metabolism that are differentially expressed in endometrial carcinoma were only found using this method. This report highlights the potential of ICA in the analysis of microarray data.

Genomotipificación de Camylobacter jejuni mediante microarrays de ADN

354 p. (Bibliogr. 271-303) - Correo electrónico de la autora: andrea.guridi@gmail.com

Estudo do polimorfismo Ser49Gly do gene do receptor beta adrenérgico 1 (β-adr 1) na população do estado do Rio de Janeiro, Brasil estratificada por cor da pele e ancestralidade genômica

As doenças cardiovasculares possuem a maior taxa de óbitos no mundo, e notavelmente nos últimos anos as pesquisas genéticas sobre as mesmas estão baseadas em estudos de associação, no qual o gene suspeito que esteja em maior frequência entre os pacientes passa a ser considerado um possível fator causal. Os polimorfismos genéticos que ocorrem no receptor beta-adrenérgico podem resultar em mudanças significativas na função do receptor, podendo acarretar fisiopatologias. Neste trabalho, o objetivo foi estimar a diversidade e a frequência do polimorfismo Ser49Gly do gene do receptor beta-adrenérgico 1 a partir de uma amostra de 188 indivíduos da população do Estado do Rio de Janeiro. As frequências também foram analisadas a partir da estratificação da amostra por critério fenotípico em função do padrão de cor da pele em (negros e não negros) ou ancestralidade genética em (afrodescendente e não afrodescendente), definida através da informação dos marcadores de ancestralidade Indels e SNP de cromossomo Y, para avaliar se os padrões de ancestralidade ou cor da pele são fundamentais para a diferenciação e distanciamento genético. Fragmentos de interesse foram amplificados por PCR (reação de cadeia de polimerase) com primers específicos para o marcador Ser49Gly e as reações de genotipagem foram realizadas com enzimas de restrição Eco0109I. Os valores da heterozigosidade variaram entre 0,25-0,50 e 0,20-0,41 nos grupos estratificados por ancestralidade e cor da pele, respectivamente. No que diz respeito à análise do equilíbrio de Hardy-Weinberg, não houve um desvio significativo na distribuição do marcador nas amostras gerais do Estado do Rio de Janeiro, ou mesmo nas amostras estratificadas. A distribuição dos alelos na amostra dos 188 indivíduos da população geral do Rio de Janeiro (AC_RJ) mostrou uma frequência de 80,30% e 19,70% para o alelo selvagem e mutado Ser49Gly, respectivamente. A comparação das análises sobre a distribuição das frequências alélicas para este marcador mostrou a ocorrência de diferenças significativas na distribuição das frequências alélicas entre negros e não negros e afrodescendentes e não afrodescendentes. A diferença significativa observada entre os negros e afrodescendentes, foi em menor grau de distanciamento. A informação obtida em relação à ancestralidade foi crucial para a obtenção dos dados sobre o aumento da variável mutada do polimorfismo Ser49Gly nas populações negras e afrodescendentes do Estado Rio de Janeiro. Tal evidência, em combinação com estudos clínicos podem contribuir para uma análise pormenorizada do padrão de susceptibilidade à doença em questão, em falhas do mecanismo deste receptor.