Deciphering the complex nature of bolting time regulation in Beta vulgaris

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 F₃ 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 ² > 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 ₁ 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.

Crop wild relative populations of Beta vulgaris allow direct mapping of agronomically important genes

Rapid identification of agronomically important genes is of pivotal interest for crop breeding. One source of such genes are crop wild relative (CWR) populations. Here we used a CWR population of <200 wild beets (B. vulgaris ssp. maritima), sampled in their natural habitat, to identify the sugar beet (Beta vulgaris ssp. vulgaris) resistance gene Rz2 with a modified version of mapping-by-sequencing (MBS). For that, we generated a draft genome sequence of the wild beet. Our results show the importance of preserving CWR in situ and demonstrate the great potential of CWR for rapid discovery of causal genes relevant for crop improvement. The candidate gene for Rz2 was identified by MBS and subsequently corroborated via RNA interference (RNAi). Rz2 encodes a CC-NB-LRR protein. Access to the DNA sequence of Rz2 opens the path to improvement of resistance towards rhizomania not only by marker-assisted breeding but also by genome editing.

Variation among wild relatives of crop plants can be used to identify genes underlying traits of agronomic importance. Here, the authors show that a modified mapping-by-sequencing approach can rapidly identify the genetic basis for viral resistance in sugar beet using wild beet populations in their natural habitat.

A Detailed Analysis of the BR1 Locus Suggests a New Mechanism for Bolting after Winter in Sugar Beet (Beta vulgaris L.)

Sugar beet (Beta vulgaris ssp. vulgaris) is a biennial, sucrose-storing plant, which is mainly cultivated as a spring crop and harvested in the vegetative stage before winter. For increasing beet yield, over-winter cultivation would be advantageous. However, bolting is induced after winter and drastically reduces yield. Thus, post-winter bolting control is essential for winter beet cultivation. To identify genetic factors controlling bolting after winter, a F₂ population was previously developed by crossing the sugar beet accessions BETA 1773 with reduced bolting tendency and 93161P with complete bolting after winter. For a mapping-by-sequencing analysis, pools of 26 bolting-resistant and 297 bolting F₂ plants were used. Thereby, a single continuous homozygous region of 103 kb was co-localized to the previously published BR₁ QTL for post-winter bolting resistance (Pfeiffer et al., 2014). The BR₁ locus was narrowed down to 11 candidate genes from which a homolog of the Arabidopsis CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 73-I (CPSF73-I) was identified as the most promising candidate. A 2 bp deletion within the BETA 1773 allele of BvCPSF73-Ia results in a truncated protein. However, the null allele of BvCPSF73-Ia might partially be compensated by a second BvCPSF73-Ib gene. This gene is located 954 bp upstream of BvCPSF73-Ia and could be responsible for the incomplete penetrance of the post-winter bolting resistance allele of BETA 1773. This result is an important milestone for breeding winter beets with complete bolting resistance after winter.

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.

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

Widespread disruption of host transcription termination in HSV-1 infection

Herpes simplex virus 1 (HSV-1) is an important human pathogen and a paradigm for virus-induced host shut-off. Here we show that global changes in transcription and RNA processing and their impact on translation can be analysed in a single experimental setting by applying 4sU-tagging of newly transcribed RNA and ribosome profiling to lytic HSV-1 infection. Unexpectedly, we find that HSV-1 triggers the disruption of transcription termination of cellular, but not viral, genes. This results in extensive transcription for tens of thousands of nucleotides beyond poly(A) sites and into downstream genes, leading to novel intergenic splicing between exons of neighbouring cellular genes. As a consequence, hundreds of cellular genes seem to be transcriptionally induced but are not translated. In contrast to previous reports, we show that HSV-1 does not inhibit co-transcriptional splicing. Our approach thus substantially advances our understanding of HSV-1 biology and establishes HSV-1 as a model system for studying transcription termination.

Herpes simplex virus 1 (HSV-1) efficiently shuts down host gene expression in infected cells. Here Rutkowski et al. analyse the genome-wide changes in transcription and translation in infected cells, and show that HSV-1 triggers an extensive disruption of transcription termination of cellular genes.

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.

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.

Comparison of the bacterial composition and structure in symptomatic and asymptomatic endodontic infections associated with root-filled teeth using pyrosequencing

Residual microorganisms and/or re-infections are a major cause for root canal therapy failure. Understanding of the bacterial content could improve treatment protocols. Fifty samples from 25 symptomatic and 25 asymptomatic previously root-filled teeth were collected from Sudanese patients with periradicular lesions. Amplified 16S rRNA gene (V1-V2) variable regions were subjected to pyrosequencing (FLX 454) to determine the bacterial profile. Obtained quality-controlled sequences from forty samples were classified into 741 operational taxonomic units (OTUs) at 3% dissimilarity, 525 at 5% dissimilarity and 297 at 10% dissimilarity, approximately corresponding to species-, genus- and class levels. The most abundant phyla were: Firmicutes (29.9%), Proteobacteria (26.1%), Actinobacteria (22.72%), Bacteroidetes (13.31%) and Fusobacteria (4.55%). Symptomatic patients had more Firmicutes and Fusobacteria than asymptomatic patients, while asymptomatic patients showed more Proteobacteria and Actinobacteria. Interaction of disease status and age was observed by two-way ANOSIM. Canonical correspondence analysis for age, tooth restoration and disease status showed a correlation of disease status with the composition and prevalence of different members of the microbial community. The pyrosequencing analysis revealed a distinctly higher diversity of the microbiota compared to earlier reports. The comparison of symptomatic and asymptomatic patients showed a clear association of the composition of the bacterial community with the presence and absence of symptoms in conjunction with the patients’ age.

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.

The Transcriptome Analysis and Comparison Explorer--T-ACE: a platform-independent, graphical tool to process large RNAseq datasets of non-model organisms

MOTIVATION

Next generation sequencing (NGS) technologies allow a rapid and cost-effective compilation of large RNA sequence datasets in model and non-model organisms. However, the storage and analysis of transcriptome information from different NGS platforms is still a significant bottleneck, leading to a delay in data dissemination and subsequent biological understanding. Especially database interfaces with transcriptome analysis modules going beyond mere read counts are missing. Here, we present the Transcriptome Analysis and Comparison Explorer (T-ACE), a tool designed for the organization and analysis of large sequence datasets, and especially suited for transcriptome projects of non-model organisms with little or no a priori sequence information. T-ACE offers a TCL-based interface, which accesses a PostgreSQL database via a php-script. Within T-ACE, information belonging to single sequences or contigs, such as annotation or read coverage, is linked to the respective sequence and immediately accessible. Sequences and assigned information can be searched via keyword- or BLAST-search. Additionally, T-ACE provides within and between transcriptome analysis modules on the level of expression, GO terms, KEGG pathways and protein domains. Results are visualized and can be easily exported for external analysis. We developed T-ACE for laboratory environments, which have only a limited amount of bioinformatics support, and for collaborative projects in which different partners work on the same dataset from different locations or platforms (Windows/Linux/MacOS). For laboratories with some experience in bioinformatics and programming, the low complexity of the database structure and open-source code provides a framework that can be customized according to the different needs of the user and transcriptome project.

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.

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.

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.

DNA sequence and analysis of human chromosome 8

The International Human Genome Sequencing Consortium (IHGSC) recently completed a sequence of the human genome. As part of this project, we have focused on chromosome 8. Although some chromosomes exhibit extreme characteristics in terms of length, gene content, repeat content and fraction segmentally duplicated, chromosome 8 is distinctly typical in character, being very close to the genome median in each of these aspects. This work describes a finished sequence and gene catalogue for the chromosome, which represents just over 5% of the euchromatic human genome. A unique feature of the chromosome is a vast region of approximately 15 megabases on distal 8p that appears to have a strikingly high mutation rate, which has accelerated in the hominids relative to other sequenced mammals. This fast-evolving region contains a number of genes related to innate immunity and the nervous system, including loci that appear to be under positive selection--these include the major defensin (DEF) gene cluster and MCPH1, a gene that may have contributed to the evolution of expanded brain size in the great apes. The data from chromosome 8 should allow a better understanding of both normal and disease biology and genome evolution.

Characterization of 8p21.3 chromosomal deletions in B-cell lymphoma: TRAIL-R1 and TRAIL-R2 as candidate dosage-dependent tumor suppressor genes

Deletions of chromosome 8p are a recurrent event in B-cell non-Hodgkin lymphoma (B-NHL), suggesting the presence of a tumor suppressor gene. We have characterized these deletions using comparative genomic hybridization to microarrays, fluorescence in situ hybridization (FISH) mapping, DNA sequencing, and functional studies. A minimal deleted region (MDR) of 600 kb was defined in chromosome 8p21.3, with one mantle cell lymphoma cell line (Z138) exhibiting monoallelic deletion of 650 kb. The MDR extended from bacterial artificial chromosome (BAC) clones RP11-382J24 and RP11-109B10 and included the tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) receptor gene loci. Sequence analysis of the individual expressed genes within the MDR and DNA sequencing of the entire MDR in Z138 did not reveal any mutation. Gene expression analysis and quantitative reverse transcriptase-polymerase chain reaction (QRT-PCR) showed down-regulation of TRAIL-R1 and TRAIL-R2 receptor genes as a consistent event in B-NHL with 8p21.3 loss. Epigenetic inactivation was excluded via promoter methylation analysis. In vitro studies showed that TRAIL-induced apoptosis was dependent on TRAIL-R1 and/or -R2 dosage in most tumors. Resistance to apoptosis of cell lines with 8p21.3 deletion was reversed by restoration of TRAIL-R1 or TRAIL-R2 expression by gene transfection. Our data suggest that TRAIL-R1 and TRAIL-R2 act as dosage-dependent tumor suppressor genes whose monoallelic deletion can impair TRAIL-induced apoptosis in B-cell lymphoma.

Polymorphic segmental duplications at 8p23.1 challenge the determination of individual defensin gene repertoires and the assembly of a contiguous human reference sequence

Background

Defensins are important components of innate immunity to combat bacterial and viral infections, and can even elicit antitumor responses. Clusters of defensin (DEF) genes are located in a 2 Mb range of the human chromosome 8p23.1. This DEF locus, however, represents one of the regions in the euchromatic part of the final human genome sequence which contains segmental duplications, and recalcitrant gaps indicating high structural dynamics.

Results

We find that inter- and intraindividual genetic variations within this locus prevent a correct automatic assembly of the human reference genome (NCBI Build 34) which currently even contains misassemblies. Manual clone-by-clone alignment and gene annotation as well as repeat and SNP/haplotype analyses result in an alternative alignment significantly improving the DEF locus representation. Our assembly better reflects the experimentally verified variability of DEF gene and DEF cluster copy numbers. It contains an additional DEF cluster which we propose to reside between two already known clusters. Furthermore, manual annotation revealed a novel DEF gene and several pseudogenes expanding the hitherto known DEF repertoire. Analyses of BAC and working draft sequences of the chimpanzee indicates that its DEF region is also complex as in humans and DEF genes and a cluster are multiplied. Comparative analysis of human and chimpanzee DEF genes identified differences affecting the protein structure. Whether this might contribute to differences in disease susceptibility between man and ape remains to be solved. For the determination of individual DEF gene repertoires we provide a molecular approach based on DEF haplotypes.

Conclusions

Complexity and variability seem to be essential genomic features of the human DEF locus at 8p23.1 and provides an ongoing challenge for the best possible representation in the human reference sequence. Dissection of paralogous sequence variations, duplicon SNPs ans multisite variations as well as haplotypes by sequencing based methods is the way for future studies of interindividual DEF locus variability and its disease association.

Comparative analysis of the Borrelia garinii genome

Three members of the genus Borrelia (B.burgdorferi, B.garinii, B.afzelii) cause tick-borne borreliosis. Depending on the Borrelia species involved, the borreliosis differs in its clinical symptoms. Comparative genomics opens up a way to elucidate the underlying differences in Borrelia species. We analysed a low redundancy whole-genome shotgun (WGS) assembly of a B.garinii strain isolated from a patient with neuroborreliosis in comparison to the B.burgdorferi genome. This analysis reveals that most of the chromosome is conserved (92.7% identity on DNA as well as on amino acid level) in the two species, and no chromosomal rearrangement or larger insertions/deletions could be observed. Furthermore, two collinear plasmids (lp54 and cp26) seem to belong to the basic genome inventory of Borrelia species. These three collinear parts of the Borrelia genome encode 861 genes, which are orthologous in the two species examined. The majority of the genetic information of the other plasmids of B.burgdorferii is also present in B.garinii although orthology is not easy to define due to a high redundancy of the plasmid fraction. Yet, we did not find counterparts of the B.burgdorferi plasmids lp36 and lp38 or their respective gene repertoire in the B.garinii genome. Thus, phenotypic differences between the two species could be attributable to the presence or absence of these two plasmids as well as to the potentially positively selected genes.

A Leishmania donovani gene that confers accelerated recovery from stationary phase growth arrest

We have isolated a gene, LdGF1, from the protozoan parasite Leishmania donovani. Overexpression of this gene confers a strong selective advantage in liquid culture after stationary phase growth arrest. We could show that recombinant L. donovani or Leishmania major, when overexpressing LdGF1, recover faster from a stationary phase growth arrest than control parasite strains. While no advantage of LdGF1 overexpression could be observed in log phase cultures or after a hydroxyurea-induced S-phase growth arrest, recovery from a cell cycle arrest due to serum deprivation was faster in LdGF1-overexpressing strains. This was found to be due to an accelerated release from a G(1) cell cycle arrest. By contrast, in a BALB/c mouse infection system, overexpression of LdGF1 in L. major resulted in reduced virulence. We conclude that increased levels of LdGF1 are beneficiary during recovery from G(1) cell cycle arrest, but pose a disadvantage inside a mammalian host. These results are discussed in the context of the observed loss of virulence during in vitro passage of Leishmania parasites.

Use of genetic complementation to identify gene(s) which specify species-specific organ tropism of Leishmania

We have employed a genetic complementation screening to identify genetic markers of heat stress tolerance and visceralisation of Leishmania infection. Leishmania major, which has a low thermotolerance and which causes cutaneous lesions, was transfected with a cosmid library of L. donovani DNA. The recombinant parasites were then screened either for thermotolerance or selected by repeated passage in BALB/c mice. Cosmids which conferred selective advantage were isolated. Several strategies were tested to identify the gene(s) within the cosmids responsible for the observed selective advantages. Of the approaches tested, the complete sequence analysis of the cosmids and subsequent screening of defined candidate ORFs proved to be the method of choice. Other approaches, such as creation of sub-libraries or transposon insertion strategies proved to be unsuccessful.

A physical map of the human genome

The human genome is by far the largest genome to be sequenced, and its size and complexity present many challenges for sequence assembly. The International Human Genome Sequencing Consortium constructed a map of the whole genome to enable the selection of clones for sequencing and for the accurate assembly of the genome sequence. Here we report the construction of the whole-genome bacterial artificial chromosome (BAC) map and its integration with previous landmark maps and information from mapping efforts focused on specific chromosomal regions. We also describe the integration of sequence data with the map.