Original Research ARTICLE
Evaluation of a High-Intensity Green Fluorescent Protein Fluorophage Method for Drug- Resistance Diagnosis in Tuberculosis for Isoniazid, Rifampin, and Streptomycin
- 1National Clinical Laboratory on Tuberculosis, Beijing Key Laboratory for Drug-Resistant Tuberculosis Research, Beijing Chest Hospital, Beijing Tuberculosis and Thoracic Tumor Institute, Capital Medical University, Beijing, China
- 2Department of Microbiology and Immunology, Albert Einstein College of Medicine, Bronx, NY, USA
- 3Division of Pulmonary, Allergy, and Critical Care Medicine, Columbia University Medical Center, New York, NY, USA
- 4Department of Epidemiology, Mailman School of Public Health, Columbia University Medical Center, New York, NY, USA
- 5Howard Hughes Medical Institute, Chevy Chase, MD, USA
A novel method for detecting drug resistance in Mycobacterium tuberculosis using mycobacteriophage Φ2GFP10 was evaluated with clinical isolates. The phage facilitates microscopic fluorescence detection due to the high expression of green fluorescence protein which also simplifies the operative protocol as well. A total of 128 clinical isolates were tested by the phage assay for isoniazid (INH), rifampin (RIF), and streptomycin (STR) resistance while conventional drug susceptibility test, by MGIT960, was used as reference. The sensitivities of Φ2GFP10 assay for INH, RIF, and STR resistance detection were 100, 98.2, and 89.3%, respectively while their specificities were 85.1, 98.6, and 95.8%, respectively. The agreement between phage and conventional assay for detecting INH, RIF, and STR resistance was 92.2, 98.4, and 93.0%, respectively. The Φ2GFP10-phage results could be available in 2 days for RIF and STR, while it takes 3 days for INH, with an estimated cost of less than $2 to test all the three antibiotics. The Φ2GFP10-phage method has the potential to be a valuable, rapid and economical screening method for detecting drug-resistant tuberculosis.
Reduction in transmission of drug resistant tuberculosis (TB) and improved patient management requires timely diagnosis of drug resistant bacilli. Presence of drug-resistant M. tuberculosis bacilli is confirmed by genotypic and phenotypic drug susceptibility testing (DST). Genotypic approaches with short turnaround time are based on identification of well-known antibiotic resistance-conferring gene mutations. Despite active research, not all TB drugs can have genotypic drug-susceptibility testing since the genetic basis of resistance may be complex or incompletely characterized (Kruuner et al., 2003; Jain et al., 2011; Cui et al., 2013; Zhang et al., 2015). Phenotypic assays observe the bacterial response to antibiotics in vitro without limitation to any particular antibiotic, allele, or working mechanism although the phenotypic DST is often time consuming.The conventional DST, based on solid medium, takes about 4–6 weeks after the isolation of M. tuberculosis, on the other hand BACTEC MGIT960 requires 10–14 days after acquiring a positive culture. Therefore, it is imperative to outline a further rapid, accurate, inexpensive DST method for the diagnosis of drug resistance M. tuberculosis.
Phage based methods have been used for drug resistance detection of M. tuberculosis since about two decades ago, such as bacteriophage D29 (Wilson et al., 1997) and luciferase reporter phages (LRPs) (Jacobs et al., 1993). A novel, high-intensity mycobacteria-specific fluorophage (Φ2GFP10) was described recently with good results in pre-clinical evaluation of drug resistant tuberculosis (Jain et al., 2012). Φ2GFP10 is an improved second generation fluorescent reporter phage which expresses gfp (green fluorescence protein) under the control of PL promoter of mycobacterium phage L5 (Guo and Ao, 2012; Jain et al., 2012) and have intensity 100 times brighter than the previous generation of fluorescent reporter phages (Piuri et al., 2009). Unlike the LRPs (Jacobs et al., 1993), the Φ2GFP10 reporter phage does not require exogenous luciferase substrate and can yield more stable and microscopically detectable intensive fluorescence. The presented study has evaluated the performance of this in-house fluorophage method to determine the rifampicin (RIF), isoniazid (INH), and streptomycin (STR) resistance in a high drug-resistant setting. Here we have evaluated the performance of Φ2GFP10 phage for detecting drug resistance from M. tuberculosis clinical isolates, and also developed a new phage assay method using fluorescent microplate reader.
Materials and Methods
A total of 128 clinical M. tuberculosis strains were isolated from patients with suspected drug-resistant TB patients visiting the Beijing chest hospital (Beijing, China) from April to June 2014. All of the isolates were identified as M. tuberculosis complex (MTBC) strains by performing a growth test on 500 μg/ml p-Nitrobenzoic Acid containing Löwenstein-Jensen medium (Tsukamura and Tsukamura, 1964). The drug susceptibility of the isolates was determined by MGIT960 SIRE kits according to the manufacturer's protocol. The critical concentrations used were as follows: INH (0.1 μg/ml), RIF(1.0 μg/ml), STR(1.0 μg/ml). Ultimately, 61, 56, and 56 were defined as resistant to INH, RIF and STR respectively by MGIT960 system. The laboratory M. tuberculosis H37Rv (ATCC 27294) strain was used in all batches as control.Technique round was demonstrated in Figure 1.
Phage Stock Preparation
The fluorophage, Φ2GFP10 was constructed by Dr William R. Jacobs' laboratory in Albert Einstein University (New York, USA). The stocks were prepared by growing M. smegmatis strain (mc2155) in the 7H9, containing 10% oleic albumin dextrose catalase (OADC), to an optical density of 1.0 detected by spectrophotometer at 600 nm wavelength. Then 300 μl cell suspension and 200 μl Φ2GFP10 were mixed and incubated at room temperature for 30 min. Subsequently, 3 ml 7H9 containing 0.75% agar was added to the tube, the contents were briefly mixed and poured onto a 100-mm Petri dish containing 7H10-OADC. After the top-agar was solidified, the dish was incubated at 30°C for 48 h. Then 3 ml MP buffer [50 mM Tris (pH 7.6), 150 mM NaCl, 10 mM MgCl2, and 2 mM CaCl2] was pipetted onto the plates and then incubated on an orbital shaker for 6 h. The buffer was then collected and passed through a 0.22 μm Millipore filter membrane. The phage titers were determined by serial dilution method and the phage density was adjusted to approximately 1010 plaque forming unit (pfu)/ml for further use. The phage stock was stable when stored at 4°C as no appreciable drop in titer occurred for at least 6 months.
Φ2GFP10 Assay by Microscopy Examination
INH, RIF, and STR were obtained from Sigma-Aldrich. INH and STR were dissolved in water; RIF was dissolved in DMSO. The stock solutions were 20 mg/ml, 2 mg/ml and 4 mg/ml for INH, RIF, and STR, respectively. All stock solutions were stored at −80°C during experimentation. Fresh colonies on Löwenstein-Jensen media were suspended and homogenized with sterile saline and the cell density was adjusted to 1 McFarland turbidity. A 1 ml aliquot was harvested by centrifugation and the sediment was re-suspended with 250 μl 7H9-OADC in the absence of Tween 80. The desired drugs would be added at this step, when needed, at the following finial concentrations: 20 μg/ml for INH, 2 μg/ml for RIF, or 4 μg/ml for STR. The tubes were incubated at 37°C for 24 h, or 48 h specifically for INH tests. From the stock, containing 1010 pfu fluorophage per ml, 100 μl volume was added to each tube to obtain a multiplicity of infection (MOI) of 100, the tubes were then incubated at 37°C for 16 h. After adding 350 μl of 4% paraformaldehyde in phosphate-buffered saline (PBS) to each tube, the reaction mixtures, left at room temperature for 90 min to ensure killing of the bacteria, were then centrifuged at 12,000 g for 15 min to remove excess phage and media. The pellet was re-suspended with 10 μl 7H9 media. Five micro liter re-suspension was spotted on a glass slide and visualized by a fluorescence microscope. The criterion for resistance with phage assay was the presence of at least one fluorescent bacillus per 50 high power fields (Rondon et al., 2011).
Φ2GFP10 Assay by Microplate Reader
To detect fluorescence of Φ2GFP10 in a 96 well format relative light unit (RLU) of fluorescence was measured by a fluorescent microplate reader (Berthold LB970, Germany). The wavelength of excitation and emitted light was 485 and 535 nm respectively and the counting time was 5 s. Each test was performed in triplicate and the mean value was used for interpretation. After removing excess phage and media by centrifuging, the pellet was re-suspended in 200 μl PBS, and then, the RLU of the plate was read. The influence of drug exposure on the fluorescence production of the bacilli was interpreted by the remaining fluorescence rate (RFR), calculated by the following formula: (reaction counts-blank counts)/(positive control counts –blank counts) × 100% (“reaction” indicates the well with drug exposure; “blank” indicates the well without drug and phage; “positive control” indicates the reaction without drug). Precision of the assay was determined with quintuple measurements of the following three strains: H37Rv which was susceptible to all the three tested drugs, strain 14,301 which was resistant to all the three tested drugs and strain 14,161 which was resistant to INH and STR but susceptible to RIF. The stability of fluorescence was determined by detecting the same samples at 0, 4, and 24 h, respectively. All the reactions were performed in triplicate.
According to MGIT960 phenotypic outcomes, receiver-operation characteristic (ROC) analysis was performed using SPSS (version 17.0) and used to define a cutoff value for phage plate reader assay. Kappa value was calculated to compare the method between Φ2GFP10 assay and MGIT960 phenotypic outcomes. Sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV) was calculated from the following URL: http://vassarstats.net. The discriminative power of fluorescent microplate reader for drug resistance was analyzed using ROC curve and area under the curve (AUC). The optimal cut-off value was defined as the one with the least (1 − sensitivity)2 + (1 − specificity)2 value (Shu et al., 2013).
Φ2GFP10 Assay by Microscopy
According to the outcomes of MGIT960, the sensitivity of Φ2GFP10 assay for detecting INH, RIF and STR resistance were 100, 98.2, and 89.3% respectively whereas the specificities were 85.1, 98.6, and 95.8% respectively. The agreement between phage assay and MGIT960 for detecting INH, RIF and STR resistance was 92.2, 98.4, and 93.0%, respectively (see Table 1).The Kappa coefficient between Φ2GFP10 assay and MGIT960 phenotypic outcomes for INH, RIF and STR were 0.85, 0.97, and 0.856, respectively.
Φ2GFP10 Assay by Fluorescent Microplate Reader
In the preliminary validation of the assay among the 3 tested isolates, the stability tests demonstrated that the fluorescence was stable for at least 24 h at 4°C (see Table 2). The coefficient of variation (CV) of RLU among the tested isolates in quintuple ranged from 1.22 to 6.73%, with the mean value at 4.67 ± 1.26%. In the validation assays using clinical isolates, in contrast to the phenotypic DST, the AUC under the ROC curve of INH, RIF and STR were 0.957, 0.960, 0.917, and the optimal cutoff value for each drug has been listed in Table 3 and Figure 2. The sensitivity of Φ2GFP10 assay by plate reader for detecting INH, RIF and STR resistance, according to the cutoff values, were 86.9, 89.3, and 83.9%, respectively, while the specificity for INH, RIF and STR resistance detection were 94.0, 90.3, and 87.5%, respectively. The Kappa coefficient between Φ2GFP10 assay and MGIT960 phenotypic outcomes for INH, RIF and STR were 0.73, 0.79, and 0.71, respectively. Furthermore, evaluation of the RFR value according to the susceptibility predictability demonstrated that when the RFR is lower than 45% the NPV for INH, RIF and STR were 98.0, 98.2, and 87.5%, respectively. Similarly, when the ratio was higher than 60%, the PPV for INH, RIF and STR were 92.6, 97.9, and 95.3%, respectively. When the ratio was between 45 and 60%, then the outcome should be interpreted as undetermined (see Table 4).
Figure 2. Receiver operating characteristic curves of fluorescent microplate reader in detecting drug resistance according to BactecMGIT 960 system. (A) INH; (B) RIF; (C) STR.
Table 4. Results and performance parameters by sectional outcomes interpretation: fluorescent microplate reader vs. MGIT960 system.
Significant efforts are being made to develop rapid, simple and accurate tests for drug- resistant M. tuberculosis. One of such technologies currently been worked at is based on using bacteriophages (Banaiee et al., 2008; Piuri et al., 2009; Rondon et al., 2011; Jain et al., 2012; Sivaramakrishnan et al., 2013). LRP have been used for identifying drug resistance by detecting the luminescence produced by LRP after their infection in M. tuberculosis. These assays classify the samples as drug-susceptible if no luminescence is detected in the luminescence drug-containing samples. A meta-analysis showed that the sensitivity and specificity of luciferase reporter phage for rapid detection of RIF resistance in M. tuberculosis ranged from 92 to 100% and 89 to 100% (Minion and Pai, 2010). Although this method has shown high sensitivity and high specificity, a relatively expensive luminometer is required. The newly constructed phage Φ2GFP10 facilitates the microscopic signal detection since it yields 100-folds higher fluorescence signal per-cell than any previously described reporter phages (Jain et al., 2011). Meanwhile, as the intensive fluorescence increased the contrast between reaction and background, tedious washing steps were not necessary for phage assay, so the Φ2GFP10 assay was easier to perform than the other phage assays (Rondon et al., 2011).
In this study, drug resistance detection by Φ2GFP10 assay yielding applicable sensitivities and specificities, especially for RIF. A recent report also obtained great consistency between Φ2GFP10 assay and GeneXpert for RIF resistance diagnosis (O'Donnell et al., 2015). Fluorescent microscopy may be a limiting factor for Φ2GFP10 assay application due to its cost, whereas less expensive and simpler fluorescent microscopes using light emitting diodes (LED) have been applied for smear test of tuberculosis in various clinical laboratories (Marzouk et al., 2013; Reza et al., 2013; Xia et al., 2013) and they could readily be adapted for Φ2GFP10 assay. The reagents for fluorophage growth and amplification are inexpensive, safe, and universally available. The total cost for the 3-drug test per strain was less than US $2. Additionally, fixation of the sample with paraformaldehyde prior to analysis overcomes the substantial biosafety concerns.
In our study, even though a longer drug exposure time of 48 h for INH, the specificity was still lower than those of RIF and STR whose drug exposure time was 24 h. Those observations were consistent with other mycobacteriophage-based assays (Rondon et al., 2011). The anti-TB mechanism of INH involves inhibiting the synthesis of mycolic acid of cell wall (Eltringham et al., 1999), which is a slow process compared with the quickly bactericidal activities of RIF and STR, therefore the bacilli can stay alive for quite a while even when they are sensitive to INH. The fluorophage could infect such bacilli and produce fluorescence signal which decreased the specificity of the assay. To enhance the bactericidal activity of INH, we increased the INH concentration for phage assay to 20 μg/ml, which was dramatically higher than those of previous studies (Xia et al., 2013; Zhang et al., 2015). However, the specificity was still lower for INH than those of RIF and STR, which indicating the INH tolerance deviation for clinical strains was universal.
First-line-drug ethambutol (EMB) also affects the cell wall construction in M. tuberculosis, and is a bacteriostatic drug rather than a bactericidal drug. We attempted to develop a Φ2GFP10 assay for EMB, but the pilot test demonstrated that even an exposure concentration of 50 μg/ml could not lead to obvious CFU loss. A previous study reported low concordance (53%) between the phage amplified biologically assay and resistance ratio method for EMB, which they attributed to the bacteriostatic role of EMB both in vitro and in a macrophage model (Eltringham et al., 1999). Therefore, the Φ2GFP10 assay might only be feasible to test the rapid-action bacteriocidal drugs.
Although we acquired favorable sensitivity and specificity for Φ2GFP10 assay by microscopy, we found that microscopic examination process is experience dependent, and a lot of time is required when handing a batch of samples. So we developed a phage assay to interpret the outcomes simply by fluorescence plate reader, which can read a 96-well plate in 10 min. We monitored the fluorescence diminishing after the bacilli were exposed to the tested drugs. When categorized the RFR lower than 45% as “susceptible,” greater than 60% as “resistant,” while between 45 and 60% as undetermined, plate reader could interpret over 80% of the test strains successfully for all the 3 drugs, with 95.24% (100/105), 98.06% (101/103), and 90.65% (97/107) consistency with MGIT960 system for INH, RIF, and STR susceptibility tests, respectively. Therefore, we suggest interpreting the outcomes in sections to increase the reliability of the assay. For undetermined samples, we suggest to recheck by microscopy to produce more reliable results.
Like other fluoromycobacteriophages, Φ2GFP10 assay could not discriminate between M. tuberculosis and other members of MTBC species (Riska et al., 1997). For the settings with high isolation rate of nontuberculous mycobacteria, one should be cautious to report drug resistance without MTBC identification.
In summary, the Φ2GFP10 assay is a rapid, simple and inexpensive DST method. In clinical laboratories with limited resources, the Φ2GFP10 based assay has the potential to be used for drug resistance diagnosis of tuberculosis.
HH designed the study. HH and XY wrote the paper. WJ modified the paper. XY and YG performed the Φ2GFP10 assay by microscopy. XY and ZS performed Φ2GFP10 assay by Microplate reader. GJ, YM, and LZ performed DST by MGIT960 kits. XY conducted statistical analysis. PJ, MO, and ML established the assay for phage stock preparation. All authors reviewed the results and approved the final version of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The work was supported by the research funding from Infectious Diseases Special Project, Minister of Health of China (2012ZX10003002-009, 2016ZX10003001-12), The capital health research and development of special (2016-4-1042). WJ is supported by the National Institute of Health (AI26170). ML is supported by NIH (1R01AI114900).
Banaiee, N., January, V., Barthus, C., Lambrick, M., Roditi, D., Behr, M. A., et al. (2008). Evaluation of a semi-automated reporter phage assay for susceptibility testing of Mycobacterium tuberculosis isolates in South Africa. Tuberculosis (Edinb) 88, 64–68. doi: 10.1016/j.tube.2007.08.006
Cui, Z., Wang, J., Lu, J., Huang, X., Zheng, R., and Hu, Z. (2013). Evaluation of methods for testing the susceptibility of clinical Mycobacterium tuberculosis isolates to pyrazinamide. J. Clin. Microbiol. 51, 1374–1380. doi: 10.1128/jcm.03197-12
Eltringham, I. J., Wilson, S. M., and Drobniewski, F. A. (1999). Evaluation of a bacteriophage-based assay (phage amplified biologically assay) as a rapid screen for resistance to isoniazid, ethambutol, streptomycin, pyrazinamide, and ciprofloxacin among clinical isolates of Mycobacterium tuberculosis. J. Clin. Microbiol. 37, 3528–3532.
Jacobs, W. R. Jr. Barletta, R. G., Udani, R., Chan, J., Kalkut, G., Sosne, G., et al. (1993). Rapid assessment of drug susceptibilities of Mycobacterium tuberculosis by means of luciferase reporter phages. Science 260, 819–822. doi: 10.1126/science.8484123
Jain, P., Hartman, T. E., Eisenberg, N., O'Donnell, M. R., Kriakov, J., Govender, K., et al. (2012). phi(2)GFP10, a high-intensity fluorophage, enables detection and rapid drug susceptibility testing of Mycobacterium tuberculosis directly from sputum samples. J. Clin. Microbiol. 50, 1362–1369. doi: 10.1128/jcm.06192-11
Jain, P., Thaler, D. S., Maiga, M., Timmins, G. S., Bishai, W. R., Hatfull, G. F., et al. (2011). Reporter phage and breath tests: emerging phenotypic assays for diagnosing active tuberculosis, antibiotic resistance, and treatment efficacy. J. Infect. Dis. 204(suppl. 4), S1142–S1150. doi: 10.1093/infdis/jir454
Kruuner, A., Jureen, P., Levina, K., Ghebremichael, S., and Hoffner, S. (2003). Discordant resistance to kanamycin and amikacin in drug-resistant Mycobacterium tuberculosis. Antimicrob. Agents Chemother. 47, 2971–2973. doi: 10.1128/AAC.47.9.2971-2973.2003
Marzouk, M., Ferjani, A., Dhaou, M., Ali, M. H., Hannachi, N., and Boukadida, J. (2013). Comparison of LED and conventional fluorescence microscopy for detection of acid-fast bacilli in an area with high tuberculosis incidence. Diagn. Microbiol. Infect. Dis. 76, 306–308. doi: 10.1016/j.diagmicrobio.2013.03.023
O'Donnell, M. R., Pym, A., Jain, P., Munsamy, V., Wolf, A., Karim, F., et al. (2015). A novel reporter phage to detect tuberculosis and rifampin resistance in a high-HIV-burden population. J. Clin. Microbiol. 53, 2188–2194. doi: 10.1128/jcm.03530-14
Piuri, M., Jacobs, W. R. Jr., and Hatfull, G. F. (2009). Fluoromycobacteriophages for rapid, specific, and sensitive antibiotic susceptibility testing of Mycobacterium tuberculosis. PLoS ONE 4:e4870. doi: 10.1371/journal.pone.0004870
Reza, L. W., Satyanarayna, S., Enarson, D. A., Kumar, A. M., Sagili, K., Kumar, S., et al. (2013). LED-fluorescence microscopy for diagnosis of pulmonary tuberculosis under programmatic conditions in India. PLoS ONE 8:e75566. doi: 10.1371/journal.pone.0075566
Riska, P. F., Jacobs, W. R. Jr. Bloom, B. R., McKitrick, J., and Chan, J. (1997). Specific identification of Mycobacterium tuberculosis with the luciferase reporter mycobacteriophage: use of p-nitro-alpha-acetylamino-beta-hydroxy propiophenone. J Clin Microbiol 35, 3225–3231.
Rondon, L., Piuri, M., Jacobs, W. R. Jr. De Waard, J., Hatfull, G. F., and Takiff, H. E. (2011). Evaluation of fluoromycobacteriophages for detecting drug resistance in Mycobacterium tuberculosis. J. Clin. Microbiol. 49, 1838–1842. doi: 10.1128/jcm.02476-10
Shu, C. C., Wu, M. F., Hsu, C. L., Huang, C. T., Wang, J. Y., Hsieh, S. L., et al. (2013). Apoptosis-associated biomarkers in tuberculosis: promising for diagnosis and prognosis prediction. BMC Infect. Dis. 13:45. doi: 10.1186/1471-2334-13-45
Sivaramakrishnan, G., Subramanyam, B., C, P., and Kumar, V. (2013). Luciferase reporter phage phAE85 for rapid detection of rifampicin resistance in clinical isolates of Mycobacterium tuberculosis. Asian Pac. J. Trop. Med. 6, 728–731. doi: 10.1016/s1995-7645(13)60127-3
Tsukamura, M., and Tsukamura, S. (1964). Differentiation of Mycobacterium tuberculosis and Mycobacterium bovis by p-nitrobenzoic acid susceptibility. Tubercle 45, 64–65. doi: 10.1016/S0041-3879(64)80091-X
Wilson, S. M., Al-Suwaidi, Z., McNerney, R., Porter, J., and Drobniewski, F. (1997). Evaluation of a new rapid bacteriophage-based method for the drug susceptibility testing of Mycobacterium tuberculosis. Nat. Med. 3, 465–468. doi: 10.1038/nm0497-465
Xia, H., Song, Y. Y., Zhao, B., Kam, K. M., O'Brien, R. J., Zhang, Z. Y., et al. (2013). Multicentre evaluation of Ziehl-Neelsen and light-emitting diode fluorescence microscopy in China. Int. J. Tuberc. Lung. Dis. 17, 107–112. doi: 10.5588/ijtld.12.0184
Zhang, X., Liu, L., Zhang, Y., Dai, G., Huang, H., and Jin, Q. (2015). Genetic determinants involved in p-aminosalicylic acid resistance in clinical isolates from tuberculosis patients in northern China from 2006 to 2012. Antimicrob. Agents Chemother. 59, 1320–1324. doi: 10.1128/aac.03695-14
Keywords: tuberculosis, mycobacteriophage, drug resistance, diagnosis, Φ2GFP10
Citation: Yu X, Gu Y, Jiang G, Ma Y, Zhao L, Sun Z, Jain P, O'Donnell M, Larsen M, Jacobs WR Jr and Huang H (2016) Evaluation of a High-Intensity Green Fluorescent Protein Fluorophage Method for Drug- Resistance Diagnosis in Tuberculosis for Isoniazid, Rifampin, and Streptomycin. Front. Microbiol. 7:922. doi: 10.3389/fmicb.2016.00922
Received: 25 August 2015; Accepted: 31 May 2016;
Published: 17 June 2016.
Edited by:Julio Alvarez, University of Minnesota, USA
Reviewed by:Kuangnan Xiong, Health Research, Inc., USA
Raju Mukherjee, Indian Institute of Science Education and Research, India
Copyright © 2016 Yu, Gu, Jiang, Ma, Zhao, Sun, Jain, O'Donnell, Larsen, Jacobs and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hairong Huang, firstname.lastname@example.org