An intragenic distribution bias of DNA uptake sequences in Pasteurellaceae and Neisseriae
© van Passel; licensee BioMed Central Ltd. 2008
Received: 21 February 2008
Accepted: 27 March 2008
Published: 27 March 2008
Most sequenced strains from Pasteurellaceae and Neisseriae contain hundreds to thousands of uptake sequence (US) motifs in their genome, which are associated with natural competence for DNA uptake. The mechanism of their recognition is still unclear, and I searched for intragenic location patterns of these motifs for clues about their distribution. In all cases, one orientation of the US has a higher occurrence in the reading frame, and in all Pasteurellaceae, the US and the reverse complement motifs are biased towards the gene termini. These findings could help design experimental set-ups to study preferential DNA uptake, thereby further unravelling the phenomenon of natural competence.
This article was reviewed by Arcady Mushegian and I. King Jordan.
Gene repertoires are not stable features of prokaryotes , and it has been suggested that most genes have been affected by horizontal gene transfer at some point during their evolutionary history . Environmental DNA not only serves as a potential nutrient for bacteria , but also represents an important substrate for the acquisition of new coding sequences via incorporation into the genome . The ability to take up DNA and thereby acquiring new genes from the environment is termed natural competence, but this trait is not shared between all bacteria to the same extent. This implies that for some bacteria, or in selected environments, natural competence is a more favourable ability.
In the Pasteurellaceae and Neisseriaea, natural competence has been linked to specific DNA motifs that are substantially overrepresented in the genome. The DNA Uptake Sequence (DUS) 5'-ATGCCGTCTGAA-3' in Neisseria species occurs almost 1,500 times in the genome [5, 6], and represents approximately 2% of the genome. The Uptake Signal Sequence (USS) in most Pasteurellaceae is 5'-AAAGTGCGGT-3' [7, 8], although a different motif is encountered in a subclade of this family . If the conspecific DNA uptake sequence (or US, as I refer to them collectively) is present in the donor DNA, it increases the transformation efficiency many fold . The exact mechanism how these motifs enhance DNA uptake remains unclear. It has been postulated that USs may be involved in some sort of transformation-barrier, as they selectively increase DNA uptake from a source containing the same US , even though this US can be shared by distantly related Pasteurellaceae species . Alternatively, these genomically overrepresented motifs were suggested to be associated with the replication and repair mechanisms due to their genome distribution , but no general orientation bias has ever been observed for these sequences in the genome sequences in which they reside. Another study proposed a role for transcription termination, due to the occurrence of inverted repeat forms of the intergenic USs, allowing the formation of hairpins . However, as a substantial fraction of the USs are not in this inverted repeat configuration, and because many USs are intragenic, the transcription terminator explanation is still contentious .
Despite research into the genomic distribution of these sequence motifs [6, 14], no analyses have focused on the distribution patterns of USs within the protein coding regions. Recently, we identified an intragenic positional preference of homopolymeric tracts in prokaryotes, which suggested that there has been intensified selection against these mutagenic repeats in the middle of genes . In the current study I test whether USs also show a skewed distribution in the protein coding genes of the Pasteurellaceae and Neisseriae.
The US in Neisseria is unrelated to the motifs found in the Pasteurellaceae, and due to its specific sequence, it hinders straightforward distribution pattern analyses. Initially, I found a strong bias towards the 5' end of protein coding genes, until I detected that dozens of annotated open reading frames in neisserial genomes actually start with the US. In contrary to the US in the Pasteurellaceae, the neisserial US motif starts with the potential start codon ATG, which could result in incorrect annotations. Eliminating the 30 genes in N. meningitidis MC58 that start with an US shows a different pattern altogether; even though many more rcUS are present in protein coding genes than US (245 vs. 123, respectively, p < 0.05), there was no similarity with the distribution bias of the USs from the Pasteurellaceae (Additional File 3). However, the USs and rcUSs motifs in Neisseriae show different, though not significant (p > 0.05), distributions, with the USs apparently underrepresented in the 5' end of genes.
In the first ever genome-wide sequence motif analysis, Smith and co-workers found that the H. influenzae US was distributed apparently randomly over the genome , similar to later observations in the genome of A. succinogenes . However, a somewhat smaller than expected number of these USs was found within the coding fraction of the H. influenzae genome . This observation was recently complemented by Davidsen and co-workers who found a biased distribution of USs towards genome maintenance genes, such as replication, recombination and DNA repair, in different representative species from the phylogenetically distinct Pasteurellaceae and Neisseria . They proposed that the observed bias of USs towards this class of genes was somehow connected to their function in these unrelated taxonomic groups, indicating convergent evolution. In this study I show that also within coding regions, the distribution of distinct competence-associated sequences is non-random. The increased incidence of these motifs in gene termini of the tested Pasteurellaceae species suggests that their presence at these locations could convey a selective benefit, as these distributions are thought to have arisen independently. The exact nature of this proposed selective benefit is unknown, but these observations could give rise to an experimental setup to examine potential differences in transformation efficiency, as genes flanked by USs in different orientations could be acquired with different efficiencies. However, experiments in Actinobacillus actinomycetemcomitans do not point in this direction .
Alternatively, the observed location biases may be caused by potential misannotations of transcriptional start sites (TSS). As TSS are still mostly assigned by bioinformatic procedures , inadequate annotations could be responsible for a skew in the distribution patterns. However, I do not expect this to negate the biased occurrence of US motifs in protein coding genes, as the observed bias is not restricted to the 5' end of genes.
Referee 1: Arcady Mushegian, Stowers Institute, Kansas City, USA.
AM1: The GC-content of both US motifs is exactly 50%, which is close to the average GC-content of a Neisseria genome (51–52%?), but is not quite close to those of Pasteurellaceae (for example, 35–38% in Haemophilus).
But more relevant are the GC-contents values of each of the 5 quintiles of the ORF lengths (i.e., what is plotted in Fig. 1). Can all or some of the biases in the US distribution be explained by the bias of the underlying nucleotide frequency distribution?
MWJvP1: I tested this for the Haemophilus influenzae Rd and Pasteurella multocida genome sequences. Considering all protein coding genes (>100 bp) in H. influenzae, I find that, on average, the GC% for the 5 quintiles are 37.8%, 39.0%, 38.9%, 38.8% and 37.3%, respectively. As the US has a GC% of 50%, one would assume to actually find lower numbers of these motifs in the first and last quintile of genes, based solely on the GC%. Instead, I find higher numbers. For P. multocida, a similar GC% progression is observed (39.5%, 41.1%, 41.1%, 41.3% and 39.8%, for each of the five quintiles, respectively).
AM2: All biases and overrepresentations should be accompanied by the appropriate statistics.
MWJvP2: My apologies for not including these in the first draft. I apply a standard Chi-square goodness of fit test (4 degrees of freedom), with a null-hypothesis that these USs are represented proportionally (i.e., equally) in each quintile. I find for H. influenzae Rd a p-value < 0.05 that this distribution represents an equally proportionate distribution, so this null-hypothesis is rejected. I expand this chi-square for all sampled genomes (Additional File 1) and find that for the rcUS in every tested Pasteurellaceae genome the observed distribution is significantly different from the expected distribution (H0 = equally proportioned). As for the US, 5 out of 8 tested genomes show no significant distribution bias (p > 0.05), although similar patterns are observed with respect to their location preference.
I conclude that the tested rcUS motifs are not proportionally distributed in the gene quintiles. Only the full 10 bp motif was tested, and not the shorter (incomplete) sequence motifs, which show less of a bias, and were merely used in Figure 1to illustrate that the observed bias is not caused by a subset of the US sequence.
As for the counts of the US and the rcUS, I find for all genomes except H. somnus that they occur in significantly non-equal amounts in the protein coding regions of the tested genomes.
For Neisseria, a significantly different number of USs and rcUSs is observed in open reading frames (p < 0.05), but there is no significant difference (p > 0.05) between the intragenic location biases of the 12 bp US and the rcUS, although there is for the 10 bp US (p < 0.05).
Referee 2: I. King Jordan, Georgia Tech University, School of Biology, Atlanta, USA.
MWJ van Passel reports on biases in the distributions of uptake sequences (US) that mediate natural competence in Pasteurellaceae and Neisseriae. Within genes, the reverse complement of the US (rcUS) is found more frequently than the US, and the rcUS is biased towards the 5' ends of genes while the US is biased towards the 3' ends of genes. This is a noteworthy discovery and I support publication of the report as a Discovery Note in Biology Direct. I have several questions related mainly to the clarity of presentation.
IKJ1: The most important point is that the biases in the distributions need to be supported by some statistical analyses. Some sort of goodness-of-fit such as chi-square with an appropriate correction for multiple tests should suffice.
MWJvP1: As reviewer 1 also asked for this addition, I would like to refer to the answer given above (MWJvP2 to AM2). I hope this adequately addresses this point.
IKJ2: According to the author, intergenic US distributions have been evaluated previously, but this is the first report on the occurrences of intragenic US. To underscore the importance of this difference, it would help the reader to known what the relative frequencies of intergenic versus intragenic US for the species analyzed are.
MWJvP2: It was Hamilton O. Smith et al. in 1995, who noted that only 65% of all USs occur in the ORFs, whereas ~86% of the H. influenzae Rd genome had been annotated as coding sequences. But Redfield et al. (2006) already counted fewer of these US (1115 instead of 1465). Therefore, we checked for the five genomes depicted in Figure 1the total counts of the (rc)US in the CDS and in the entire genome. For the genomes of H. influenzae Rd, H. somnus, P. multocida, M. succiniciproducens and A. sucinogenes, the percentages of (rc)USs in the coding regions are 59%, 70%, 55%, 55% and 47%, respectively, whereas their coding fractions represent 84%, 88%, 88%, 89% and 87%, respectively (Additional File 4depict these data). This suggests a non-random distribution of DNA uptake sequence motifs with respect to the coding fractions, similar as to what Hamilton O. Smith et al. found.
IKJ3: In the discussion, the distribution of US across different functional classes of H. influenzae genes is mentioned. This is a very interesting point and a similar analysis for Pasteurellaceae and Neisseriae would strengthen this manuscript. Along these same lines, I was wondering if US sequences are over-represented among genes that have been horizontally transferred and/or are strain-specific.
MWJvP3: This distribution of (rc)USs has been analyzed thoroughly by Davidsen et al. (2004), which was in fact performed for both Pasteurellaceae and Neisseriae (I hope to have clarified this now in the text). They concluded that for both Pasteurellaceae and Neisseriae, the motifs are more common in the functional group of genome maintenance genes (GMGs). I quote "These results imply that the high frequency of DUS in genome maintenance genes is conserved among phylogenetically divergent species and thus are of significant biological importance" (Davidsen et al., 2004). There is however some debate whether these genes are frequently transferred: no aberrant nucleotide composition has been detected, but then again, it has been suggested that these US promote conspecific DNA uptake. Therefore, I did not carry out an analysis for US occurrence in horizontally transferred genes in this project.
IKJ4: In Figure 1, it would helpful to indicate which are the US and which are the rcUS.
MWJvP4: A reference to which lanes are depict the USs and which depict the rcUSs is now incorporated in the legend of Figure 1.
MWJvP5: I am sorry that this wasn't mentioned properly in the legend of Figure 1. In this figure, four lanes per genome are visible: the two on the left are the US with increasing stringency (the motif is depicted below the lane, 8 bp and 10 bp), and the two on the right are the rcUS with increasing stringency. I hope the legend is now clearer.
IKJ6: The intragenic distributions of US and rcUS in H. influenzae are reported as 'data not shown'. Given that Biology Direct is an online journal, it is preferable to report these distributions as Supplementary Information.
MWJvP6: As there are several different strains of H. influenzae sequenced, I thought it better not to include these strains next to the H. influenzae rd KW20, as they do not represent independent measurements. However, I did expect a reviewer to ask whether these patterns are also found in the other strains, which they are (although not as pervasive with respect to the US distribution). I have included Additional File 1with the counts of the motifs in all the Pasteurellaceae genomes I tested, so readers also have access to the individual counts per genome.
The author gratefully acknowledges Rosie Redfield, Kevin Vogel and Rogier van Doorn for helpful discussions concerning earlier versions of this manuscript.
- Lerat E, Daubin V, Ochman H, Moran NA: Evolutionary origins of genomic repertoires in bacteria. PLoS Biol. 2005, 3 (5): e130-10.1371/journal.pbio.0030130.PubMedPubMed CentralView ArticleGoogle Scholar
- Dagan T, Martin W: Ancestral genome sizes specify the minimum rate of lateral gene transfer during prokaryote evolution. Proc Natl Acad Sci USA. 2007, 104 (3): 870-875. 10.1073/pnas.0606318104.PubMedPubMed CentralView ArticleGoogle Scholar
- Redfield RJ: Genes for breakfast: the have-your-cake-and-eat-it-too of bacterial transformation. J Hered. 1993, 84 (5): 400-404.PubMedGoogle Scholar
- Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405 (6784): 299-304. 10.1038/35012500.PubMedView ArticleGoogle Scholar
- Goodman SD, Scocca JJ: Identification and arrangement of the DNA sequence recognized in specific transformation of Neisseria gonorrhoeae. Proc Natl Acad Sci USA. 1988, 85 (18): 6982-6986. 10.1073/pnas.85.18.6982.PubMedPubMed CentralView ArticleGoogle Scholar
- Ambur OH, Frye SA, Tonjum T: New functional identity for the DNA uptake sequence in transformation and its presence in transcriptional terminators. J Bacteriol. 2007, 189 (5): 2077-2085. 10.1128/JB.01408-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Danner DB, Deich RA, Sisco KL, Smith HO: An eleven-base-pair sequence determines the specificity of DNA uptake in Haemophilus transformation. Gene. 1980, 11 (3–4): 311-318. 10.1016/0378-1119(80)90071-2.PubMedView ArticleGoogle Scholar
- Redfield RJ, Findlay WA, Bosse J, Kroll JS, Cameron AD, Nash JH: Evolution of competence and DNA uptake specificity in the Pasteurellaceae. BMC Evol Biol. 2006, 6: 82-10.1186/1471-2148-6-82.PubMedPubMed CentralView ArticleGoogle Scholar
- Sisco KL, Smith HO: Sequence-specific DNA uptake in Haemophilus transformation. Proc Natl Acad Sci USA. 1979, 76 (2): 972-976. 10.1073/pnas.76.2.972.PubMedPubMed CentralView ArticleGoogle Scholar
- Scocca JJ, Poland RL, Zoon KC: Specificity in deoxyribonucleic acid uptake by transformable Haemophilus influenzae. J Bacteriol. 1974, 118 (2): 369-373.PubMedPubMed CentralGoogle Scholar
- Bakkali M, Chen TY, Lee HC, Redfield RJ: Evolutionary stability of DNA uptake signal sequences in the Pasteurellaceae. Proc Natl Acad Sci USA. 2004, 101 (13): 4513-4518. 10.1073/pnas.0306366101.PubMedPubMed CentralView ArticleGoogle Scholar
- Karlin S, Mrazek J, Campbell AM: Frequent oligonucleotides and peptides of the Haemophilus influenzae genome. Nucleic Acids Res. 1996, 24 (21): 4263-4272. 10.1093/nar/24.21.4263.PubMedPubMed CentralView ArticleGoogle Scholar
- Bakkali M: Genome dynamics of short oligonucleotides: the example of bacterial DNA uptake enhancing sequences. PLoS ONE. 2007, 2 (1): e741-10.1371/journal.pone.0000741.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith HO, Tomb JF, Dougherty BA, Fleischmann RD, Venter JC: Frequency and distribution of DNA uptake signal sequences in the Haemophilus influenzae Rd genome. Science. 1995, 269 (5223): 538-540. 10.1126/science.7542802.PubMedView ArticleGoogle Scholar
- van Passel MW, Ochman H: Selection on the genic location of disruptive elements. Trends Genet. 2007, 23 (12): 601-604. 10.1016/j.tig.2007.08.017.PubMedView ArticleGoogle Scholar
- Wang Y, Orvis J, Dyer D, Chen C: Genomic distribution and functions of uptake signal sequences in Actinobacillus actinomycetemcomitans. Microbiology. 2006, 152 (Pt 11): 3319-3325. 10.1099/mic.0.29018-0.PubMedView ArticleGoogle Scholar
- Davidsen T, Rodland EA, Lagesen K, Seeberg E, Rognes T, Tonjum T: Biased distribution of DNA uptake sequences towards genome maintenance genes. Nucleic Acids Res. 2004, 32 (3): 1050-1058. 10.1093/nar/gkh255.PubMedPubMed CentralView ArticleGoogle Scholar
- Rison SC, Mattow J, Jungblut PR, Stoker NG: Experimental determination of translational starts using peptide mass mapping and tandem mass spectrometry within the proteome of Mycobacterium tuberculosis. Microbiology. 2007, 153 (Pt 2): 521-528. 10.1099/mic.0.2006/001537-0.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.