16S rRNA序列分析法在医学微生物鉴定中的应用

１６Ｓ ｒＲＮＡ序列分析作为微生物系统分类的主要依据已得到了广泛认同,随着微生物核糖体数据库的日益完善,该技术成为细菌分类和鉴定的一个有力工具。本文概述了 １６５ ｒＲＮＡ序列分析法的技术步骤以及该技术在医学微生物研究中的应用,总结了目前文献报导的各种致病微生物种属特异性 １６５ ｒＲＮＡ引物和探针序列,同时分析了该技术在应用中存在的一些问题。     

Questionable 16S ribosomal RNA gene annotations are frequent in completed microbial genomes

According to recent reports, many ribosomal RNA gene annotations are still questionable, and the use of inappropriate tools for annotation has been blamed. However, we believe that the abundant 16S rRNA partial sequence in the databases, mainly created by culture-independent PCR methods, is another main cause of the ambiguous annotations of 16S rRNA. To examine the current status of 16S rRNA gene annotations in complete microbial genomes, we used as a criterion the conserved anti-SD sequence, located at the 3′ end of the 16S rRNA gene, which is commonly overlooked by culture-independent PCR methods. In our large survey, 859 16S rRNA gene sequences from 252 different species of the microbial complete genomes were inspected. 67 species (234 genes) were detected with ambiguous annotations. The common anti-SD sequence and other conserved 16S rRNA sequence features could be detected in the downstream-intergenic regions for almost every questionable sequence, indicating that many of the 16S rRNA genes were annotated incorrectly. Furthermore, we found that more than 91.5% of the 93,716 sequences of the available 16S rRNA in the main databases are partial sequences. We also performed BLAST analysis for every questionable rRNA sequence, and most of the best hits in the analysis were rRNA partial sequences. This result indicates that partial sequences are prevalent in the databases, and that these sequences have significantly affected the accuracy of microbial genomic annotation. We suggest that the annotation of 16S rRNA genes in newly complete microbial genomes must be done in more detail, and that revision of questionable rRNA annotations should commence as soon as possible.     

小鲵科线粒体16S rRNA基因序列分析及其系统发育

To study the phylogeny of Hynobiidae, we amplified DNA fragments of 470 bp 16S ribosomal RNA (16S rRNA) gene on mitochondrial DNA from Ranodon sibiricus and Ranodon tsinpaensis. PCR products were cloned into PMD18 T vector after purification. These sequences were determined and deposited in the GenBank (accession numbers: AY373459 for Ranodon sibiricus, AY372534 for Ranodon tsinpaensis). By comparing the nucleotide differences of 16S ribosomal RNA sequences among Liua shihi, Pseudohynobius flavomaculatus and Batrachuperus genus from GenBank database, we analyzed the divergences and base substitution among these sequences with the MEGA software. The molecular results support that B. tibetanus, B. pinchonii and B. karlschmidti are classified into three valid species. Liua shihi has closer phylogenetic relationships to Ranodon tsinpaensis than to other species. More our results reveal that Pseudohynobius flavomaculatus is not a synonym of Ranodon tsinpaensis. [Acta Zoologica Sinica 50 (3) : 464 - 469,2004].     

A novel mutation 3090 G>A of the mitochondrial 16S ribosomal RNA associated with myopathy

We describe a young woman who presented with a progressive myopathy since the age of 9. Spectrophotometric analysis of the respiratory chain in muscle tissue revealed combined and profound complex I, III, II+III, and IV deficiency ranging from 60% to 95% associated with morphological and histochemical abnormalities of the muscle. An exhaustive screening of mitochondrial transfer and ribosomal RNAs showed a novel G>A substitution at nucleotide position 3090 which was detected only in urine sediment and muscle of the patient and was not found in her mother's blood cells and urine sample. We suggest that this novel de novo mutation in the 16S ribosomal RNA, a nucleotide which is highly conserved in different species, would impair mitochondrial protein synthesis and would cause a severe myopathy.     

基于18S rRNA和线粒体16S rRNA基因序列的柳蚕进化分析

为探讨柳蚕Actias selene Hübner与鳞翅目昆虫的系统发育关系,本研究利用PCR扩增获得了柳蚕核糖体18S rRNA和线粒体16S rRNA基因的部分序列,长度分别为391bp和428bp。并采用邻近距离法(NJ)、最大简约法(MP)、类平均聚类法(UPGMA)构建系统进化树。结果表明,柳蚕线粒体16SrRNA基因序列与大蚕蛾科昆虫的16SrRNA基因序列均表现出偏好于碱基AT的倾向。柳蚕与所研究的其它蚕类的遗传距离介于0.016至0.140之间,其中与温带柞蚕Antheraea roylii的遗传距离最小,与野桑蚕Bombyx mandarina的遗传距离最大。而基于鳞翅目昆虫18S rRNA基因部分序列的进化分析显示,柳蚕与柞蚕Antheraea pernyi之间的遗传距离最小(0.010),与蓖麻蚕Samia ricini的遗传距离最大(0.017)。     

Phylogenic homogeneity of Coxiella burnetii strains as determinated by 16S ribosomal RNA sequencing

Abstract DNA coding for the 16S rRNA of six strains of the obligate intracellular bacterium Coxiella burnetti was directly amplified from lysed host cells using the polymerase chain reaction. The amplification product was sequenced using a linear-PCR procedure and compared with other published 16S rRNA sequences. The results of this analysis confirm the position of C. burnetii in the gamma subgroup of the proteobacteria. The data show that all of the C. burnetii strains are highly related (> 99%) on the basis of 16S rRNA sequences although they had different geographic origins and phenotypic characteristics. The data support a phylogenetic homogeneity of the genus Coxiella with only one species which is C. burnetii .     

Role of N-terminal helix in interaction of ribosomal protein S15 with 16S rRNA

The position and conformation of the N-terminal helix of free ribosomal protein S15 was earlier found to be modified under various conditions. This variability was supposed to provide the recognition by the protein of its specific site on 16S rRNA. To test this hypothesis, we substituted some amino acid residues in this helix and assessed effects of these substitutions on the affinity of the protein for 16S rRNA. The crystal structure of the complex of one of these mutants (Thr3Cys S15) with the 16S rRNA fragment was determined, and a computer model of the complex containing another mutant (Gln8Met S15) was designed. The available and new information was analyzed in detail, and the N-terminal helix was concluded to play no significant role in the specific binding of the S15 protein to its target on 16S rRNA.     

Assembly of the 5' and 3' minor domains of 16S ribosomal RNA as monitored by tethered probing from ribosomal protein S20

The ribosomal protein (r-protein) S20 is a primary binding protein. As such, it interacts directly and independently with the 5′ domain as well as the 3′ minor domain of 16S ribosomal RNA (rRNA) in minimal particles and the fully assembled 30S subunit. The interactions observed between r-protein S20 and the 5′ domain of 16S rRNA are quite extensive, while those between r-protein S20 and the 3′ minor domain are significantly more limited. In this study, directed hydroxyl radical probing mediated by Fe(II)-derivatized S20 proteins was used to monitor the folding of 16S rRNA during r-protein association and 30S subunit assembly. An analysis of the cleavage patterns in the minimal complexes [16S rRNA and Fe(II)-S20] and the fully assembled 30S subunit containing the same Fe(II)-derivatized proteins shows intriguing similarities and differences. These results suggest that the two domains, 5′ and 3′ minor, are organized relative to S20 at different stages of assembly. The 5′ domain acquires, in a less complex ribonucleoprotein particle than the 3′ minor domain, the same architecture as observed in mature subunits. These results are similar to what would be predicted of subunit assembly by the 5′-to-3′ direction assembly model.     

Selection of primers for optimal taxonomic classification of environmental 16S rRNA gene sequences

Microbial community profiling using 16S rRNA gene sequences requires accurate taxonomy assignments. ‘Universal'' primers target conserved sequences and amplify sequences from many taxa, but they provide variable coverage of different environments, and regions of the rRNA gene differ in taxonomic informativeness—especially when high-throughput short-read sequencing technologies (for example, 454 and Illumina) are used. We introduce a new evaluation procedure that provides an improved measure of expected taxonomic precision when classifying environmental sequence reads from a given primer. Applying this measure to thousands of combinations of primers and read lengths, simulating single-ended and paired-end sequencing, reveals that these choices greatly affect taxonomic informativeness. The most informative sequence region may differ by environment, partly due to variable coverage of different environments in reference databases. Using our Rtax method of classifying paired-end reads, we found that paired-end sequencing provides substantial benefit in some environments including human gut, but not in others. Optimal primer choice for short reads totaling 96 nt provides 82–100% of the confident genus classifications available from longer reads.     

Assembly of the central domain of the 30S ribosomal subunit: roles for the primary binding ribosomal proteins S15 and S8

Assembly of the 30S ribosomal subunit occurs in a highly ordered and sequential manner. The ordered addition of ribosomal proteins to the growing ribonucleoprotein particle is initiated by the association of primary binding proteins. These proteins bind specifically and independently to 16S ribosomal RNA (rRNA). Two primary binding proteins, S8 and S15, interact exclusively with the central domain of 16S rRNA. Binding of S15 to the central domain results in a conformational change in the RNA and is followed by the ordered assembly of the S6/S18 dimer, S11 and finally S21 to form the platform of the 30S subunit. In contrast, S8 is not part of this major platform assembly branch. Of the remaining central domain binding proteins, only S21 association is slightly dependent on S8. Thus, although S8 is a primary binding protein that extensively contacts the central domain, its role in assembly of this domain remains unclear. Here, we used directed hydroxyl radical probing from four unique positions on S15 to assess organization of the central domain of 16S rRNA as a consequence of S8 association. Hydroxyl radical probing of Fe(II)-S15/16S rRNA and Fe(II)-S15/S8/16S rRNA ribonucleoprotein particles reveal changes in the 16S rRNA environment of S15 upon addition of S8. These changes occur predominantly in helices 24 and 26 near previously identified S8 binding sites. These S8-dependent conformational changes are consistent with 16S rRNA folding in complete 30S subunits. Thus, while S8 binding is not absolutely required for assembly of the platform, it appears to affect significantly the 16S rRNA environment of S15 by influencing central domain organization.     

Rapid identification of marine bioluminescent bacteria by amplified 16S ribosomal RNA gene restriction analysis

To rapidly identify natural isolates of marine bioluminescent bacteria, we developed amplified ribosomal DNA restriction analysis (ARDRA) methods. ARDRA, which is based on the restriction patterns of 16S rRNA gene digested with five enzymes (EcoRI, DdeI, HhaI, HinfI, RsaI), clearly distinguished the 14 species of marine bioluminescent bacteria currently known, which belong to the genera Vibrio, Photobacterium, and Shewanella. When we applied ARDRA to 129 natural isolates from two cruises in Sagami Bay, Japan, 127 were grouped into six ARDRA types with distinctive restriction patterns; these isolates represented the bioluminescent species, P. angustum, P. leiognathi, P. phosphoreum, S. woodyi, V. fischeri, and V. harveyi. The other two isolates showing unexpected ARDRA patterns turned out to have 16S rRNA gene sequences similar to P. leiognathi and P. phosphoreum. Nevertheless, ARDRA provides a simple and fairly robust means for rapid identification of the natural isolates of marine bioluminescent bacteria, and is therefore useful in studying their diversity.     

A comparison of the 16S ribosomal RNAs from mesophilic and thermophilic bacilli: Some modifications in the sanger method for RNA sequencing

Summary Two modifications in the Sanger two dimensional electrophoretic procedure for RNA analysis are reported. One increases resolution on the primary fingerprint to the point that digests of large RNAs, of the size 1500–3000 nucleotides yield well resolved fingerprint patterns. The other is a novel endonucleolytic procedure that proves useful in determining sequences of the large oligonucleotides produced by T1 ribonuclease.These modifications have been used in determining the catalogs of oligomers produced by T1 ribonuclease digestion of 16S rRNAs from three related organisms,Bacillus subtilis, B.pumilus andB.stearothermophilus. The possible effects of adaptation to a thermophilic niche on ribosomal RNA primary structure and the phylogenetic relatedness of the two mesophilic Bacilli are discussed.This is contribution No.6 in a series on procaryote phylogeny.     

A Comparison Between Elastic Network Interpolation and MD Simulation of 16S Ribosomal RNA

Abstract

In this paper a coarse-grained method called elastic network interpolation (ENI) is used to generate feasible transition pathways between two given conformations of the core central domain of 16S Ribosomal RNA (16S rRNA). The two given conformations are the extremes generated by a molecular dynamics (MD) simulation, which differ from each other by 10Å in root-mean-square deviation (RMSD). It takes only several hours to build an ENI pathway on a 1.5GHz Pentium with 512 MB memory, while the MD takes several weeks on high-performance multi-processor servers such as the SGI ORIGIN 2000/2100. It is shown that multiple ENI pathways capture the essential anharmonic motions of millions of timesteps in a particular MD simulation. A coarse-grained normal mode analysis (NMA) is performed on each intermediate ENI conformation, and the lowest 1% of the normal modes (representing about 40 degrees of freedom (DOF)) are used to parameterize fluctuations. This combined ENI/NMA method captures all intermediate conformations in the MD run with 1.5Å RMSD on average. In addition, if we restrict attention to the time interval of the MD run between the two extreme conformations, the RMSD between the closest ENI/NMA pathway and the MD results is about 1Å. These results may serve as a paradigm for reducedDOF dynamic simulations of large biological macromolecules as well as a method for the reduced-parameter interpretation of massive amounts of MD data.     

Identification of Escherichia coli through analysis of 16S rRNA and 16S-23S rRNA internal transcribed spacer region sequences

A bacterial strain, designated BzDS03 was isolated from water sample, collected from Dal Lake Srinagar. The strain was characterized by using 16S ribosomal RNA gene and 16S-23S rRNA internal transcribed spacer region sequences. Phylogenetic analysis showed that 16S rRNA sequence of the isolate formed a monophyletic clade with genera Escherichia. The closest phylogenetic relative was Escherichia coli with 99% 16S rRNA gene sequence similarity. The result of Ribosomal database project's classifier tool revealed that the strain BzDS03 belongs to genera Escherichia.16S rRNA sequence of isolate was deposited in GenBank with accession number FJ961336. Further analysis of 16S-23S rRNA sequence of isolate confirms that the identified strain BzDS03 be assigned as the type strain of Escherichia coli with 98% 16S-23S rRNA sequence similarity. The GenBank accession number allotted for 16S-23S rRNA intergenic spacer sequence of isolate is FJ961337.     

A phylogenetic analysis of the genus Catellatospora based on 16S ribosomal DNA sequences, including transfer of Catellatospora matsumotoense to the genus Micromonospora as Micromonospora matsumotoense comb. nov.

Phylogenetic studies based on the 16S ribosomal gene sequences showed that members of the genus Catellatospora revealed phylogenetic heterogeneity within the family Micromonosporaceae as well as a heterogeneous menaquinone composition. Among them, Catellatospora matsumotoense was closely related to members of the genus Micromonospora, indicating that this organism should be excluded from the genus Catellatospora. On the basis of classical taxonomic characteristics and phylogenetic evidence, Catellatospora matsumotoense is proposed to be transferred to the genus Micromonospora as M. matsumotoense comb. nov.     

Functional studies of the 900 tetraloop capping helix 27 of 16S ribosomal RNA

The 900 tetraloop (positions 898-901) of Escherichia coli 16S rRNA caps helix 27, which is involved in a conformational switch crucial for the decoding function of the ribosome. This tetraloop forms a GNRA motif involved in intramolecular RNA-RNA interactions with its receptor in helix 24 of 16S rRNA. It is involved also in an intersubunit bridge, via an interaction with helix 67 in domain IV of 23S rRNA. Using a specialized ribosome system and an instant-evolution procedure, the four nucleotides of this loop were randomized and 15 functional mutants were selected in vivo. Positions 899 and 900, responsible for most of the tetraloop/receptor interactions, were found to be the most critical for ribosome activity. Functional studies showed that mutations in the 900 tetraloop impair subunit association and decrease translational fidelity. Computer modeling of the mutations allows correlation of the effect of mutations with perturbations of the tetraloop/receptor interactions.     

基于16S rDNA部分序列探讨中国近海30种石斑鱼类的分子系统进化关系

石斑鱼因其种类繁多、分布广泛及缺乏显著的形体特征,使其系统分类的研究颇为困难。为探讨中国近海石斑鱼类的系统进化关系,通过PCR扩增获得了石斑鱼亚科(Epinephelinae)6属30个种类的线粒体16SrDNA基因片段序列。采用多个生物软件对序列变异和碱基组成进行分析,计算了Kimura2parameter遗传距离、转换/颠换比等遗传信息指数,并结合GenBank石斑鱼属的同源序列,以多纹长尾(Pronotogrammusmultifasciatus)和皮氏叫姑鱼(Johniusbelengerii)为外群构建NJ、MP和ML系统树。根据所得分子依据并结合形态学特征,推论如下:1)在本研究的30种石斑鱼中,鳃棘鲈属(Plectropomus)最先分化,并呈明显单系性;九棘鲈属(Cephalopholis)是一个单系群,并且较石斑鱼属(Epinephelus)原始;侧牙鲈属(Variola)的进化地位介于鳃棘鲈属与九棘鲈属之间;2)宽额鲈(Promicropslanceolatus)可以归入石斑鱼属,而驼背鲈(Cromileptesaltivelis)也与石斑鱼属有很近的亲缘关系,甚至可能是石斑鱼属内的特化类群;3)石斑鱼属内部存在两个平行进化的姐妹分支,分支内部的种类组成与地理分布无关,暗示了石斑鱼属早期的分化模式。     

Ribosomal protein-dependent orientation of the 16 S rRNA environment of S15

Ribosomal protein S15 binds specifically to the central domain of 16 S ribosomal RNA (16 S rRNA) and directs the assembly of four additional proteins to this domain. The central domain of 16 S rRNA along with these five proteins form the platform of the 30 S subunit. Previously, directed hydroxyl radical probing from Fe(II)-S15 in small ribonucleoprotein complexes was used to study assembly of the central domain of 16 S rRNA. Here, this same approach was used to understand the 16 S rRNA environment of Fe(II)-S15 in 30 S subunits and to determine the ribosomal proteins that are involved in forming the mature S15-16 S rRNA environment. We have identified additional sites of Fe(II)-S15-directed cleavage in 30S subunits compared to the binary complex of Fe(II)-S15/16 S rRNA. Along with novel targets in the central domain, sites within the 5' and 3' minor domains are also cleaved. This suggests that during the course of 30S subunit assembly these elements are positioned in the vicinity of S15. Besides the previously determined role for S8, roles for S5, S6+S18, and S16 in altering the 16 S rRNA environment of S15 were established. These studies reveal that ribosomal proteins can alter the assembly of regions of the 30 S subunit from a considerable distance and influence the overall conformation of this ribonucleoprotein particle.