Users Online: 2621
Home Print this page Email this page
Home About us Editorial board Search Browse articles Submit article Ahead of Print Instructions Subscribe Contacts Special issues Login 

Previous article Browse articles Next article 
Adv Biomed Res 2015,  4:152

Molecular typing of Iranian mycobacteria isolates by polymerase chain reaction-restriction fragment length polymorphism analysis of 360-bp rpoB gene

1 Department of Microbiology, School of Medicine, Isfahan University of Medical Sciences, Isfahan, Iran
2 Department of Health, Isfahan Provincial Health Office, Isfahan University of Medical Sciences, Isfahan, Iran

Date of Submission08-Feb-2014
Date of Acceptance15-Feb-2015
Date of Web Publication27-Jul-2015

Correspondence Address:
Dr. Bahram Nasr Esfahani
Department of Microbiology, School of Medicine, Isfahan University of Medical Sciences, Isfahan
Login to access the Email id

Source of Support: None, Conflict of Interest: None

DOI: 10.4103/2277-9175.161579

Rights and Permissions

Background: Diagnosis and typing of Mycobacterium genus provides basic tools for investigating the epidemiology and pathogenesis of this group of bacteria. Polymerase chain reaction (PCR)-restriction fragment length polymorphism analysis (PRA) is an accurate method providing diagnosis and typing of species of mycobacteria. The present study is conducted by the purpose of determining restriction fragment profiles of common types of mycobacteria by PRA method of rpoB gene in this geographical region.
Materials and Methods: Totally 60 clinical and environmental isolates from February to October, 2013 were collected and subcultured and identified by phenotypic methods. A 360 bp fragment of the rpoB gene amplified by PCR and products were digested by MspI and HaeIII enzymes.
Results: In the present study, of all mycobacteria isolates identified by PRA method, 13 isolates (21.66%) were Mycobacterium tuberculosis, 34 isolates (56.66%) were rapidly growing Nontuberculosis Mycobacteria (NTM) that including 26 clinical isolates (43.33%) and 8 environmental isolates (13.33%), 11 isolates (18.33%) were clinical slowly growing NTM. among the clinical NTM isolates, Mycobacterium fortuitum Type I with the frequency of 57.77% was the most prevalent type isolates. Furthermore, an unrecorded of the PRA pattern of Mycobacterium conceptionense (HeaIII: 120/90/80, MspI: 120/105/80) was found. This study demonstrated that the PRA method was high discriminatory power for identification and typing of mycobacteria species and was able to identify 96.6% of all isolates.
Conclusion: Based on the result of this study, rpoB gene could be a potentially useful tool for identification and investigation of molecular epidemiology of mycobacterial species.

Keywords: Molecular typing, mycobacteria, polymerase chain reaction-restriction fragment length polymorphism analysis, rpoB gene

How to cite this article:
Hadifar S, Moghim S, Fazeli H, GhasemianSafaei H, Havaei SA, Farid F, Esfahani BN. Molecular typing of Iranian mycobacteria isolates by polymerase chain reaction-restriction fragment length polymorphism analysis of 360-bp rpoB gene. Adv Biomed Res 2015;4:152

How to cite this URL:
Hadifar S, Moghim S, Fazeli H, GhasemianSafaei H, Havaei SA, Farid F, Esfahani BN. Molecular typing of Iranian mycobacteria isolates by polymerase chain reaction-restriction fragment length polymorphism analysis of 360-bp rpoB gene. Adv Biomed Res [serial online] 2015 [cited 2023 Sep 26];4:152. Available from:

  Introduction Top

Mycobacterium genus includes more than 160 distinct species that are mostly involved in human infections. [1] Mycobacteria are responsible for considerable human morbidity and mortality worldwide. Some species, such as Mycobacterium tuberculosis are the primary cause of mortality in individual adult's infections. [2] This successful human bacterial parasite has infected one-third of the world population and each year causes more than 8.8 million new cases of infection and kills closely 2 million people. [3],[4] The World Health Organization also estimates that between 2000 and 2020, 35 million people will die from tuberculosis (TB). [2],[5]

Nontuberculous Mycobacteria (NTM) were only known as environmental bacteria. However, they have recently been identified as important pathogens due to the increased prevalence of immunodeficiency diseases. NTM, have more than 100 species but only around 15 being cause of a wide variety of human infections. [6],[7] They are one of the cause for nosocomial and occupational infections and induce pneumonitis, hypersensitivity, asthma and bronchitis in workers exposed to liquid water, metal and aerosols. [8] This group of mycobacteria also causes opportunistic infections, particularly in patients who are immune compromised. [9] With regarding the widespread prevalence of AIDS, infection with NTM as well as TB has increased in many parts of the world more than last decades. [10] Hence, clinical diagnosis and treatment of NTM infections have become a challenge for clinicians. [11]

Diagnosis and typing of Mycobacterium genus provides a basis for investigating their epidemiology and pathogenesis. Regarding to the prevalence of mycobacteria infection in Iran and because of the neighborhood of Iran with countries that are among 22 high-burden countries and the countries with high prevalence of multidrug-resistant TB, increasing attention to mycobacteria diseases and introducing molecular epidemiology of mycobacteria seems would be necessary to deal with this challenge. In other hands, because there are different treatments strategies for various related disease of mycobacteria species; notice the frequency of species in each region would be useful in different geographic regions in adopting a promising method of control and treatment by physicians. [12]

Mycobacteria are routinely identified to the species level by the phenotypic approaches, but these phenotypic-based methods have some limitations in the diagnosis including being time-consuming, requirement of living organism, confusing and sometimes incorrect results. [12],[13] Today, molecular approaches like sequencing, INNO-LiPA mycobacteria v2, and AccuProbe are used for precise and rapid identification of the species, however, high cost of these methods and also the requirement for special facilities impose limitations on their application in developing countries. [14],[15],[16] Common typing methods like bacteriophage typing and serotyping are replaced by ribotyping, polymerase chain reaction (PCR)-based methods, fingerprinting plasmid, and analysis of restriction fragments of chromosomal DNA by pulse-field gel electrophoresis in the last decades. [17] Among molecular approaches for typing, analysis PCR-restriction fragment length polymorphism (RFLP) is an accurate and inexpensive method providing diagnosis and typing of species and subspecies of mycobacteria. [14],[18] In order to typing mycobacteria by PCR-RFLP analysis (PRA) method, different genes such as 16SrRNA, rpoB, dnaj, and hsp65 were targeted. [9],[14] The variable regions of rpoB gene that is encoding the β subunit of RNA polymerase enzyme are proposed as proper gene for phylogenic analysis, determination of inter-species diversity and diagnosis of bacteria specially ones with close relativity. [18],[19] In this study, the PRA profiles and inter-species diversities of Iranian clinical and environmental mycobacteria isolates was determined and their profiles was compared with the global patterns and relevant information about the epidemiology and genetic basis of Iran isolates was provided.

  Materials and methods Top

Mycobacteria strains and species

Totally 60 mycobacteria isolates were obtained from February to October, 2013 at the Mycobacterium collection of Department of Microbiology, Isfahan University of Medical Sciences and Tuberculosis Center of Isfahan. Of these 60 specimens, 40 were respiratory specimens, including sputum, bronchoalveolar lavage, and bronchial wash specimens, 11 were nonrespiratory and 9 specimens were environmental sample such as water due to break up the mucin, each respiratory specimen was treated with a same volume of 3.5% NaOH in a 50-ml centrifuge tube and, followed by vortexing for 30 s. specimens decontaminated after the tubes were incubated at room temperature for 15 min. Sterile phosphate buffer was added to stop the digestion decontamination process, and then the tube were mixed by inversion, and the tubes were centrifuged at 3000 × g for 15 min. The supernatant was discarded, and the remaining pellet was resuspended in 3.0 ml of sterile phosphate-buffered saline. Nonrepirotry after the homogenies processed in the same way sampling of water was done using the grab sampling method. To 2 L' sterile Erlenmeyer flasks, sodium thio-sulfate as antichlor and 0.04% Cetyl pridinium chloride as an antimicrobial agent were added. Five hundred milliliters of samples were passed from 0.45 μm filters. The filters were transferred directly onto 7H10 Middlebrook solid media, include 15% oleic acid, albumin, dextrose, catalase.

Primary identification of isolates by conventional methods

The mycobacteria species were subcultured on LJ media and Middlebrook 7H9. The isolates were identified by primary conventional methods including acid-fast staining, colony characteristics, growth at 25, 37 and 42°C, pigment production, semi-quantitative catalase test, Tween 80 hydrolysis, arylsulfatase test (3 and 14 days), heat-stable catalase (pH 7, 68°C), pyrazinamidase (4 and 7 days), urease, nitrate reduction test and colony morphology. The reference strains used in this study were M. tuberculosis H37Rv (ATCC 27294) and Mycobacterium fortuitum (ATCC 49403).

Polymerase chain reaction-restriction fragment length polymorphism analysis

Preparation of genomic DNA

Fresh cultures of bacteria were prepared on 7H10 medium and then DNA was extracted using cetyltrimethylammonium bromide method. [20] Purified DNA was stored at − 70°C for subsequent experiments.

Polymerase chain reaction amplification of rpoB

The reaction was performed in a thermal cycler from PCR Hybid (Omnigene). Partial rpoB gene (360 bp) was amplified using primers (RPO5') 5-TCAAGGAGAAGCGCTACGA-3' and (RPO3') 5'- GGATGTTGATCAGGGTCTGC-3'. [11] Each 25 μl PCR mixture contained 2 μl DNA supernatant (5 ng genomic DNA(, 1 μl of each primer (10 pmol/ml), 1.25 μl MgCl 2 (1.5 mM), 0.5 μl dNTP (200 mM), 0.25 μl Taq polymerase (500U), 2.5 μl × 10 Buffer. The PCR program was performed in a thermocycler (Eppendorf) including: 94°C for 5 min, 35 cycle of 94°C for 1 min, 60°C for 1 min, and 72°C for 1 min, and a final step at 72°C for 7 min. The PCR products were run on 1.5% agarose gel.

Restriction fragment length polymorphism

The amplified products of rpoB gene regions were digested with two restriction enzymes of MspI and Hea III according to recommendations of the manufacturers. Fourteen μl of PCR products were digested in a 20 μl reaction volume containing 2 μl of × 10 reaction buffer supplied by manufacture, 3.5 μl dH 2 O and 5U MspI enzyme and were incubated at 37°C for 4 h. A similar condition was used for the Hea III enzyme. Finally, samples were run on 4% metaphor agarose gel, and fragment sizes were determined based on the restriction patterns first described by Lee et al. [11] For undetermined isolates based on the reference pattern or methods also sequence analysis was performed.


Of 60 isolates, out of 13 isolates (21.66%) identified as M. tuberculosis and 47 isolates (78.33%) as NTM by phenotypic method. Among of 47 NTM isolate, 38 isolates (80.85%) were clinical and 9 isolates (19.14%) were environmental [Table 1].
Table 1: List of mycobacteria strains and their sources identified by phenotypic tests and fragment sizes of mycobacterial 360 bp rpoB PCR products after digestion by HaeIII and MspI

Click here to view

Of 60 mycobacteria isolates identified by PRA method, 13 isolates (21.66%) were M. tuberculosis, 34 isolates (56.66%) were rapidly growing NTM that including 26 clinical isolates (43.33%) and 8enviromental isolates (13.33%), 11 isolates (18.33%) were slowly growing NTM which were clinical isolates.

Based on PRA of rpoB gene, among 34 rapidly growing isolates, 32 isolates were M. fortuitum Type I; 1 isolate Mycobacterium smegmatis and 1isolate Mycobacterium conceptionense. An unrecorded pattern profile of M. conceptionense was identified [Figure 1]a.
Figure 1: (a) Result of polymerase chain reaction-restriction fragment length polymorphism analysis (PRA) studied mycobacteria lanes: M, size marker (50-bp); 1– 7, Mycobacterium fortuitum Type I; 8, unknown 9, Mycobacterium kansasil Type L; 10, Mycobacterium smegmatis, 11, Mycobacterium avium, 12– 13, Mycobacterium gordonae Type II. (b) Result of PRA studied mycobacteria lanes: M, size marker (50-bp); 1– 6, Mycobacterium fortuitum Type I; 7, Mycobacterium intracelluare; 8, Mycobacterium conceptionese, 9, Mycobacterium kansasii Type 1; 10, Mycobacterium fortuitum Type I, 11– 13, Mycobacterium gordonae Type I

Click here to view

Among 11 slowly growing isolates, Mycobacterium avium, 1isolate; Mycobacterium intracellular subtype I, 1 isolate; Mycobacterium kansasii Type I, 2 isolates; Mycobacterium gordonae Type I, 4 isolates, Type II, 2 isolates and Type IV, 1 isolate were identified. Two isolates that were identified as M. gordonae by phenotypic methods were not identified by PRA because their profiles were not recorded in any of the last reported studies [Figure 1]b. After sequencing of 360 bp rpoB gene, they were identified as Actinosynnema mirum. The PCR-RFLP profiles of all isolates are listed in [Table 1]. M. fortuitum (ATCC 49403) shown its specific pattern (HaeIII: 70/100/12O, MspI: 80/90/175).

  Discussion Top

In the present study, based on PRA of rpoB gene, among the clinical NTM isolates, M. fortuitum Type I with the frequency of 57.77% was the most prevalent type isolates. Other dominant clinical isolates were M. gordonae Type I with the frequency of 8.88% and M. gordonae Type II, M. kansasii Type I, both with the frequency of 4.44%. In 2013, Choi et al., reported M. intracellular Type I, M. avium, M. abscessus with the frequency of 26.7%, 25% and19.5% respectively, as the species with higher frequency. [21] In 2012, Ong et al. reported M. fortuitum Type I, M. chelonae/M. abscessus, M. avium, M. fortuitum Type II and M. gordonae Type I with 30%, 22.85%, 10%, 11.42%, and 8.57%, respectively, as NTM species with higher frequency. [15] The PRA pattern in M. Intracellular subtype, I in the present study, was MspI: 175/105/80 while in other studies reported MspI: 175/100/80 and 170/95/80. [15],[22] These slight differences are ignorable as sometimes in different studies; the 10-12-bp variations are due to the not exactly identical procedures and standards. [11],[21],[22] PRA patterns of M. kansasii Type I (HeaIII: 90/205, MspI: 30/40/60/175) and M. avium (HeaIII: 270, MspI: 40/80/105) in this study are identical to some studies and different from the profile of some others. [15],[22] For example indifferent pattern, Macheras et al., showed that profiles of M. kansasii Type I and M. avium are MspI: 105/70/60/45/40; HaeIII: 210/80/60 and MspI: 180/95/40; HaeIII: 180/90, respectively. [23] In the present study, an unrecorded pattern of the profile of M. conceptionense (HeaIII: 120/90/80, MspI: 120/105/80) was identified which is not reported in other recorded studies. This result represents a wide inter-species diversity between mycobacteria in different geographical areas.

In this study, the applied method for identification and typing of isolates showed high discriminatory power and was able to identify 96.6% of all isolates. In 2012, Kazumia reported that the sequence of the gene rpoB improves diagnostic system and precise identification of mycobacteria species. [12] In 2010, Ong et al. applied PRA method for rpoB and hsp65 genes for typing simultaneously and showed that PRA of the rpoB gene was able to identify 85.6% of all isolates and is more efficient than hsp65 gene for identification and typing of rapidly growing mycobacteria. [15] Based on the results of Lee et al. in 2000, mycobacteria can be detected by PRA of rpoB gene to the level of species and subspecies. They reported that M. gordonae, M. kansasii, M. fortuitum, and M. gastrican be identified by this method too while 16SrRNA has some limitations in these cases. [11] In some studies, PRA for 16SrRNA, hsp65, and dnaj genes were also conducted but need other restriction enzymes, which made the interpretation of the patterns very difficult. Additionally, some studies confirm that rpoB is a powerful identification tool for identifying and typing of NTMs specially rapid growing mycobacteria and a powerful tool for investigating the inter and intra-species relativity. [22],[23],[24],[25],[26],[27],[28] In cases of strains that cannot be identified completely by one or two enzymes, a third enzyme can be helpful. [11] With regard to the fact that There are fewer restriction patterns documented in the rpoB PRA [11] this study showed that can be useful to improving the small database by introduced the new pattern in conclusion, based on result of this study, rpoB gene could be a potentially useful tool for rapid identification and investigation of molecular epidemiology of mycobacterial species.

  Acknowledgments Top

This assay was supported by grant 191123 from Isfahan University of Medical Sciences. The authors are grateful to Manuchehr Homaee and Maryam Mokhtari; staff of Tuberculosis Center, Isfahan Provincial Health Office, Isfahan University of Medical Sciences, and I.R. Iran.

  References Top

de Zwaan R, van Ingen J, van Soolingen D. Utility of rpoB gene sequencing for identification of nontuberculous mycobacteria in The Netherlands. J Clin Microbiol 2014;52:2544-51.  Back to cited text no. 1
Neonakis IK, Gitti Z, Krambovitis E, Spandidos DA. Molecular diagnostic tools in mycobacteriology. J Microbiol Methods 2008;75:1-11.  Back to cited text no. 2
Mathema B, Kurepina NE, Bifani PJ, Kreiswirth BN. Molecular epidemiology of tuberculosis: Current insights. Clin Microbiol Rev 2006;19:658-85.  Back to cited text no. 3
Okada M, Kobayashi K. Article in Japanese Recent progress in mycobacteriology. Kekkaku 2007;82:783-99.  Back to cited text no. 4
Tan Y, Hu Z, Zhao Y, Cai X, Luo C, Zou C, et al. The beginning of the rpoB gene in addition to the rifampin resistance determination region might be needed for identifying rifampin/rifabutin cross-resistance in multidrug-resistant Mycobacterium tuberculosis isolates from Southern China. J Clin Microbiol 2012;50:81-5.  Back to cited text no. 5
Moore JE, Kruijshaar ME, Ormerod LP, Drobniewski F, Abubakar I. Increasing reports of non-tuberculous mycobacteria in England, Wales and Northern Ireland, 1995-2006. BMC Public Health 2010;10:612.  Back to cited text no. 6
Griffith DE, Aksamit T, Brown-Elliott BA, Catanzaro A, Daley C, Gordin F, et al. An official ATS/IDSA statement: Diagnosis, treatment, and prevention of nontuberculous mycobacterial diseases. Am J Respir Crit Care Med 2007;175:367-416.  Back to cited text no. 7
Khan IU, Yadav JS. Development of a single-tube, cell lysis-based, genus-specific PCR method for rapid identification of mycobacteria: Optimization of cell lysis, PCR primers and conditions, and restriction pattern analysis. J Clin Microbiol 2004;42:453-7.  Back to cited text no. 8
Covert TC, Rodgers MR, Reyes AL, Stelma GN Jr. Occurrence of nontuberculous mycobacteria in environmental samples. Appl Environ Microbiol 1999;65:2492-6.  Back to cited text no. 9
Lee H, Bang HE, Bai GH, Cho SN. Novel polymorphic region of the rpoB gene containing Mycobacterium species-specific sequences and its use in identification of mycobacteria. J Clin Microbiol 2003;41:2213-8.  Back to cited text no. 10
Lee H, Park HJ, Cho SN, Bai GH, Kim SJ. Species identification of mycobacteria by PCR-restriction fragment length polymorphism of the rpoB gene. J Clin Microbiol 2000;38:2966-71.  Back to cited text no. 11
Kazumia SM. The evaluation of an identification algorithm for Mycobacterium species using the 16S rRNA coding gene and rpoB. Int J Mycobacteriol 2012;1:21-8.  Back to cited text no. 12
Nasr Esfahani B, Rezaei Yazdi H, Moghim S, Ghasemian Safaei H, Zarkesh Esfahani H. Rapid and accurate identification of Mycobacterium tuberculosis complex and common non-tuberculous mycobacteria by multiplex real-time PCR targeting different housekeeping genes. Curr Microbiol 2012;65:493-9.  Back to cited text no. 13
Varma-Basil M, Garima K, Pathak R, Dwivedi SK, Narang A, Bhatnagar A, et al. Development of a novel PCR restriction analysis of the hsp65 gene as a rapid method to screen for the Mycobacterium tuberculosis complex and nontuberculous mycobacteria in high-burden countries. J Clin Microbiol 2013;51:1165-70.  Back to cited text no. 14
Ong CS, Ngeow YF, Yap SF, Tay ST. Evaluation of PCR-RFLP analysis targeting hsp65 and rpoB genes for the typing of mycobacterial isolates in Malaysia. J Med Microbiol 2010;59:1311-6.  Back to cited text no. 15
Pourahmad F, Thompson KD, Adams A, Richards RH. Comparative evaluation of polymerase chain reaction-restriction enzyme analysis (PRA) and sequencing of heat shock protein 65 (hsp65) gene for identification of aquatic mycobacteria. J Microbiol Methods 2009;76:128-35.  Back to cited text no. 16
Tenover FC, Arbeit RD, Goering RV, Mickelsen PA, Murray BE, Persing DH, et al. Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: Criteria for bacterial strain typing. J Clin Microbiol 1995;33:2233-9.  Back to cited text no. 17
Adékambi T, Drancourt M, Raoult D. The rpoB gene as a tool for clinical microbiologists. Trends Microbiol 2009;17:37-45.  Back to cited text no. 18
Dvorska L, Bartos M, Martin G, Erler W, Pavlik I. Strategies for differentiation, identification and typing of medically important species of mycobacteria by molecular methods. Vet Med Czech 2001;46:309-28.  Back to cited text no. 19
Clarke JD. Cetyltrimethyl ammonium bromide (CTAB) DNA miniprep for plant DNA isolation. Cold Spring Harb Protoc 2009;2009:pdb.prot5177.  Back to cited text no. 20
Choi SR, Kang MJ, Park GH, Kim DH, Jeong DW, Seo EH, et al. Species identification of nontuberculous mycobacteria (NTM) by PCR-restriction fragment length polymorphism (PRA) of the rpoB Gene from Three Hospitals of Busan-Kyeongnam Area. Korean J Clin Lab Sci 2013;45:48-53.  Back to cited text no. 21
Whang J, Lee BS, Choi GE, Cho SN, Kil PY, Collins MT, et al. Polymerase chain reaction-restriction fragment length polymorphism of the rpoB gene for identification of Mycobacterium avium subsp. Paratuberculosis and differentiation of Mycobacterium avium subspecies. Diagn Microbiol Infect Dis 2011;70:65-71.  Back to cited text no. 22
Macheras E, Roux AL, Bastian S, Leão SC, Palaci M, Sivadon-Tardy V, et al. Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (sensu lato) strains. J Clin Microbiol 2011;49:491-9.  Back to cited text no. 23
Simmon KE, Low YY, Brown-Elliott BA, Wallace RJ Jr, Petti CA. Phylogenetic analysis of Mycobacterium aurum and Mycobacterium neoaurum with redescription of M. aurum culture collection strains. Int J Syst Evol Microbiol 2009;59:1371-5.  Back to cited text no. 24
Adékambi T, Colson P, Drancourt M. rpoB-based identification of nonpigmented and late-pigmenting rapidly growing mycobacteria. J Clin Microbiol 2003;41:5699-708.  Back to cited text no. 25
Kim BJ, Lee SH, Lyu MA, Kim SJ, Bai GH, Chae GT, et al. Identification of mycobacterial species by comparative sequence analysis of the RNA polymerase gene (rpoB). J Clin Microbiol 1999;37:1714-20.  Back to cited text no. 26
Devulder G, Pérouse de Montclos M, Flandrois JP. A multigene approach to phylogenetic analysis using the genus Mycobacterium as a model. Int J Syst Evol Microbiol 2005;55:293-302.  Back to cited text no. 27
Adékambi T, Berger P, Raoult D, Drancourt M. rpoB gene sequence-based characterization of emerging non-tuberculous mycobacteria with descriptions of Mycobacterium bolletii sp. nov. Mycobacterium phocaicum sp. nov. and Mycobacterium aubagnense sp. nov. Int J Syst Evol Microbiol 2006;56:133-43.  Back to cited text no. 28


  [Figure 1]

  [Table 1]

This article has been cited by
1 Rapid molecular diagnosis of live Mycobacterium tuberculosis on an integrated microfluidic system
Chih-Hung Wang, Jia-Ru Chang, Shang-Cheng Hung, Horng-Yunn Dou, Gwo-Bin Lee
Sensors and Actuators B: Chemical. 2022; : 131968
[Pubmed] | [DOI]


Previous article  Next article
Similar in PUBMED
   Search Pubmed for
   Search in Google Scholar for
 Related articles
Access Statistics
Email Alert *
Add to My List *
* Registration required (free)

  In this article
Materials and me...
Article Figures
Article Tables

 Article Access Statistics
    PDF Downloaded221    
    Comments [Add]    
    Cited by others 1    

Recommend this journal