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.

Mutations as molecular tools : The metabolic-rate dependent molecular clock and DNA barcoding of allied species

Mutation and recombination are the fundamental processes leading to genetic variation in natural populations. This variation forms the raw material for evolution through natural selection and drift. Therefore, studying mutation rates may reveal information about evolutionary histories as well as phylogenetic interrelationships of organisms. In this thesis two molecular tools, DNA barcoding and the molecular clock were examined. In the first part, the efficiency of mutations to delineate closely related species was tested and the implications for conservation practices were assessed. The second part investigated the proposition that a constant mutation rate exists within invertebrates, in form of a metabolic-rate dependent molecular clock, which can be applied to accurately date speciation events. DNA barcoding aspires to be an efficient technique to not only distinguish between species but also reveal population-level variation solely relying on mutations found on a short stretch of a single gene. In this thesis barcoding was applied to discriminate between Hylochares populations from Russian Karelia and new Hylochares findings from the greater Helsinki region in Finland. Although barcoding failed to delineate the two reproductively isolated groups, their distinct morphological features and differing life-history traits led to their classification as two closely related, although separate species. The lack of genetic differentiation appears to be due to a recent divergence event not yet reflected in the beetles molecular make-up. Thus, the Russian Hylochares was described as a new species. The Finnish species, previously considered as locally extinct, was recognized as endangered. Even if, due to their identical genetic make-up, the populations had been regarded as conspecific, conservation strategies based on prior knowledge from Russia would not have guaranteed the survival of the Finnish beetle. Therefore, new conservation actions based on detailed studies of the biology and life-history of the Finnish Hylochares were conducted to protect this endemic rarity in Finland. The idea behind the strict molecular clock is that mutation rates are constant over evolutionary time and may thus be used to infer species divergence dates. However, one of the most recent theories argues that a strict clock does not tick per unit of time but that it has a constant substitution rate per unit of mass-specific metabolic energy. Therefore, according to this hypothesis, molecular clocks have to be recalibrated taking body size and temperature into account. This thesis tested the temperature effect on mutation rates in equally sized invertebrates. For the first dataset (family Eucnemidae, Coleoptera) the phylogenetic interrelationships and evolutionary history of the genus Arrhipis had to be inferred before the influence of temperature on substitution rates could be studied. Further, a second, larger invertebrate dataset (family Syrphidae, Diptera) was employed. Several methodological approaches, a number of genes and multiple molecular clock models revealed that there was no consistent relationship between temperature and mutation rate for the taxa under study. Thus, the body size effect, observed in vertebrates but controversial for invertebrates, rather than temperature may be the underlying driving force behind the metabolic-rate dependent molecular clock. Therefore, the metabolic-rate dependent molecular clock does not hold for the here studied invertebrate groups. This thesis emphasizes that molecular techniques relying on mutation rates have to be applied with caution. Whereas they may work satisfactorily under certain conditions for specific taxa, they may fail for others. The molecular clock as well as DNA barcoding should incorporate all the information and data available to obtain comprehensive estimations of the existing biodiversity and its evolutionary history.

Mutations of mitochondrial DNA polymerase gamma: an important cause of neurological disorders

Thesis focuses on mutations of POLG1 gene encoding catalytic subunit polγ-α of mitochondrial DNA polymerase gamma holoenzyme (polG) and the association of mutations with different clinical phenotypes. In addition, particular defective mutant variants of the protein were characterized biochemically in vitro. PolG-holoenzyme is the sole DNA polymerase found in mitochondria. It is involved in replication and repair of the mitochondrial genome, mtDNA. Holoenzyme also includes the accessory subunit polγ-β, which is required for the enhanced processivity of polγ-α. Defective polγ-α causes accumulation of secondary mutations on mtDNA, which leads to a defective oxidative phosphorylation system. The clinical consequences of such mutations are variable, affecting nervous system, skeletal muscles, liver and other post-mitotic tissues. The aims of the studies included: 1) Determination of the role of POLG1 mutations in neurological syndromes with features of mitochondrial dysfunction and an unknown molecular cause. 2) Development and set up of diagnostic tests for routine clinical purposes. 3) Biochemical characterization of the functional consequences of the identified polγ-α variants. Studies describe new neurological phenotypes in addition to PEO caused by POLG1 mutations, including parkinsonism, premature amenorrhea, ataxia and Parkinson s disease (PD). POLG1 mutations and polymorphisms are both common and/or potential genetic risk factors at least among the Finnish population. The major findings and applications reported here are: 1) POLG1 mutations cause parkinsonism and premature menopause in PEO families in either a recessive or a dominant manner. 2) A common recessive POLG1 mutations (A467T and W748S) in the homozygous state causes severe adult or juvenile-onset ataxia without muscular symptoms or histological or mtDNA abnormalities in muscles. 3) A common recessive pathogenic change A467T can also cause a mild dominant disease in heterozygote carriers. 4) The A467T variant shows reduced polymerase activity due to defective template binding. 5) Rare polyglutamine tract length variants of POLG1 are significantly enriched in Finnish idiopathic Parkinson s disease patients. 6) Dominant mutations are clearly restricted to the highly conserved polymerase domain motifs, whereas recessive ones are more evenly distributed along the protein. The present results highlight and confirm the new role of mitochondria in parkinsonism/Parkinson s disease and describe a new mitochondrial ataxia. Based on these results, a POLG1 diagnostic routine has been set up in Helsinki University Central Hospital (HUSLAB).

MECHANISMS AND EFFECTS OF MITOCHONDRIAL DNA INSTABILITY AND COPY NUMBER MANIPULATION

Defects in mitochondrial DNA (mtDNA) maintenance cause a range of human diseases, including autosomal dominant progressive external ophthalmoplegia (adPEO). This study aimed to clarify the molecular background of adPEO. We discovered that deoxynucleoside triphosphate (dNTP) metabolism plays a crucial in mtDNA maintenance and were thus prompted to search for therapeutic strategies based on the modulation of cellular dNTP pools or mtDNA copy number. Human mtDNA is a 16.6 kb circular molecule present in hundreds to thousands of copies per cell. mtDNA is compacted into nucleoprotein clusters called nucleoids. mtDNA maintenance diseases result from defects in nuclear encoded proteins that maintain the mtDNA. These syndromes typically afflict highly differentiated, post-mitotic tissues such as muscle and nerve, but virtually any organ can be affected. adPEO is a disease where mtDNA molecules with large-scale deletions accumulate in patients tissues, particularly in skeletal muscle. Mutations in five nuclear genes, encoding the proteins ANT1, Twinkle, POLG, POLG2 and OPA1, have previously been shown to cause adPEO. Here, we studied a large North American pedigree with adPEO, and identified a novel heterozygous mutation in the gene RRM2B, which encodes the p53R2 subunit of the enzyme ribonucleotide reductase (RNR). RNR is the rate-limiting enzyme in dNTP biosynthesis, and is required both for nuclear and mitochondrial DNA replication. The mutation results in the expression of a truncated form of p53R2, which is likely to compete with the wild-type allele. A change in enzyme function leads to defective mtDNA replication due to altered dNTP pools. Therefore, RRM2B is a novel adPEO disease gene. The importance of adequate dNTP pools and RNR function for mtDNA maintenance has been established in many organisms. In yeast, induction of RNR has previously been shown to increase mtDNA copy number, and to rescue the phenotype caused by mutations in the yeast mtDNA polymerase. To further study the role of RNR in mammalian mtDNA maintenance, we used mice that broadly overexpress the RNR subunits Rrm1, Rrm2 or p53R2. Active RNR is a heterotetramer consisting of two large subunits (Rrm1) and two small subunits (either Rrm2 or p53R2). We also created bitransgenic mice that overexpress Rrm1 together with either Rrm2 or p53R2. In contrast to the previous findings in yeast, bitransgenic RNR overexpression led to mtDNA depletion in mouse skeletal muscle, without mtDNA deletions or point mutations. The mtDNA depletion was associated with imbalanced dNTP pools. Furthermore, the mRNA expression levels of Rrm1 and p53R2 were found to correlate with mtDNA copy number in two independent mouse models, suggesting nuclear-mitochondrial cross talk with regard to mtDNA copy number. We conclude that tight regulation of RNR is needed to prevent harmful alterations in the dNTP pool balance, which can lead to disordered mtDNA maintenance. Increasing the copy number of wild-type mtDNA has been suggested as a strategy for treating PEO and other mitochondrial diseases. Only two proteins are known to cause a robust increase in mtDNA copy number when overexpressed in mice; the mitochondrial transcription factor A (TFAM), and the mitochondrial replicative helicase Twinkle. We studied the mechanisms by which Twinkle and TFAM elevate mtDNA levels, and showed that Twinkle specifically implements mtDNA synthesis. Furthermore, both Twinkle and TFAM were found to increase mtDNA content per nucleoid. Increased mtDNA content in mouse tissues correlated with an age-related accumulation of mtDNA deletions, depletion of mitochondrial transcripts, and progressive respiratory dysfunction. Simultaneous overexpression of Twinkle and TFAM led to a further increase in the mtDNA content of nucleoids, and aggravated the respiratory deficiency. These results suggested that high mtDNA levels have detrimental long-term effects in mice. These data have to be considered when developing and evaluating treatment strategies for elevating mtDNA copy number.

Uniparental DNA markers and forensic genetics in Finland

Over the years, a wide range of methods to verify identity have been developed. Molecular markers have been used for identification since the 1920s, commencing with blood types and culminating with the advent of DNA techniques in the 1980s. Identification is required by authorities in many occasions, e.g. in disputed paternity cases, identification of deceased, or crime investigation. To clarify maternal and paternal lineages, uniparental DNA markers in mtDNA and Y-chromosome can be utilized. These markers have several advantages: male specific Y-chromosome can be used to identify a male from a mixture of male and female cells, e.g. in rape cases. MtDNA is durable and has a high copy number, allowing analyses even from old or degraded samples. However, both markers are lineage-specific, not individualizing, and susceptible to genetic drift. Prior to the application of any DNA marker in forensic casework, it is of utmost importance to investigate its qualities and peculiarities in the target population. Earlier studies on the Finnish population have shown reduced variation in the Y-chromosome, but in mtDNA results have been ambiguous. The obtained results confirmed the low diversity in Y-chromosome in Finland. Detailed population analysis revealed large regional differences, and extremely reduced diversity especially in East Finland. Analysis of the qualities affecting Y-chromosomal short tandem repeat (Y-STR) variation and mutation frequencies, and search of new polymorphic markers resulted a set of Y-STRs with especially high diversity in Finland. Contrary to Y-chromosome, neither reduced diversity nor regional differences were found in mtDNA within Finland. In fact, mtDNA diversity was found similar to other European populations. The revealed peculiarities in the uniparental markers are a legacy of the Finnish population history. The obtained results challenge the traditional explanation which emphasizes relatively recent founder effects creating the observed east-west patterns. Uniparentally inherited markers, both mtDNA and Y-chromosome, are applicable for identification purposes in Finland. By adjusting the analysed Y marker set to meet the characteristics of Finnish population, Y-chromosomal diversity increases and the regional differentiation decreases, resulting increase in discrimination power and thus usefulness of Y-chromosomal analysis in forensic casework.

Integration of miRNA and mRNA expression with DNA copy number of Ewing sarcoma cell lines

Ewing sarcoma is an aggressive and poorly differentiated malignancy of bone and soft tissue. It primarily affects children, adolescents, and young adults, with a slight male predominance. It is characterized by a translocation between chromosomes 11 and 22 resulting in the EWSR1-FLI1fusion transcription factor. The aim of this study is to identify putative Ewing sarcoma target genes through an integrative analysis of three microarray data sets. Array comparative genomic hybridization is used to measure changes in DNA copy number, and analyzed to detect common chromosomal aberrations. mRNA and miRNA microarrays are used to measure expression of protein-coding and miRNA genes, and these results integrated with the copy number data. Chromosomal aberrations typically contain also bystanders in addition to the driving tumor suppressor and oncogenes, and integration with expression helps to identify the true targets. Correlation between expression of miRNAs and their predicted target mRNAs is also evaluated to assess the results of post-transcriptional miRNA regulation on mRNA levels. The highest frequencies of copy number gains were identified in chromosome 8, 1q, and X. Losses were most frequent in 9p21.3, which also showed an enrichment of copy number breakpoints relative to the rest of the genome. Copy number losses in 9p21.3 were found have a statistically significant effect on the expression of MTAP, but not on CDKN2A, which is a known tumor-suppressor in the same locus. MTAP was also down-regulated in the Ewing sarcoma cell lines compared to mesenchymal stem cells. Genes exhibiting elevated expression in association with copy number gains and up-regulation compared to the reference samples included DCAF7, ENO2, MTCP1, andSTK40. Differentially expressed miRNAs were detected by comparing Ewing sarcoma cell lines against mesenchymal stem cells. 21 up-regulated and 32 down-regulated miRNAs were identified, includingmiR-145, which has been previously linked to Ewing sarcoma. The EWSR1-FLI1 fusion gene represses miR-145, which in turn targets FLI1 forming a mutually repressive feedback loop. In addition higher expression linked to copy number gains and compared to mesenchymal stem cells, STK40 was also found to be a target of four different miRNAs that were all down-regulated in Ewing sarcoma cell lines compared to the reference samples. SLCO5A1 was identified as the only up-regulated gene within a frequently gained region in chromosome 8. This region was gained in over 90 % of the cell lines, and also with a higher frequency than the neighboring regions. In addition, SLCO5A1 was found to be a target of three miRNAs that were down-regulated compared to the mesenchymal stem cells.

DNA damage recognition in the normal epithelium of human prostate and seminal vesicles

Prostate cancer is one of the most prevalent cancer types in men. The development of prostate tumors is known to require androgen exposure, and several pathways governing cell growth are deregulated in prostate tumorigenesis. Recent genetic studies have revealed that complex gene fusions and copy - number alterations are frequent in prostate cancer, a unique feature among solid tumors. These chromosomal aberrations are though to arise as a consequence of faulty repair of DNA double strand breaks (DSB). Most repair mechanisms have been studied in detail in cancer cell lines, but how DNA damage is detected and repaired in normal differentiated human cells has not been widely addressed. The events leading to the gene fusions in prostate cancer are under rigorous studies, as they not only shed light on the basic pathobiologic mechanisms but may also produce molecular targets for prostate cancer treatment and prevention. Prostate and seminal vesicles are part of the male reproductive system. They share similar structure and function but differ dramatically in their cancer incidence. Approximately fifty primary seminal vesicle carcinomas have been reported worldwide. Surprisingly, only little is known on why seminal vesicles are resistant to neoplastic changes. As both tissues are androgen dependent, it is a mystery that androgen signaling would only lead to tumors in prostate tissue. In this work, we set up novel ex vivo human tissue culture models of prostate and seminal vesicles, and used them to study how DNA damage is recognized in normal epithelium. One of the major DNA - damage inducible pathways, mediated by the ATM kinase, was robustly activated in all main cell types of both tissues. Interestingly, we discovered that secretory epithelial cells had less histone variant H2A.X and after DNA damage lower levels of H2AX were phosphorylated on serine 139 (γH2AX) than in basal or stromal cells. γH2AX has been considered essential for efficient DSB repair, but as there were no significant differences in the γH2AX levels between the two tissues, it seems more likely that the role of γH2AX is less important in postmitotic cells. We also gained insight into the regulation of p53, an important transcription factor that protects genomic integrity via multiple mechanisms, in human tissues. DSBs did not lead to a pronounced activation of p53, but treatments causing transcriptional stress, on the other hand, were able to launch a notable p53 response in both tissue types. In general, ex vivo culturing of human tissues provided unique means to study differentiated cells in their relevant tissue context, and is suited for testing novel therapeutic drugs before clinical trials. In order to study how prostate and seminal vesicle epithelial cells are able to activate DNA damage induced cell cycle checkpoints, we used primary cultures of prostate and seminal vesicle epithelial cells. To our knowledge, we are the first to report isolation of human primary seminal vesicle cells. Surprisingly, human prostate epithelial cells did not activate cell cycle checkpoints after DSBs in part due to low levels of Wee1A, a kinase regulating CDK activity, while primary seminal vesicle epithelial cells possessed proficient cell cycle checkpoints and expressed high levels of Wee1A. Similarly, seminal vesicle cells showed a distinct activation of the p53 - pathway after DSBs that did not occur in prostate epithelial cells. This indicates that p53 protein function is under different control mechanisms in the two cell types, which together with proficient cell cycle checkpoints may be crucial in protecting seminal vesicles from endogenous and exogenous DNA damaging factors and, as a consequence, from carcinogenesis. These data indicate that two very similar organs of male reproductive system do not respond to DNA damage similarly. The differentiated, non - replicating cells of both tissues were able to recognize DSBs, but under proliferation human prostate epithelial cells had deficient activation of the DNA damage response. This suggests that prostate epithelium is most vulnerable to accumulating genomic aberrations under conditions where it needs to proliferate, for example after inflammatory cellular damage.

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.