Dataset Information

Consistency of gene starts among Burkholderia genomes.

ABSTRACT: BACKGROUND: Evolutionary divergence in the position of the translational start site among orthologous genes can have significant functional impacts. Divergence can alter the translation rate, degradation rate, subcellular location, and function of the encoded proteins. RESULTS: Existing Genbank gene maps for Burkholderia genomes suggest that extensive divergence has occurred--53% of ortholog sets based on Genbank gene maps had inconsistent gene start sites. However, most of these inconsistencies appear to be gene-calling errors. Evolutionary divergence was the most plausible explanation for only 17% of the ortholog sets. Correcting probable errors in the Genbank gene maps decreased the percentage of ortholog sets with inconsistent starts by 68%, increased the percentage of ortholog sets with extractable upstream intergenic regions by 32%, increased the sequence similarity of intergenic regions and predicted proteins, and increased the number of proteins with identifiable signal peptides. CONCLUSIONS: Our findings highlight an emerging problem in comparative genomics: single-digit percent errors in gene predictions can lead to double-digit percentages of inconsistent ortholog sets. The work demonstrates a simple approach to evaluate and improve the quality of gene maps.

SUBMITTER: Dunbar J

PROVIDER: S-EPMC3049151 | biostudies-literature | 2011

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Consistency of gene starts among Burkholderia genomes.

Dunbar John J Cohn Judith D JD Wall Michael E ME

BMC genomics 20110222

<h4>Background</h4>Evolutionary divergence in the position of the translational start site among orthologous genes can have significant functional impacts. Divergence can alter the translation rate, degradation rate, subcellular location, and function of the encoded proteins.<h4>Results</h4>Existing Genbank gene maps for Burkholderia genomes suggest that extensive divergence has occurred--53% of ortholog sets based on Genbank gene maps had inconsistent gene start sites. However, most of these in ...[more]

PMID: 21342528

Similar Datasets

Project description:BACKGROUND:Structural variations in human genomes, such as insertions, deletion, or rearrangements, play an important role in cancer development. Next-Generation Sequencing technologies have been central in providing ways to detect such variations. Most existing methods however are limited to the analysis of a single genome, and it is only recently that the comparison of closely related genomes has been considered. In particular, a few recent works considered the analysis of data sets obtained by sequencing both tumor and healthy tissues of the same cancer patient. In that context, the goal is to detect variations that are specific to exactly one of the genomes, for example to differentiate between patient-specific and tumor-specific variations. This is a difficult task, especially when facing the additional challenge of the possible contamination of healthy tissues by tumor cells and conversely. RESULTS:In the current work, we analyzed a data set of paired-end short-reads, obtained by sequencing tumor tissues and healthy tissues, both from the same cancer patient. Based on a combinatorial notion of conflict between deletions, we show that in the tumor data, more deletions are predicted than there could actually be in a diploid genome. In contrast, the predictions for the data from normal tissues are almost conflict-free. We designed and applied a method, specific to the analysis of such pooled and contaminated data sets, to detect potential tumor-specific deletions. Our method takes the deletion calls from both data sets and assigns reads from the mixed tumor/normal data to the normal one with the goal to minimize the number of reads that need to be discarded to obtain a set of conflict-free deletion clusters. We observed that, on the specific data set we analyze, only a very small fraction of the reads needs to be discarded to obtain a set of consistent deletions. CONCLUSIONS:We present a framework based on a rigorous definition of consistency between deletions and the assumption that the tumor sample also contains normal cells. A combined analysis of both data sets based on this model allowed a consistent explanation of almost all data, providing a detailed picture of candidate patient- and tumor-specific deletions.

Project description:The flagellar gene fliO of Salmonella typhimurium can be translated from an AUG codon that overlaps the termination codon of fliN (K. Ohnishi et al., J. Bacteriol. 179:6092-6099, 1997). However, it had been concluded on the basis of complementation analysis that in Escherichia coli a second start codon 60 bp downstream was the authentic one (J. Malakooti et al., J. Bacteriol. 176:189-197, 1994). This raised the possibility of tandem translational starts, such as occur for the chemotaxis gene cheA; this possibility was increased by the existence of a stem-loop sequence covering the second start, a feature also found with cheA. Protein translated from the first start codon was detected regardless of whether the second start codon was present; it was also detected when the stem-loop structure was disrupted or deleted. Translation from the second start codon, either as the natural one (GUG) or as AUG, was not detected when the first start and intervening sequence were intact. Nor was it detected when the first codon was attenuated (by conversion of AUGAUG to AUAAUA; in S. typhimurium there is a second, adjacent, AUG) or eliminated (by conversion to CGCCGC); disruption of the stem-loop structure still did not yield detectable translation from the second start. When the entire sequence up to the second start was deleted, translation from the second start was detected provided the natural codon GUG had been converted to AUG. A fliO null mutant could be fully complemented in swarm assays whenever the first start and intervening sequence were present, regardless of the state of the second start. Reasonably good complementation occurred when the first start and intervening sequence were absent provided the second start was intact, either as AUG or as GUG; thus translation from the GUG codon must have been occurring even though protein levels were too low to be detected. The translated intervening sequence is rather divergent between S. typhimurium and E. coli and corresponds to a substantial cytoplasmic domain prior to the sole transmembrane segment, which is highly conserved; the sequence following the second start begins immediately prior to that transmembrane segment. The significance of the data for FliO is discussed and compared to the equivalent data for CheA. Attention is also drawn to the fact that given an optimal ribosome binding site, AUA can serve as a fairly efficient start codon even though it seldom if ever appears to be used in nature.

Dataset Information

Consistency of gene starts among Burkholderia genomes.

Publications

Consistency of gene starts among Burkholderia genomes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure