Unfinished Genome Sequence Research Articles

Non-model organisms, a species endangered by proteogenomics

Previously, large-scale proteomics was possible only for organisms whose genomes were sequenced, meaning the most common model organisms. The use of next-generation sequencers is now changing the deal. With "proteogenomics", the use of experimental proteomics data to refine genome annotations, a higher integration of omics data is gaining ground. By extension, combining genomic and proteomic data is becoming routine in many research projects. "Proteogenomic"-flavored approaches are currently expanding, enabling the molecular studies of non-model organisms at an unprecedented depth. Today draft genomes can be obtained using next-generation sequencers in a rather straightforward way and at a reasonable cost for any organism. Unfinished genome sequences can be used to interpret tandem mass spectrometry proteomics data without the need for time-consuming genome annotation, and the use of RNA-seq to establish nucleotide sequences that are directly translated into protein sequences appears promising. There are, however, certain drawbacks that deserve further attention for RNA-seq to become more efficient. Here, we discuss the opportunities of working with non-model organisms, the proteomic methods that have been used until now, and the dramatic improvements proffered by proteogenomics. These put the distinction between model and non-model organisms in great danger, at least in terms of proteomics! Model organisms have been crucial for in-depth analysis of cellular and molecular processes of life. Focusing the efforts of thousands of researchers on the Escherichia coli bacterium, Saccharomyces cerevisiae yeast, Arabidopsis thaliana plant, Danio rerio fish and other models for which genetic manipulation was possible was certainly worthwhile in terms of fundamental and invaluable biological insights. Until recently, proteomics of non-model organisms was limited to tedious, homology-based techniques, but today draft genomes or RNA-seq data can be straightforwardly obtained using next-generation sequencers, allowing the establishment of a draft protein database for any organism. Thus, proteogenomics opens new perspectives for molecular studies of non-model organisms, although they are still difficult experimental organisms. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Proteogenomics for environmental microbiology

Proteogenomics sensu stricto refers to the use of proteomic data to refine the annotation of genomes from model organisms. Because of the limitations of automatic annotation pipelines, a relatively high number of errors occur during the structural annotation of genes coding for proteins. Whether putative orphan sequences or short genes encoding low-molecular-weight proteins really exist is still frequently a mystery. Whether start codons are well defined is also an open debate. These problems are exacerbated for genomes of microorganisms belonging to poorly documented genera, as related sequences are not always available for homology-guided annotation. The functional annotation of a significant proportion of genes is also another well-known issue when annotating environmental microorganisms. High-throughput shotgun proteomics has recently greatly evolved, allowing the exploration of the proteome from any microorganism at an unprecedented depth. The structural and functional annotation process may be usefully complemented with experimental data. Indeed, proteogenomic mapping has been successfully performed for a wide variety of organisms. Specific approaches devoted to systematically establishing the N-termini of a large set of proteins are being developed. N-terminomics is giving rise to datasets of experimentally proven translational start codons as well as validated peptide signals for secreted proteins. By extension, combining genomic and proteomic data is becoming routine in many research projects. The proteomic analysis of organisms with unfinished genome sequences, the so-called composite proteomics, and the search for microbial biomarkers by bottom-up and top-down combined approaches are some examples of proteogenomic-flavored studies. They illustrate the advent of a new era of environmental microbiology where proteomics and genomics are intimately integrated to answer key biological questions.

High Throughput Sequencing and Proteomics to Identify Immunogenic Proteins of a New Pathogen: The Dirty Genome Approach

BackgroundWith the availability of new generation sequencing technologies, bacterial genome projects have undergone a major boost. Still, chromosome completion needs a costly and time-consuming gap closure, especially when containing highly repetitive elements. However, incomplete genome data may be sufficiently informative to derive the pursued information. For emerging pathogens, i.e. newly identified pathogens, lack of release of genome data during gap closure stage is clearly medically counterproductive.Methods/Principal FindingsWe thus investigated the feasibility of a dirty genome approach, i.e. the release of unfinished genome sequences to develop serological diagnostic tools. We showed that almost the whole genome sequence of the emerging pathogen Parachlamydia acanthamoebae was retrieved even with relatively short reads from Genome Sequencer 20 and Solexa. The bacterial proteome was analyzed to select immunogenic proteins, which were then expressed and used to elaborate the first steps of an ELISA.Conclusions/SignificanceThis work constitutes the proof of principle for a dirty genome approach, i.e. the use of unfinished genome sequences of pathogenic bacteria, coupled with proteomics to rapidly identify new immunogenic proteins useful to develop in the future specific diagnostic tests such as ELISA, immunohistochemistry and direct antigen detection. Although applied here to an emerging pathogen, this combined dirty genome sequencing/proteomic approach may be used for any pathogen for which better diagnostics are needed. These genome sequences may also be very useful to develop DNA based diagnostic tests. All these diagnostic tools will allow further evaluations of the pathogenic potential of this obligate intracellular bacterium.

Mining the Moraxella catarrhalis Genome: Identification of Potential Vaccine Antigens Expressed during Human Infection

Moraxella catarrhalis is an important cause of respiratory infections in adults and otitis media in children. Developing an effective vaccine would reduce the morbidity, mortality, and costs associated with such infections. An unfinished genome sequence of a strain of M. catarrhalis available in the GenBank database was analyzed, and open reading frames predicted to encode potential vaccine candidates were identified. Three genes encoding proteins having molecular masses of approximately 22, 75, and 78 kDa (designated Msp [Moraxella surface proteins]) (msp22, msp75, and msp78, respectively) were determined to be conserved by competitive hybridization using a microarray, PCR, and sequencing of the genes in clinical isolates of M. catarrhalis. The genes were transcribed when M. catarrhalis was grown in vitro. These genes were amplified by PCR and cloned into Escherichia coli expression vectors. Recombinant proteins were generated and then studied using enzyme-linked immunosorbent assays with preacquisition and postclearance serum and sputum samples from 31 adults with chronic obstructive pulmonary disease (COPD) who acquired and cleared M. catarrhalis. New antibody responses to the three proteins were observed for a small proportion of the patients with COPD, indicating that these proteins were expressed during human infection. These studies indicate that the Msp22, Msp75, and Msp78 proteins, whose genes were discovered using genome mining, are highly conserved among strains, are expressed during human infection with M. catarrhalis, and represent potential vaccine antigens.

Identification, isolation, and analysis of a gene cluster involved in iron acquisition by Pseudomonas mendocina ymp

Microbial acquisition of iron from natural sources in aerobic environments is a little-studied process that may lead to mineral instability and trace metal mobilization. Pseudomonas mendocina ymp was isolated from the Yucca Mountain Site for long-term nuclear waste storage. Its ability to solubilize a variety of Fe-containing minerals under aerobic conditions has been previously investigated but its molecular and genetic potential remained uncharacterized. Here, we have shown that the organism produces a hydroxamate and not a catecholate-based siderophore that is synthesized via non-ribosomal peptide synthetases. Gene clustering patterns observed in other Pseudomonads suggested that hybridizing multiple probes to the same library could allow for the identification of one or more clusters of syntenic siderophore-associated genes. Using this approach, two independent clusters were identified. An unfinished draft genome sequence of P. mendocina ymp indicated that these mapped to two independent contigs. The sequenced clusters were investigated informatically and shown to contain respectively a potentially complete set of genes responsible for siderophore biosynthesis, uptake, and regulation, and an incomplete set of genes with low individual homology to siderophore-associated genes. A mutation in the cluster's pvdA homolog (pmhA) resulted in a siderophore-null phenotype, which could be reversed by complementation. The organism likely produces one siderophore with possibly different isoforms and a peptide backbone structure containing seven residues (predicted sequence: Acyl-Asp-Dab-Ser-fOHOrn-Ser-fOHorn). A similar approach could be applied for discovery of Fe- and siderophore-associated genes in unsequenced or poorly annotated organisms.

A protein database constructed from low-coverage genomic sequence of Bacillus megaterium and its use for accelerated proteomic analysis

Peptide mass fingerprint (PMF) matching is a high-throughput method used for protein spot identification in connection with two-dimensional gel electrophoresis (2DE). However, the success of PMF matching largely depends on whether the proteins to be identified exist in the database searched. Consequently, it is often necessary to apply other more sophisticated but also time-consuming technologies to generate sequence-tags for definitive protein identification. On the other hand, modern sequencing technologies are generating a large quantity of DNA sequences, first in unfinished form or with low genome coverage due to the time-consuming and thus limiting steps of finishing and annotation. We recently started to sequence the genome of Bacillus megaterium DSM 319, a bacterium of industrial interest. In this study, we demonstrate that a protein database generated from merely three-fold coverage, unfinished genomic sequences of this bacterium allows a fast and reliable protein spot identification solely based on PMF from high-throughput MALDI-TOF MS analysis. We further show that the strain-specific protein database from low coverage genomic sequence greatly outperforms the commonly used cross-species databases constructed from 13 completely sequenced Bacillus strains for protein spot identification via PMF.

Large-scale enhancer detection in the zebrafish genome

Murine retroviral vectors carrying an enhancer detection cassette were used to generate 95 transgenic lines of fish in which reporter expression is observed in distinct patterns during embryonic development. We mapped 65 insertion sites to the as yet unfinished zebrafish genome sequence. Many integrations map close to previously known developmental genes, including transcription factors of the Pax, Hox, Sox, Pou, Otx, Emx, zinc-finger and bHLH gene families. In most cases, the activated provirus is located in, or within a 15 kb interval around, the corresponding transcriptional unit. The exceptions include four insertions into a gene desert on chromosome 20 upstream of sox11b, and an insertion upstream of otx1. In these cases, the activated insertions are found at a distance of between 32 kb and 132 kb from the coding region. These as well as seven other insertions described here identify genes that have recently been associated with ultra conserved non-coding elements found in all vertebrate genomes.

Projector 2: contig mapping for efficient gap-closure of prokaryotic genome sequence assemblies

With genome sequencing efforts increasing exponentially, valuable information accumulates on genomic content of the various organisms sequenced. Projector 2 uses (un)finished genomic sequences of an organism as a template to infer linkage information for a genome sequence assembly of a related organism being sequenced. The remaining gaps between contigs for which no linkage information is present can subsequently be closed with direct PCR strategies. Compared with other implementations, Projector 2 has several distinctive features: a user-friendly web interface, automatic removal of repetitive elements (repeat-masking) and automated primer design for gap-closure purposes. Moreover, when using multiple fragments of a template genome, primers for multiplex PCR strategies can also be designed. Primer design takes into account that, in many cases, contig ends contain unreliable DNA sequences and repetitive sequences. Closing the remaining gaps in prokaryotic genome sequence assemblies is thereby made very efficient and virtually effortless. We demonstrate that the use of single or multiple fragments of a template genome (i.e. unfinished genome sequences) in combination with repeat-masking results in mapping success rates close to 100%. The web interface is freely accessible at .

The bioinformatics resource for oral pathogens

Complete genomic sequences of several oral pathogens have been deciphered and multiple sources of independently annotated data are available for the same genomes. Different gene identification schemes and functional annotation methods used in these databases present a challenge for cross-referencing and the efficient use of the data. The Bioinformatics Resource for Oral Pathogens (BROP) aims to integrate bioinformatics data from multiple sources for easy comparison, analysis and data-mining through specially designed software interfaces. Currently, databases and tools provided by BROP include: (i) a graphical genome viewer (Genome Viewer) that allows side-by-side visual comparison of independently annotated datasets for the same genome; (ii) a pipeline of automatic data-mining algorithms to keep the genome annotation always up-to-date; (iii) comparative genomic tools such as Genome-wide ORF Alignment (GOAL); and (iv) the Oral Pathogen Microarray Database. BROP can also handle unfinished genomic sequences and provides secure yet flexible control over data access. The concept of providing an integrated source of genomic data, as well as the data-mining model used in BROP can be applied to other organisms. BROP can be publicly accessed at .

Identification of nine human-specific frameshift mutations by comparative analysis of the human and the chimpanzee genome sequences

The recent release of the draft sequence of the chimpanzee genome is an invaluable resource for finding genome-wide genetic differences that might explain phenotypic differences between humans and chimpanzees. In this paper, we describe a simple procedure to identify potential human-specific frameshift mutations that occurred after the divergence of human and chimpanzee. The procedure involves collecting human coding exons bearing insertions or deletions compared with the chimpanzee genome and identification of homologs from other species, in support of the mutations being human-specific. Using this procedure, we identified nine genes, BASE, DNAJB3, FLJ33674, HEJ1, NTSR2, RPL13AP, SCGB1D4, WBSCR27 and ZCCHC13, that show human-specific alterations including truncations of the C-terminus. In some cases, the frameshift mutation results in gene inactivation or decay. In other cases, the altered protein seems to be functional. This study demonstrates that even the unfinished chimpanzee genome sequence can be useful in identifying modification of genes that are specific to the human lineage and, therefore, could potentially be relevant to the study of the acquisition of human-specific traits.

Identification and preliminary characterization of cell-wall-anchored proteins of Staphylococcus epidermidis

Staphylococcus epidermidis is a ubiquitous human skin commensal that has emerged as a major cause of foreign-body infections. Eleven genes encoding putative cell-wall-anchored proteins were identified by computer analysis of the publicly available S. epidermidis unfinished genomic sequence. Four genes encode previously described proteins (Aap, Bhp, SdrF and SdrG), while the remaining seven have not been characterized. Analysis of primary sequences of the Staphylococcus epidermidis surface (Ses) proteins indicates that they have a structural organization similar to the previously described cell-wall-anchored proteins from S. aureus and other Gram-positive cocci. However, not all of the Ses proteins are direct homologues of the S. aureus proteins. Secondary and tertiary structure predictions suggest that most of the Ses proteins are composed of several contiguous subdomains, and that the majority of these predicted subdomains are folded into beta-rich structures. PCR analysis indicates that certain genes may be found more frequently in disease isolates compared to strains isolated from healthy skin. Patients recovering from S. epidermidis infections had higher antibody titres against some Ses proteins, implying that these proteins are expressed during human infection. Western blot analyses of early-logarithmic and late-stationary in vitro cultures suggest that different regulatory mechanisms control the expression of the Ses proteins.

Comparison of the Yeast Proteome to Other Fungal Genomes to Find Core Fungal Genes

The purpose of this research was to search for evolutionarily conserved fungal sequences to test the hypothesis that fungi have a set of core genes that are not found in other organisms, as these genes may indicate what makes fungi different from other organisms. By comparing 6355 predicted or known yeast (Saccharomyces cerevisiae) genes to the genomes of 13 other fungi using Standalone TBLASTN at an e-value <1E-5, a list of 3340 yeast genes was obtained with homologs present in at least 12 of 14 fungal genomes. By comparing these common fungal genes to complete genomes of animals (Fugu rubripes, Caenorhabditis elegans), plants (Arabidopsis thaliana, Oryza sativa), and bacteria (Agrobacterium tumefaciens, Xylella fastidiosa), a list of common fungal genes with homologs in these plants, animals, and bacteria was produced (938 genes), as well as a list of exclusively fungal genes without homologs in these other genomes (60 genes). To ensure that the 60 genes were exclusively fungal, these were compared using TBLASTN to the major sequence databases at GenBank: NR (nonredundant), EST (expressed sequence tags), GSS (genome survey sequences), and HTGS (unfinished high-throughput genome sequences). This resulted in 17 yeast genes with homologs in other fungal genomes, but without known homologs in other organisms. These 17 core, fungal genes were not found to differ from other yeast genes in GC content or codon usage patterns. More intensive study is required of these 17 genes and other common fungal genes to discover unique features of fungi compared to other organisms.

Multilocus sequence typing of Listeria monocytogenes by use of hypervariable genes reveals clonal and recombination histories of three lineages.

In an attempt to develop a method to discriminate among isolates of Listeria monocytogenes, the sequences of all of the annotated genes from the fully sequenced strain L. monocytogenes EGD-e (serotype 1/2a) were compared by BLASTn to a file of the unfinished genomic sequence of L. monocytogenes ATCC 19115 (serotype 4b). Approximately 7% of the matching genes demonstrated 90% or lower identity between the two strains, and the lowest observed identity was 80%. Nine genes (hisJ, cbiE, truB, ribC, comEA, purM, aroE, hisC, and addB) in the 80 to 90% identity group and two genes (gyrB and rnhB) with approximately 97% identity were selected for multilocus sequence analysis in two sets of L. monocytogenes isolates (a 15-strain diversity set and a set of 19 isolates from a single food-processing plant). Based on concatenated sequences, a total of 33 allotypes were differentiated among the 34 isolates tested. Population genetics analyses revealed three lineages of L. monocytogenes that differed in their history of apparent recombination. Lineage I appeared to be completely clonal, whereas representatives of the other lineages demonstrated evidence of horizontal gene transfer and recombination. Although most of the gene sequences for lineage II strains were distinct from those of lineage I, a few strains with the majority of genes characteristic of lineage II had some that were characteristic of lineage I. Genes from lineage III organisms were mostly similar to lineage I genes, with instances of genes appearing to be mosaics with lineage II genes. Even though lineage I and lineage II generally demonstrated very distinct sequences, the sequences for the 11 selected genes demonstrated little discriminatory power within each lineage. In the L. monocytogenes isolate set obtained from one food-processing plant, lineage I and lineage II were found to be almost equally prevalent. While it appears that different lineages of L. monocytogenes can share habitats, they appear to differ in their histories of horizontal gene transfer.

Methylotrophy in Methylobacterium extorquens AM1 from a genomic point of view.

Methylotrophy is defined as the ability to “grow at the expense of reduced carbon compounds containing one or more carbon atoms but containing no carbon-carbon bonds” (3). It is an intriguing example of microbial metabolic agility, with the use of a class of chemicals disregarded by the majority of organisms. Even though the ability to grow methylotrophically was first discovered in the early 1900s (cited in reference 3), it was not until the 1960s to 1970s that an understanding of the biochemical nature of this capability started to emerge. Fascination with methylotrophy in those years was fueled by the commercial interest in single-cell protein production, and as a result, the specific details of the biochemistry of methylotrophy began to be revealed. Enzymes for the primary oxidation of C1 substrates such as methanol dehydrogenase and methylamine dehydrogenase were characterized, and distinct modes of C1 assimilation, such as the ribulose monophosphate cycle and the serine cycle were discovered. The biochemical processes involved in methylotrophy that were known by 1982 are described in detail in the now classic book Biochemistry of Methylotrophs by Christopher Anthony (3). In the 20 years following the publication of Biochemistry of Methylotrophs, a few additional methylotrophy biochemical pathways have been discovered, such as the pathway for C1 transfer linked to methanopterin and methanofuran, which solved the long-standing mystery of formaldehyde oxidation in many methylotrophs (15, 53), and novel pathways for primary C1 oxidation, such as the pathways for degradation of chlorinated methanes and methanesulfonic acid (21, 50). The knowledge concerning the biochemistry and physiology of methylotrophic organisms accumulated over the past three decades suggests a new framework for understanding methylotrophy as a novel metabolic mode. In this framework, methylotrophy is envisioned as a set of specific metabolic functional modules, with different combinations of such modules being present in different methylotrophs (Fig. ​(Fig.11 for the methylotrophic metabolic modules in Methylobacterium extorquens AM1). However, until recently, a number of important details of these modules were missing, and so the picture remained incomplete. The availability of two unfinished genome sequences for the important model organisms M. extorquens AM1 (http://www.integratedgenomics.com/genomereleases.html#list6)and Methylococcus capsulatus Bath (http://tigrblast.tigr.org/ufmg/) is transforming our understanding of methylotrophy. Annotation of these two genomes combined with functional analysis will delineate the set of genes and functions that is both sufficient and necessary to define a methylotroph. Expanding genomic analyses to include other groups of methylotrophs will in turn provide clues to the origins of methylotrophy and the evolution of various methylotrophic pathways. In this publication, we summarize the existing knowledge of the genes involved in methylotrophic pathways in M. extorquens AM1, analyze its yet unfinished genome with respect to location and clustering of methylotrophy genes, and present a comprehensive list of methylotrophy genes and enzymes known at this time in M. extorquens AM1 (Table ​(Table11). FIG. 1. Methylotrophy metabolic modules in M. extorquens AM1. Known genes are in italic. For simplicity, redox reactions in the assimilatory pathways are not indicated. For details, refer to the references given in Table ​Table11. TABLE 1. Methylotrophy genes in M. extorquens AM1

Apical location of a novel EGF-like domain-containing protein of Plasmodium falciparum

Using bioinformatics analyses of the unfinished malaria genome sequence, we have identified a novel protein of Plasmodium falciparum that contains two epidermal growth factor (EGF)-like domains near the C-terminus of the protein. The sequence contains a single open reading frame of 1572bp with the potential to encode a protein of 524 residues containing hydrophobic regions at the extreme N- and C-termini which appear to represent signal peptide and glycosylphosphatidylinositol (GPI)-attachment sites, respectively. RT-PCR analysis has confirmed that the novel gene is transcribed in asexual stages of P. falciparum. Antibodies to the EGF-like domains of the novel protein are highly specific and do not cross-react with the EGF-like domains of MSP1, MSP4, MSP5 or MSP8 expressed as GST fusion proteins. Antisera to the C-terminal fragments react with two bands of 80 and 36kDa in P. falciparum parasite lysates whereas antisera to the most N-terminal fusion protein only recognises the 80kDa band, suggesting that the novel protein may undergo processing in a similar way to MSP1 and MSP8, but with fewer cleavage events. Immunoblot analysis of stage-specific parasite samples reveals that the protein is present in trophozoites, schizonts and in isolated merozoites. The protein partitions in the detergent-enriched phase after Triton X-114 fractionation and is localised to the surfaces of trophozoites, schizonts and free merozoites in an apical distribution. Based on the accepted nomenclature in the field we now designate this protein MSP10. We have shown that the MSP10 fusion proteins are in a conformation that can be recognised by human immune sera and that there is very limited sequence diversity in an approximately lkb region of MSP10, encompassing the two EGF-like domains. A sequence similar to MSP10 can be identified in the available P. yoelii genomic sequence, offering the possibility of ascertaining whether this novel protein can induce host protective responses in an in vivo model.

Genomic BLAST: custom-defined virtual databases for complete and unfinished genomes.

BLAST (Basic Local Alignment Search Tool) searches against DNA and protein sequence databases have become an indispensable tool for biomedical research. The proliferation of the genome sequencing projects is steadily increasing the fraction of genome-derived sequences in the public databases and their importance as a public resource. We report here the availability of Genomic BLAST, a novel graphical tool for simplifying BLAST searches against complete and unfinished genome sequences. This tool allows the user to compare the query sequence against a virtual database of DNA and/or protein sequences from a selected group of organisms with finished or unfinished genomes. The organisms for such a database can be selected using either a graphic taxonomy-based tree or an alphabetical list of organism-specific sequences. The first option is designed to help explore the evolutionary relationships among organisms within a certain taxonomy group when performing BLAST searches. The use of an alphabetical list allows the user to perform a more elaborate set of selections, assembling any given number of organism-specific databases from unfinished or complete genomes. This tool, available at the NCBI web site http://www.ncbi.nlm.nih.gov/cgi-bin/Entrez/genom_table_cgi, currently provides access to over 170 bacterial and archaeal genomes and over 40 eukaryotic genomes.

Genomic BLAST: custom-defined virtual databases for complete and unfinished genomes

BLAST (Basic Local Alignment Search Tool) searches against DNA and protein sequence databases have become an indispensable tool for biomedical research. The proliferation of the genome sequencing projects is steadily increasing the fraction of genome-derived sequences in the public databases and their importance as a public resource. We report here the availability of Genomic BLAST, a novel graphical tool for simplifying BLAST searches against complete and unfinished genome sequences. This tool allows the user to compare the query sequence against a virtual database of DNA and/or protein sequences from a selected group of organisms with finished or unfinished genomes. The organisms for such a database can be selected using either a graphic taxonomy-based tree or an alphabetical list of organism-specific sequences. The first option is designed to help explore the evolutionary relationships among organisms within a certain taxonomy group when performing BLAST searches. The use of an alphabetical list allows the user to perform a more elaborate set of selections, assembling any given number of organism-specific databases from unfinished or complete genomes. This tool, available at the NCBI web site http://www.ncbi.nlm.nih.gov/cgi-bin/Entrez/genom_table_cgi, currently provides access to over 170 bacterial and archaeal genomes and over 40 eukaryotic genomes.

Molecular analysis of the dormancy response in Mycobacterium smegmatis: expression analysis of genes encoding the DevR–DevS two-component system, Rv3134c and chaperone α-crystallin homologues

Mycobacteria adapt to a decrease in oxygen tension by entry into a non-replicative persistent phase. It was shown earlier that the two-component system, DevR–DevS, was induced in Mycobacterium tuberculosis and Mycobacterium bovis BCG cultures during hypoxia, suggesting that it may play a regulatory role in their adaptation to oxygen limitation. The presence of a homologous genetic system in Mycobacterium smegmatis was predicted by scanning its unfinished genome sequence with devR and devS genes of M. tuberculosis. Rv3134c, which is cotranscribed with devR–devS in M. tuberculosis, was also present in M. smegmatis at a similar location upstream from devR. The expression of all three genes was induced at the RNA and protein levels in M. smegmatis cultures grown under microaerobic and anaerobic conditions. The M. smegmatis genome also contained the hspX gene, encoding chaperone α-crystallin, Acr, that was induced during hypoxia. The similarity in sequences and hypoxia-responsive behaviour of devR–devS, Rv3134c and hspX genes in M. smegmatis and M. tuberculosis suggests that the molecular mechanisms involved in the dormancy response are likely conserved in these two species. M. smegmatis could therefore serve as a useful model for the delineation of the hypoxia response in general and DevR–DevS regulated pathways in particular.

Molecular analysis of the dormancy response in Mycobacterium smegmatis: expression analysis of genes encoding the DevR-DevS two-component system, Rv3134c and chaperone alpha-crystallin homologues.

Mycobacteria adapt to a decrease in oxygen tension by entry into a non-replicative persistent phase. It was shown earlier that the two-component system, DevR-DevS, was induced in Mycobacterium tuberculosis and Mycobacterium bovis BCG cultures during hypoxia, suggesting that it may play a regulatory role in their adaptation to oxygen limitation. The presence of a homologous genetic system in Mycobacterium smegmatis was predicted by scanning its unfinished genome sequence with devR and devS genes of M. tuberculosis. Rv3134c, which is cotranscribed with devR-devS in M. tuberculosis, was also present in M. smegmatis at a similar location upstream from devR. The expression of all three genes was induced at the RNA and protein levels in M. smegmatis cultures grown under microaerobic and anaerobic conditions. The M. smegmatis genome also contained the hspX gene, encoding chaperone alpha-crystallin, Acr, that was induced during hypoxia. The similarity in sequences and hypoxia-responsive behaviour of devR-devS, Rv3134c and hspX genes in M. smegmatis and M. tuberculosis suggests that the molecular mechanisms involved in the dormancy response are likely conserved in these two species. M. smegmatis could therefore serve as a useful model for the delineation of the hypoxia response in general and DevR-DevS regulated pathways in particular.

Identification of six novel genes by experimental validation of GeneMachine predicted genes

In silico gene identification from finished and unfinished human genome sequence has become critically important in many projects seeking to gain insights into the gene content of genomic regions implicated in diseases. To establish limitations and criteria for in silico gene identification, and to identify novel genes of potential relevance to human prostate cancer and melanoma, 3 Mb of chromosome 1 sequence have been analyzed using GeneMachine. This program is a software suite comprising of sequence similarity programs and four gene identification programs. A total of 49 potential transcripts were selected and 37 of them were selected for experimental validation. We verified 16 of the predicted genes by experimental analysis. The comparison of the predicted transcripts with their cloned forms helped to refine predicted gene models as well as to identify splice variants for several of them. Although sequences matching with ten of our verified genes have been recently deposited in the GenBank, six of them remain novel. Our studies support the feasibility of identifying novel genes from regions of interest using draft human genome sequence.