Theoretical basis for reducing time-lines to the determination of positive Mycobacterium tuberculosis cultures using thymidylate kinase (TMK) assays
© Wayengera; licensee BioMed Central Ltd. 2009
Received: 19 December 2008
Accepted: 18 March 2009
Published: 18 March 2009
In vitro culture of pathogens on growth media forms a "pillar" for both infectious disease diagnosis and drug sensitivity profiling. Conventional cultures of Mycobacterium tuberculosis (M.tb) on Lowenstein Jensen (LJ) medium, however, take over two months to yield observable growth, thereby delaying diagnosis and appropriate intervention. Since DNA duplication during interphase precedes microbial division, "para-DNA synthesis assays" could be used to predict impending microbial growth. Mycobacterial thymidylate kinase (TMKmyc) is a phosphotransferase critical for the synthesis of the thymidine triphosphate precursor necessary for M.tb DNA synthesis. Assays based on high-affinity detection of secretory TMKmyc levels in culture using specific antibodies are considered. The aim of this study was to define algorithms for predicting positive TB cultures using antibody-based assays of TMKmyc levels in vitro.
Methods and results
Systems and chemical biology were used to derive parallel correlation of "M.tb growth curves" with "TMKmyc curves" theoretically in four different scenarios, showing that changes in TMKmyc levels in culture would in each case be predictive of M.tb growth through a simple quadratic curvature, |tmk| = at2+ bt + c, consistent with the "S" pattern of microbial growth curves. Two drug resistance profiling scenarios are offered: isoniazid (INH) resistance and sensitivity. In the INH resistance scenario, it is shown that despite the presence of optimal doses of INH in LJ to stop M.tb proliferation, bacilli grow and the resulting phenotypic growth changes in colonies/units are predictable through the TMKmyc assay. According to our current model, the areas under TMKmyc curves (AUC, calculated as the integral ∫(at2+ bt + c)dt or ~1/3 at3+ 1/2 bt2+ct) could directly reveal the extent of prevailing drug resistance and thereby aid decisions about the usefulness of a resisted drug in devising "salvage combinations" within resource-limited settings, where second line TB chemotherapy options are limited.
TMKmyc assays may be useful for reducing the time-lines to positive identification of Mycobacterium tuberculosis (M.tb) cultures, thereby accelerating disease diagnosis and drug resistance profiling. Incorporating "chemiluminiscent or fluorescent" strategies may enable "photo-detection of TMKmyc changes" and hence automation of the entire assay.
Infection with Mycobacterium tuberculosis (M.tb), the causative agent of tuberculosis (TB), is one of the leading global health challenges [1, 2]. An estimated 8–10 million persons acquire tuberculosis annually, 2 million of whom die [1, 2]. The global TB epidemic has been complicated by the human immunodeficiency virus (HIV) co-epidemic . Interaction between HIV and TB: (i) is associated with a higher risk of progression to active M. tuberculosis infection (ATBI) among persons with existing latent infection (LTBI); (ii) leads to increased susceptibility to new infection with M.tb; (iii) renders diagnosis difficult and treatment/cure rates slow; (iv) results in a higher incidence of relapses; and (v) favors the evolution of drug resistance [1, 4–6]. Amongst these, the emergence of drug resistance forms the deadliest challenge to controlling the TB epidemic. In HIV/TB high burden areas such as South Africa, drug resistance has extended from the first line anti-TB drugs to include the second line drugs spared for multi-drug resistant (MDR) TB, XDR . Early detection of TB has a crucial role to play in controlling the epidemic here. In the past, diagnosis has been based on prediction of prior infection using the purified protein derivative (PPD)-tuberculin skin test (TST), detecting active TB by sputum smear staining using the Zeihl Neilson Stain (ZN) or culturing the organism on Lowenstein Jensen (LJ) medium; and radiographic imaging for TB-associated pathology [8, 9]. More recently, newer molecular assays for TB have emerged based on serology (detecting the 38 kDa antigen, Early Secretory Antigen or ESAT-6, CFP10 and other secretory antigens), nucleic acid amplification (NAATS), phage amplification and line probe assays for MDR . Despite the advent of these modern TB assays, largely because there are no inexpensive technology platforms exploiting related biomarkers, approaches based on in vitro cultures still form the most reliable and readily affordable method for diagnosing TB in many resource-limited settings . In general, in vitro culture of pathogens on appropriate growth media forms a "pillar" for both infectious disease diagnosis and drug sensitivity profiling .
Conventional cultures of the slowly-growing tubercle bacilli (M.tb) on LJ medium, though highly sensitive for TB , are nevertheless time-consuming since they take about two months to yield observable growth, thereby delaying disease diagnosis and appropriate intervention. Although reading of culture results is still widely based on physical observation of the formation of bacterial colonies, modified assays that predict changes in 'in vitro M.tb growth" by monitoring turbidity or metabolism of labeled metabolites are occasionally used to confirm the presence of actively growing tubercle bacilli [9–13] and to characterize drug resistance [9, 14–16]. Specifically, some of these modified culture assays such as phage amplification assays, line probe assays and colorimetric redox-indicator methods can predict positive M.tb cultures (or drug resistance) in shorter times than the stipulated 2–3 months of culture on unmodified LJ medium. Most of these emerging TB assays are, however, still not readily available for routine use in most resource-limited settings where the TB burden is concurrently high . We hypothesized that metabolomic-based assays of "para-DNA duplicative" changes during interphase could predict impending M.tb growth before actual growth occurs. Moreover, easier platforms such as chemilumniscence or fluorescence photodetection may be incorporated into these bioassays, thereby making them more affordable than the aforementioned existing methods for TB detection . That hypothesis is derived from and based on several biological principles of the cell cycle, which are discussed below.
First, prokaryotic growth and proliferation, unlike that of eukaryotes, is predominantly achieved through division of single cells . The period between the "birth" of two daughter cells from the parent cell and their own division is termed interphase. During interphase, cells undergo a resting (S) phase during which, though little physical activity is noted microscopically, DNA duplication occurs . Since DNA synthesis precedes actual microbial division and proliferation, monitoring "DNA synthesis" may predict microbial growth changes well before the visible manifestation of actual microbial growth in culture.
Second, phosphorylation of deoxythymidine monophosphate (dTMP) is a critical step in the pathway leading to the synthesis of thymidylate triphosphate (TTP), a necessary precursor for DNA synthesis . This process is catalyzed by enzymes called kinases or phosphotransferases . Among mammalian cells, two phosphotransferases exist: thymidine kinases (types 1 and 2) and thymidylate kinase . Sherley and Kelly have previously reported a positive correlation between levels of human TK and DNA duplicative changes in the cell cycle . Within mycobacterial proteomes, however, no thymidine kinase has been detected . Instead, mycobacterial thymidylate kinase (TMKmyc) carries out this critical phosphorylative function . From this knowledge base, we proposed that levels of TMKmyc may similarly parallel DNA synthesis during the M.tb cell cycle [18, 19]. We recently focused on research and development of an antibody-based biomarker for TMKmyc. To date, we have identified highly conserved surface linear peptides of TMKmyc as potential B-cell immunogens for in situ engineering of TMK-specific monoclonal antibodies (MAbs). Such TMKmyc-specific MAb(s) could be used to monitor of levels of TMKmyc in cultures of M.tb. In this paper, we attempt to provide a theoretical framework upon which we propose reduction of time-lines for positive identification and drug resistance profiling among Mycobacterium tuberculosis (M. tb ) cultures using thymidylate kinase (TMK) assays.
Methods and results
1. Arithmetical derivations of the systems biology of the tubercle bacillus growth curve
(a) Theoretical partition of the tubercle bacillus cell cycle
(b) Tubercle bacillus growth curve patterns and biometric algorithms
2. Geometric correlation of TMKmyc levels with the M.tb growth curve
We employed systems and chemical biology models to parallel "M.tb growth curves" theoretically with "TMKmyc curves" in four different scenarios. First, we assumed that the level of TMKmyc within an M.tb culture is equal to the total TMKmyc produced by all the functionally viable cellular units within all the colonies in vitro. This assumption was aimed at minimizing the impact of heterogeneity within the culture, which may contain an admixture (quasi) of dividing and non-dividing M.tb. Basically, as long as overall proliferation occurs, the cumulative levels of TMK would rise in proportion. Second, assuming that DNA synthesis occurs 2, 4, 6 and 8 days after inoculation (and that these are the equivalent time-lines for the G1 phase, tG1); then a further 58, 56, 54 and 52 days would be required for actual division of the tubercle bacilli to occur. Therefore, M.tb DNA synthesis may parallel the actual growth of tubercle bacilli but occur at a much earlier time. Hence the total secreted levels of enzymes such as TMKmyc synthesized during the cell cycle may parallel the "S" pattern of the microbial growth curve (and equally obey the quadratic equation); except that these changes are bound to occur much earlier.
3. Prediction of drug sensitivity profiles based on TMKmyc levels rather than observation of M.tb growth changes (colonies per unit)
This paper provides the first ever theoretical and modeling framework to support the view that assays of enzymes involved in the synthesis of DNA precursors may be applied to shorten the time-lines to positive identification of microbial cultures (and possibly of other cell lines including cancers). Specifically: if variations in TMKmyc during the S phase of the tubercle bacillus cell cycle are independent of the time taken by the tubercle in other stages of the cycle, then TMKmyc assays are a potential predictive biomarker of in vitro M.tb growth that may be applied to reducing the time-lines to positive identification of Mycobacterium tuberculosis (M.tb) cultures. There are several existing methods for detecting TB within cultures, but they mostly lack inexpensive platforms for routine use within resource-limited settings where the TB burden is highest . Moreover, many of these emerging TB assays  are not specific for detecting duplicative changes in M.tb DNA. We therefore felt it necessary to establish a highly specific M.tb DNA duplicative assay that may be mounted on inexpensive technology platforms that are currently widely used within resource-limited settings, such as lateral flow immunochromatography.
Overall, tmk-based assays are therefore likely to be cheap but still serve the purpose of accelerating diagnosis and drug resistance profiling of M.tb and possibly other infectious diseases in culture. First, we show here why variations in the durations of the G1, S and G2 phases of the M.tb cell cycle may account for a longer 'culture time' for M.tb (see Fig. 1 for illustration) than for other microbes. In principle, the generation time of any cell is determined by the time it spends in G1 and G2 . Some cells, however, are known to go into deeper states of latency or dormancy (designated G0 or even G00), only reverting to the G1 stage when a need for reproduction, proliferation or regeneration arises . Various environmental factors are known to initiate such deeper states of dormancy in microbes and eukaryotes alike including low temperature, absence of nutrients and low oxygenation states ; though the full range of factors that determine why tubercle bacilli go into related states of latency in vivo are yet to be fully established [1, 21]. It is therefore unclear whether these deeper states of dormancy explain the slowly-growing nature of M.tb (especially since host immune factors also seem to play a significant role in arresting M.tb proliferation in vivo ). What is evident from our theoretical models, though, is that experiments aimed at establishing the actual values of the durations tG1, ts and tG2 for various microbes including M.tb may yield significant knowledge about the rate-limiting steps in their "growth" activity when cultured. For instance, using our proposed TMKmyc antigen capture assays, together with radio-labeled phosphorylated (dTMP r ) nucleotide bases as metabolites, the time ts required for DNA duplication by the tubercle bacillus may be estimated.
Note that in equations I and II representing tmk variation in sensitive and resistant M.tb isolate scenarios respectively; the |tmk| levels within the test isolate culture are inversely correlated to the numerical value of the wayengera-joloba lag or "wj" window, also denoted "WJ".
The above Area Under Curve (AUC) method for drug resistance profiling offers a means of quantifying drug resistance. Such ability to quantify resistance is especially clinically significant within resource-limited settings where second-line TB chemotherapy options are limited and drugs with minimal evident resistance are useful in combination with other drugs to which the isolate being treated is not resistant. In other words, TMK assays would enable the derivation of "salvage" regimens from the limited available drug options. Supposing that tmk levels follow a quadratic curve as above, parallel to the "S" growth curve of microbes, then this AUC can be calculated as the integral (∫) of the equation: |tmk| = at2+ bt + c or simply 1/3 at3+1/2 bt2+ct. From this equation, it can also be derived that the drug sensitivity of tubercle bacilli is an inverse factor of the AUC while existing resistance profiles are directly correlated to the same Area. However, in order to compare the prevailing resistance profiles between two isolates, one would need first to define the time over which to make the above calculations (say t1-t2); although we recommend that future universal algorithms be based on the AUC between t = 0 and the t-value when the tmk levels first reach the stationery stage. Note also that TMKmyc assays have the advantage of being applicable to all drug scenarios without the need for modifications such as those required for colorimetric redox-indicator methods or nitrate reductase assays, which have been found to be highly sensitive only for rifampicin and isoniazid resistance testing . Moreover, targeted detection of other related microbial enzymes such as HIV thymidine kinase may enable drug resistance profiling for HAART to be conducted within resource-limited settings, where inexpensive platforms for phenotyping or genotyping for drug sensitivity have remained mostly lacking . In the case of the M.tb scenario, whether or not the proposed TMKmyc-based drug profiling assay would yield consistent results relative to several existing rapid drug resistance tests such as phage amplification, or line probe assays such as INNO-LipA Rif.TB(LiPA) or Genotype MTBDR, is a subject that requires future comparative evaluations once TMKmyc assays are in practical use .
Several limitations are evident in our theoretical model, which must be dealt with prior to the actualization of TMKmyc assays as a predictive biomarker for in vitro tubercle bacillus growth. First, it is currently unclear whether secreted TMKmyc is present in the culture medium at levels detectable by MAb(s). Specifically, existing data provide contradictory views about the possible secretory nature of TMKmyc. Munier-Lehmman et al.  have previously shown that TMKmt forms over 30% of all colony filtrate proteins in a strain of E. coli genetically engineered to express TMKmt. However, other studies on host sera using in vivo-induced antigen technology (IVIAT) [23–25] and crossed immunoelectrophoresis (CIE) , though identifying over 11 genes involved in M.tb metabolism, do not categorically list TMKmt among them. Note, however, that the absence to date of methods for specific monitoring of TMKmyc (such as the one we propose) among the latter studies [[23–25] and ] relative to work by Munier-Lehmann and colleagues  may possibly explain these discrepancies. Second, it is assumed that TB cultures constitute a homogenous mixture of actively proliferating cells. Within wild type M.tb cultures, however, this may not be the case as active and dormant tubercle may co-exist. Such heterogeneity of growth among the M.tb in culture may also affect TMKmyc secretory levels. We hold that since qualitative rather than quantitative measures of TMKmyc are required to determine M.tb culture positivity, TMKmyc-based assays are still valid for the purpose of reducing time-lines to the reading of TB cultures. What may be affected is drug resistance profiling, since this requires measurement of the levels of TMK in the culture. Therefore, for purposes of drug resistance profiling, more algorithms may need to be integrated into the TMKmyc assays to render them more standardized and reliable in the face of heterogeneity of M.tb proliferation in vitro. Third, in the absence of in vitro evidence to support the view that TMKmyc may be adequately and specifically detected by monoclonal antibodies, the proposed assays for TMKmyc levels in culture are still speculative. Whether TMKmyc-based assays will perform better than, worse than or the same as existing TB tests  therefore remains elusive and requires future comparative clinical trials. More work needs to be done to affirm the binding affinity and specificity of the antigen capture assays for TMKmt. This work is, however, complicated by the findings by Steingart and colleagues, who recently reported inconsistencies in most commercially available TB serodiagnostics . Fourth, although it is clear that DNA duplication precedes actual microbial division , the duration of this process among various microbes including M.tb is largely unclear. Studies such as that by Sherley and Kelly  are therefore needed for M.tb and other infectious pathogens alike. Overall, if the above-listed potential shortcomings are overcome, the proposed antibody-based TMKmyc assays provide the flexibility to be mounted on a cheap technology platform for use within resource-limited settings. For instance, incorporation of "chemiluminiscent or fluorescent" strategies may enable "photo-detection of TMKmyc changes in culture" and thereby automation of the entire assay. Specifically, if the initiation of light emission is tailored to antigen (TMKmyc) capture, then chemiluminiscent or fluorescent TMKmyc-specific antibodies (CAbs and FAbs respectively) may be integrated within the LJ medium so that visual observation of color change or photometric measurements may be used to detect the levels of TMKmyc secreted. This simple automation could remove the need to conduct laborious quantitative ELISA or antigen capture assays to monitor TMKmyc levels, and thereby allow use by persons without extensive specialized training. Lastly, as TMKmyc is a specific antigen of M.tb , targeted detection of TMKmyc in sputum may provide an alternative strategy for diagnosing symptomatic TB to that based on staining for AFBs.
TMKmyc assays may be applicable to reducing time-lines to determining positivity of Mycobacterium tuberculosis (M.tb). Moreover, TMKmyc assays may also serve to hasten drug resistance profiling in vitro. The inherent flexibility of the proposed antigen capture assay allows the potential incorporation of simple technologies such as chemiluminiscent or fluorescent strategies and could enable easy-to-use platforms to be built such as those based on "photo-detection of TMKmyc changes in culture", as well as automation of the entire assay.
I thank Dr Joloba L. Moses at the Division of Molecular Biology, Dept of Medical Microbiology, College of Health Sciences, Makerere University, for helpful insights. Dr Joloba L. Moses, Dr Okwera Alfonse (Director of the NTLP Treatment Reference Center, Old Mulago) and Prof Byarugaba Wilson (Director, Postgraduate Research, Kampala International University Western Campus) are acknowledged co-investigators in on-going proof of concept experiments. A Sida/Sarec Faculty small grant (College of Health Sciences, Makerere University) was received towards this work.
- Snider DEJ, Raviglione M, Kochi A: Tuberculosis: Pathogenesis, Protection, and Control. Edited by: Bloom BR. 1996, Washington DC: Am. Soc. Microbiol, 2-11.Google Scholar
- Barry CE, Mdluli K: Drug sensitivity and environmental adaptation of mycobacterial cell wall components. Trends Microbiol. 1996, 4: 275-281. 10.1016/0966-842X(96)10031-7.View ArticlePubMedGoogle Scholar
- World Health Organization: Preventive therapy against tuberculosis in people living with HIV. Wkly Epidemiol Rec. 1999, 74: 385-398.Google Scholar
- Wheeler PR, Ratledge C: Tuberculosis: Pathogenesis, Protection, and Control. Edited by: Bloom BR. 1994, Washington DC; Am. Soc. Microbiol, 353-385.Google Scholar
- Tanaka G, Shojima J, Matsushita I, Nagai H, Kurashima A, Nakata K, Toyota E, Kobayashi N, Kudo K, Keicho N: Pulmonary Mycobacterium avium complex infection: association with NRAMP1 polymorphisms. Eur Respir. 2007, 30: 90-96. 10.1183/09031936.00042506.View ArticleGoogle Scholar
- Prince DS, Peterson DD, Steiner RM, Gottlieb JE, Scott R, Israel HL, Figueroa WG, Fish JE: Infection with Mycobacterium avium complex in patients without predisposing conditions. N Engl J Med. 1989, 321 (13): 863-868.View ArticlePubMedGoogle Scholar
- Andrews J, Basu S, Scales D, Maru DSR, Subbaraman R: XDR-TB in South Africa: Theory and Practice. PLoS Med. 2007, 4 (4): e163-10.1371/journal.pmed.0040163.PubMed CentralView ArticlePubMedGoogle Scholar
- Whalen CC: Diagnosis of latent tuberculosis infection: measure for measure. JAMA. 2005, 293: 2785-2787. 10.1001/jama.293.22.2785.View ArticlePubMedGoogle Scholar
- Pai M, Ramsay A, O'Brien R: Evidence-Based Tuberculosis Diagnosis. PLoS Med. 2008, 5 (7): e156-10.1371/journal.pmed.0050156.PubMed CentralView ArticlePubMedGoogle Scholar
- Weldingh K, Andersen P: ESAT-6/CFP10 Skin Test Predicts Disease in M. tuberculosis-Infected Guinea Pigs. PLoS ONE. 2008, 3 (4): e1978-10.1371/journal.pone.0001978.PubMed CentralView ArticlePubMedGoogle Scholar
- Morgan M, Kalantri S, Flores L, Pai M: A commercial line probe assay for rapid detection of rifampicin resistance in Mycobacterium tuberculosis: A systemic review and meta-analysis. BMC. 2005, 5: 62-10.1186/1471-2210-5-S1-P62.View ArticleGoogle Scholar
- Ling DI, Zwerking A, Pai M: Genotype MTBDR assays for diagnosis of multi-drug resistant tuberculosis: A meta-analysis. Eur Respir J. 2008, 32 (5): 1165-1174. 10.1183/09031936.00061808.View ArticlePubMedGoogle Scholar
- Pai M, Kalantri S, Pascopella L, Riley LW, Reingold AL: Bacteriophage based assays for rapid detection of rifampicin resistance in mycobacteria tuberculosis: a meta-analysis. J Infect. 2005, 51: 175-187. 10.1016/j.jinf.2005.05.017.View ArticlePubMedGoogle Scholar
- Martin A, Portaels F, Palomino JC: Colorimetric redox-reductase methods for rapi detection of multidrug resistance in Mycobacteria tuberculosis: a metanalysis and meta-regression. PLoS ONE. 2008, 3: e1536-10.1371/journal.pone.0001536.View ArticleGoogle Scholar
- Flores LL, Pai M, Colford JM, Riley LW: In-house nucleic acid amplification tests for the detection of Mycobacterium tuberculosis in sputum specimens: meta-analysis and Meta-regression. BMC Microbiol. 2005, 5: 55-10.1186/1471-2180-5-55.PubMed CentralView ArticlePubMedGoogle Scholar
- Steingart KR, Henry M, Laal S, Hopewell PC, Ramsay A, Menzies D, Cunningham J, Weldingh K, Pai M: Commercial serological antibody detection tests for the diagnosis of pulmonary tuberculosis: a systematic review. PLoS Med. 2007, 4 (6): e202-10.1371/journal.pmed.0040202.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper NG, Berg P: The Eukaryotic Cell cycle-Understanding inheritance. The Human Genome Project: Deciphering the Blueprint of heredity. 1994, University Science Books, 9-1Google Scholar
- Munier-lehmann H, Chaffotte A, Pochet S, Labesse G: Thymidylate kinase of Mycobacterium tuberculosis: A chimera sharing properties common to eukaryotic andbacterial enzymes. Protein Science. 2001, 10: 1195-1205. 10.1110/ps.45701.PubMed CentralView ArticlePubMedGoogle Scholar
- Sherley JLTJ: Regulation of human thymidine kinase during the cell cycle. J Biol Chem. 1988, 263 (17): 8350-8358.PubMedGoogle Scholar
- Moors MA, Portnoy DA: Identification of bacterial genes that contribute to survival and growth in an intracellular environment. Trends Microbiol. 1995, 3: 83-85. 10.1016/S0966-842X(00)88885-X.View ArticlePubMedGoogle Scholar
- Zahrt TC: Molecular mechanisms regulating persistent Mycobacterium tuberculosis infection. Microbes and Infection. 2003, 5 (2): 159-167. 10.1016/S1286-4579(02)00083-7.View ArticlePubMedGoogle Scholar
- Wayengera M, Kajumbula H, Byarugaba W: Harnessing pharmacogenomics to tackle resistance to the "Nucleoside Reverse Transcriptase Inhibitor" backbone of highly active antiretroviral therapy in resource limited settings. Open AIDS Journal. 2008, 2: 78-81. 10.2174/1874613600802010078.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang ZD, Li ZH, DU BP, Jia HY, Liu ZQ, Chen X, Huang HR, Xing AY, Gu SX, Ma Y: [Screening and analysis of in vivo induced genes of Mycobacterium tuberculosis]. Zhonghua Yi Xue Za Zhi. 2008, 88 (3): 189-193.PubMedGoogle Scholar
- Deb DK, Dahiya P, Srivastava KK, Srivastava R, Srivastava BS: Selective identification of new therapeutic targets of Mycobacterium tuberculosis by IVIAT approach. Tuberculosis (Edinb). 2002, 82 (4–5): 175-182. 10.1054/tube.2002.0337.View ArticleGoogle Scholar
- Singh KK, Zhang X, Patibandla AS, Chien P, Laal S: Antigens of Mycobacterium tuberculosis expressed during preclinical tuberculosis: serological immunodominance of proteins with repetitive amino acid sequences. Infect Immun. 2001, 9 (6): 4185-4191. 10.1128/IAI.69.6.4185-4191.2001.View ArticleGoogle Scholar
- Wiker HG, Harboe M, Bennedsen J, Closs O: The antigens of Mycobacterium tuberculosis, H37Rv, studied by crossed immunoelectrophoresis. Comparison with a reference system for Mycobacterium bovis, BCG. Scand J Immunol. 1988, 27 (2): 223-239. 10.1111/j.1365-3083.1988.tb02342.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.