Here we describe the SweGen data set, a comprehensive map of genetic variation in the Swedish population. These data represent a basic resource for clinical genetics laboratories as well as for sequencing-based association studies by providing information on genetic variant frequencies in a cohort that is well matched to national patient cohorts. To select samples for this study, we first examined the genetic structure of the Swedish population using high-density SNP-array data from a nation-wide cohort of over 10 000 Swedish-born individuals included in the Swedish Twin Registry. A total of 1000 individuals, reflecting a cross-section of the population and capturing the main genetic structure, were selected for whole-genome sequencing. Analysis pipelines were developed for automated alignment, variant calling and quality control of the sequencing data. This resulted in a genome-wide collection of aggregated variant frequencies in the Swedish population that we have made available to the scientific community through the website https://swefreq.nbis.se. A total of 29.2 million single-nucleotide variants and 3.8 million indels were detected in the 1000 samples, with 9.9 million of these variants not present in current databases. Each sample contributed with an average of 7199 individual-specific variants. In addition, an average of 8645 larger structural variants (SVs) were detected per individual, and we demonstrate that the population frequencies of these SVs can be used for efficient filtering analyses. Finally, our results show that the genetic diversity within Sweden is substantial compared with the diversity among continental European populations, underscoring the relevance of establishing a local reference data set.

Bartonella henselae is a natural intracellular colonizer of cats, and is transferred by blood-sucking insect vectors. It is also an opportunistic human pathogen. Two strains of B. henselae, thought to be representative of the diversity of the species, were selected for low-coverage 454 sequencing. The comparison of these two strains to the published Houston-1 reveals very high nucleotide identity and low substitution and recombination, with the remarkable exception of phages and host-interaction genes such as type IV and V secretion systems. Among the few variable genes of unknown function, BH14680, an alpha-Proteobacteria-specific gene, shows faster evolution in Bartonella compared to other alpha-Proteobacteria. Its 5’ end, which is likely coding for a domain exposed extracellularly, is under positive or very relaxed selection, and might be involved in host-interaction processes. Finally, we show that a simple genome coverage analysis reveal major genomic events such as duplications and unusual replication modes, such as the run-off replication. The latter, combined with a gene transfer agent, is thought to be a novel way to increase substitution and recombination frequencies. An extensive analysis of all bacterial pyrosequencing projects showed that it is probably Bartonella-specific.

Gene transfer agents (GTAs) randomly transfer short fragments of a bacterial genome. A novel putative GTA was recently discovered in the mouse-infecting bacterium Bartonella grahamii. Although GTAs are widespread in phylogenetically diverse bacteria, their role in evolution is largely unknown. Here, we present a comparative analysis of 16 Bartonella genomes ranging from 1.4 to 2.6 Mb in size, including six novel genomes from Bartonella isolated from a cow, two moose, two dogs, and a kangaroo. A phylogenetic tree inferred from 428 orthologous core genes indicates that the deadly human pathogen B. bacilliformis is related to the ruminant-adapted clade, rather than being the earliest diverging species in the genus as previously thought. A gene flux analysis identified 12 genes for a GTA and a phage-derived origin of replication as the most conserved innovations. These are located in a region of a few hundred kb that also contains 8 insertions of gene clusters for type III, IV, and V secretion systems, and genes for putatively secreted molecules such as cholera-like toxins. The phylogenies indicate a recent transfer of seven genes in the virB gene cluster for a type IV secretion system from a catadapted B. henselae to a dog-adapted B. vinsonii strain. We show that the B. henselae GTA is functional and can transfer genes in vitro. We suggest that the maintenance of the GTA is driven by selection to increase the likelihood of horizontal gene transfer and argue that this process is beneficial at the population level, by facilitating adaptive evolution of the host-adaptation systems and thereby expansion of the host range size. The process counters gene loss and forces all cells to contribute to the production of the GTA and the secreted molecules. The results advance our understanding of the role that GTAs play for the evolution of bacterial genomes.

Background: The polytene nuclei of the dipteran Chironomus tentans (Ch. tentans) with their Balbiani ring (BR) genes constitute an exceptional model system for studies of the expression of endogenous eukaryotic genes. Here, we report the first draft genome of Ch. tentans and characterize its gene expression machineries and genomic architecture of the BR genes. Results: The genome of Ch. tentans is approximately 200 Mb in size, and has a low GC content (31%) and a low repeat fraction (15%) compared to other Dipteran species. Phylogenetic inference revealed that Ch. tentans is a sister clade to mosquitoes, with a split 150-250 million years ago. To characterize the Ch. tentans gene expression machineries, we identified potential orthologus sequences to more than 600 Drosophila melanogaster (D. melanogaster) proteins involved in the expression of protein-coding genes. We report novel data on the organization of the BR gene loci, including a novel putative BR gene, and we present a model for the organization of chromatin bundles in the BR2 puff based on genic and intergenic in situ hybridizations. Conclusions: We show that the molecular machineries operating in gene expression are largely conserved between Ch. tentans and D. melanogaster, and we provide enhanced insight into the organization and expression of the BR genes. Our data strengthen the generality of the BR genes as a unique model system and provide essential background for in-depth studies of the biogenesis of messenger ribonucleoprotein complexes.

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.

Many bacteria live in close association with other organisms such as plants and animals, with important implications for both health and disease. This thesis investigates bacteria that are well adapted to live inside an animal host, and describes the molecular evolutionary processes underlying host-adaptation, based on bacterial genome comparisons.

Insect-transmitted bacteria of the genus Bartonella infect the red blood cells of mammals, and we investigate host adaptation and genome evolution in this genus. In Bartonella, many host-interaction systems are encoded in a highly variable chromosomal segment previously shown to be amplified and packaged into bacteriophage particles. Among all genes imported into the Bartonella ancestor, we identify the short gene cluster encoding these phage particles as the most evolutionary conserved, indicating a strong selective advantage and a role in niche adaptation. We also provide an overview of the remarkable evolutionary dynamics of type IV and type V secretion systems, including a detailed analysis of the type IV secretion system trw. Our results highlight the importance of recombination and gene conversion in the evolution of host-adaptation systems, and reveal how these mutational mechanisms result in strikingly different outcomes depending on the selective constraints.

In the insect endosymbionts Buchnera and Blochmannia, we show that genes frameshifted at poly(A) tracts can remain functional due to transcriptional slippage. Selection against poly(A) tracts is very inefficient in these genomes compared to other bacteria, and we discuss why this can lead to increased rates of gene loss. Using the human pathogen Helicobacter pylori as a model, we provide a deeper understanding of why highly expressed genes evolve slowly.

This thesis emphasizes the power of using complete genome sequences to study evolutionary processes. In particular, we argue that knowledge about the complex evolution of duplicated gene segments is crucial to understand host adaptation in bacteria.

We have studied the evolution of a type TV secretion system (T4SS), in Bartonella, which is thought to have changed function from conjugation to erythrocyte adherence following a recent horizontal gene transfer event. The system, called Trw, is unique among T4SSs in that genes encoding both exo- and intracellular components are located within the same duplicated fragment. This provides an opportunity to study the influence of selection on proteins involved in host-pathogen interactions. We sequenced the trw locus from several strains of Bartonella henselae and investigated its evolutionary history by comparisons to other Bartonella species. Several instances of recombination and gene conversion events where detected in the 2- to 5-fold duplicated gene fragments encompassing trwJIH, explaining the homogenization of the anchoring protein TrwI and the divergence of the minor pilus protein TrwJ. A phylogenetic analysis of the 7- to 8-fold duplicated gene coding for the major pilus protein TrwL displayed 2 distinct clades, likely representing a subfunctionalization event. The analyses of the B. henselae strains also identified a recent horizontal transfer event of almost the complete trwL region. We suggest that the switch in function of the T4SS was mediated by the duplication of the genes encoding pilus components and their diversification by combinatorial sequence shuffling within and among genomes. We suggest that the pilus proteins have evolved by diversifying selection to match a divergent set of erythrocyte surface structures, consistent with the trench warfare coevolutionary model.

Bacteria of the genus Bartonella are facultative intracellular bacteria infecting the red blood cells of mammals. Bartonella isolates have now been reported from a wide range of mammalian host species, including humans, domestic animals such as pets and livestock, as well as many wild animals such as deer, moose, kangaroo, and whales. Here, we present the first major genus-wide investigation of host-adaptation systems in Bartonella, using 5 published and 5 draft genome sequences. The sampling includes both clinical and natural isolates, and represent well the major phylogenetic diversity of the genus. Our study reveals four distinct protein families of Type V Secretion Systems (T5SS) shared by all sequenced members of the genus. We also show that a recently identified gene transfer agent (GTA) consisting of a defective phage is, surprisingly, the most conserved gene cluster among all Bartonella-specific or imported genes, strongly emphasizing the functional importance of this system for the life-style and evolution of Bartonella.

Vascular and haematopoietic cells organize into specialized tissues during early embryogenesis to supply essential nutrients to all organs and thus play critical roles in development and disease. At the top of the haemato-vascular specification cascade lies cloche, a gene that when mutated in zebrafish leads to the striking phenotype of loss of most endothelial and haematopoietic cells(1-4) and a significant increase in cardiomyocyte numbers(5). Although this mutant has been analysed extensively to investigate mesoderm diversification and differentiation(1-7) and continues to be broadly used as a unique avascular model, the isolation of the cloche gene has been challenging due to its telomeric location. Here we used a deletion allele of cloche to identify several new cloche candidate genes within this genomic region, and systematically genome-edited each candidate. Through this comprehensive interrogation, we succeeded in isolating the cloche gene and discovered that it encodes a PAS-domain-containing bHLH transcription factor, and that it is expressed in a highly specific spatiotemporal pattern starting during late gastrulation. Gain-of-function experiments show that it can potently induce endothelial gene expression. Epistasis experiments reveal that it functions upstream of etv2 and tal1, the earliest expressed endothelial and haematopoietic transcription factor genes identified to date. A mammalian cloche orthologue can also rescue blood vessel formation in zebrafish cloche mutants, indicating a highly conserved role in vertebrate vasculogenesis and haematopoiesis. The identification of this master regulator of endothelial and haematopoietic fate enhances our understanding of early mesoderm diversification and may lead to improved protocols for the generation of endothelial and haematopoietic cells in vivo and in vitro.

More than 20 α-proteobacterial genomes are currently available. These range in size from 1–9 Mb and represent excellent model systems for evolutionary studies of the organizational features of bacterial genomes. Computational inferences have shown that genome reductions have occurred independently in lineages such as Rickettsia and Bartonella that are associated with intracellular lifestyles. Analyses of these reduced genomes have provided insights into the evolution of vector-borne transmission pathways. Further research into the population biology of bacteria, arthropods and vertebrate hosts will help to refine the biology of host–pathogen interactions and will facilitate the design of vaccines and vector-control programs.

More than 20 alpha-proteobacterial genomes are currently available. These range in size from 1–9 Mb and represent excellent model systems for evolutionary studies of the organizational features of bacterial genomes. Computational inferences have shown that genome reductions have occurred independently in lineages such as Rickettsia and Bartonella that are associated with intracellular lifestyles. Analyses of these reduced genomes have provided insights into the evolution of vector-borne transmission pathways. Further research into the population biology of bacteria, arthropods and vertebrate hosts will help to refine the biology of host–pathogen interactions and will facilitate the design of vaccines and vector-control programs.

In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites.

In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites.

This study presents development and behaviour of a feedback turnover model that mimics asymmetric circadian oscillations of body temperature, blood pressure and heart rate in rats.The study also includes an application to drug-induced hypothermia, tolerance and handling effects. Data were collected inn normotensive Sprague-Dawley rats, housed at 25 degrees C with a 12:12 hr light dark cycle (light on at 06:00 am) and with free access of food and water. The model consisted of two intertwined parallel compartments which captured a free-running rhythm with a period close to but not exactly 24 hrs. The free-running rhythm was synchronised to exactly 24 hrs by the environmental timekeeper (12:12 hr light on/off cycle) in experimental settings. The baseline model was fitted to a standardised 24-hr period derived from mean data of six animals over a period of nine consecutive days. The first-order rate constants related to the turnover of the baseline temperature, alpha and beta, were 0.026 min(-1) (+/-5%) and 0.0037 min(-1) (+/-3%). The alpha and beta parameters are approximately 2/transition time between day and night and 2/night time, respectively. The day:night timekeeper g(t), reference point T(ref) and amplitude were 0.053(+/-2%),37.3(+/-0.02%) and 3.3% (+/-2%), respectively. Simulations with the baseline model revealed stable oscillations (free-running rhythm) in the absence of the timekeeper. This temperature-time profile was then symmetric and had a smaller amplitude, with a slightly shorter period and less pronounced temperature shift as compared to the profile in the presence of an external Timekeeper. Fitting the model to 96 hr mean profiles of blood pressure and heart rate from 10 control animals demonstrated the usefulness of the model.Simulations of the integrated temperature model succeeded in mimicking other modes of administration such as oral dosing.

The aim of the present investigation was to develop a pharmacokinetic-pharmacodynamic model for the characterization of clomethiazole (CMZ)-induced hypothermia and the rapid development of long-lasting tolerance in rats while taking into account circadian rhythm in baseline and the influence of handling. CMZ-induced hypothermia and tolerance was measured using body temperature telemetry in male Sprague-Dawley rats, which were given s.c. bolus injections of 0, 15, 150, 300, and 600 micromol kg(-1) and 24-h s.c. continuous infusions of 0, 20, and 40 micromol kg(-1) h(-1) using osmotic pumps. The duration of tolerance was studied by repeated injections of 300 micromol kg(-1) at 3- to 32-day intervals. Plasma exposure to CMZ was obtained in satellite groups of catheterized rats. Fitted population concentration-time profiles served as input for the pharmacodynamic analysis. The asymmetric circadian rhythm in baseline body temperature was successfully described by a novel negative feedback model incorporating external light-dark conditions. An empirical function characterized the transient increase in temperature upon handling of the animal. A feedback model for temperature regulation and tolerance development allowed estimation of CMZ potency at 30 +/- 1 microM. The delay in onset of tolerance was estimated via a series of four transit compartments at 7.6 +/- 2 h. The long-lasting tolerance was assumed to be caused by inactivation of a mediator with an estimated turnover time of 46 +/- 3 days. This multicomponent turnover model was able to quantify the CMZ-induced hypothermia, circadian rhythm in baseline, and rapid onset of a long-lasting tolerance to CMZ in rats.

Complete and accurate genome assembly and annotation is a crucial foundation for comparative and functional genomics. Despite this, few complete eukaryotic genomes are available, and genome annotation remains a major challenge. Here, we present a complete genome assembly of the skin commensal yeast Malassezia sympodialis and demonstrate how proteogenomics can substantially improve gene an-notation. Through long-read DNA sequencing, we obtained a gap-free genome assembly for M. sympodi-alis (ATCC 42132), comprising eight nuclear and one mitochondrial chromosome. We also sequenced and assembled four M. sympodialis clinical isolates, and showed their value for understanding Malassezia reproduction by confirming four alternative allele combinations at the two mating-type loci. Importantly, we demonstrated how proteomics data could be readily integrated with transcriptomics data in standard annotation tools. This increased the number of annotated protein-coding genes by 14% (from 3612 to 4113), compared to using transcriptomics evidence alone. Manual curation further increased the number of protein-coding genes by 9% (to 4493). All of these genes have RNA-seq evidence and 87% were confirmed by proteomics. The M. sympodialis genome assembly and annotation presented here is at a quality yet achieved only for a few eukaryotic organisms, and constitutes an important reference for future host-microbe interaction studies.