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.

Analysis of body composition in individuals with high bone mass reveals a marked increase in fat mass in women but not men

Context: High bone mass (HBM), detected in 0.2% of dual-energy x-ray absorptiometry (DXA) scans, is characterized by raised body mass index, the basis for which is unclear. Objective: To investigate why body mass index is elevated in individuals with HBM, we characterized body composition and examined whether differences could be explained by bone phenotypes, eg, bone mass and/or bone turnover. Design, Setting, and Participants: We conducted a case-control study of 153 cases with unexplained HBM recruited from 4 UK centers by screening 219 088 DXA scans. Atotal of 138 first-degree relatives (of whom 51 had HBM) and 39 spouses were also recruited. Unaffected individuals served as controls. Main Outcome Measures: We measured fat mass, by DXA, and bone turnover markers. Results: Amongwomen, fat mass was inversely related to age in controls (P<.01), but not in HBM cases (P<.96) in whom mean fat mass was 8.9 [95% CI 4.7, 13.0] kg higher compared with controls (fully adjusted mean difference, P<.001). Increased fat mass in male HBM cases was less marked (gender interaction P = .03). Compared with controls, lean mass was also increased in female HBM cases (by 3.3 [1.2, 5.4] kg; P<.002); however, lean mass increases wereless marked than fat mass increases, resulting in 4.5% lower percentage lean mass in HBM cases (P<.001). Osteocalcin was also lower in female HBM cases compared with controls (by 2.8 [0.1, 5.5]μg/L; P = .04). Differences in fat mass were fully attenuated after hip bone mineral density (BMD) adjustment (P = .52) but unchanged after adjustment for bone turnover (P < .001), whereas the greater hip BMD in female HBM cases was minimally attenuated by fat mass adjustment (P<.001). Conclusions: HBM is characterized by a marked increase in fat mass in females, statistically explained by their greater BMD, but not by markers of bone turnover. Copyright © 2013 by The Endocrine Society.

Elevated circulating sclerostin concentrations in individuals with high bone mass, with and without LRP5 mutations

CONTEXT: The role and importance of circulating sclerostin is poorly understood. High bone mass (HBM) caused by activating LRP5 mutations has been reported to be associated with increased plasma sclerostin concentrations; whether the same applies to HBM due to other causes is unknown. OBJECTIVE: Our objective was to determine circulating sclerostin concentrations in HBM. DESIGN AND PARTICIPANTS: In this case-control study, 406 HBM index cases were identified by screening dual-energy x-ray absorptiometry (DXA) databases from 4 United Kingdom centers (n = 219 088), excluding significant osteoarthritis/artifact. Controls comprised unaffected relatives and spouses. MAIN MEASURES: Plasma sclerostin; lumbar spine L1, total hip, and total body DXA; and radial and tibial peripheral quantitative computed tomography (subgroup only) were evaluated. RESULTS: Sclerostin concentrations were significantly higher in both LRP5 HBM and non-LRP5 HBM cases compared with controls: mean (SD) 130.1 (61.7) and 88.0 (39.3) vs 66.4 (32.3) pmol/L (both P < .001, which persisted after adjustment for a priori confounders). In combined adjusted analyses of cases and controls, sclerostin concentrations were positively related to all bone parameters found to be increased in HBM cases (ie, L1, total hip, and total body DXA bone mineral density and radial/tibial cortical area, cortical bone mineral density, and trabecular density). Although these relationships were broadly equivalent in HBM cases and controls, there was some evidence that associations between sclerostin and trabecular phenotypes were stronger in HBM cases, particularly for radial trabecular density (interaction P < .01). CONCLUSIONS: Circulating plasma sclerostin concentrations are increased in both LRP5 and non-LRP5 HBM compared with controls. In addition to the general positive relationship between sclerostin and DXA/peripheral quantitative computed tomography parameters, genetic factors predisposing to HBM may contribute to increased sclerostin levels.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNAs' Target Sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)-, and femoral neck (FN)-bone mineral density (BMD). In stage I, 41,102 poly-miRTSs were meta-analyzed in 7 cohorts with a genome-wide significance (GWS) α=0.05/41,102=1.22×10-6. By applying α=5×10-5 (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P-value=7.67×10-6 and 1.58×10-5) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P-value=5.08×10-3) at α=0.10/11=9.09×10-3. PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P-value=7.55×10-6) at α=0.05/2=0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P-value=8.87×10-12). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation. © The Author 2015. Published by Oxford University Press. All rights reserved.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.

Mutations in known monogenic high bone mass loci only explain a small proportion of high bone mass cases

High bone mass (HBM) can be an incidental clinical finding; however, monogenic HBM disorders (eg, LRP5 or SOST mutations) are rare. We aimed to determine to what extent HBM is explained by mutations in known HBM genes. A total of 258 unrelated HBM cases were identified from a review of 335,115 DXA scans from 13 UK centers. Cases were assessed clinically and underwent sequencing of known anabolic HBM loci: LRP5 (exons 2, 3, 4), LRP4 (exons 25, 26), SOST (exons 1, 2, and the van Buchem's disease [VBD] 52-kb intronic deletion 3'). Family members were assessed for HBM segregation with identified variants. Three-dimensional protein models were constructed for identified variants. Two novel missense LRP5 HBM mutations ([c.518C>T; p.Thr173Met], [c.796C>T; p.Arg266Cys]) were identified, plus three previously reported missense LRP5 mutations ([c.593A>G; p.Asn198Ser], [c.724G>A; p.Ala242Thr], [c.266A>G; p.Gln89Arg]), associated with HBM in 11 adults from seven families. Individuals with LRP5 HBM ( approximately prevalence 5/100,000) displayed a variable phenotype of skeletal dysplasia with increased trabecular BMD and cortical thickness on HRpQCT, and gynoid fat mass accumulation on DXA, compared with both non-LRP5 HBM and controls. One mostly asymptomatic woman carried a novel heterozygous nonsense SOST mutation (c.530C>A; p.Ser177X) predicted to prematurely truncate sclerostin. Protein modeling suggests the severity of the LRP5-HBM phenotype corresponds to the degree of protein disruption and the consequent effect on SOST-LRP5 binding. We predict p.Asn198Ser and p.Ala242Thr directly disrupt SOST binding; both correspond to severe HBM phenotypes (BMD Z-scores +3.1 to +12.2, inability to float). Less disruptive structural alterations predicted from p.Arg266Cys, p.Thr173Met, and p.Gln89Arg were associated with less severe phenotypes (Z-scores +2.4 to +6.2, ability to float). In conclusion, although mutations in known HBM loci may be asymptomatic, they only account for a very small proportion ( approximately 3%) of HBM individuals, suggesting the great majority are explained by either unknown monogenic causes or polygenic inheritance.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Nephrin-associated molecules in human glomerular diseases

The glomerular epithelial cells and their intercellular junctions, termed slit diaphragms, are essential components of the filtration barrier in the kidney glomerulus. Nephrin is a transmembrane adhesion protein of the slit diaphragm and a signalling molecule regulating podocyte physiology. In congenital nephrotic syndrome of the Finnish type, mutation of nephrin leads to disruption of the permeability barrier and leakage of plasma proteins into the urine. This doctoral thesis hypothesises that novel nephrin-associated molecules are involved in the function of the filtration barrier in health and disease. Bioinformatics tools were utilized to identify novel nephrin-like molecules in genomic databases, and their distribution in the kidney and other tissues was investigated. Filtrin, a novel nephrin homologue, is expressed in the glomerular podocytes and, according to immunoelectron microscopy, localizes at the slit diaphragm. Interestingly, the nephrin and filtrin genes, NPHS1 and KIRREL2, locate in a head-to-head orientation on chromosome 19q13.12. Another nephrin-like molecule, Nphs1as was cloned in mouse, however, no expression was detected in the kidney but instead in the brain and lymphoid tissue. Notably, Nphs1as is transcribed from the nephrin locus in an antisense orientation. The glomerular mRNA and protein levels of filtrin were measured in kidney biopsies of patients with proteinuric diseases, and marked reduction of filtrin mRNA levels was detected in the proteinuric samples as compared to controls. In addition, altered distribution of filtrin in injured glomeruli was observed, with the most prominent decrease of the expression in focal segmental glomerulosclerosis. The role of the slit diaphragm-associated genes for the development of diabetic nephropathy was investigated by analysing single nucleotide polymorphisms. The genes encoding filtrin, densin-180, NEPH1, podocin, and alpha-actinin-4 were analysed, and polymorphisms at the alpha-actinin-4 gene were associated with diabetic nephropathy in a gender-dependent manner. Filtrin is a novel podocyte-expressed protein with localization at the slit diaphragm, and the downregulation of filtrin seems to be characteristic for human proteinuric diseases. In the context of the crucial role of nephrin for the glomerular filter, filtrin appears to be a potential candidate molecule for proteinuria. Although not expressed in the kidney, the nephrin antisense Nphs1as may regulate the expression of nephrin in extrarenal tissues. The genetic association analysis suggested that the alpha-actinin-4 gene, encoding an actin-filament cross-linking protein of the podocytes, may contribute to susceptibility for diabetic nephropathy.

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.

Applications of coralline hydroxyapatite with bioabsorbable containment and reinforcement as bone graft substitute : An experimental study

The aim of the present experimental study was to find out if the applications of coralline hydroxyapatite (HA) can be improved by using bioabsorbable containment or binding substance with particulate HA in mandibular contour augmentation and by using bioabsorbable fibre-reinforced HA blocks in filling bone defects and in anterior lumbar interbody fusion. The use of a separate curved polyglycolide (PGA) containment alone or together with a fast resorbing polyglycolide/polylactide (PGA/PLA) binding substance were compared to the conventional non-contained method in ridge augmentation in sheep. The contained methods decreased HA migration, but the augmentations did not differ significantly. The use of the containment caused a risk for wound dehiscence and infection. Histologically there was a rapid connective tissue ingrowth into the HA graft and it was more abundant with the PGA containment compared to the non-contained augmentation and even additionally rich when the HA particles were bound with PGA/PLA copolymer. However, the bone ingrowth was best in the non-contained augmentation exceeding 10-12 % of the total graft area at 24 weeks. Negligible or no bone ingrowth was seen in the cases where the polymer composite was added to the HA particles and, related to that, foreign-body type cells were seen at the interface between the HA and host bone. The PGA and poly-dl/l-lactide (PDLLA) fibre-reinforced coralline HA blocks were studied in the metaphyseal and in the diaphyseal defects in rabbits. A rapid bone ingrowth was seen inside the both types of implants. Both PGA and PDLLA fibres induced an inflammatory fibrous reaction around themselves but it did not hinder the bone ingrowth. The bone ingrowth pattern was directed according to the loading conditions so that the load-carrying cortical ends of the implants as well as the implants sited in the diaphyseal defects were the most ossified. The fibre-reinforced coralline HA implants were further studied as stand-alone grafts in the lumbar anterior interbody implantation in pigs. The strength of the HA implants proved not to be adequate, the implants fractured in six weeks and the disc space was gradually lost similarly to that of the discectomized spaces. Histologically, small quantities of bone ingrowth was seen in some of the PGA and PDLLA reinforced coralline implants while no bone formation was identified in any of the PDLLA reinforced synthetic porous HA implants. While fragmented, the inner structure of the implants was lost, the bone ingrowth was minimal, and the disc was replaced by the fibrous connective tissue. When evaluated radiologically the grade of ossification was assessed as better than histologically, and, when related to the histologic findings, CT was more dependable than the plain films to show ossification of the implanted disc space. Local kyphosis was a frequent finding along with anterior bone bridging and ligament ossification as a consequence of instability of the implanted segment.