- Research article
- Open Access
An in silico analysis of T-box regulated genes and T-box evolution in prokaryotes, with emphasis on prediction of substrate specificity of transporters
BMC Genomics volume 9, Article number: 330 (2008)
T-box anti-termination is an elegant and sensitive mechanism by which many bacteria maintain constant levels of amino acid-charged tRNAs. The amino acid specificity of the regulatory element is related to a so-called specifier codon and can in principle be used to guide the functional annotation of the genes controlled via the T-box anti-termination mechanism.
Hidden Markov Models were defined to search the T-box regulatory element and were applied to all completed prokaryotic genomes. The vast majority of the genes found downstream of the retrieved elements encoded functionalities related to transport and synthesis of amino acids and the charging of tRNA. This is completely in line with findings reported in literature and with the proposed biological role of the regulatory element. For several species, the functional annotation of a large number of genes encoding proteins involved in amino acid transport could be improved significantly on basis of the amino acid specificity of the identified T-boxes. In addition, these annotations could be extrapolated to a larger number of orthologous systems in other species. Analysis of T-box distribution confirmed that the element is restricted predominantly to species of the phylum Firmicutes. Furthermore, it appeared that the distribution was highly species specific and that in the case of amino acid transport some boxes seemed to "pop-up" only recently.
We have demonstrated that the identification of the molecular specificity of a regulatory element can be of great help in solving notoriously difficult annotation issues, e.g. by defining the substrate specificity of genes encoding amino acid transporters on basis of the amino acid specificity of the regulatory T-box. Furthermore, our analysis of the species-dependency of the occurrence of specific T-boxes indicated that these regulatory elements propagate in a semi-independent way from the genes that they control.
Transcription anti-termination is a regulatory mechanism commonly encountered in all lineages within the bacterial kingdom (see e.g. [1–3]). In transcription anti-termination, the regulation of transcription occurs after the initiation of RNA synthesis, but before transcription of the coding region. The mechanism of anti-termination involves a structural change in the RNA transcript that is dependent on the interaction of the transcript with, for instance, a regulatory protein , a tRNA  or a metabolite . The structural elements that compose these anti-terminators are encoded by conserved sequences on the DNA and can be found by searches for the related sequence motifs in upstream regions of regulated genes .
A well-studied anti-termination element is the so-called T-box. T-box anti-termination is an elegant and sensitive mechanism by which many bacteria maintain constant levels of tRNA charged with amino acids [2, 8]. When there is a sufficient supply of charged tRNA in a cell, the T-box folds into a terminator structure, thereby blocking further transcription. Transcription can only proceed upon conversion into an anti-terminator structure, which is induced by binding of a highly conserved 5'-NCCA-3' of the uncharged tRNA with a conserved '5-UGGN-3' sequence in the T-box [9–11]. Although anti-terminator formation involves contacts between many nucleotides, the specificity of the interaction seems largely dependent on the interaction of a tri-nucleotide (anti-anti)-codon in the so-called specifier loop of the T-box with the anti-codon of an amino acid-specific tRNA [12–15]. The structural and kinetic details of this interaction have been well-studied [16–22]. The appropriate assignment of the specifier codon has been used previously to improve the functional annotation of various genes located downstream of the T-box [3, 23–31]. The T-box controlled genes identified thus far encode functionalities that reflect perfectly the pivotal role of uncharged tRNAs in the regulatory mechanism. These functionalities include not only tRNA ligation, but also amino acid biosynthesis and transport [11, 17, 24, 26, 28, 32–34]. The encoded proteins are involved in modulation of the level of uncharged tRNA in the cell, either directly by charging the corresponding tRNA with its cognate amino acid or indirectly by controlling the intracellular concentration of the specific amino acid.
To date, T-boxes have been identified predominantly in the genomes of bacterial species of the phyla Firmicutes (including Mollicutes) and Actinobacteria , although anti-termination systems have been argued to be among the oldest regulatory systems in bacteria because of their independence of regulatory proteins . To investigate this further, we have explored the occurrence of T-boxes in all sequenced prokaryotic genomes. To circumvent potential differences between T-box systems in different bacterial lineages, an iterative HMM-based identification search was performed using the best conserved region of the T-box sequence. Species- and amino acid-specific T-box regulation networks were reconstructed. Most importantly, the acquired knowledge on amino acid specificity could be used to propose an improved functional annotation for many T-box controlled genes and to shed light on the evolution of the regulatory element itself.
Results and discussion
I) A comprehensive collection of T-boxes
The analysis of the taxonomic and functional distribution of T-boxes was started by de novo identification of T-box motif characteristics. Conserved nucleotide sequence motifs upstream of tRNA-ligase encoding genes in species of the phylum Firmicutes were recovered and used to identify T-boxes located at other positions in the same genome as well as in the genomes of other species (see methods). These searches showed that a T-box could be specified best by a 30 nt motif that is extremely well-conserved and positioned in the 3'-region of the terminator/anti-terminator loop (motif 1 in Figure 1). In fact, this motif is known as 'the T-box sequence' since its discovery . Later it was recognized that this conserved region belongs to a larger conserved RNA structure known as the T-box element . This element contains four other highly conserved regions (see Figure 1 motifs 2–5).
The initial search showed prominent variations in the number of T-boxes per genome between different classes of the phylum Firmicutes and between different phyla. Therefore, additional searches with phylum-specific and class-specific T-box HMMs were performed, but generally did not yield novel hits. Only in the case of the Clostridia a limited number of 10 additional T-boxes were identified. Further iterations did not expand the dataset. Visual inspection of the upstream regions of all genes encoding a t-RNA ligase in the Firmicutes indicated that indeed all those regions that contain the distinctive T-box motifs were identified by our algorithm. A comparison of the number of T-boxes identified by us for a representative set of organisms with the number obtained using the Rfam T-box model  proved that our recovery procedure was very efficient (see methods for details).
Identification of the specifier codon and amino acid specificity
Although T-boxes were readily identified, it was more difficult to define their amino acid specificity. To that end, the phylogeny of homologous genes preceded by a T-box from different species was determined and the upstream regions corresponding to each orthologous group were aligned. In all the cases of T-boxes for which the specifier codon had been identified experimentally [3, 5, 9, 11, 12, 14–18, 21, 28, 29, 33], we observed that the specifier codon aligned perfectly within the related orthologous sequences. In fact, this was true for almost all orthologous groups of sequences. Moreover, most of the alignments could easily be clustered by eye into larger groups for which the specifier codon remained directly apparent from the alignment. The resulting alignments and the annotation of the specifier codon can be found at . Nevertheless, there remained a few less clear cases. For some (~5%) a secondary structure prediction could be used to provide the additional information required to define the specifier codon in the specifier loop . Taken together, a specifier codon could be identified directly for over 90% of the identified T-boxes.
Codon usage in the specifier codon
Most amino acids are encoded by multiple codons. Leu for instance, is encoded by six different codons (CUA, CUU, CUG, CUC, UUA and UUG). Remarkably, the T-boxes had a conserved preference for certain codons within as well as between species (Additional file 1). Evaluation of these preferences showed that they complied almost perfectly with the rules observed by Elf et al. for the codon usage by E. coli . In an elegant study these authors analyzed the dependence of the charging of various codon-specific tRNAs on the use of various codons in particular proteins. They concluded that: "when codon reading is part of a control loop that regulates synthesis of missing amino acid, the translation rate of the selected codon should be as sensitive as possible to starvation" . And, in their paper they showed which codons are the most sensitive in E. coli. We found that for all but one of the most predominantly used specifier codons in T-boxes, the corresponding tRNA is among the highest in sensing shortage of that specific amino acid in E. coli as reported by . The only exception was the T-box codon for Ala (GCU). Therefore, assuming the conclusions by Elf et al. are also valid for Gram-positive bacteria, our findings suggest that the codons that are sensitive to depletion are preferentially used in T-box regulation.
Functionalities controlled by a T-box
As expected, the proposed regulatory role of the T-box elements appeared to be perfectly reflected by the genes under their control. The majority of the T-boxes (62%) were found to precede genes encoding tRNA ligases, while most others were found upstream of genes encoding proteins involved in amino acid transport (12%) or amino acid biosynthesis (18%). The remaining T-boxes (8%; 71 genes in total) were found upstream of genes encoding proteins with unknown function (54 genes), or a function that lacks an apparent relation to amino acid metabolism (17 genes). A complete and species-specific subdivision of T-boxes based on function prediction of the proteins encoded downstream and a list of genes with no apparent relation to amino acid metabolism is provided in the supplementary material, in Additional files 1 and 2.
II) The use of regulator specificity to improve annotation of molecular function and biological role
We made two important observations: i) In all cases, the T-boxes identified upstream of the genes encoding a tRNA ligase contained a specifier codon that corresponded with the amino acid specificity of the ligase; and ii) in all other cases where the function of the protein encoded by the gene downstream of the T-box had experimentally been verified, the specifier-codon corresponded to the established functionality of the gene. These observations implied that the employed method for the identification of T-box specificity was reliable and, consequently, that predicted T-box specificities could be extrapolated to the molecular function of the protein encoded by the gene located downstream, as had occasionally been done before. Many of the genes preceded by a T-box had not been specifically annotated to date in the sense that, although the functional category was often evident (e.g. proton symport, ABC transport family, etc.), a specific molecular function had not been attributed. In fact, more than two-third of the non-tRNA ligase genes preceded by a T-box lacked such a specific annotation of molecular function. As importantly, the functional annotation of the genes could be extended to a different level entirely by using the knowledge on T-box (regulator) specificity, as this knowledge discloses (in part) under which conditions the regulated genes will play their biological role (for the distinction between molecular function and biological role see Francke et al. ).
a) T-box regulation of amino acid transport
Many of the genes encoding amino acid transporters were found to be preceded by a T-box, especially in the genomes of the Lactobacilli and Bacilli of the Bacillus cereus-group. The transporters controlled by T-boxes belonged to no less than seven distinct transporter families (MFS, 2.A.1; APC, 2.A.3; NSS, 2.A.22; DAACS, 2.A.23; LIVCS, 2.A.26; NhaC, 2.A.35; and ABC-cassette, 3.A.1; Transporter Classification described by Saier ). Table 1 gives an overview of the distribution of the transport systems regulated by a T-box over the various Firmicutes species.
Overall, for more than 85% of the T-box regulated transporters the functional annotation (molecular function and/or biological role) could be improved as compared to the entries in the reference database of NCBI. A full list can be found in Additional file 1. We have limited the substrate specificity definition in our annotation to putatively dominant substrates based on the amino acid specificity of the T-box. However, broader substrate specificity is probably more common for transporters. Especially transport systems consisting of only a permease are expected (and have been shown) to display broader substrate specificity (see  and  for examples), whereas systems that require prior substrate-binding (like in ABC transport) will be more specific. We discuss the T-box based functional annotation of some transport systems in more detail in the following paragraphs and in Additional file 3.
The ABC family
T-box regulation of ABC transport systems was found in most lineages of the Firmicutes but not in the Bacilli. The T-box regulated ABC transporters could be sub-divided into four sub-families, based on the specificity of the substrate-binding protein and the permease. A striking use of extensive T-box regulation in ABC transport was observed in L. plantarum. It appears that in the absence of methionine, L. plantarum uses a single mechanism to switch on not only transport of the amino acid itself, but also of the precursors and co-factors needed for its biosynthesis (see Figure 2 and Additional file 3).
The APC family
A T-box was identified in front of an APC-family protein encoding gene in all the studied Bacillus genomes. In B. subtilis and B. licheniformis the gene ybvW is preceded by a Leu T-box, whereas such a box is lacking upstream of the orthologous genes, which are found in E. faecalis, G. kaustophilus and in L. lactis (co-orthologs: yibG and ysjA) (Figure 3). The Leu T-box suggests that the YbvW protein is a Leucine transporter, in line with the general functionality of transporters of the APC family (family characteristics described in [44, 45]). Surprisingly, in the members of the Bacillus cereus-group another APC family gene is preceded by a T-box, specific for threonine. Although an orthologous gene is present in most of the Firmicutes genomes, e.g. ykbA in B. subtilis, it is regulated by a T-box only in the species of the Bacillus cereus-group (Figure 3). The protein encoded by ykbA in B. subtilis has recently been shown to be a Ser/Thr exchanger and was consequently renamed SteT . A similar functionality of the protein ortholog in the members of the Bacillus cereus group is supported by the codon identification of the T-box.
The LIVCS family
The Bacilli of the Bacillus cereus group, the Lactobacilli and the Clostridia contain several branched-chain amino acid cation symporters of the LIVCS-family , some of which are T-box regulated (Figure 3). Since the three branched-chain amino acids share very similar molecular properties (e.g. size and hydrophobicity) we expect that these transporters are not highly specific despite their proposed amino acid specific control, but merely that expression of the "multi-specific" system has been brought under the control of the individual amino acids. Indeed, the orthologous transporters that have been characterized in L. delbrueckii (BrnQ; ), C. glutamicum (BrnQ; [47, 48]) and P. aeruginosa (BraZ; ) displayed transport of all three branched-chain amino acids.
The NSS family
Finally, the LIVCS family branched-chain amino acid transporter BraZ of P. aeruginosa, was shown to have a clear preference for isoleucine and valine over leucine . In this respect it is noteworthy that in the Bacilli of the Bacillus cereus group the expression of one of the homologs of the NSS family (neurotransmitter:sodium symport) is controlled by a Leu T-box, whereas such a box is lacking for the LIVCS homologs in those species (Figure 4). Besides the Leu T-box regulated NSS transporter, the Bacilli of the Bacillus cereus group contain three other homologs of the same family, two of which are controlled by a Trp and a Phe T-box, respectively (Figure 3). These two amino acids agree well with the experimentally determined tryptophan transport functionality of the NSS homolog TnaT in S. thermophilum . The presence of a regulatory T-box ranging from Leu to Trp, and Phe suggest that the members of the NSS transporter family may display a rather broad amino acid specificity.
b) Hypothetical proteins controlled by a T-box
Another class of proteins to which significant functional information could be added (when compared to the NCBI-annotation) using the specificity of the detected T-box is that of the so-called hypothetical proteins or unknown function proteins (data accumulated in Table 2). Obviously, when orthologous proteins in related species were also of unknown function, specifier codon information clearly improved the annotation. Examples of new annotations related to amino acid biosynthesis or transport are enzymes (methionine synthase, cystathionine gamma synthase, chorismate mutase, anthranilate synthase), transporters (Leu-, Lys- and His-specific permeases), tRNA-ligase related functions, and regulation (anti-TRAP protein).
III) Taxonomic variation and T-box evolution
The comprehensive list of T-boxes that was generated for all sequenced genomes (see Table 3 for the phylogenetic distribution) confirmed the previous attribution that T-boxes are predominantly encountered in species of the phylum Firmicutes (>95% of the hits). Our analyses uncovered many previously unidentified T-boxes. In species of the class Mollicutes two T-box elements were found, but only in the subclass Endoplasmatales , in Proteobacteria (i.e. in Geobacter sulfurreducens and Pelobacter carbinolicus) a typical T-box element was identified upstream of the leuA gene (2-isopropylmalate synthase) in both species. Deinococcus radiodurans contained two T-boxes (related to ile and gly t-RNA ligase) whereas species of the phylum of Chloroflexi (Dehalococcoides CBDB1 and Dehalococcoides ethenogenes 195) contained three T-box elements (one related to an ile tRNA ligase and two related to tryptophan biosynthesis (trpE and trpB-like). Earlier analysis of riboswitches in Actinobacteria showed that some species belonging to this phylum contain a T-box upstream of ileS . However, Symbiobacterium thermophilum contained not less than eighteen T-boxes, comparable to species of the Firmicutes. This finding is in line with the conclusion of  that S. thermophilum is probably more closely related to Firmicutes than to Actinobacteria.
To evaluate the phylogenetic distribution of T-boxes in more detail, the correlation between the presence of T-box regulatory elements and the regulated genes was analyzed for the Firmicutes and will be described shortly in the next sections. Furthermore, the scattered appearance of these regulatory boxes as observed for the various transporter families will be discussed.
T-box regulation of genes encoding tRNA ligases in the Firmicutes
It appeared (Figure 4) that regulation by T-boxes is conserved in almost all (at least 29 out of 34) Firmicutes for several tRNA ligases (ileS, a laS, serS and thrS), whereas some tRNA ligases (lysS, asnS and gltX; the latter gene encodes a Glu-tRNA ligase (see also the legend of Figure 4)) appeared to be controlled by a T-box in only a few species. The genes encoding the tRNA ligases for cysteine and asparagine were often found as the second gene in a putative operon that was T-box regulated. In several organisms, multiple copies of amino-acid specific tRNA-ligase encoding genes are found and in more than half of the cases (58%) only one of them is subject to T-box regulation.
A clear phylogenetic effect is observed when the four major orders within the Firmicutes (Bacillales, Clostridia, Lactobacillales and Mollicutes) are compared. This is true for both the number and the type of tRNA-ligase encoding genes regulated by a T-box. Most T-box regulated tRNA-ligase encoding genes are found in the Bacillales and especially in species of the Bacillus cereus group. Within the Lactobacillales, it appears that Streptococci have far less tRNA- ligase encoding genes regulated by T-boxes than Lactobacillus species. The lowest number of tRNA ligases regulated by T-boxes was found for S. thermophilus, S. pneumoniae and some strains of S. pyogenes. The relatively low amount of T-box regulation in these species could be the result of regressive evolution, a process that was suggested to be the underlying mechanism for the large loss of functionally active genes in S. thermophilus .
T-box regulation of genes involved in amino acid biosynthesis in the Firmicutes
T-box regulation of genes related to amino acid biosynthesis has been described previously for various amino acids [3, 23, 24]. Like in the case of the t-RNA ligases, T-box control of amino acid biosynthesis displayed clear phylogenetic patterns (Figure 5). For instance, the biosynthesis of Branched Chain Amino Acids (BCA: isoleucine, leucine, valine) was found to be T-box regulated in Bacillales and Clostridia, whereas several families within the Bacillales (e.g. Staphylococci and Listeria) as well as several Streptococci consistently lack T-box control of BCA biosynthesis. Similarly, we found that the species of the B. cereus group contain a T-box in the upstream region of the tyrosine biosynthesis operon consistent with the experimental data that showed that tyrosine biosynthesis from shikimate is T-box regulated in B. anthracis . The B. cereus group representatives are the only organisms in our study that encode a phenylalanine-4-hydroxylase ortholog (converting phenylalanine into tyrosine). We found this gene also to be T-box regulated in all members of the B. cereus group.
The evolution and propagation of T-boxes
An important observation related to T-box evolution was made by Grundy et al. [9, 12]. They showed that a single nucleotide change in the specifier codon of the Tyr T-box of tyrS in B. subtilis was enough to change the amino acid specificity. In addition, while analyzing the distribution of T-boxes over the various transporter families we were struck by the fact that although some of these families are very large, there was only one (or a few) family-member(s) found to be regulated by a T-box and only in a restricted number of species (see Table 1 and e.g. Figure 3). In our opinion, the only likely scenario to explain the phylogenetically limited occurrence of the transporter T-box associations that would not imply massive loss of the T-box regulation was that of acquisition of the regulatory element by the transporter encoding gene in a specific lineage. Moreover, the results of Grundy et al. [9, 12] imply that in principle the T-Box can change specificity easily.
To analyze this further, the T-boxes preceding the genes that encode the t-RNA-ligases for the branched-chain amino acids (ileS, leuS and valS) were examined. These were chosen because T-box regulation of these genes is most wide-spread (Figure 4) and because the proximal genes themselves, ileS, leuS and valS form a separate tRNA-ligase sub-family with a very clear evolutionary lineage (Figure 6A). In sharp contrast, the NJ-tree of the aligned complete T-boxes (approximately 200 – 300 nt) appeared extremely unreliable (Figure 6B). This finding implies that amino acid specificity of a T-box is not prominent on the overall sequence level and that the apparent overall sequence variability in time of the T-box is relatively high. Simultaneously, in case T-box sequences do cluster in reliable clusters (high bootstrap support) in a NJ-tree they must therefore be closely related in time.
To limit possible obscuring effects of comparing sequences between species, we collected and compared the T-box sequences within species and for clarity restricted the comparison to the T-boxes that accompany transport systems and the related tRNA ligases. For all three analyzed species: B. anthracis, L. acidophilus and L. plantarum. similar phenomena were observed. The multiple sequence alignment of the analyzed T-box sequences and the associated NJ-tree (depicted in Figure 6C for B. anthracis) are strongly suggestive of a close evolutionary relationship between several of the T-boxes. For example in B. anthracis, the Thr T-boxes found in front of BA4970 (transport systems of LIVCS-type) and the Thr-tRNA ligase (BA4820) were highly similar and the same was observed for the Ile T-boxes found in front of another LIVCS homolog and the Ile-tRNA ligase. As the LIVCS homologs appear closely related in time (apparent duplication in the B. cereus group ancestor, see Figure 3), the data imply that the regulatory T-box was not inherited in a similar way but 'acquired' independently. In fact, this explanation fits the observed scattered appearance of the T-boxes for the various transporter families perfectly.
The results presented in Figure 6 are also suggestive of another way in which the T-boxes have evolved. The NJ-tree relates the Phe T-box found in front of one of the NSS family transporters to the Tyr T-box associated with the Tyr-tRNA ligase (Figure 6C). It thus seems that the Tyr T-box of the tRNA ligase was duplicated -as this T-box is present in various Firmicutes species- and has diverged/adapted to control a Phe transporter in the Bacilli of the B. cereus group. In fact, the similarity between the T-box upstream of Tyr-tRNA ligase and the Phe T-box in front of the transporter is higher than between the Tyr T-box and the Tyr T-box preceding BA4353 (NhaC family transporter). This is consistent with the fact that: the Tyr T-box acquisition of the NhaC ortholog should have occurred earlier in history, as the Tyr T-box control of the NhaC ortholog is present in several Firmicutes, and the sequences thus had more time to diverge.
The sequence signature of a T-box is very specific and as a result T-boxes can be readily identified. Using specific T-box HMMs, we identified a large number of the T-boxes and their amino acid specificity in sequenced prokaryote genomes.
An important aspect of this work is that we show that the prediction of the amino acid specificity of the various T-boxes can be used to improve the functional annotation of a large number of genes. In particular, the functional annotation of genes related to amino acid transport and genes with unknown substrate specificity, genes for which it is normally quite difficult to find functional attributes, could be improved significantly. In our opinion, the procedure of improving annotation through knowledge of the regulatory signals can be generalized and should be used on a much broader scale than currently is being done.
Riboswitches have been argued to be among the oldest regulatory systems in bacteria because of their independence of regulatory proteins and widespread biological distribution . One might therefore have expected that T-boxes are abundantly present among all different lineages of bacteria. This clearly can not be concluded from our results and those presented in other studies [52, 55]. In fact, these regulatory elements can only be found in a few bacterial phyla and only abundantly in the phylum Firmicutes. This implies that either Firmicutes developed T-box regulation after their branching off from the other bacteria or that the other bacteria lost the system soon after the branching off of the Firmicutes to evolve more complex regulatory systems. Which of the two scenarios is most likely remains unclear.
Nevertheless our data do allow some extrapolation of the propagation of T-boxes within the phylum Firmicutes. We conclude on basis of our observations that the T-boxes have evolved in four clearly distinct ways: i) by co-evolution with the regulated gene or operon; ii) by co-evolution and divergence with the regulated gene or operon to adopt a new specificity; by iii) duplication and insertion of the regulatory element in front of a gene or operon that encodes functions related to the T-box-specified amino acid; and finally iv) by duplication and divergence toward a new amino acid specificity after duplication. In short, this means that every T-box regulatory element acts as a connected yet independent "functional module". The fact that the isoleucine specific T-box connected to the ile-tRNA ligase encoding gene is the only box present in all T-box containing bacteria suggests this box could very well compose the archetype T-box.
Sequence information and tools
Genome sequences and annotation files of completely sequenced bacterial and archaeal genomes were downloaded from the NCBI repository . In addition, for the analysis of the molecular functions encoded by T-box controlled genes, genomic information was obtained from the ERGO genome analysis and discovery system . Multiple sequence alignments were created with MUSCLE  and bootstrapped Neighbor-Joining trees with CLUSTALX  (corrected for multiple substitutions). HMMs were constructed from a multiple sequence alignment using HMMER 2.3.2 . Conserved nucleotide sequence motifs were recovered via MEME  (settings: maximally 5 motifs, modus ZOOPs, minimal width 10 nt, maximal width 30 nt) and visualized using Weblogo . Hidden Markov Models (HMMs) were constructed on basis of the recovered motifs to perform genome-wide searches. Although HMMsearch  yielded essentially identical results to the search-tool MAST , in our hands the procedure was much faster and the output was easier to interpret (including motif e-values rather than p-values). RNA secondary structure predictions were performed by Mfold  using default settings.
Identification of a general T-box sequence
The T-boxes reported so far in literature were mainly located upstream of tRNA-ligase encoding genes in species of the phylum Firmicutes [7, 26]. We therefore initiated the search for a general T-box sequence motif by collecting the nucleotide sequences (300 nt) preceding all tRNA ligases (n = 910) found in sequenced genomes of the Firmicutes. Five characteristic consecutive motifs were recovered in the nucleotide sequences (displayed in Figure 1). They were part of a generic T-box sequence that spans a length of about 250 nt, as described by . The best-conserved motif (Figure 1, motif 1, E-value = 1.2e-2841), which is located closest to the translation start, had a length of 30 nt and was found in about 61% (n = 553) of the nucleotide sequences. A genome-wide HMMsearch based upon the best-conserved motif yielded 374 new hits. To probe eventual in-homogeneity induced by the fact that only T-boxes associated with tRNA ligases were taken to create the initial T-box HMM, a new HMM based on the non-tRNA ligase associated T-boxes was made and the genomes were searched anew. We found that all t-RNA ligase associated T-boxes were recovered (with an e-value < 1) with this new HMM. This finding implied that the set of recovered T-boxes was relatively homogeneous. The upstream 500 nt of every putative best-conserved T-box motif was then checked for the presence of at least one of the four other conserved T-box specific motifs recovered using MEME. In 40 cases (less than 4% of total) none of the other motifs was detected. Without exception, these were not located in the proximity of a coding sequence and thus unlikely to be related to transcription attenuation. The hits were therefore considered false-positives and were removed. To potentially increase the recovery rate, the remaining 887 T-box sequences were subdivided on basis of the taxonomy of the species and taxonomic class-specific HMMs were built as before. This procedure yielded only 10 additional T-boxes. Finally, the position of all boxes with respect to the start-codon of the proximal gene was determined. It appeared that In 839 cases, motif 1 was found to be located within 300 nucleotides upstream of a predicted gene start. The remaining 48 T-boxes were located 300 to 480 nt upstream of a start-codon. Manual inspection of those T-boxes revealed an overrepresentation of the proximal genes thrZ (11×), trpE (11×), branched amino acid transferase (5×), chorismate mutase (5×) and a sodium symporter (4×). In fact, for both thrZ and trpE it was previously shown that they are preceeded by multiple (2 or 3) adjacently located T-boxes in different bacilli [3, 29].
Functional classification of the genes regulated by T-boxes
The genes downstream of the recovered T-boxes were divided into four different classes on basis of the gene annotation information of the proximal gene: 'tRNA ligation', 'amino acid biosynthesis, 'amino acid transport' and 'other'. In those cases the T-box preceded an operon (45%) it appeared that most of these operons (>75%) contained genes of only one functional class.
Comparison with the Rfam T-box model
The Rfam database  is an excellent reference database to start the identification of specific RNA-motifs. To evaluate the performance of our recovery procedure, we compared our results with predictions based on the Rfam T-box HMM. Unfortunately, a direct comparison with the contents of the database proved not informative as the Rfam database includes all nucleotide sequences present in the TREMBL database, which includes many genome fragments or partially sequenced genomes, whereas our analysis was restricted to completed genomes.
Therefore, we used the Rfam model to search for T-boxes against the same set of complete genome sequences that was used for our analysis. A striking difference between the Rfam-search and our analysis was that the Rfam-search was computationally far more expensive (more than three weeks on an 8 node,16 core, linux cluster compared to 16 hours on a 2 core linux system). For the selected species we predicted 883 T-boxes characterized by the presence of all the characteristic motifs, of which 835 (95%) were within the first 300 nt upstream of a gene start. When using a cut-off of 53.000 bits (described by Rfam as reliable) only 501 (60%) of the 883 T-boxes were predicted using Rfam. In addition to these 501 shared T-boxes Rfam identified 8 additional genuine T-boxes (~0.9 gain) within the first 300 nt. upstream of a gene start. At a lower cutoff-value (25.000 bits), chosen such that ~81% (683) of the boxes identified by us were recovered, the prediction with our HMM and the Rfam results were more alike albeit that Rfam now also yielded a considerable number of false positive identifications: only ~90% of the total number of T-boxes was found within the first 300 nt. upstream of a gene start. Using this cutoff, 6 additional genuine T-boxes were found (~0.7% gain) to be present within the first 300 nt. upstream of a gene start. When applying no cutoff, the Rfam-search retrieved 89% (745) of the T-boxes that were found in our analysis. However, in this case the number of putative false-positives (not located within 300 nt upstream of a gene start) had increased to 65% of the total.
Henkin TM, Yanofsky C: Regulation by transcription attenuation in bacteria: how RNA provides instructions for transcription termination/antitermination decisions. Bioessays. 2002, 24 (8): 700-707. 10.1002/bies.10125.
Merino E, Yanofsky C: Transcription attenuation: a highly conserved regulatory strategy used by bacteria. Trends Genet. 2005, 21 (5): 260-264. 10.1016/j.tig.2005.03.002.
Gutierrez-Preciado A, Jensen RA, Yanofsky C, Merino E: New insights into regulation of the tryptophan biosynthetic operon in Gram-positive bacteria. Trends Genet. 2005, 21 (8): 432-436. 10.1016/j.tig.2005.06.001.
Schilling O, Langbein I, Muller M, Schmalisch MH, Stulke J: A protein-dependent riboswitch controlling ptsGHI operon expression in Bacillus subtilis: RNA structure rather than sequence provides interaction specificity. Nucleic Acids Res. 2004, 32 (9): 2853-2864. 10.1093/nar/gkh611.
Grundy FJ, Yousef MR, Henkin TM: Monitoring uncharged tRNA during transcription of the Bacillus subtilis glyQS gene. J Mol Biol. 2005, 346 (1): 73-81. 10.1016/j.jmb.2004.11.051.
Rodionov DA, Vitreschak AG, Mironov AA, Gelfand MS: Comparative genomics of thiamin biosynthesis in procaryotes. New genes and regulatory mechanisms. J Biol Chem. 2002, 277 (50): 48949-48959. 10.1074/jbc.M208965200.
Abreu-Goodger C, Ontiveros-Palacios N, Ciria R, Merino E: Conserved regulatory motifs in bacteria: riboswitches and beyond. Trends Genet. 2004, 20 (10): 475-479. 10.1016/j.tig.2004.08.003.
Grundy FJ, Henkin TM: The T box and S box transcription termination control systems. Front Biosci. 2003, 8: d20-31. 10.2741/908.
Grundy FJ, Henkin TM: tRNA as a positive regulator of transcription antitermination in B. subtilis. Cell. 1993, 74 (3): 475-482. 10.1016/0092-8674(93)80049-K.
Putzer H, Laalami S, Brakhage AA, Condon C, Grunberg-Manago M: Aminoacyl-tRNA synthetase gene regulation in Bacillus subtilis: induction, repression and growth-rate regulation. Mol Microbiol. 1995, 16 (4): 709-718. 10.1111/j.1365-2958.1995.tb02432.x.
van de Guchte M, Ehrlich DS, Chopin A: tRNATrp as a key element of antitermination in the Lactococcus lactis trp operon. Mol Microbiol. 1998, 29 (1): 61-74. 10.1046/j.1365-2958.1998.00903.x.
Grundy FJ, Hodil SE, Rollins SM, Henkin TM: Specificity of tRNA-mRNA interactions in Bacillus subtilis tyrS antitermination. J Bacteriol. 1997, 179 (8): 2587-2594.
Grundy FJ, Rollins SM, Henkin TM: Interaction between the acceptor end of tRNA and the T box stimulates antitermination in the Bacillus subtilis tyrS gene: a new role for the discriminator base. J Bacteriol. 1994, 176 (15): 4518-4526.
Marta PT, Ladner RD, Grandoni JA: A CUC triplet confers leucine-dependent regulation of the Bacillus subtilis ilv-leu operon. J Bacteriol. 1996, 178 (7): 2150-2153.
Grundy FJ, Winkler WC, Henkin TM: tRNA-mediated transcription antitermination in vitro: codon-anticodon pairing independent of the ribosome. Proc Natl Acad Sci U S A. 2002, 99 (17): 11121-11126. 10.1073/pnas.162366799.
Luo D, Condon C, Grunberg-Manago M, Putzer H: In vitro and in vivo secondary structure probing of the thrS leader in Bacillus subtilis. Nucleic Acids Res. 1998, 26 (23): 5379-5387. 10.1093/nar/26.23.5379.
Pelchat M, Lapointe J: In vivo and in vitro processing of the Bacillus subtilis transcript coding for glutamyl-tRNA synthetase, serine acetyltransferase, and cysteinyl-tRNA synthetase. Rna. 1999, 5 (2): 281-289. 10.1017/S1355838299980858.
Grundy FJ, Henkin TM: Kinetic analysis of tRNA-directed transcription antitermination of the Bacillus subtilis glyQS gene in vitro. J Bacteriol. 2004, 186 (16): 5392-5399. 10.1128/JB.186.16.5392-5399.2004.
Grundy FJ, Collins JA, Rollins SM, Henkin TM: tRNA determinants for transcription antitermination of the Bacillus subtilis tyrS gene. Rna. 2000, 6 (8): 1131-1141. 10.1017/S1355838200992100.
van de Guchte M, Ehrlich SD, Chopin A: Identity elements in tRNA-mediated transcription antitermination: implication of tRNA D- and T-arms in mRNA recognition. Microbiology. 2001, 147 (Pt 5): 1223-1233.
Yousef MR, Grundy FJ, Henkin TM: Structural transitions induced by the interaction between tRNA(Gly) and the Bacillus subtilis glyQS T box leader RNA. J Mol Biol. 2005, 349 (2): 273-287. 10.1016/j.jmb.2005.03.061.
Fauzi H, Jack KD, Hines JV: In vitro selection to identify determinants in tRNA for Bacillus subtilis tyrS T box antiterminator mRNA binding. Nucleic Acids Res. 2005, 33 (8): 2595-2602. 10.1093/nar/gki546.
Rodionov DA, Vitreschak AG, Mironov AA, Gelfand MS: Comparative genomics of the methionine metabolism in Gram-positive bacteria: a variety of regulatory systems. Nucleic Acids Res. 2004, 32 (11): 3340-3353. 10.1093/nar/gkh659.
Panina EM, Vitreschak AG, Mironov AA, Gelfand MS: Regulation of biosynthesis and transport of aromatic amino acids in low-GC Gram-positive bacteria. FEMS Microbiol Lett. 2003, 222 (2): 211-220.
Henkin TM, Glass BL, Grundy FJ: Analysis of the Bacillus subtilis tyrS gene: conservation of a regulatory sequence in multiple tRNA synthetase genes. J Bacteriol. 1992, 174 (4): 1299-1306.
Grundy FJ, Henkin TM: Conservation of a transcription antitermination mechanism in aminoacyl-tRNA synthetase and amino acid biosynthesis genes in gram-positive bacteria. J Mol Biol. 1994, 235 (2): 798-804. 10.1006/jmbi.1994.1038.
Grundy FJ, Moir TR, Haldeman MT, Henkin TM: Sequence requirements for terminators and antiterminators in the T box transcription antitermination system: disparity between conservation and functional requirements. Nucleic Acids Res. 2002, 30 (7): 1646-1655. 10.1093/nar/30.7.1646.
Delorme C, Ehrlich SD, Renault P: Regulation of expression of the Lactococcus lactis histidine operon. J Bacteriol. 1999, 181 (7): 2026-2037.
Putzer H, Gendron N, Grunberg-Manago M: Co-ordinate expression of the two threonyl-tRNA synthetase genes in Bacillus subtilis: control by transcriptional antitermination involving a conserved regulatory sequence. Embo J. 1992, 11 (8): 3117-3127.
Chopin A, Biaudet V, Ehrlich SD: Analysis of the Bacillus subtilis genome sequence reveals nine new T-box leaders. Mol Microbiol. 1998, 29 (2): 662-664. 10.1046/j.1365-2958.1998.00912.x.
Mwangi MM, Siggia ED: Genome wide identification of regulatory motifs in Bacillus subtilis. BMC Bioinformatics. 2003, 4: 18-10.1186/1471-2105-4-18.
Henkin TM: tRNA-directed transcription antitermination. Mol Microbiol. 1994, 13 (3): 381-387. 10.1111/j.1365-2958.1994.tb00432.x.
Grundy FJ, Haldeman MT, Hornblow GM, Ward JM, Chalker AF, Henkin TM: The Staphylococcus aureus ileS gene, encoding isoleucyl-tRNA synthetase, is a member of the T-box family. J Bacteriol. 1997, 179 (11): 3767-3772.
Sarsero JP, Merino E, Yanofsky C: A Bacillus subtilis operon containing genes of unknown function senses tRNATrp charging and regulates expression of the genes of tryptophan biosynthesis. Proc Natl Acad Sci U S A. 2000, 97 (6): 2656-2661. 10.1073/pnas.050578997.
Vitreschak AG, Rodionov DA, Mironov AA, Gelfand MS: Riboswitches: the oldest mechanism for the regulation of gene expression?. Trends Genet. 2004, 20 (1): 44-50. 10.1016/j.tig.2003.11.008.
Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Eddy SR, Bateman A: Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2005, 33 (Database issue): D121-4. 10.1093/nar/gki081.
T-box data. [http://www.cmbi.ru.nl/T_box_analysis]
Zuker M, Jacobson AB: Using reliability information to annotate RNA secondary structures. Rna. 1998, 4 (6): 669-679. 10.1017/S1355838298980116.
Elf J, Nilsson D, Tenson T, Ehrenberg M: Selective charging of tRNA isoacceptors explains patterns of codon usage. Science. 2003, 300 (5626): 1718-1722. 10.1126/science.1083811.
Francke C, Siezen RJ, Teusink B: Reconstructing the metabolic network of a bacterium from its genome. Trends Microbiol. 2005, 13 (11): 550-558. 10.1016/j.tim.2005.09.001.
Saier MH: A functional-phylogenetic classification system for transmembrane solute transporters. Microbiol Mol Biol Rev. 2000, 64 (2): 354-411. 10.1128/MMBR.64.2.354-411.2000.
Bandell M, Ansanay V, Rachidi N, Dequin S, Lolkema JS: Membrane potential-generating malate (MleP) and citrate (CitP) transporters of lactic acid bacteria are homologous proteins. Substrate specificity of the 2-hydroxycarboxylate transporter family. J Biol Chem. 1997, 272 (29): 18140-18146. 10.1074/jbc.272.29.18140.
Stucky K, Hagting A, Klein JR, Matern H, Henrich B, Konings WN, Plapp R: Cloning and characterization of brnQ, a gene encoding a low-affinity, branched-chain amino acid carrier in Lactobacillus delbruckii subsp. lactis DSM7290. Mol Gen Genet. 1995, 249 (6): 682-690. 10.1007/BF00418038.
Jack DL, Paulsen IT, Saier MH: The amino acid/polyamine/organocation (APC) superfamily of transporters specific for amino acids, polyamines and organocations. Microbiology. 2000, 146 ( Pt 8): 1797-1814.
Reig N, del Rio C, Casagrande F, Ratera M, Gelpi JL, Torrents D, Henderson PJ, Xie H, Baldwin SA, Zorzano A, Fotiadis D, Palacin M: Functional and structural characterization of the first prokaryotic member of the L-amino acid transporter (LAT) family: a model for APC transporters. J Biol Chem. 2007, 282 (18): 13270-13281. 10.1074/jbc.M610695200.
Ren Q, Kang KH, Paulsen IT: TransportDB: a relational database of cellular membrane transport systems. Nucleic Acids Res. 2004, 32 (Database issue): D284-8. 10.1093/nar/gkh016.
Ebbighausen H, Weil B, Kramer R: Transport of branched-chain amino acids in Corynebacterium glutamicum. Arch Microbiol. 1989, 151 (3): 238-244. 10.1007/BF00413136.
Tauch A, Hermann T, Burkovski A, Kramer R, Puhler A, Kalinowski J: Isoleucine uptake in Corynebacterium glutamicum ATCC 13032 is directed by the brnQ gene product. Arch Microbiol. 1998, 169 (4): 303-312. 10.1007/s002030050576.
Hoshino T, Kose-Terai K, Uratani Y: Isolation of the braZ gene encoding the carrier for a novel branched-chain amino acid transport system in Pseudomonas aeruginosa PAO. J Bacteriol. 1991, 173 (6): 1855-1861.
Androutsellis-Theotokis A, Goldberg NR, Ueda K, Beppu T, Beckman ML, Das S, Javitch JA, Rudnick G: Characterization of a functional bacterial homologue of sodium-dependent neurotransmitter transporters. J Biol Chem. 2003, 278 (15): 12703-12709. 10.1074/jbc.M206563200.
Vasconcelos AT, Ferreira HB, Bizarro CV, Bonatto SL, Carvalho MO, Pinto PM, Almeida DF, Almeida LG, Almeida R, Alves-Filho L, Assuncao EN, Azevedo VA, Bogo MR, Brigido MM, Brocchi M, Burity HA, Camargo AA, Camargo SS, Carepo MS, Carraro DM, de Mattos Cascardo JC, Castro LA, Cavalcanti G, Chemale G, Collevatti RG, Cunha CW, Dallagiovanna B, Dambros BP, Dellagostin OA, Falcao C, Fantinatti-Garboggini F, Felipe MS, Fiorentin L, Franco GR, Freitas NS, Frias D, Grangeiro TB, Grisard EC, Guimaraes CT, Hungria M, Jardim SN, Krieger MA, Laurino JP, Lima LF, Lopes MI, Loreto EL, Madeira HM, Manfio GP, Maranhao AQ, Martinkovics CT, Medeiros SR, Moreira MA, Neiva M, Ramalho-Neto CE, Nicolas MF, Oliveira SC, Paixao RF, Pedrosa FO, Pena SD, Pereira M, Pereira-Ferrari L, Piffer I, Pinto LS, Potrich DP, Salim AC, Santos FR, Schmitt R, Schneider MP, Schrank A, Schrank IS, Schuck AF, Seuanez HN, Silva DW, Silva R, Silva SC, Soares CM, Souza KR, Souza RC, Staats CC, Steffens MB, Teixeira SM, Urmenyi TP, Vainstein MH, Zuccherato LW, Simpson AJ, Zaha A: Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae. J Bacteriol. 2005, 187 (16): 5568-5577. 10.1128/JB.187.16.5568-5577.2005.
Seliverstov AV, Putzer H, Gelfand MS, Lyubetsky VA: Comparative analysis of RNA regulatory elements of amino acid metabolism genes in Actinobacteria. BMC Microbiol. 2005, 5: 54-10.1186/1471-2180-5-54.
Ueda K, Yamashita A, Ishikawa J, Shimada M, Watsuji TO, Morimura K, Ikeda H, Hattori M, Beppu T: Genome sequence of Symbiobacterium thermophilum, an uncultivable bacterium that depends on microbial commensalism. Nucleic Acids Res. 2004, 32 (16): 4937-4944. 10.1093/nar/gkh830.
Bolotin A, Quinquis B, Renault P, Sorokin A, Ehrlich SD, Kulakauskas S, Lapidus A, Goltsman E, Mazur M, Pusch GD, Fonstein M, Overbeek R, Kyprides N, Purnelle B, Prozzi D, Ngui K, Masuy D, Hancy F, Burteau S, Boutry M, Delcour J, Goffeau A, Hols P: Complete sequence and comparative genome analysis of the dairy bacterium Streptococcus thermophilus. Nat Biotechnol. 2004, 22 (12): 1554-1558. 10.1038/nbt1034.
Henkin TM: Control of transcription termination in prokaryotes. Annu Rev Genet. 1996, 30: 35-57. 10.1146/annurev.genet.30.1.35.
NCBI repository. [ftp://ftp.ncbi.nlm.nih.gov/genomes/Bacteria/]
Overbeek R, Larsen N, Walunas T, D'Souza M, Pusch G, Selkov E, Liolios K, Joukov V, Kaznadzey D, Anderson I, Bhattacharyya A, Burd H, Gardner W, Hanke P, Kapatral V, Mikhailova N, Vasieva O, Osterman A, Vonstein V, Fonstein M, Ivanova N, Kyrpides N: The ERGO genome analysis and discovery system. Nucleic Acids Res. 2003, 31 (1): 164-171. 10.1093/nar/gkg148.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.
Durbin R ES: Biological Sequence Analysis: Probabilistic Models of Proteins and Nucleic Acids. 1998, Cambridge , Cambridge University Press
Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.
Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.
Bailey TL, Gribskov M: Combining evidence using p-values: application to sequence homology searches. Bioinformatics. 1998, 14 (1): 48-54. 10.1093/bioinformatics/14.1.48.
Schneider TD, Stephens RM: Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 1990, 18 (20): 6097-6100. 10.1093/nar/18.20.6097.
Teusink B, Wiersma A, Molenaar D, Francke C, de Vos WM, Siezen RJ, Smid EJ: Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model. J Biol Chem. 2006, 281 (52): 40041-40048. 10.1074/jbc.M606263200.
Strauch MA, Zalkin H, Aronson AI: Characterization of the glutamyl-tRNA(Gln)-to-glutaminyl-tRNA(Gln) amidotransferase reaction of Bacillus subtilis. J Bacteriol. 1988, 170 (2): 916-920.
Curnow AW, Ibba M, Soll D: tRNA-dependent asparagine formation. Nature. 1996, 382 (6592): 589-590. 10.1038/382589b0.
NCBI taxonomy webpage. [http://www.ncbi.nlm.nih.gov/Taxonomy/]
We would like to thank Victor de Jager for his assistance in applying the Rfam calculations, Jos Boekhorst and Richard Notebaart for insightful discussions. Christof Francke gratefully acknowledges the support of NBIC/the Netherlands Genomics Initiative via the Kluyver Centre for Genomics of Industrial Fermentations and the BioRange program. For our work we used the life-science grid, which is part of the Dutch e-Science Grid (BIG GRID) and which is a collaborative effort of Netherlands National Computing Facilities foundation (NCF), the Netherlands Bioinformatics Centre (NBIC) and the National Institute for Nuclear Physics and High Energy Physics (NIKHEF).
MW conceived, designed and carried out the tRNA-ligase and amino acid biosynthesis analysis, Rfam comparison and drafted and revised the manuscript. TGK carried out the T-box detection and specificity analysis, RJS and MK helped in the design and coordination of the research and drafting and revising the manuscript. CF helped in the design and coordination of the research, carried out the transporter analysis and helped drafting and revising the manuscript. All authors have read and approved the final manuscript
Electronic supplementary material
Additional file 1: T-box location organized per genome. The occurrence of T-boxes in all analyzed genomes is shown in detail. Position, direction, e-value and specifier codon are shown for the T-box, together with position, name and function of the first gene located downstream and the size of the operon located downstream. All T-boxes are color coded according to the function of the genes located downstream: red; tRNA synthesis; green; amino acid biosynthesis; blue; amino acid transport and purple: other/unknown. (XLS 204 KB)
Additional file 2: T-box distribution among different bacterial species. For all species the number of T-boxes is displayed, divided over four different categories (tRNA synthesis, amino acid transport, amino acid biosynthesis and other). Between brackets the number of regulated genes is shown. This number is based on the operon structure of the genes downstream of the T-box. In case multiple strains of a specific species were sequenced, these are only shown when differences between the strains were observed. (XLS 21 KB)
Authors’ original submitted files for images
About this article
Cite this article
Wels, M., Kormelink, T.G., Kleerebezem, M. et al. An in silico analysis of T-box regulated genes and T-box evolution in prokaryotes, with emphasis on prediction of substrate specificity of transporters. BMC Genomics 9, 330 (2008). https://doi.org/10.1186/1471-2164-9-330
- Amino Acid Transport
- Amino Acid Specificity
- Amino Acid Biosynthesis
- Gene Start
- Proximal Gene