Filling annotation gaps in yeast genomes using genome-wide contact maps.

MOTIVATIONS: De novo sequencing of genomes is followed by annotation analyses aiming at identifying functional genomic features such as genes, non-coding RNAs or regulatory sequences, taking advantage of diverse datasets. These steps sometimes fail at detecting non-coding functional sequences: for example, origins of replication, centromeres and rDNA positions have proven difficult to annotate with high confidence. Here, we demonstrate an unconventional application of Chromosome Conformation Capture (3C) technique, which typically aims at deciphering the average 3D organization of genomes, by showing how functional information about the sequence can be extracted solely from the chromosome contact map. RESULTS: Specifically, we describe a combined experimental and bioinformatic procedure that determines the genomic positions of centromeres and ribosomal DNA clusters in yeasts, including species where classical computational approaches fail. For instance, we determined the centromere positions in Naumovozyma castellii, where these coordinates could not be obtained previously. Although computed centromere positions were characterized by conserved synteny with neighboring species, no consensus sequences could be found, suggesting that centromeric binding proteins or mechanisms have significantly diverged. We also used our approach to refine centromere positions in Kuraishia capsulata and to identify rDNA positions in Debaryomyces hansenii. Our study demonstrates how 3C data can be used to complete the functional annotation of eukaryotic genomes. AVAILABILITY AND IMPLEMENTATION: The source code is provided in the Supplementary Material. This includes a zipped file with the Python code and a contact matrix of Saccharomyces cerevisiae. CONTACT: romain.koszul@pasteur.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Selection for unequal densities of sigma70 promoter-like signals in different regions of large bacterial genomes.

The evolutionary processes operating in the DNA regions that participate in the regulation of gene expression are poorly understood. In Escherichia coli, we have established a sequence pattern that distinguishes regulatory from nonregulatory regions. The density of promoter-like sequences, that could be recognizable by RNA polymerase and may function as potential promoters, is high within regulatory regions, in contrast to coding regions and regions located between convergently transcribed genes. Moreover, functional promoter sites identified experimentally are often found in the subregions of highest density of promoter-like signals, even when individual sites with higher binding affinity for RNA polymerase exist elsewhere within the regulatory region. In order to see the generality of this pattern, we have analyzed 43 additional genomes belonging to most established bacterial phyla. Differential densities between regulatory and nonregulatory regions are detectable in most of the analyzed genomes, with the exception of those that have evolved toward extreme genome reduction. Thus, presence of this pattern follows that of genes and other genomic features that require weak selection to be effective in order to persist. On this basis, we suggest that the loss of differential densities in the reduced genomes of host-restricted pathogens and symbionts is an outcome of the process of genome degradation resulting from the decreased efficiency of purifying selection in highly structured small populations. This implies that the differential distribution of promoter-like signals between regulatory and nonregulatory regions detected in large bacterial genomes confers a significant, although small, fitness advantage. This study paves the way for further identification of the specific types of selective constraints that affect the organization of regulatory regions and the overall distribution of promoter-like signals through more detailed comparative analyses among closely related bacterial genomes.

A novel mycobacterial antigen relevant to cellular immunity belongs to a family of secreted lipoproteins.

The gene sequence of a novel 24.1 kDa Mycobacterium tuberculosis protein was identified within the Sanger Centre (UK) M. tuberculosis genome database (cosmid MTCY24G1) by searching with a 126 bp DNA sequence isolated from a genomic M. leprae lambda gt11 library with M. leprae reactive human T cell clones as probes. The 24.1 kDa antigen is common to the vaccine strain Mycobacterium bovis BCG, as well as Mycobacterium leprae. The 699 bp open reading frame encodes a 233 amino acid long precursor protein with a signal peptide sequence for secretion and a consensus motif for lipid conjugation, which suggests that the mature protein is an exported lipoprotein antigen. The molecular mass of the mature protein antigen from M. leprae sonicate was shown to correspond to the deduced size of the M. tuberculosis protein by T cell Western analysis. Homology searches revealed two other similarly sized hypothetical secreted mycobacterial lipoproteins within the M. tuberculosis genome database.

Interactions of OxyR with the promoter region of the oxyR and ahpC genes from Mycobacterium leprae and Mycobacterium tuberculosis.

In contrast to the intact oxyR gene (a homolog of the central regulator of peroxide stress response in enteric bacteria) in Mycobacterium leprae, this gene is inactive in all strains of M. tuberculosis. In both species, oxyR is divergently transcribed from ahpC, which encodes a homolog of alkyl hydroperoxide reductase. To initiate investigations of the regulation of oxidative stress in mycobacteria and consequences of the elimination of oxyR in M. tuberculosis, in this work we tested the hypothesis that mycobacterial OxyR acts as a DNA binding protein and analyzed its interactions with the oxyR and ahpC promoters. M. leprae OxyR was overproduced and purified, and its binding to the oxyR-ahpC intergenic region of M. leprae was demonstrated. By using a sequential series of overlapping DNA fragments, the minimal OxyR binding site was delimited to a 30-bp DNA segment which included a palindromic sequence conforming with the established rules for the LysR family of regulators. A consensus sequence for the mycobacterial OxyR recognition site (cTTATCggc-N3-gccGATAAg) was deduced based on its conservation in different mycobacteria. A variance in two potentially critical nucleotides within this site was observed in M. tuberculosis, in keeping with its reduced affinity for OxyR. Transcription of plasmid-borne M. leprae oxyR and ahpC was investigated in M. smegmatis and M. bovis BCG by S1 nuclease protection and transcriptional fusion analyses. Two mRNA 5' ends were detected in each direction: (i) P1oxyR and P2oxyR and (ii) P1ahpC and P2ahpC. The binding site for OxyR overlapped P1oxyR, reminiscent of the autoregulatory loops controlling expression of oxyR in enteric bacteria and characteristic of the LysR superfamily in general. This site was also centered 65 bp upstream of P1ahpC, matching the usual position of LysR-type recognition sequences in relationship to positively controlled promoters. Superimposed on these features was the less orthodox presence of multiple transcripts and their unique arrangement, including a region of complementarity at the 5' ends of the P2ahpC and P2oxyR mRNAs, suggesting the existence of complex regulatory relationships controlling oxyR and ahpC expression in mycobacteria.

Genetic analysis of the Mycobacterium smegmatis rpsL promoter.

The DNA sequence of the promoter region of the Mycobacterium smegmatis rpsL gene, which encodes the S12 ribosomal protein, was determined. Primer extension analysis and S1 nuclease protection experiments identified the 5' end of the rpsL mRNA to be 199 bp upstream of the translation initiation codon. The rpsL promoter contained sequences upstream of this start point for transcription that were similar to the canonical hexamers found at the -10 and -35 regions of promoters recognized by Esigma70, the major form of RNA polymerase in Escherichia coli. To define the promoter of the rpsL gene, DNA fragments containing progressive deletions of the upstream region of the rpsL gene were inserted into a plasmid vector containing a promoterless xylE gene. These insertions revealed that the 200 bp of DNA sequence immediately upstream from the translation initiation codon was not essential for promoter function. In addition, 5' deletions removing all but 34 bp upstream of the transcription start point retained greater than 90% promoter activity, suggesting that the -35 hexamer was not essential for promoter activity. To determine which nucleotides were critical for promoter function, oligonucleotide-directed mutagenesis and mutagenic PCR amplification were used to produce point mutations in the region upstream of the start point of transcription. Single base substitutions in the -10 hexamer, but not in the -35 hexamer, severely reduced rpsL promoter activity in vivo. Within the -10 hexamer, nucleotide substitutions causing divergence from the E. Coli sigma70 consensus reduced promoter activity. The DNA sequence immediately upstream from the - 10 hexamer contained the TGn motif described as an extended -10 region in prokaryotic promoters. Mutations in this motif, in combination with a transition at either the -38 or -37 position within the -35 hexamer, severely reduced promoter activity, indicating that in the absence of a functional -35 region, the rpsL promoter is dependent on the TGn sequence upstream from the -10 hexamer. Comparison of the nucleotide sequence of the rpsL promoter region of M. smegmatis with the homologous sequences from Mycobacterium leprae, Mycobacterium bovis, and Mycobacterium tuberculosis showed the presence in these slowly growing mycobacterial species of conserved promoter elements a similar distance upstream of the translation initiation codon of the rpsL gene, but these other mycobacterial promoters did not contain the extended -10 motif.

Activity of mycobacterial promoters during intracellular and extracellular growth.

pUS933, a bifunctional Mycobacterium-Escherichia coli translational fusion vector containing an amino-terminally truncated E. coli lacZ reporter gene, was constructed. Derivatives of pUS933, containing the promoter, RBS and start codon of the Mycobacterium bovis BCG hsp60 gene, the Mycobacterium leprae 28 kDa gene and the M. leprae 18 kDa gene were constructed and introduced into E. coli, Mycobacterium smegmatis and M. bovis BCG. beta-Galactosidase activity was measured for mycobacteria grown in liquid culture. Primer-extension analysis was used to determine the transcriptional start point for the 18 kDa promoter in M. smegmatis. Murine macrophages were infected with recombinant BCG containing the pUS933 derivatives and expression levels were examined, by fluorescence microscopy and fluorometry, during intracellular growth of BCG. Both the BCG hsp60 gene promoter and the M. leprae 28 kDa gene promoter gave high levels of beta-galactosidase expression in all situations examined. In contrast, the M. leprae 18 kDa promoter fragment gave very low levels of expression in M. smegmatis and BCG grown in liquid culture, but in BCG growing within macrophages it was induced to levels almost as high as the other promoters. This indicated that the 18 kDa gene is specifically activated during intracellular growth and may therefore be involved in survival of M. leprae within macrophages. This pattern of regulation may be useful for controlling expression of foreign genes in recombinant BCG strains.

The Mycobacterium leprae genome: systematic sequence analysis identifies key catabolic enzymes, ATP-dependent transport systems and a novel polA locus associated with genomic variability.

In the framework of the mycobacterial genome sequencing project, a continuous 37,049 bp sequence from the Mycobacterium leprae chromosome has been determined. Computer analysis revealed 10 complete open reading frames, and nine of their products show similarity to known proteins. Seven of these were identified as the enzyme isocitrate lyase, two P-type ATPase cation transporters, two AMP-binding proteins, the ribosomal protein S1, and DNA polymerase I. Interestingly, the polA gene, encoding DNA polymerase, is flanked by two inverted copies of a new class of the M. leprae specific repetitive sequence, RLEP, and this structure resembles a transposable element. A second copy of this element was found at another locus in the genome, but the two copies were not present in equal amounts and could not be found in all isolates of M. leprae. This is the first evidence for genomic variability in the leprosy bacillus and might ultimately be useful for developing a molecular test capable of distinguishing between strains of M. leprae.

Isolation and characterization of the aspartokinase and aspartate semialdehyde dehydrogenase operon from mycobacteria.

Diaminopimelic acid (DAP) is a major component of the peptidoglycan layer of the mycobacterial cell wall. The mycobacterial cell wall has been implicated as a potential virulence factor and is highly immunogenic. The pathway for biosynthesis of DAP may serve as a target in the design of antimycobacterial agents and construction of in vivo selection systems. Despite its significance, this biosynthetic pathway is poorly understood in mycobacteria. In order to develop a better understanding of mycobacterial DAP biosynthesis, the aspartate semialdehyde dehydrogenase (asd) genes of Mycobacterium smegmatis, bacille Calmette-Guerin (BCG), Mycobacterium avium, Mycobacterium leprae, and Mycobacterium tuberculosis were isolated. The M. smegmatis asd gene was isolated by complementation in Escherichia coli. This gene was then used to isolate the asd genes from other mycobacteria. The asd-complementing fragments from BCG and M. smegmatis were sequenced. An open reading frame upstream of the mycobacterial asd gene was identified as the mycobacterial aspartokinase gene (ask). Primer extension analysis revealed that the only transcriptional start in this region is found 5' of the ask gene. This observation indicates that the mycobacterial ask and asd genes are in an operon.

Isolation and characterization of an autonomously replicating sequence (ARSD) from the marine yeast Debaryomyces hansenii.

The marine yeast Debaryomyces hansenii is known to tolerate salinities ranging from 0 to 24%. As a first step toward the molecular analysis of halotolerance in this organism, we report the isolation of an autonomously replicating sequence (ARS) and its use in the construction of a shuttle vector. The ARS from D. hansenii (ARSD) is 0.4 kbp long, and the function rests in 0.13 kbp of the sequence. Sequence analysis of ARSD shows strong homology to ARS from other organisms, including a 12-bp consensus sequence common to all ARS functional in Saccharomyces cerevisiae.