Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Type I signal peptidase from Leishmania is a target of the immune response in human cutaneous and visceral leishmaniasis.

The gene encoding type I signal peptidase (Lmjsp) has been cloned from Leishmania major. Lmjsp encodes a protein of 180 amino residues with a predicted molecular mass of 20.5 kDa. Comparison of the protein sequence with those of known type I signal peptidases indicates homology in five conserved domains A-E which are known to be important, or essential, for catalytic activity. Southern blot hybridisation analysis indicates that there is a single copy of the Lmjsp gene. A recombinant SPase protein and a synthetic peptide of the L. major signal peptidase were used to examine the presence of specific antibodies in sera from either recovered or active individuals of both cutaneous and visceral leishmaniasis. This evaluation demonstrated that sera from cutaneous and visceral forms of leishmaniasis are highly reactive to both the recombinant and synthetic signal peptidase antigens. Therefore, the Leishmania signal peptidase, albeit localised intracellularly, is a significant target of the Leishmania specific immune response and highlights its potential use for serodiagnosis of cutaneous and visceral leishmaniasis.

Sense and antisense transcripts in the histone H1 (HIS-1) locus of Leishmania major.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by two genes, HIS-1.1 and HIS-1.2. These genes are separated by approximately 20 kb of sequence and are located on the same DNA strand of chromosome 27. When Northern blots of parasite RNA were probed with HIS-1 strand-specific riboprobes, we detected sense and antisense transcripts that were polyadenylated and developmentally regulated. When the HIS-1.2 coding region was replaced with the coding region of the neomycin phosphotransferase gene, antisense transcription of this gene was unaffected, indicating that the regulatory elements controlling antisense transcription were located outside of the HIS-1.2 gene, and that transcription in Leishmania can occur from both DNA strands even in the presence of transcription of a selectable marker in the complementary strand. A search for other antisense transcripts within the HIS-1 locus identified an additional transcript (SC-1) within the intervening HIS-1 sequence, downstream of adenine and thymine-rich sequences. These results show that gene expression in Leishmania is not only regulated polycistronically from the sense strand of genomic DNA, but that the complementary strand of DNA also contains sequences that could drive expression of open reading frames from the antisense strand of DNA. These findings suggest that the parasite has evolved in such a way as to maximise the transcription of its genome, a mechanism that might be important for it to maintain virulence.

Superantigen-reactive CD4+ T cells are required to stimulate B cells after infection with mouse mammary tumor virus.

Superantigens are defined by their ability to stimulate a large fraction of T cells via interaction with the T cell receptor (TCR) V beta domain. Endogenous superantigens, classically termed minor lymphocyte-stimulating (Mls) antigens, were recently identified as products of open reading frames (ORF) in integrated proviral copies of mouse mammary tumor virus (MMTV). We have described an infectious MMTV homologue of the classical endogenous superantigen Mls-1a (Mtv-7). The ORF molecules of both the endogenous Mtv-7 and the infectious MMTV(SW) interact with T cells expressing the TCR V beta 6, 7, 8.1, and 9 domains. Furthermore, the COOH termini of their ORF molecules, thought to confer TCR specificity, are very similar. Since successful transport of MMTV from the site of infection in the gut to the mammary gland depends on a functional immune system, we were interested in determining the early events after and requirements for MMTV infection. We show that MMTV(SW) infection induces a massive response of V beta 6+ CDC4+ T cells, which interact with the viral ORF. Concomitantly, we observed a B cell response and differentiation that depends on both the presence and stimulation of the superantigen-reactive T cells. Furthermore, we show that B cells are the main target of the initial MMTV infection as judged by the presence of the reverse-transcribed viral genome and ORF transcripts. Thus, we suggest that MMTV infection of B cells leads to ORF-mediated B-T cell interaction, which maintains and possibly amplifies viral infection.

Clonal deletion of V beta 14-bearing T cells in mice transgenic for mammary tumour virus.

Autoreactive T lymphocytes are clonally deleted during maturation in the thymus. Deletion of T cells expressing particular receptor V beta elements is controlled by poorly defined autosomal dominant genes. A gene has now been identified by expression of transgenes in mice which causes deletion of V beta 14+ T cells. The gene lies in the open reading frame of the long terminal repeat of the mouse mammary tumour virus.

Autoinduction of 2,4-diacetylphloroglucinol biosynthesis in the biocontrol agent Pseudomonas fluorescens CHA0 and repression by the bacterial metabolites salicylate and pyoluteorin.

The antimicrobial metabolite 2,4-diacetylphloroglucinol (2,4-DAPG) contributes to the capacity of Pseudomonas fluorescens strain CHA0 to control plant diseases caused by soilborne pathogens. A 2, 4-DAPG-negative Tn5 insertion mutant of strain CHA0 was isolated, and the nucleotide sequence of the 4-kb genomic DNA region adjacent to the Tn5 insertion site was determined. Four open reading frames were identified, two of which were homologous to phlA, the first gene of the 2,4-DAPG biosynthetic operon, and to the phlF gene encoding a pathway-specific transcriptional repressor. The Tn5 insertion was located in an open reading frame, tentatively named phlH, which is not related to known phl genes. In wild-type CHA0, 2, 4-DAPG production paralleled expression of a phlA'-'lacZ translational fusion, reaching a maximum in the late exponential growth phase. Thereafter, the compound appeared to be degraded to monoacetylphloroglucinol by the bacterium. 2,4-DAPG was identified as the active compound in extracts from culture supernatants of strain CHA0 specifically inducing phlA'-'lacZ expression about sixfold during exponential growth. Induction by exogenous 2,4-DAPG was most conspicuous in a phlA mutant, which was unable to produce 2, 4-DAPG. In a phlF mutant, 2,4-DAPG production was enhanced severalfold and phlA'-'lacZ was expressed at a level corresponding to that in the wild type with 2,4-DAPG added. The phlF mutant was insensitive to 2,4-DAPG addition. A transcriptional phlA-lacZ fusion was used to demonstrate that the repressor PhlF acts at the level of transcription. Expression of phlA'-'lacZ and 2,4-DAPG synthesis in strain CHA0 was strongly repressed by the bacterial extracellular metabolites salicylate and pyoluteorin as well as by fusaric acid, a toxin produced by the pythopathogenic fungus Fusarium. In the phlF mutant, these compounds did not affect phlA'-'lacZ expression and 2, 4-DAPG production. PhlF-mediated induction by 2,4-DAPG and repression by salicylate of phlA'-'lacZ expression was confirmed by using Escherichia coli as a heterologous host. In conclusion, our results show that autoinduction of 2,4-DAPG biosynthesis can be countered by certain bacterial (and fungal) metabolites. This mechanism, which depends on phlF function, may help P. fluorescens to produce homeostatically balanced amounts of extracellular metabolites.

Cloning and characterization of Helicobacter pylori succinyl CoA:acetoacetate CoA-transferase, a novel prokaryotic member of the CoA-transferase family.

Sequencing of a fragment of Helicobacter pylori genome led to the identification of two open reading frames showing striking homology with Coenzyme A (CoA) transferases, enzymes catalyzing the reversible transfer of CoA from one carboxylic acid to another. The genes were present in all H. pylori strains tested by polymerase chain reaction or slot blotting but not in Campylobacter jejuni. Genes for the putative A and B subunits of H. pylori CoA-transferase were introduced into the bacterial expression vector pKK223-3 and expressed in Escherichia coli JM105 cells. Amino acid sequence comparisons, combined with measurements of enzyme activities using different CoA donors and acceptors, identified the H. pylori CoA-transferase as a succinyl CoA:acetoacetate CoA-transferase. This activity was consistently observed in different H. pylori strains. Antibodies raised against either recombinant A or B subunits recognized two distinct subunits of Mr approximately 26,000 and 24, 000 that are both necessary for H. pylori CoA-transferase function. The lack of alpha-ketoglutarate dehydrogenase and of succinyl CoA synthetase activities indicates that the generation of succinyl CoA is not mediated by the tricarboxylic acid cycle in H. pylori. We postulate the existence of an alternative pathway where the CoA-transferase is essential for energy metabolism.

Loss-of-function mutations of an inhibitory upstream ORF in the human hairless transcript cause Marie Unna hereditary hypotrichosis.

Marie Unna hereditary hypotrichosis (MUHH) is an autosomal dominant form of genetic hair loss. In a large Chinese family carrying MUHH, we identified a pathogenic initiation codon mutation in U2HR, an inhibitory upstream ORF in the 5' UTR of the gene encoding the human hairless homolog (HR). U2HR is predicted to encode a 34-amino acid peptide that is highly conserved among mammals. In 18 more families from different ancestral groups, we identified a range of defects in U2HR, including loss of initiation, delayed termination codon and nonsense and missense mutations. Functional analysis showed that these classes of mutations all resulted in increased translation of the main HR physiological ORF. Our results establish the link between MUHH and U2HR, show that fine-tuning of HR protein levels is important in control of hair growth, and identify a potential mechanism for preventing hair loss or promoting hair removal.

Prominent use of distal 5' transcription start sites and discovery of a large number of additional exons in ENCODE regions.

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

The small RNA PhrS stimulates synthesis of the Pseudomonas aeruginosa quinolone signal.

Quorum sensing, a cell-to-cell communication system based on small signal molecules, is employed by the human pathogen Pseudomonas aeruginosa to regulate virulence and biofilm development. Moreover, regulation by small trans-encoded RNAs has become a focal issue in studies of virulence gene expression of bacterial pathogens. In this study, we have identified the small RNA PhrS as an activator of PqsR synthesis, one of the key quorum-sensing regulators in P. aeruginosa. Genetic studies revealed a novel mode of regulation by a sRNA, whereby PhrS uses a base-pairing mechanism to activate a short upstream open reading frame to which the pqsR gene is translationally coupled. Expression of phrS requires the oxygen-responsive regulator ANR. Thus, PhrS is the first bacterial sRNA that provides a regulatory link between oxygen availability and quorum sensing, which may impact on oxygen-limited growth in P. aeruginosa biofilms.

New infectious mammary tumor virus superantigen with V beta-specificity identical to staphylococcal enterotoxin B (SEB).

Only few infectious mouse mammary tumor viruses (MMTV) have been characterized which induce a potent superantigen response in vivo. Here we describe the characterization of an MMTV which was isolated from milk of the highly mammary tumor-prone SHN mouse strain. Exposure of newborn mice to milk-borne MMTV (SHN) results in a very slow deletion of V beta 7, 8.1, 8.2 and 8.3 expressing peripheral T cells. Subcutaneous injection of adult mice with this virus induces a rapid and strong stimulation of all four affected V beta-subsets in vivo. Besides the strong T cell effect we observed an early proliferation and activation of the local B cell pool leading to the initial secretion of IgM followed by preferential secretion of IgG2a by day 6. Sequence comparison of the polymorphic C terminus with known open reading frames revealed high homology to the endogenous provirus Mtv-RCS. This is the first report of a virus having a complete overlap in V beta-specificity with a bacterial superantigen stimulating as many as 35% of the whole CD4+ T cell repertoire including V beta 8.2.

Identification of Genes Affecting Vacuole Membrane Fragmentation in Saccharomyces cerevisiae.

The equilibrium of membrane fusion and fission influences the volume and copy number of organelles. Fusion of yeast vacuoles has been well characterized but their fission and the mechanisms determining vacuole size and abundance remain poorly understood. We therefore attempted to systematically characterize factors necessary for vacuole fission. Here, we present results of an in vivo screening for deficiencies in vacuolar fragmentation activity of an ordered collection deletion mutants, representing 4881 non-essential genes of the yeast Saccharomyces cerevisiae. The screen identified 133 mutants with strong defects in vacuole fragmentation. These comprise numerous known fragmentation factors, such as the Fab1p complex, Tor1p, Sit4p and the V-ATPase, thus validating the approach. The screen identified many novel factors promoting vacuole fragmentation. Among those are 22 open reading frames of unknown function and three conspicuous clusters of proteins with known function. The clusters concern the ESCRT machinery, adaptins, and lipases, which influence the production of diacylglycerol and phosphatidic acid. A common feature of these factors of known function is their capacity to change membrane curvature, suggesting that they might promote vacuole fragmentation via this property.

Tn5-directed cloning of pqq genes from Pseudomonas fluorescens CHA0: mutational inactivation of the genes results in overproduction of the antibiotic pyoluteorin.

Pseudomonas fluorescens CHA0 produces several secondary metabolites, e.g., the antibiotics pyoluteorin (Plt) and 2,4-diacetylphloroglucinol (Phl), which are important for the suppression of root diseases caused by soil-borne fungal pathogens. A Tn5 insertion mutant of strain CHA0, CHA625, does not produce Phl, shows enhanced Plt production on malt agar, and has lost part of the ability to suppress black root rot in tobacco plants and take-all in wheat. We used a rapid, two-step cloning-out procedure for isolating the wild-type genes corresponding to those inactivated by the Tn5 insertion in strain CHA625. This cloning method should be widely applicable to bacterial genes tagged with Tn5. The region cloned from P. fluorescens contained three complete open reading frames. The deduced gene products, designated PqqFAB, showed extensive similarities to proteins involved in the biosynthesis of pyrroloquinoline quinone (PQQ) in Klebsiella pneumoniae, Acinetobacter calcoaceticus, and Methylobacterium extorquens. PQQ-negative mutants of strain CHA0 were constructed by gene replacement. They lacked glucose dehydrogenase activity, could not utilize ethanol as a carbon source, and showed a strongly enhanced production of Plt on malt agar. These effects were all reversed by complementation with pqq+ recombinant plasmids. The growth of a pqqF mutant on ethanol and normal Plt production were restored by the addition of 16 nM PQQ. However, the Phl- phenotype of strain CHA625 was due not to the pqq defect but presumably to a secondary mutation. In conclusion, a lack of PQQ markedly stimulates the production of Plt in P. fluorescens.

Genome-wide association study identifies eight loci associated with blood pressure.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N ≤ 71,225 European ancestry, N ≤ 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10(-24)), CYP1A2 (P = 1 × 10(-23)), FGF5 (P = 1 × 10(-21)), SH2B3 (P = 3 × 10(-18)), MTHFR (P = 2 × 10(-13)), c10orf107 (P = 1 × 10(-9)), ZNF652 (P = 5 × 10(-9)) and PLCD3 (P = 1 × 10(-8)) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.