αr35 is a family of bacterial small non-coding RNAs with representatives in a reduced group of Alphaproteobacteria from the order Hyphomicrobiales. The first member of this family (Smr35B) was found in a Sinorhizobium meliloti 1021 locus located in the symbiotic plasmid B (pSymB). Further homology and structure conservation analysis have identified full-length SmrB35 homologs in other legume symbionts (i.e. Rhizobium leguminosarum bv.viciae, R. leguminosarum bv. trifolii and R. etli), as well as in the human and plant pathogens Brucella anthropi and Agrobacterium tumefaciens, respectively. αr35 RNA species are 139-142 nt long (Table 1) and share a common secondary structure consisting of two stem loops and a well conserved rho independent terminator (Figure 1, 2, 3). Most of the αr35 transcripts can be catalogued as trans-acting sRNAs expressed from well-defined promoter regions of independent transcription units within intergenic regions of the Alphaproteobacterial genomes (Figure 5).
Smr35B sRNA was firstly described by del Val et al.,[1] as a result of a computational comparative genomic approach in the intergenic regions (IGRs) of the reference S. meliloti 1021 strain. Northern hybridization experiments confirmed that the predicted smr35B locus did express a single transcript of the expected size, which accumulated differentially in free-living and endosymbiotic bacteria. TAP-based 5’-RACE experiments mapped the transcription start site (TSS) of the full-length Smr35B transcript to the 577,730 nt position in the S. meliloti 1021 genome (http://iant.toulouse.inra.fr/bacteria/annotation/cgi/rhime.cgi) whereas the 3’-end was assumed to be located at the 577,868 nt position matching the last residue of the consecutive stretch of Us of a bona fide Rho-independent terminator (Figure 5). Recent deep sequencing-based characterization of the small RNA fraction (50-350 nt) of S. meliloti further confirmed the expression of Smr35B (here referred to as SmelB053), and mapped the 5’- and 3´-ends of the molecule to the positions proposed earlier.[2]
The nucleotide sequence of Smr35B was initially used as query to search against the Rfam database. This homology search rendered no matches to known bacterial sRNA in this database. Smr35B was next BLASTed with default parameters against all the currently available bacterial genomes (1,615 sequences at 20 April 2011; https://www.ncbi.nlm.nih.gov;). The regions exhibiting significant homology to the query sequence (78-89% similarity) were extracted to create a Covariance Model (CM) from a seed alignment using Infernal (version1.0)[3] (Figure 2).
This CM was used in a further search for new members of the αr35 family in the existing bacterial genomic databases.
CM model | Name | GI accession number | begin | end | strand | %GC | length | Organism | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
class=cellinside | αr35 | class=cellinside | Smr35B | class=cellinside | gi|16263748|ref|NC_003078.1| | class=cellinside | 577730 | class=cellinside | 577868 | class=cellinside | + | class=cellinside | 52 | class=cellinside | 139 | class=cellinside | Sinorhizobium meliloti 1021 plasmid pSymB | |
class=cellinside | αr35 | class=cellinside | Atr35C | class=cellinside | gi|159185562|ref|NC_003063.2| | class=cellinside | 132595 | class=cellinside | 132733 | class=cellinside | + | class=cellinside | 48 | class=cellinside | 139 | class=cellinside | Agrobacterium tumefaciens str. C58 chromosome linear | |
class=cellinside | αr35 | class=cellinside | Rlvr35C | class=cellinside | gi|116249766|ref|NC_008380.1| | class=cellinside | 2256716 | class=cellinside | 2256853 | class=cellinside | + | class=cellinside | 55 | class=cellinside | 138 | class=cellinside | Rhizobium leguminosarum bv. viciae 3841 | |
class=cellinside | αr35 | class=cellinside | Rlt1325r35p04 | class=cellinside | gi|241258599|ref|NC_012852.1| | class=cellinside | 114247 | class=cellinside | 114385 | class=cellinside | - | class=cellinside | 56 | class=cellinside | 139 | class=cellinside | Rhizobium leguminosarum bv. trifolii WSM1325 plasmid pR132504 | |
class=cellinside | αr35 | class=cellinside | Rlt1325r35p02 | class=cellinside | gi|241666492|ref|NC_012858.1| | class=cellinside | 466255 | class=cellinside | 466394 | class=cellinside | - | class=cellinside | ? | class=cellinside | 140 | class=cellinside | Rhizobium leguminosarum bv. trifolii WSM1325 plasmid pR132502 | |
class=cellinside | αr35 | class=cellinside | ReCFNr35f | class=cellinside | gi|86360734|ref|NC_007766.1| | class=cellinside | 136368 | class=cellinside | 136508 | class=cellinside | + | class=cellinside | 57 | class=cellinside | 141 | class=cellinside | Rhizobium etli CFN 42 plasmid p42f | |
class=cellinside | αr35 | class=cellinside | Oar35CII | class=cellinside | gi|153010078|ref|NC_009668.1| | class=cellinside | 1587138 | class=cellinside | 1587279 | class=cellinside | - | class=cellinside | 52 | class=cellinside | 142 | class=cellinside | Brucella anthropi ATCC 49188 chromosome 2 |
The results were manually inspected to deduce a consensus secondary structure for the family (Figure 1 and Figure 2). The consensus structure was also independently predicted with the program locARNATE[4] with very similar predictions. The manual inspection of the 84 sequences found with the CM using Infernal allowed finding seven true homolog sequences: two copies in Rhizobium leguminosarum bv. viciae (chromosome and plasmid pRL11), two copies in Rhizobium leguminosarumbv. trifolii WSM1325 (plasmid pR132504 and plasmid pR132502), in Rhizobium etli CFN 42 plasmid p42f and in the chromosomes of Agrobacterium tumefaciens and Brucella anthropi. All these sequences showed significant Infernal E-values (1.38e-33 – 1.05e-11) and bit-scores. In the case of S. meliloti a second copy was identified in the symbiotic plasmid pSymB (574630-574766) with a significant E-value (3.73e-07) but no expression has been detected under any of the tested conditions (unpublished data). The rest of the sequences found with the model showed high E-values between (8.76e-12 and 1.e-3) but very low bit-scores, which usually is a sign of a remote homologue. However, a manual inspection of these cases showed that the rho independent terminator and the second stem were the only conserved regions, failing the first stem. This two stem arregment construction was largely extended in all the Alphaproteobacteria, being specially conserved in Brucella species.
Smr35B expression was first assessed by del Val et al. in S. meliloti 1021 under different biological conditions; i.e. bacterial growth in TY, minimal medium (MM) and luteolin-MM broth and endosymbiotic bacteria (i.e. mature symbiotic alfalfa nodules). Expression of Smr35B in free-living bacteria was found to be growth-dependent, being the gene down-regulated when bacteria entered the stationary phase. Supplementation of MM with luteolin, the plant flavone that specifically induces transcription of the S. meliloti nodulation genes, stimulated the expression of Smr35B by ~4 fold. In contrast, the Smr35B transcript was not detected in mature nodule tissues. Schlüter et al. further described up-regulation of Smr35B upon an osmotic upshift.
All αr35 loci have recognizable σ70-dependent promoters showing a -35/-10 consensus motif CTTAGAC-n17-CTATAT previously shown to be widely conserved among several other genera in the Alphaproteobacteria.[5] To identify binding sites for other known transcription factors we used the fasta sequences provided by RegPredict[6] (http://regpredict.lbl.gov/regpredict/help.html), and used those position weight matrices (PSWM) provided by RegulonDB[7] (http://regulondb.ccg.unam.mx). We built PSWM for each transcription factor from the RegPredict sequences using the Consensus/Patser program, choosing the best final matrix for motif lengths between 14 and 30 if the corresponding length had not been previously specified (see "Consensus matrices" threshold (average E-value < 10E-10) for each matrix was established (see "Thresholded consensus" in http://gps-tools2.its.yale.edu). Moreover, we searched for conserved unknown motifs using MEME[8] (http://meme.sdsc.edu/meme4_6_1/intro.html) and used relaxed regular expressions (i.e. pattern matching) over all Smr35B homologs promoters. Only an inverted repeat structure built around the motif T-N11-A was found 55 nt upstream of the transcription start site of SmrB35 in S. meliloti which is a degenerated motif of the known conserved nod boxes (Figure 4). This characteristic sequence has been proposed as the specific binding site for the LysR-type proteins.[9] All promoter regions of the seed SmrB35 homologs presented the motif as well.
Most of the members of the αr35 family are trans-encoded sRNAs transcribed from independent promoters in the IGRs of the rhizobial megaplasmids. Exceptions are SmrB35 homologs of R. leguminosarum bv. viciae (Rlvr35C),and R. etli CFN 42 plasmid p42f (ReCFNr35f), which are encoded in the opposite strand of annotated genes, partially overlapping ORFs. The predicted protein products of these overlapping ORFs could not be assigned to any functional category on the basis of the amino acid sequence homology.[10] [11] [12] Thus, these αr35 members are putative cis-encoded antisense sRNAs.The genomic regions of the trans-encoded αr35 sRNAs exhibit partial conservation mainly limited to the sRNA-coding sequence and one flanking gene. Most of the flanking genes of the αr35 loci encode transcription factors and proteins related to nitrogen and glutamine metabolism.
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | SM_b20551 | ALIGN=LEFT | R | ALIGN=RIGHT | 576952 | ALIGN=RIGHT | 577398 | ALIGN=LEFT | NP_437070.1 | ALIGN=LEFT | proteolysis | ALIGN=LEFT | Sinorhizobium meliloti 1021 plasmid pSymB (NC_003078) | |
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Smr35B | ALIGN=LEFT | D | ALIGN=RIGHT | 577730 | ALIGN=RIGHT | 577868 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Sinorhizobium meliloti 1021 plasmid pSymB (NC_003078) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | SM_b20552 | ALIGN=LEFT | D | ALIGN=RIGHT | 578150 | ALIGN=RIGHT | 578881 | ALIGN=LEFT | NP_437071.1 | ALIGN=LEFT | nitrogen compound metabolic process | ALIGN=LEFT | Sinorhizobium meliloti 1021 plasmid pSymB (NC_003078) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Oant_4157 | ALIGN=LEFT | D | ALIGN=RIGHT | 1586007 | ALIGN=RIGHT | 1587065 | ALIGN=LEFT | YP_001372686.1 | ALIGN=LEFT | nitrogen compound metabolic process | ALIGN=LEFT | Brucella anthropi ATCC 49188 chromosome 2 (NC_009668) | |
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Oar35CII | ALIGN=LEFT | R | ALIGN=RIGHT | 1587138 | ALIGN=RIGHT | 1587279 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Brucella anthropi ATCC 49188 chromosome 2 (NC_009668) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Oant_4158 | ALIGN=LEFT | R | ALIGN=RIGHT | 1587338 | ALIGN=RIGHT | 1587724 | ALIGN=LEFT | YP_001372687.1 | ALIGN=LEFT | proteolysis | ALIGN=LEFT | Brucella anthropi ATCC 49188 chromosome 2 (NC_009668) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | RHE_PF00127 | ALIGN=LEFT | R | ALIGN=RIGHT | 133963 | ALIGN=RIGHT | 134406 | ALIGN=LEFT | YP_472745.1 | ALIGN=LEFT | hypothetical protein | ALIGN=LEFT | Rhizobium etli CFN 42 plasmid p42f (NC_007766) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | RHE_PF00128 | ALIGN=LEFT | D | ALIGN=RIGHT | 136269 | ALIGN=RIGHT | 136700 | ALIGN=LEFT | YP_472746.1 | ALIGN=LEFT | hypothetical protein | ALIGN=LEFT | Rhizobium etli CFN 42 plasmid p42f (NC_007766) | |
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | ReCFNr35f | ALIGN=LEFT | D | ALIGN=RIGHT | 136368 | ALIGN=RIGHT | 136508 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Rhizobium etli CFN 42 plasmid p42f (NC_007766) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | RHE_PF00129 | ALIGN=LEFT | D | ALIGN=RIGHT | 137962 | ALIGN=RIGHT | 138264 | ALIGN=LEFT | YP_472747.1 | ALIGN=LEFT | membrane protein | ALIGN=LEFT | Rhizobium etli CFN 42 plasmid p42f (NC_007766) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Atu3124 | ALIGN=LEFT | D | ALIGN=RIGHT | 132103 | ALIGN=RIGHT | 132318 | ALIGN=LEFT | NP_357476.1 | ALIGN=LEFT | ALIGN=LEFT | Agrobacterium tumefaciens str. C58 chromosome linear (NC_003063) | ||
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Atr35C | ALIGN=LEFT | D | ALIGN=RIGHT | 132595 | ALIGN=RIGHT | 132733 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Agrobacterium tumefaciens str. C58 chromosome linear (NC_003063) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Atu3126 | ALIGN=LEFT | D | ALIGN=RIGHT | 133057 | ALIGN=RIGHT | 133344 | ALIGN=LEFT | NP_357475.1 | ALIGN=LEFT | nitrogen compound metabolic process | ALIGN=LEFT | Agrobacterium tumefaciens str. C58 chromosome linear (NC_003063) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | RL2133 | ALIGN=LEFT | D | ALIGN=RIGHT | 2256297 | ALIGN=RIGHT | 2256500 | ALIGN=LEFT | YP_767731.1 | ALIGN=LEFT | hypothetical protein | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 (NC_008380) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | RL2134 | ALIGN=LEFT | R | ALIGN=RIGHT | 2256617 | ALIGN=RIGHT | 2256982 | ALIGN=LEFT | YP_767732.1 | ALIGN=LEFT | hyphotetical protein | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 (NC_008380) | |
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Rlvr35C | ALIGN=LEFT | D | ALIGN=RIGHT | 2256716 | ALIGN=RIGHT | 2256853 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 (NC_008380) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | RL2135 | ALIGN=LEFT | D | ALIGN=RIGHT | 2256994 | ALIGN=RIGHT | 2257383 | ALIGN=LEFT | YP_767733.1 | ALIGN=LEFT | transposase-related protein | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 (NC_008380) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Rleg_6079 | ALIGN=LEFT | D | ALIGN=RIGHT | 113829 | ALIGN=RIGHT | 114197 | ALIGN=LEFT | YP_002978585.1 | ALIGN=LEFT | membrane proteiin | ALIGN=LEFT | Rhizobium leguminosarum trifolii WSM1325 plasmid pR132502 (NC_012852) | |
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Rlt132504r35p04 | ALIGN=LEFT | R | ALIGN=RIGHT | 114247 | ALIGN=RIGHT | 114385 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Rhizobium leguminosarum trifolii WSM1325 plasmid pR132502 (NC_012852) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Rleg_6080 | ALIGN=LEFT | R | ALIGN=RIGHT | 114489 | ALIGN=RIGHT | 115121 | ALIGN=LEFT | YP_002978586.1 | ALIGN=LEFT | endonuclease | ALIGN=LEFT | Rhizobium leguminosarum trifolii WSM1325 plasmid pR132502 (NC_012852) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Rleg_7049 | ALIGN=LEFT | D | ALIGN=RIGHT | 465959 | ALIGN=RIGHT | 466222 | ALIGN=LEFT | YP_002985022.1 | ALIGN=LEFT | ALIGN=LEFT | Rhizobium leguminosarum trifolii WSM1325 plasmid pR132504 (NC_012858) | ||
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Rlt132502r35p02 | ALIGN=LEFT | R | ALIGN=RIGHT | 466255 | ALIGN=RIGHT | 466394 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Rhizobium leguminosarum trifolii WSM1325 plasmid pR132504 (NC_012858) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | Rleg_7050 | ALIGN=LEFT | R | ALIGN=RIGHT | 466934 | ALIGN=RIGHT | 467824 | ALIGN=LEFT | YP_002985023.1 | ALIGN=LEFT | transcription regulator | ALIGN=LEFT | Rhizobium leguminosarum trifolii WSM1325 plasmid pR132504 (NC_012858) | |
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | pRL110105 | ALIGN=LEFT | D | ALIGN=RIGHT | 122566 | ALIGN=RIGHT | 123456 | ALIGN=LEFT | YP_771137.1 | ALIGN=LEFT | transcription regulator | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 plasmid pRL11 (NC_008384) | |
αr35 | ALIGN=LEFT | sRNA | ALIGN=LEFT | Rlvr35p11 | ALIGN=LEFT | D | ALIGN=RIGHT | 124030 | ALIGN=RIGHT | 124162 | ALIGN=LEFT | ALIGN=LEFT | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 plasmid pRL11 (NC_008384) | |||
αr35 | ALIGN=LEFT | gene | ALIGN=LEFT | pRL110106 | ALIGN=LEFT | R | ALIGN=RIGHT | 124229 | ALIGN=RIGHT | 124447 | ALIGN=LEFT | YP_771138.1 | ALIGN=LEFT | hyphotetical protein | ALIGN=LEFT | Rhizobium leguminosarum bv. viciae 3841 plasmid pRL11 (NC_008384) |