MINCR is a MYC-induced lncRNA able to modulate MYC's transcriptional network in Burkitt lymphoma cells.

Despite the established role of the transcription factor MYC in cancer, little is known about the impact of a new class of transcriptional regulators, the long noncoding RNAs (lncRNAs), on MYC ability to influence the cellular transcriptome. Here, we have intersected RNA-sequencing data from two MYC-inducible cell lines and a cohort of 91 B-cell lymphomas with or without genetic variants resulting in MYC overexpression. We identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them, we focused on a lncRNA that we named MYC-induced long noncoding RNA (MINCR), showing a strong correlation with MYC expression in MYC-positive lymphomas. To understand its cellular role, we performed RNAi and found that MINCR knockdown is associated with an impairment in cell cycle progression. Differential gene expression analysis after RNAi showed a significant enrichment of cell cycle genes among the genes down-regulated after MINCR knockdown. Interestingly, these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of the MYC transcriptional program. Accordingly, MINCR knockdown was associated with a reduction in MYC binding to the promoters of selected cell cycle genes. Finally, we show that down-regulation of Aurora kinases A and B and chromatin licensing and DNA replication factor 1 may explain the reduction in cellular proliferation observed on MINCR knockdown. We, therefore, suggest that MINCR is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

Truncating and missense mutations in IGHMBP2 cause Charcot-Marie Tooth disease type 2.

Using a combination of exome sequencing and linkage analysis, we investigated an English family with two affected siblings in their 40s with recessive Charcot-Marie Tooth disease type 2 (CMT2). Compound heterozygous mutations in the immunoglobulin-helicase-µ-binding protein 2 (IGHMBP2) gene were identified. Further sequencing revealed a total of 11 CMT2 families with recessively inherited IGHMBP2 gene mutations. IGHMBP2 mutations usually lead to spinal muscular atrophy with respiratory distress type 1 (SMARD1), where most infants die before 1 year of age. The individuals with CMT2 described here, have slowly progressive weakness, wasting and sensory loss, with an axonal neuropathy typical of CMT2, but no significant respiratory compromise. Segregating IGHMBP2 mutations in CMT2 were mainly loss-of-function nonsense in the 5' region of the gene in combination with a truncating frameshift, missense, or homozygous frameshift mutations in the last exon. Mutations in CMT2 were predicted to be less aggressive as compared to those in SMARD1, and fibroblast and lymphoblast studies indicate that the IGHMBP2 protein levels are significantly higher in CMT2 than SMARD1, but lower than controls, suggesting that the clinical phenotype differences are related to the IGHMBP2 protein levels.

Deciphering the complex nature of bolting time regulation in Beta vulgaris.

KEY MESSAGE: Only few genetic loci are sufficient to increase the variation of bolting time in Beta vulgaris dramatically, regarding vernalization requirement, seasonal bolting time and reproduction type. Beta species show a wide variation of bolting time regarding the year of first reproduction, seasonal bolting time and the number of reproduction cycles. To elucidate the genetics of bolting time control, we used three F3 mapping populations that were produced by crossing a semelparous, annual sugar beet with iteroparous, vernalization-requiring wild beet genotypes. The semelparous plants died after reproduction, whereas iteroparous plants reproduced at least twice. All populations segregated for vernalization requirement, seasonal bolting time and the number of reproduction cycles. We found that vernalization requirement co-segregated with the bolting locus B on chromosome 2 and was inherited independently from semel- or iteroparous reproduction. Furthermore, we found that seasonal bolting time is a highly heritable trait (h 2 > 0.84), which is primarily controlled by two major QTL located on chromosome 4 and 9. Late bolting alleles of both loci act in a partially recessive manner and were identified in both iteroparous pollinators. We observed an additive interaction of both loci for bolting delay. The QTL region on chromosome 4 encompasses the floral promoter gene BvFT2, whereas the QTL on chromosome 9 co-localizes with the BR 1 locus, which controls post-winter bolting resistance. Our findings are applicable for marker-assisted sugar beet breeding regarding early bolting to accelerate generation cycles and late bolting to develop bolting-resistant spring and winter beets. Unexpectedly, one population segregated also for dwarf growth that was found to be controlled by a single locus on chromosome 9.

Alterations of microRNA and microRNA-regulated messenger RNA expression in germinal center B-cell lymphomas determined by integrative sequencing analysis.

MicroRNA are well-established players in post-transcriptional gene regulation. However, information on the effects of microRNA deregulation mainly relies on bioinformatic prediction of potential targets, whereas proof of the direct physical microRNA/target messenger RNA interaction is mostly lacking. Within the International Cancer Genome Consortium Project "Determining Molecular Mechanisms in Malignant Lymphoma by Sequencing", we performed miRnome sequencing from 16 Burkitt lymphomas, 19 diffuse large B-cell lymphomas, and 21 follicular lymphomas. Twenty-two miRNA separated Burkitt lymphomas from diffuse large B-cell lymphomas/follicular lymphomas, of which 13 have shown regulation by MYC. Moreover, we found expression of three hitherto unreported microRNA. Additionally, we detected recurrent mutations of hsa-miR-142 in diffuse large B-cell lymphomas and follicular lymphomas, and editing of the hsa-miR-376 cluster, providing evidence for microRNA editing in lymphomagenesis. To interrogate the direct physical interactions of microRNA with messenger RNA, we performed Argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation experiments. MicroRNA directly targeted 208 messsenger RNA in the Burkitt lymphomas and 328 messenger RNA in the non-Burkitt lymphoma models. This integrative analysis discovered several regulatory pathways of relevance in lymphomagenesis including Ras, PI3K-Akt and MAPK signaling pathways, also recurrently deregulated in lymphomas by mutations. Our dataset reveals that messenger RNA deregulation through microRNA is a highly relevant mechanism in lymphomagenesis.

The PCBP1 gene encoding poly(rC) binding protein I is recurrently mutated in Burkitt lymphoma.

The genetic hallmark of Burkitt lymphoma is the translocation t(8;14)(q24;q32), or one of its light chain variants, resulting in IG-MYC juxtaposition. However, these translocations alone are insufficient to drive lymphomagenesis, which requires additional genetic changes for malignant transformation. Recent studies of Burkitt lymphoma using next generation sequencing approaches have identified various recurrently mutated genes including ID3, TCF3, CCND3, and TP53. Here, by using similar approaches, we show that PCBP1 is a recurrently mutated gene in Burkitt lymphoma. By whole-genome sequencing, we identified somatic mutations in PCBP1 in 3/17 (18%) Burkitt lymphomas. We confirmed the recurrence of PCBP1 mutations by Sanger sequencing in an independent validation cohort, finding mutations in 3/28 (11%) Burkitt lymphomas and in 6/16 (38%) Burkitt lymphoma cell lines. PCBP1 is an intron-less gene encoding the 356 amino acid poly(rC) binding protein 1, which contains three K-Homology (KH) domains and two nuclear localization signals. The mutations predominantly (10/12, 83%) affect the KH III domain, either by complete domain loss or amino acid changes. Thus, these changes are predicted to alter the various functions of PCBP1, including nuclear trafficking and pre-mRNA splicing. Remarkably, all six primary Burkitt lymphomas with a PCBP1 mutation expressed MUM1/IRF4, which is otherwise detected in around 20-40% of Burkitt lymphomas. We conclude that PCBP1 mutations are recurrent in Burkitt lymphomas and might contribute, in cooperation with other mutations, to its pathogenesis.

Genetic differentiation of hypothalamus parentally biased transcripts in populations of the house mouse implicate the Prader-Willi syndrome imprinted region as a possible source of behavioral divergence.

Parentally biased expression of transcripts (genomic imprinting) in adult tissues, including the brain, can influence and possibly drive the evolution of behavioral traits. We have previously found that paternally determined cues are involved in population-specific mate choice decisions between two populations of the Western house mouse (Mus musculus domesticus). Here, we ask whether this could be mediated by genomically imprinted transcripts that are subject to fast differentiation between these populations. We focus on three organs that are of special relevance for mate choice and behavior: The vomeronasal organ (VNO), the hypothalamus, and the liver. To first identify candidate transcripts at a genome-wide scale, we used reciprocal crosses between M. m. domesticus and M. m. musculus inbred strains and RNA sequencing of the respective tissues. Using a false discovery cutoff derived from mock reciprocal cross comparisons, we find a total of 66 imprinted transcripts, 13 of which have previously not been described as imprinted. The largest number of imprinted transcripts were found in the hypothalamus; fewer were found in the VNO, and the least were found in the liver. To assess molecular differentiation and imprinting in the wild-derived M. m. domesticus populations, we sequenced the RNA of the hypothalamus from individuals of these populations. This confirmed the presence of the above identified transcripts also in wild populations and allowed us to search for those that show a high genetic differentiation between these populations. Our results identify the Ube3a-Snrpn imprinted region on chromosome 7 as a region that encompasses the largest number of previously not described transcripts with paternal expression bias, several of which are at the same time highly differentiated. For four of these, we confirmed their imprinting status via single nucleotide polymorphism-specific pyrosequencing assays with RNA from reciprocal crosses. In addition, we find the paternally expressed Peg13 transcript within the Trappc9 gene region on chromosome 15 to be highly differentiated. Interestingly, both regions have been implicated in Prader-Willi nervous system disorder phenotypes in humans. We suggest that these genomically imprinted regions are candidates for influencing the population-specific mate-choice in mice.

Ultrashort and progressive 4sU-tagging reveals key characteristics of RNA processing at nucleotide resolution.

RNA synthesis and decay rates determine the steady-state levels of cellular RNAs. Metabolic tagging of newly transcribed RNA by 4-thiouridine (4sU) can reveal the relative contributions of RNA synthesis and decay rates. The kinetics of RNA processing, however, had so far remained unresolved. Here, we show that ultrashort 4sU-tagging not only provides snapshot pictures of eukaryotic gene expression but, when combined with progressive 4sU-tagging and RNA-seq, reveals global RNA processing kinetics at nucleotide resolution. Using this method, we identified classes of rapidly and slowly spliced/degraded introns. Interestingly, each class of splicing kinetics was characterized by a distinct association with intron length, gene length, and splice site strength. For a large group of introns, we also observed long lasting retention in the primary transcript, but efficient secondary splicing or degradation at later time points. Finally, we show that processing of most, but not all small nucleolar (sno)RNA-containing introns is remarkably inefficient with the majority of introns being spliced and degraded rather than processed into mature snoRNAs. In summary, our study yields unparalleled insights into the kinetics of RNA processing and provides the tools to study molecular mechanisms of RNA processing and their contribution to the regulation of gene expression.

From next-generation sequencing alignments to accurate comparison and validation of single-nucleotide variants: the pibase software.

Scientists working with single-nucleotide variants (SNVs), inferred by next-generation sequencing software, often need further information regarding true variants, artifacts and sequence coverage gaps. In clinical diagnostics, e.g. SNVs must usually be validated by visual inspection or several independent SNV-callers. We here demonstrate that 0.5-60% of relevant SNVs might not be detected due to coverage gaps, or might be misidentified. Even low error rates can overwhelm the true biological signal, especially in clinical diagnostics, in research comparing healthy with affected cells, in archaeogenetic dating or in forensics. For these reasons, we have developed a package called pibase, which is applicable to diploid and haploid genome, exome or targeted enrichment data. pibase extracts details on nucleotides from alignment files at user-specified coordinates and identifies reproducible genotypes, if present. In test cases pibase identifies genotypes at 99.98% specificity, 10-fold better than other tools. pibase also provides pair-wise comparisons between healthy and affected cells using nucleotide signals (10-fold more accurately than a genotype-based approach, as we show in our case study of monozygotic twins). This comparison tool also solves the problem of detecting allelic imbalance within heterozygous SNVs in copy number variation loci, or in heterogeneous tumor sequences.

Pyrosequencing of supra- and subgingival biofilms from inflamed peri-implant and periodontal sites.

BACKGROUND: To investigate the microbial composition of biofilms at inflamed peri-implant and periodontal tissues in the same subject, using 16S rRNA sequencing. METHODS: Supra- and submucosal, and supra- and subgingival plaque samples were collected from 7 subjects suffering from diseased peri-implant and periodontal tissues. Bacterial DNA was isolated and 16S rRNA genes were amplified, sequenced and aligned for the identification of bacterial genera. RESULTS: 43734 chimera-depleted, denoised sequences were identified, corresponding to 1 phylum, 8 classes, 10 orders, 44 families and 150 genera. The most abundant families or genera found in supramucosal or supragingival plaque were Streptoccocaceae, Rothia and Porphyromonas. In submucosal plaque, the most abundant family or genera found were Rothia, Streptococcaceae and Porphyromonas on implants. The most abundant subgingival bacteria on teeth were Prevotella, Streptococcaceae, and TG5. The number of sequences found for the genera Tannerella and Aggregatibacter on implants differed significantly between supra- and submucosal locations before multiple testing. The analyses demonstrated no significant differences between microbiomes on implants and teeth in supra- or submucosal and supra- or subgingival biofilms. CONCLUSION: Diseased peri-implant and periodontal tissues in the same subject share similiar bacterial genera and based on the analysis of taxa on a genus level biofilm compositions may not account for the potentially distinct pathologies at implants or teeth.

Genome sequence of Comamonas testosteroni ATCC 11996, a representative strain involved in steroid degradation.

Comamonas testosteroni strains belong to the family of Comamonadaceae and are known for their ability to utilize steroid compounds as carbon source. Here, we present the draft genome sequence of strain ATCC 11996, with a G+C content of 61.48%.

Accurate variant detection across non-amplified and whole genome amplified DNA using targeted next generation sequencing.

BACKGROUND: Many hypothesis-driven genetic studies require the ability to comprehensively and efficiently target specific regions of the genome to detect sequence variations. Often, sample availability is limited requiring the use of whole genome amplification (WGA). We evaluated a high-throughput microdroplet-based PCR approach in combination with next generation sequencing (NGS) to target 384 discrete exons from 373 genes involved in cancer. In our evaluation, we compared the performance of six non-amplified gDNA samples from two HapMap family trios. Three of these samples were also preamplified by WGA and evaluated. We tested sample pooling or multiplexing strategies at different stages of the tested targeted NGS (T-NGS) workflow. RESULTS: The results demonstrated comparable sequence performance between non-amplified and preamplified samples and between different indexing strategies [sequence specificity of 66.0% ± 3.4%, uniformity (coverage at 0.2× of the mean) of 85.6% ± 0.6%]. The average genotype concordance maintained across all the samples was 99.5% ± 0.4%, regardless of sample type or pooling strategy. We did not detect any errors in the Mendelian patterns of inheritance of genotypes between the parents and offspring within each trio. We also demonstrated the ability to detect minor allele frequencies within the pooled samples that conform to predicted models. CONCLUSION: Our described PCR-based sample multiplex approach and the ability to use WGA material for NGS may enable researchers to perform deep resequencing studies and explore variants at very low frequencies and cost.

A tissue-specific landscape of sense/antisense transcription in the mouse intestine.

BACKGROUND: The intestinal mucosa is characterized by complex metabolic and immunological processes driven highly dynamic gene expression programs. With the advent of next generation sequencing and its utilization for the analysis of the RNA sequence space, the level of detail on the global architecture of the transcriptome reached a new order of magnitude compared to microarrays. RESULTS: We report the ultra-deep characterization of the polyadenylated transcriptome in two closely related, yet distinct regions of the mouse intestinal tract (small intestine and colon). We assessed tissue-specific transcriptomal architecture and the presence of novel transcriptionally active regions (nTARs). In the first step, signatures of 20,541 NCBI RefSeq transcripts could be identified in the intestine (74.1% of annotated genes), thereof 16,742 are common in both tissues. Although the majority of reads could be linked to annotated genes, 27,543 nTARs not consistent with current gene annotations in RefSeq or ENSEMBL were identified. By use of a second independent strand-specific RNA-Seq protocol, 20,966 of these nTARs were confirmed, most of them in vicinity of known genes. We further categorized our findings by their relative adjacency to described exonic elements and investigated regional differences of novel transcribed elements in small intestine and colon. CONCLUSIONS: The current study demonstrates the complexity of an archetypal mammalian intestinal mRNA transcriptome in high resolution and identifies novel transcriptionally active regions at strand-specific, single base resolution. Our analysis for the first time shows a strand-specific comparative picture of nTARs in two tissues and represents a resource for further investigating the transcriptional processes that contribute to tissue identity.

Comprehensive analysis of candidate genes for photosensitivity using a complementary bioinformatic and experimental approach.

Photoparoxysmal response (PPR) is a highly heritable electroencephalographic trait characterized by an increased sensitivity to photic stimulation. It may serve as an endophenotype for idiopathic generalized epilepsy. Family linkage studies identified susceptibility loci for PPR on chromosomes 5q35.3, 8q21.13, and 16p13.3. This study aimed to identify key candidate genes within these loci. We used bioinformatics tools for gene prioritization integrating information on biologic function, sequence data, gene expression, and others. The prime candidate gene from this analysis was sequenced in 48 photopositive probands. Presumed functional implications of identified polymorphisms were investigated using bioinformatics methods. The glutamate receptor subunit gene GRIN2A was identified as a prime candidate gene. Sequence analysis revealed various new polymorphisms. None of the identified variants was predicted to be functionally relevant. We objectified the selection of candidate genes for PPR without an a priori hypothesis. Particularly among the various ion channel genes in the linkage regions, GRIN2A was identified as the prime candidate gene. GRIN2A mutations have recently been identified in various epilepsies. Even though our mutation analysis failed to demonstrate direct involvement of GRIN2A in photosensitivity, in silico gene prioritization may provide a useful tool for the identification of candidate genes within large genomic regions.

SNP discovery performance of two second-generation sequencing platforms in the NOD2 gene region.

A potentially important application of second generation sequencing technologies is to identify disease-associated variation. For comparison of the performance in SNP detection, the Crohn's disease (CD)-associated NOD2 gene was subjected to targeted resequencing using two different second-generation sequencing technologies. Eleven CD patients were selected based on their haplotype background at the NOD2 locus. The 40-kb large NOD2 gene region was amplified using long-range PCR (LR-PCR), and sequenced with the Roche 454/FLX system, an Applied Biosystems SOLiD mate-pair library (2 x 25 bp), and a SOLiD fragment (50 bp) library. The entire NOD2 region was also sequenced using conventional Sanger technology. Four-hundred forty-two single nucleotide polymorphisms (SNPs) were discovered with the SOLiD mate-pair library, 454 with the fragment library, and 441 with the 454/FLX. For the homozygous SNPs, 98% were confirmed by Sanger for the mate-pair library, 100% for the fragment library and 99% for the 454/FLX. Ninety-six percent of the heterozygous SNPs detected with the SOLiD mate-pair library, 91% with the fragment library and 96% with the 454/FLX were confirmed by Sanger. In a simulation, the SNP detection performance fell rapidly when the achieved coverage was below 40 x. Due to uneven representation of the target region when using LR-PCR, oversequencing of other regions is necessary.

Recent transfer of an iron-regulated gene from the plastid to the nuclear genome in an oceanic diatom adapted to chronic iron limitation.

BACKGROUND: Although the importance and widespread occurrence of iron limitation in the contemporary ocean is well documented, we still know relatively little about genetic adaptation of phytoplankton to these environments. Compared to its coastal relative Thalassiosira pseudonana, the oceanic diatom Thalassiosira oceanica is highly tolerant to iron limitation. The adaptation to low-iron conditions in T. oceanica has been attributed to a decrease in the photosynthetic components that are rich in iron. Genomic information on T. oceanica may shed light on the genetic basis of the physiological differences between the two species. RESULTS: The complete 141790 bp sequence of the T. oceanica chloroplast genome [GenBank: GU323224], assembled from massively parallel pyrosequencing (454) shotgun reads, revealed that the petF gene encoding for ferredoxin, which is localized in the chloroplast genome in T. pseudonana and other diatoms, has been transferred to the nucleus in T. oceanica. The iron-sulfur protein ferredoxin, a key element of the chloroplast electron transport chain, can be replaced by the iron-free flavodoxin under iron-limited growth conditions thereby contributing to a reduction in the cellular iron requirements. From a comparison to the genomic context of the T. pseudonana petF gene, the T. oceanica ortholog can be traced back to its chloroplast origin. The coding potential of the T. oceanica chloroplast genome is comparable to that of T. pseudonana and Phaeodactylum tricornutum, though a novel expressed ORF appears in the genomic region that has been subjected to rearrangements linked to the petF gene transfer event. CONCLUSIONS: The transfer of the petF from the cp to the nuclear genome in T. oceanica represents a major difference between the two closely related species. The ability of T. oceanica to tolerate iron limitation suggests that the transfer of petF from the chloroplast to the nuclear genome might have contributed to the ecological success of this species.

Nucleotide divergence vs. gene expression differentiation: comparative transcriptome sequencing in natural isolates from the carrion crow and its hybrid zone with the hooded crow.

Recent advances in sequencing technology promise to provide new strategies for studying population differentiation and speciation phenomena in their earliest phases. We focus here on the black carrion crow (Corvus [corone] corone), which forms a zone of hybridization and overlap with the grey coated hooded crow (Corvus [corone] cornix). However, although these semispecies are taxonomically distinct, previous analyses based on several types of genetic markers did not reveal significant molecular differentiation between them. We here corroborate this result with sequence data obtained from a set of 25 nuclear intronic loci. Thus, the system represents a case of a very early phase of species divergence that requires new molecular approaches for its description. We have therefore generated RNAseq expression profiles using barcoded massively parallel pyrosequencing of brain mRNA from six individuals of the carrion crow and five individuals from a hybrid zone with the hooded crow. We obtained 856 675 reads from two runs, with average read length of 270 nt and coverage of 8.44. Reads were assembled de novo into 19 552 contigs, 70% of which could be assigned to annotated genes in chicken and zebra finch. This resulted in a total of 7637 orthologous genes and a core set of 1301 genes that could be compared across all individuals. We find a clear clustering of expression profiles for the pure carrion crow animals and disperse profiles for the animals from the hybrid zone. These results suggest that gene expression differences may indeed be a sensitive indicator of initial species divergence.

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma.

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

GenColors: annotation and comparative genomics of prokaryotes made easy.

GenColors (gencolors.fli-leibniz.de) is a new web-based software/database system aimed at an improved and accelerated annotation of prokaryotic genomes considering information on related genomes and making extensive use of genome comparison. It offers a seamless integration of data from ongoing sequencing projects and annotated genomic sequences obtained from GenBank. A variety of export/import filters manages an effective data flow from sequence assembly and manipulation programs (e.g., GAP4) to GenColors and back as well as to standard GenBank file(s). The genome comparison tools include best bidirectional hits, gene conservation, syntenies, and gene core sets. Precomputed UniProt matches allow annotation and analysis in an effective manner. In addition to these analysis options, base-specific quality data (coverage and confidence) can also be handled if available. The GenColors system can be used both for annotation purposes in ongoing genome projects and as an analysis tool for finished genomes. GenColors comes in two types, as dedicated genome browsers and as the Jena Prokaryotic Genome Viewer (JPGV). Dedicated genome browsers contain genomic information on a set of related genomes and offer a large number of options for genome comparison. The system has been efficiently used in the genomic sequencing of Borrelia garinii and is currently applied to various ongoing genome projects on Borrelia, Legionella, Escherichia, and Pseudomonas genomes. One of these dedicated browsers, the Spirochetes Genome Browser (sgb.fli-leibniz.de) with Borrelia, Leptospira, and Treponema genomes, is freely accessible. The others will be released after finalization of the corresponding genome projects. JPGV (jpgv.fli-leibniz.de) offers information on almost all finished bacterial genomes, as compared to the dedicated browsers with reduced genome comparison functionality, however. As of January 2006, this viewer includes 632 genomic elements (e.g., chromosomes and plasmids) of 293 species. The system provides versatile quick and advanced search options for all currently known prokaryotic genomes and generates circular and linear genome plots. Gene information sheets contain basic gene information, database search options, and links to external databases. GenColors is also available on request for local installation.

Comparative genome analysis: selection pressure on the Borrelia vls cassettes is essential for infectivity.

BACKGROUND: At least three species of Borrelia burgdorferi sensu lato (Bbsl) cause tick-borne Lyme disease. Previous work including the genome analysis of B. burgdorferi B31 and B. garinii PBi suggested a highly variable plasmid part. The frequent occurrence of duplicated sequence stretches, the observed plasmid redundancy, as well as the mainly unknown function and variability of plasmid encoded genes rendered the relationships between plasmids within and between species largely unresolvable. RESULTS: To gain further insight into Borreliae genome properties we completed the plasmid sequences of B. garinii PBi, added the genome of a further species, B. afzelii PKo, to our analysis, and compared for both species the genomes of pathogenic and apathogenic strains. The core of all Bbsl genomes consists of the chromosome and two plasmids collinear between all species. We also found additional groups of plasmids, which share large parts of their sequences. This makes it very likely that these plasmids are relatively stable and share common ancestors before the diversification of Borrelia species. The analysis of the differences between B. garinii PBi and B. afzelii PKo genomes of low and high passages revealed that the loss of infectivity is accompanied in both species by a loss of similar genetic material. Whereas B. garinii PBi suffered only from the break-off of a plasmid end, B. afzelii PKo lost more material, probably an entire plasmid. In both cases the vls gene locus encoding for variable surface proteins is affected. CONCLUSION: The complete genome sequences of a B. garinii and a B. afzelii strain facilitate further comparative studies within the genus Borrellia. Our study shows that loss of infectivity can be traced back to only one single event in B. garinii PBi: the loss of the vls cassettes possibly due to error prone gene conversion. Similar albeit extended losses in B. afzelii PKo support the hypothesis that infectivity of Borrelia species depends heavily on the evasion from the host response.