Separation and Mapping of Chromosomes of Parasitic Protozoa

Many protozoan parasites represent an important group of human pathogens. Pulsed Field Gradient Gel Electrophoresis (PFGE) analysis has been an important tool for fundamental genetic studies of parasites like Trypanosoma, Leishmania, Giardia or the human malaria parasite Plasmodium falciparum. We present PFGE conditions allowing a high resolution separation of chromosomes ranging from 500 to 4000 kb within a two day electrophoresis run. In addition, we present conditions for separating large chromosomes (2000-6000 kb) within 36 hr. We demontrate that the application of two dimentional PFGE (2D-PFGE) technique to parasite karyotypes is a very useful method for the analysis of dispersed gene families and comparative studies of the intrachomosomal genome organization

Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Characterization of the hcnABC gene cluster encoding hydrogen cyanide synthase and anaerobic regulation by ANR in the strictly aerobic biocontrol agent Pseudomonas fluorescens CHA0.

The secondary metabolite hydrogen cyanide (HCN) is produced by Pseudomonas fluorescens from glycine, essentially under microaerophilic conditions. The genetic basis of HCN synthesis in P. fluorescens CHA0 was investigated. The contiguous structural genes hcnABC encoding HCN synthase were expressed from the T7 promoter in Escherichia coli, resulting in HCN production in this bacterium. Analysis of the nucleotide sequence of the hcnABC genes showed that each HCN synthase subunit was similar to known enzymes involved in hydrogen transfer, i.e., to formate dehydrogenase (for HcnA) or amino acid oxidases (for HcnB and HcnC). These similarities and the presence of flavin adenine dinucleotide- or NAD(P)-binding motifs in HcnB and HcnC suggest that HCN synthase may act as a dehydrogenase in the reaction leading from glycine to HCN and CO2. The hcnA promoter was mapped by primer extension; the -40 sequence (TTGGC ... ATCAA) resembled the consensus FNR (fumarate and nitrate reductase regulator) binding sequence (TTGAT ... ATCAA). The gene encoding the FNR-like protein ANR (anaerobic regulator) was cloned from P. fluorescens CHA0 and sequenced. ANR of strain CHA0 was most similar to ANR of P. aeruginosa and CydR of Azotobacter vinelandii. An anr mutant of P. fluorescens (CHA21) produced little HCN and was unable to express an hcnA-lacZ translational fusion, whereas in wild-type strain CHA0, microaerophilic conditions strongly favored the expression of the hcnA-lacZ fusion. Mutant CHA21 as well as an hcn deletion mutant were impaired in their capacity to suppress black root rot of tobacco, a disease caused by Thielaviopsis basicola, under gnotobiotic conditions. This effect was most pronounced in water-saturated artificial soil, where the anr mutant had lost about 30% of disease suppression ability, compared with wild-type strain CHA0. These results show that the anaerobic regulator ANR is required for cyanide synthesis in the strictly aerobic strain CHA0 and suggest that ANR-mediated cyanogenesis contributes to the suppression of black root rot.

Re-mapping the molecular features of the human immunodeficiency virus type 1 and human T-cell lymphotropic virus type 1 Brazilian sequences using a bioinformatics unit established in Salvador, Bahia, Brazil, to give support to the viral epidemiology studies

The analysis of genetic data for human immunodeficiency virus type 1 (HIV-1) and human T-cell lymphotropic virus type 1 (HTLV-1) is essential to improve treatment and public health strategies as well as to select strains for vaccine programs. However, the analysis of large quantities of genetic data requires collaborative efforts in bioinformatics, computer biology, molecular biology, evolution, and medical science. The objective of this study was to review and improve the molecular epidemiology of HIV-1 and HTLV-1 viruses isolated in Brazil using bioinformatic tools available in the Laboratório Avançado de Sáude Pública (Lasp) bioinformatics unit. The analysis of HIV-1 isolates confirmed a heterogeneous distribution of the viral genotypes circulating in the country. The Brazilian HIV-1 epidemic is characterized by the presence of multiple subtypes (B, F1, C) and B/F1 recombinant virus while, on the other hand, most of the HTLV-1 sequences were classified as Transcontinental subgroup of the Cosmopolitan subtype. Despite the high variation among HIV-1 subtypes, protein glycosylation and phosphorylation domains were conserved in the pol, gag, and env genes of the Brazilian HIV-1 strains suggesting constraints in the HIV-1 evolution process. As expected, the functional protein sites were highly conservative in the HTLV-1 env gene sequences. Furthermore, the presence of these functional sites in HIV-1 and HTLV-1 strains could help in the development of vaccines that pre-empt the viral escape process.

Multiple adenosine 3',5'-cyclic [corrected] monophosphate response element DNA-binding proteins generated by gene diversification and alternative exon splicing.

The protein sequence deduced from the open reading frame of a human placental cDNA encoding a cAMP-responsive enhancer (CRE)-binding protein (CREB-327) has structural features characteristic of several other transcriptional transactivator proteins including jun, fos, C/EBP, myc, and CRE-BP1. Results of Southwestern analysis of nuclear extracts from several different cell lines show that there are multiple CRE-binding proteins, which vary in size in cell lines derived from different tissues and animal species. To examine the molecular diversity of CREB-327 and related proteins at the nucleic acid level, we used labeled cDNAs from human placenta that encode two different CRE-binding proteins (CREB-327 and CRE-BP1) to probe Northern and Southern blots. Both probes hybridized to multiple fragments on Southern blots of genomic DNA from various species. Alternatively, when a human placental c-jun probe was hybridized to the same blot, a single fragment was detected in most cases, consistent with the intronless nature of the human c-jun gene. The CREB-327 probe hybridized to multiple mRNAs, derived from human placenta, ranging in size from 2-9 kilobases. In contrast, the CRE-BP1 probe identified a single 4-kilobase mRNA. Sequence analyses of several overlapping human genomic cosmid clones containing CREB-327 sequences in conjunction with polymerase chain reaction indicates that the CREB-327/341 cDNAs are composed of at least eight or nine exons, and analyses of human placental cDNAs provide direct evidence for at least one alternatively spliced exon. Analyses of mouse/hamster-human hybridoma DNAs by Southern blotting and polymerase chain reaction localizes the CREB-327/341 gene to human chromosome 2. The results indicate that there is a dichotomy of CREB-like proteins, those that are related by overall structure and DNA-binding specificity as well as those that are related by close similarities of primary sequences.

Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

The effect of translocation-induced nuclear reorganization on gene expression.

Translocations are known to affect the expression of genes at the breakpoints and, in the case of unbalanced translocations, alter the gene copy number. However, a comprehensive understanding of the functional impact of this class of variation is lacking. Here, we have studied the effect of balanced chromosomal rearrangements on gene expression by comparing the transcriptomes of cell lines from controls and individuals with the t(11;22)(q23;q11) translocation. The number of differentially expressed transcripts between translocation-carrying and control cohorts is significantly higher than that observed between control samples alone, suggesting that balanced rearrangements have a greater effect on gene expression than normal variation. Many of the affected genes are located along the length of the derived chromosome 11. We show that this chromosome is concomitantly altered in its spatial organization, occupying a more central position in the nucleus than its nonrearranged counterpart. Derivative 22-mapping chromosome 22 genes, on the other hand, remain in their usual environment. Our results are consistent with recent studies that experimentally altered nuclear organization, and indicated that nuclear position plays a functional role in regulating the expression of some genes in mammalian cells. Our study suggests that chromosomal translocations can result in hitherto unforeseen, large-scale changes in gene expression that are the consequence of alterations in normal chromosome territory positioning. This has consequences for the patterns of gene expression change seen during tumorigenesis-associated genome instability and during the karyotype changes that lead to speciation.

The S. pombe cdc15 gene is a key element in the reorganization of F-actin at mitosis.

The S. pombe cdc15 gene is essential for cell division. cdc15ts mutants do not form a septum, but growth and nuclear division continue, leading to formation of multinucleate cells. The earliest step in septum formation and cytokinesis, rearrangement of actin to the center of the cell, is associated with appearance of hypophosphorylated cdc15p and formation of a cdc15p ring, which colocalizes with actin. Loss of cdc15p function impairs formation of the actin ring. The abundance of cdc15 mRNA varies through the cell division cycle, peaking in early mitosis before septation. Expression of cdc15 in G2-arrested cells induces actin rearrangement to the center of the cell. These data implicate cdc15p as a key element in mediating the cytoskeletal rearrangements required for cytokinesis.

Occurrence of blaOXA-23 gene in imipenem-susceptible Acinetobacter baumannii

The aim of the current study was to describe the occurrence of the blaOXA-23 gene and the ISAba1 element in imipenem-susceptible Acinetobacter baumannii strains. By performing the polymerase chain reaction mapping using combinations of ISAba1 forward primers and the blaOXA-23-like gene reverse primers, we demonstrated that the ISAba1 element did not occur upstream of the blaOXA-23 gene in five of 31 isolates, which explained the lack of resistance to imipenem despite the presence of the blaOXA-23 gene. All of the blaOXA-23-positive isolates were susceptible to imipenem and meropenem with minimal inhibitory concentration < 4 µg/mL. Pulsed-field gel electrophoresis analysis revealed four genotypes among the five blaOXA-23-positive isolates. The current report of the blaOXA-23 gene in imipenem-susceptible isolates provided evidence that this gene may be silently spread in a hospital environment and highlighted the threat of undetected reservoirs of carbapenemase genes.

IL2RA/CD25 gene polymorphisms: uneven association with multiple sclerosis (MS) and type 1 diabetes (T1D).

BACKGROUND: IL-2 receptor (IL2R) alpha is the specific component of the high affinity IL2R system involved in the immune response and in the control of autoimmunity. METHODS AND RESULTS: Here we perform a replication and fine mapping of the IL2RA gene region analyzing 3 SNPs previously associated with multiple sclerosis (MS) and 5 SNPs associated with type 1 diabetes (T1D) in a collection of 798 MS patients and 927 matched Caucasian controls from the south of Spain. We observed association with MS in 6 of 8 SNPs. The rs1570538, at the 3'- UTR extreme of the gene, previously reported to have a weak association with MS, is replicated here (P = 0.032). The most associated T1D SNP (rs41295061) was not associated with MS in the present study. However, the rs35285258, belonging to another independent group of SNPs associated with T1D, showed the maximal association in this study but different risk allele. We replicated the association of only one (rs2104286) of the two IL2RA SNPs identified in the recently performed genome-wide association study of MS. CONCLUSIONS: These findings confirm and extend the association of this gene with MS and reveal a genetic heterogeneity of the associated polymorphisms and risk alleles between MS and T1D suggesting different immunopathological roles of IL2RA in these two diseases.

Describing ancient horizontal gene transfers at the nucleotide and gene levels by comparative pathogenicity island genometrics.

MOTIVATION: Lateral gene transfer is a major mechanism contributing to bacterial genome dynamics and pathovar emergence via pathogenicity island (PAI) spreading. However, since few of these genomic exchanges are experimentally reproducible, it is difficult to establish evolutionary scenarios for the successive PAI transmissions between bacterial genera. Methods initially developed at the gene and/or nucleotide level for genomics, i.e. comparisons of concatenated sequences, ortholog frequency, gene order or dinucleotide usage, were combined and applied here to homologous PAIs: we call this approach comparative PAI genometrics. RESULTS: YAPI, a Yersinia PAI, and related islands were compared with measure evolutionary relationships between related modules. Through use of our genometric approach designed for tracking codon usage adaptation and gene phylogeny, an ancient inter-genus PAI transfer was oriented for the first time by characterizing the genomic environment in which the ancestral island emerged and its subsequent transfers to other bacterial genera.

New gene functions in megakaryopoiesis and platelet formation.

Platelets are the second most abundant cell type in blood and are essential for maintaining haemostasis. Their count and volume are tightly controlled within narrow physiological ranges, but there is only limited understanding of the molecular processes controlling both traits. Here we carried out a high-powered meta-analysis of genome-wide association studies (GWAS) in up to 66,867 individuals of European ancestry, followed by extensive biological and functional assessment. We identified 68 genomic loci reliably associated with platelet count and volume mapping to established and putative novel regulators of megakaryopoiesis and platelet formation. These genes show megakaryocyte-specific gene expression patterns and extensive network connectivity. Using gene silencing in Danio rerio and Drosophila melanogaster, we identified 11 of the genes as novel regulators of blood cell formation. Taken together, our findings advance understanding of novel gene functions controlling fate-determining events during megakaryopoiesis and platelet formation, providing a new example of successful translation of GWAS to function.

Cross-species and cross-platform gene expression studies with the Bioconductor-compliant R package 'annotationTools'.

Background: The variety of DNA microarray formats and datasets presently available offers an unprecedented opportunity to perform insightful comparisons of heterogeneous data. Cross-species studies, in particular, have the power of identifying conserved, functionally important molecular processes. Validation of discoveries can now often be performed in readily available public data which frequently requires cross-platform studies.Cross-platform and cross-species analyses require matching probes on different microarray formats. This can be achieved using the information in microarray annotations and additional molecular biology databases, such as orthology databases. Although annotations and other biological information are stored using modern database models ( e. g. relational), they are very often distributed and shared as tables in text files, i.e. flat file databases. This common flat database format thus provides a simple and robust solution to flexibly integrate various sources of information and a basis for the combined analysis of heterogeneous gene expression profiles.Results: We provide annotationTools, a Bioconductor-compliant R package to annotate microarray experiments and integrate heterogeneous gene expression profiles using annotation and other molecular biology information available as flat file databases. First, annotationTools contains a specialized set of functions for mining this widely used database format in a systematic manner. It thus offers a straightforward solution for annotating microarray experiments. Second, building on these basic functions and relying on the combination of information from several databases, it provides tools to easily perform cross-species analyses of gene expression data.Here, we present two example applications of annotationTools that are of direct relevance for the analysis of heterogeneous gene expression profiles, namely a cross-platform mapping of probes and a cross-species mapping of orthologous probes using different orthology databases. We also show how to perform an explorative comparison of disease-related transcriptional changes in human patients and in a genetic mouse model.Conclusion: The R package annotationTools provides a simple solution to handle microarray annotation and orthology tables, as well as other flat molecular biology databases. Thereby, it allows easy integration and analysis of heterogeneous microarray experiments across different technological platforms or species.

The gene for the ligand binding chain of the human interferon gamma receptor.

In order to characterize the gene encoding the ligand binding (1(st); alpha) chain of the human IFN-gamma receptor, two overlapping cosmid clones were analyzed. The gene spans over 25 kilobases (kb) of the genomic DNA and has seven exons. The extracellular domain is encoded by exons 1 to 5 and by part of exon 6. The transmembrane region is also encoded by exon 6. Exon 7 encodes the intracellular domain and the 3' untranslated portion. The gene was located on chromosome 6q23.1, as determined by in situ hybridization. The 4 kb region upstream (5') of the gene was sequenced and analyzed for promoter activity. No consensus-matching TATA or CAAT boxes in the 5' region were found. Potential binding sites for Sp1, AP-1, AP-2, and CREB nuclear factors were identified. Compatible with the presence of the Sp1/AP-2 sites and the lack of TATA box, S1-nuclease mapping experiments showed multiple transcription initiation sites. Promoter activity of the 5' flanking region was analyzed with two different reporter genes: the Escherichia coli chloramphenicol acetyltransferase and human growth hormone. The smallest 5' region of the gene that still had full promoter activity was 692 base pairs in length. In addition, we found sequences belonging to the oldest family of Alu repeats, 2 - 3 kb upstream of the gene, which could be useful for genetic studies.

The TPTE gene family: cellular expression, subcellular localization and alternative splicing.

The human TPTE (Transmembrane Phosphatase with TEnsin homology) gene family encodes a PTEN-related tyrosine phosphatase with four potential transmembrane domains. Chromosomal mapping revealed multiple copies of the TPTE gene on chromosomes 13, 15, 21, 22 and Y. Human chromosomes 13 and 21 copies encode two functional proteins, TPIP (TPTE and PTEN homologous Inositol lipid Phosphatase) and TPTE, respectively, whereas only one copy of the gene exists in the mouse genome. In the present study, we show that TPTE and TPIP proteins are expressed in secondary spermatocytes and/or prespermatids. In addition, we report the existence of several novel alternatively spliced isoforms of these two proteins with variable number of transmembrane domains. The latter has no influence on the subcellular localization of these different peptides as shown by co-immunofluorescence experiments. Finally, we identify another expressed TPTE copy, mapping to human chromosome 22, whose transcription appears to be under the control of the LTR of human endogenous retrovirus RTVL-H3.