
ORIGINAL ARTICLE 



Year : 2005  Volume
: 42
 Issue : 2  Page : 7077 

Variations in clinical estimates of tumor volume regression parameters and time factor during external radiotherapy in cancer cervix: Does it mimic the linearquadratic model of cell survival?
NR Datta, A Rajkumar, R Basu
Department of Radiotherapy, Sanjay Gandhi Postgraduate Institute of Medical Sciences, Lucknow, India
Correspondence Address: N R Datta Department of Radiotherapy, Sanjay Gandhi Postgraduate Institute of Medical Sciences, Lucknow India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/0019509X.16695
BACKGROUND: Tumor regression parameters and time factor during external radiotherapy (EXTRT) are of paramount importance. AIMS: To quantify the parameters of tumor regression and time factor during EXTRT in cancer cervix. SETTINGS AND DESIGN: Patients, treated solely with radiotherapy and enrolled for other prospective studies having weekly tumor regressions recorded were considered. MATERIALS AND METHODS: Seventyseven patients received 50Gy of EXTRT followed by intracavitary brachytherapy. Locoregional regressions were assessed clinically and regression fraction (RF) was represented as RF = c + a_{1}D + a_{2}D^{2} a_{3}T, with c, D and T as constant, cumulative EXTRT dose and treatment time respectively. STATISTICAL ANALYSIS USED: Step wise linear regression was performed for RF. Scatter plots were fitted using linearquadratic fit. RESULTS: Coefficients of parameters D, D^{2} and T were computed for various dose intervals, namely 020 Gy, 030 Gy, 040 Gy and 050 Gy. At 020 Gy and 030 Gy, only the coefficient of D2 was significant (P < 0.001), while both D^{2} and T turned significant (P < 0.001) at 040 Gy. For the entire range of 050 Gy, all the coefficients of D, D^{2} and T showed significance, leading to an estimate of 26 Gy for a_{1}/a_{2} and 0.96 Gy/day for a_{3}/a_{1}. CONCLUSIONS: As with a/β and g/a of postirradiation cell survival curves, a_{1}/a_{2} and a_{3}/a_{1} represents the cumulative effect of various radiobiological factors influencing clinical regression of tumor during the course of EXTRT. The dynamic changes in the coefficients of D, D^{2}sub and T, indicate their relative importance during various phases of EXTRT.
Keywords: Cervical cancer, Radiotherapy, Linearquadratic, Alpha/beta, Time factor
How to cite this article: Datta N R, Rajkumar A, Basu R. Variations in clinical estimates of tumor volume regression parameters and time factor during external radiotherapy in cancer cervix: Does it mimic the linearquadratic model of cell survival?. Indian J Cancer 2005;42:707 
How to cite this URL: Datta N R, Rajkumar A, Basu R. Variations in clinical estimates of tumor volume regression parameters and time factor during external radiotherapy in cancer cervix: Does it mimic the linearquadratic model of cell survival?. Indian J Cancer [serial online] 2005 [cited 2020 Jun 4];42:707. Available from: http://www.indianjcancer.com/text.asp?2005/42/2/70/16695 
» Introduction   
Linearquadratic (LQ) model has provided a satisfactory mathematical description to the mechanism of radiation induced cell kill.[1] The ratio of linear (α) and quadratic (β) components of cell kill  α/β, have been used to design various therapeutic strategies. [1],[2],[3],[4] Values of α/β, especially for tumor, can be best regarded as an extension of experimentally derived in vivo or in vitro estimates. Clinical estimates of α/β are usually carried out through Fe plots between reciprocal of isoeffective dose (D) and dose per fraction (d) or logisticdose response model from fractionated radiotherapy where the coefficients of D and Dd, are considered to play a role similar to α and β.[5] Although the respective coefficients of D and Dd are not the true estimates of α and β but their ratios continues to be preserved as estimates of α/β.[6] However, in both the situations, groups of patients need to be irradiated to different dosefractionation schedule to obtain isoeffects for specific end points which could be unethical in most clinical situations.
Tumor regression during clinical fractionated radiotherapy can be regarded as a reflection of the conglomeration of the various radiobiological events. External radiotherapy (EXTRT) being one of the important prognosticators of long term locoregional disease control in cancer cervix, patterns of tumor regression during the course of EXTRT assumes special significance. [7],[8],[9]
The present study has been carried out to estimate the tumor regression parameters and time factor from doseresponse relationship in patients of cancer cervix undergoing EXTRT to interpret the curvilinear LQ curves as obtained from in vivo and in vitro studies.
» Materials and methods   
Seventyseven patients of carcinoma of the cervix were selected from the database of cervical cancer treated at our institution from January 1996 to June 2001. Ethical clearance for this retrospective study was not required as the patients were subjected to a standard treatment as per the Department's policy, and the data was obtained from the review of their treatment records. All patient records maintained in our departmental data base were sorted based on the diagnosis and the years of the study. The data was retrieved by all the authors, individual files reviewed, parameters entered and taken up for systematic analysis. Out of the database, only those patients of cancer cervix in whom the treatment was delivered as per the timedosefractionation schedule of teletherapy stated below and those in which the tumor sizes and their weekly assessment were documented were included for this analysis.
Before the start of therapy, patients were subjected to routine clinical, hematological, biochemical and radiological investigations and staged clinically as per FIGO staging guidelines. Only those patients treated with radical radiotherapy to an EXTRT dose of 50 Gy at 2 Gy per fraction with 5 fractions per week followed by intracavitary brachytherapy were considered for this analysis. Teletherapy was delivered by either telecobalt (Theratron 780C, AECL, Canada) or 6 MV / 10MV photons (ML 20DX, Mitsubishi Corporation, Japan). Following 2 weeks of EXTRT, all these patients received high dose rate intracavitary brachytherapy of 18 Gy in 3 fractions at weekly intervals to point A (Microselectron HDR, Nucletron, The Netherlands or Ralstron 20B, Shimadzu Corporation, Japan). In all patients, periodic tumor responses at corresponding doses were available during the course of EXTRT. The details of treatment and response assessments have been described in a recent publication.[8] Locoregional response at the primary and parametria were recorded as percentage regression from the pretreatment size. These were based on clinical evaluation namely, per speculum, per vaginum and rectovaginal examinations carried out during weekly monitoring of patients during the course of treatment. Even though, this forms a retrospective collection of cases, these patients had been part of other prospective studies, where weekly recording of tumor regression was mandatory. Radiological assessment of response using contrast enhanced computed tomography (CECT) or magnetic resonance imaging (MRI) were not a part of any of the prospective studies, as subjecting patients to weekly imaging studies during the course of EXTRT was not practically feasible. Thus, even though imaging assessment using CECT or MRI could have refined the quantum of tumor regression, as this was not practically possible for all patients due to waiting times on CT / MRI machine and cost implications, clinical regression was uniformly advocated for all cases.
The residual gross tumor at various dose intervals was expressed in terms of residual fraction (RF) of the pretreatment tumor. Thus, starting with a RF of 1, at commencement of radiation, RF was recorded at subsequent EXTRT doses and expressed as fraction of the residual disease. The observed RF was plotted along the Yaxis against the corresponding cumulative EXTRT doses in Xaxis. The scatter plots were fitted using quadratic fit.
Statistical analysis
Step wise forward linear regression for RF with EXTRT dose, was carried out as per the model, RF = c + a_{1}nd + a_{2}nd^{2}  a_{3}T, with c as constant; n, number of fractions; d, dose per fraction and T, the corresponding treatment time. Variables with probability of F equal to 0.05 were entered into the model and removed at 0.10. Only variables found significant at p < 0.05 were considered for further estimations. However, in view of the poor fit with the above expression, regression was repeated using dose D, (D = nd), D^{2} and T as per the model, RF = c + a_{1}D + a_{2}D^{2}  a_{3}T. As mentioned earlier, the coefficients a_{1}, a_{2} and a_{3} were assumed to play a role similar to coefficients α, β and γ of the LQ doseeffect relationship. The coefficients  a_{1}, a_{2} and a_{3} were therefore not the actual estimates of α, β and γ respectively, but the ratios a_{1}/a_{2} and a_{3}/a_{1} could be considered to represent α/β and time factor (γ/α) respectively.[6]
The cumulative EXTRT dose range of 50 Gy was divided into different segments of increments of 10 Gy each i.e., 010 Gy, 020 Gy, 030 Gy, 040 Gy and 050 Gy over which the cumulative RF was evaluated. Linear regression of RF against the corresponding cumulative EXTRT dose and treatment time were carried out to examine if there were changes in the coefficients of the total dose D, D^{2} and T at these dose intervals. The regression model r^{2}, its significance and the 95% confidence intervals of the various coefficients were computed. Only those models with r^{2} of at least 0.5 and only those coefficients which were significant were considered for further analysis and interpretation which included estimation of a_{1}/a_{2} and a_{3}/a_{1} ratios. All the statistical calculations were performed using SPSS software package for windows, version 9.0 (SPSS Inc. Chicago, IL, USA).
» Results   
The detail demographic profiles of the patients are summarized in [Table  1] with majority of patients in stages II and III. All patients had received the planned dose of 50 Gy of EXTRT followed by 18 Gy of intracavitary brachytherapy. The overall treatment time for EXTRT ranged from 33 to 52 days (median: 37 days). At the completion of teletherapy, 40 patients (51.9%) had no gross primary locoregional tumor and following one month of completion of intracavitary brachytherapy, 53 (68.8%) patients had complete locoregional clearance of their pelvic disease.
The response fraction and corresponding EXTRT doseresponse curve showed a typical curvilinear plot on a linear scale and was fitted using a quadratic fit (R^{2} = 0.838) [Figure  1]a). On plotting the same on a loglinear scale, the shape of the curve had a strong resemblance to the classic curvilinear plots of the cell survival curves. This included an initial slope followed by shoulder and a terminal straightline (R^{2} = 0.739) [Figure  1]b).
From the entire 77 patients who had received 50 Gy of teletherapy, a step wise forward linear regression for RF was computed initially using "nd", "nd2" and T using the expression RF = c + a_{1}(nd) + a_{2}(nd^{2})  a_{3}T. This failed to provide significant coefficients for "nd", thereby any further estimates for a_{1}/a_{2} and a_{3}/a_{1}, could not be carried out [Table  2]. Substituting D^{2} in place of nd^{2} in the expression RF = c + a_{1}D + a_{2}D^{2}  a_{3}T, estimations were repeated and this yielded coefficients which were significant for D (P = 0.019), D^{2} (P < 0.001) and T (P = 0.029). Thus, a_{1}/a_{2} and a_{3}/a_{1} were estimated at 26 Gy and 0.96 Gy/day respectively [Table  2].
The cumulative dose range of 50 Gy were divided into different segments of increments of 10 Gy each i.e., 010 Gy, 020 Gy, 030 Gy, 040 Gy and 050 Gy. The regression models for RF at these dose intervals were considered only if the model R^{2} was at least 0.5. As stated, the regression coefficients for D, D^{2} and T at each of these intervals were considered only if these were significant at p < 0.05. As evident from [Table  3], the model R^{2} for dose interval 010 Gy was unsatisfactory at 0.258, while the models for the remaining dose intervals were beyond 0.5. It was apparent from the models for different dose levels that D^{2} had an impact for all dose levels beyond 10 Gy (P < 0.001). The EXTRT treatment time influenced RF for dose level beyond 30 Gy (P < 0.001 for dose interval 040 Gy) while for the entire range of 50 Gy, all the coefficients of D, D^{2} and T showed significant influence (P = 0.019 for D, P < 0.001 for D^{2} and P = 0.029 for T) [Table  3], [Figure  2].
» Discussion   
The LQ model, based on in vivo and in vitro studies has inherent limitations of its use for doses beyond 8  10Gy due to the continuous bending shape of the cell survival curve at higher doses. It is assumed that the effect on cell survival for each of the multiple fractions would be identical and thus LQ model could be regarded as a low dose approximation to equations that become straight exponentials at higher doses.[1] Estimates of α/β are the uncertain parameters used in LQ model and the application of applying LQ model for multifractional studies, is based on the assumption that every equal radiation fraction will have the same biological effect on the target, independent of its position in the series.[1],[5] This led to the expression for total effect, E represented with time factor as E = αnd + βnd^{2}  γT.
The observation that clinical RF could not be satisfactorily fitted using the parameters "nd" and "nd^{2}" indicates that for fractionated clinical radiotherapy, the effect of each fraction of dose on the RF may not be equal. This is in contrast to the basic assumption of the LQ model for fractionated radiotherapy. On the contrary, it appeared that RF was a second order function of the total dose delivered and for multifractionation schedules, it could be best represented in the form, RF = c + αD + βD^{2}  γT.
This brings to the fundamental issue on the kinetics of the repair of damage caused by each fraction. Are these completely repaired before the next radiation dose delivered after 24 hours or are the cells likely to harbor residual damage at the time of each subsequent fraction of radiation? The repair of sublethal damage after each fraction is usually considered to be monoexponential with a specific halftime (T½) of around 0.2 to 4 hours, or to consist of a fast and a slow component with T½ of around 0.150.3 hours and 23 hours respectively.[5],[6] This would ensure that most of the repair gets almost completed before the next dose of radiation, say after 24 hours. Thus based on these two basic assumptions, the cell survival curve for multiple fractions has been considered to have a negative exponential straight line with a slope similar to that of the initial slope, a, caused by a dose per fraction, d, rather than the typical curvilinear plot on a loglinear scale.[1]
Over the years, the concept of a monoexponential repair having a single exponential half time has been debated.[10] The monoexponential repair concept of a near complete repair for interfraction intervals around 6 to 8 hours has been debated.[11],[12] Fowler had proposed that the repair kinetics could be a second order whereby at the first half time τ, the time constant for the unrepaired proportion would fall to ½, but at 2 τ it would fall to 1/3, 3 τ to ¼, at 4 τ to 1/5 and so on. This would indicate that in contrast to the repair being monoexponential, considerable unrepaired damage exists even at 24 hours, which may get compounded with each subsequent radiotherapy fractions.
Based on the above derivations of incomplete repair of radiation induced damage during interfraction period, the assumption of equal effect for identical fractions of radiation dose in multifractionated radiotherapy may be questionable. This would alter the biological effect E, proposed by the LQ model from the expression, E = n(αd + βd^{2}) to E = αD + βD^{2}. Incorporating the treatment time, T, E would be better represented as E= αD + βD^{2}  γT in clinical situations. The observations from this study also demonstrate the validity of the latter expression [Table  2] and [Table  3].
Variation in coefficients of D, D^{2} and T during the various cumulative EXTRT doses indicates their relative dominance during the course of EXTRT [Table  3]. The results of our study could be interpreted as an analogy to the cell survival curves and the RF curve could be considered basically in 3 parts  an initial straight exponential fall, followed by a shoulder and then a straight line exponential fall [[Figure  2]]. At low doses (010 Gy), the linear component of dose would be expected to dominate and cause lethal damage to the DNA, but the results of our study has not been able to show this with certainty [Figure  3]a). This could be due to the difficulty in quantification of the early tumor regression at doses till 10 Gy. However, with a cumulative dose of 020 Gy and beyond, cell kill was predominantly a function of quadratic component of cumulative EXTRT dose. This could be represented as lethal DNA damage due to two separate hits close to each other [Figure  3]b). With time, as in clinical fractionated EXTRT, some of these could get repaired [Figure  3]c), leading even to complete repair of the damage to the DNA [Figure  3]d). Thus the impact of treatment time starts to show at around 3rd week of radiotherapy. By the 5th week, the terminal part of the RF curve tends to become straight, possibly indicating the resurgence of the linear component of cell kill with a difference. In this situation some of the residual sublethal or potentially lethal damages caused by either of the two hits of quadratic component would need only a single hit to be inactivated [Figure  3]e). A similar change in time factor related to the tumor repopulation occurring after 2 weeks of EXTRT in headandneck cancer has been recently reported.[13] The impact of treatment time, are usually ascribed to repopulation, especially when the trials are conducted with different treatment times and other split course RT schedules. However, since in patients of this study, teletherapy was executed without any gap or changing time periods, repair kinetics could be a major contributory factor.
The population of shrinking tumors under the influence of reoxygenation, repopulation, tumor and tissue debris etc. are not expected to behave in a manner similar to cell survival cures from in vitro conditions. Thus, the clinical estimates of regression parameters and time factor for carcinoma of the cervix would not be expected to be identical to those of α/β and γ/α of that obtained through cell survival studies. The estimated a_{1}/a_{2} and a_{3}/a_{1} represent the complex radiobiological phenomena taking place during fractionated radiotherapy incorporating combined influence of cell loss, changing tumor kinetics, clonogen doubling, repair, reoxygenation, blood flow and clearance of the dead and necrotic tissue / debris and thus could be more appropriate in clinical radiotherapy.[1],[14],[15],[16],[17],[18],[19],[20],[21]
Radiobiological alterations during the course of fractioned radiotherapy are a complex interplay of the "4 Rs" of radiobiology and attempts to represent them in LQ equation would not only make the LQ model lose its innocence, but also introduce several parameters with uncertainty in their estimates.[14],[22],[23],[24] This study has tried to represent estimate α/β in cervical tumors from actual doseresponse relations. To the best of our knowledge, the only estimate of γ/α available in literature for cancer cervix has been reported to be more than 13.9 Gy.[2]
The estimates derived in this study could be applicable to cervical cancer treated with teletherapy using 2 Gy per fractions over a 5 week period with week end gaps. Extrapolating these to cervical tumours treated with different fractionation schedules may not be appropriate since the doseresponse curves could be different for different dose per fractions. The values obtained by this study indicate that it may not be always proper to calculate biological equivalent doses for various dose fractionation studies by taking a standard α/β value of 10 Gy and a time factor of 0.60 Gy/day in all cases. Inherent limitations of these to represent the complex radiobiological phenomena during fractionated radiotherapy should always be kept in mind. Thus, mathematical models used to represent complex biological phenomena should be based on careful clinical observation to guide treatment.[25]
» References   
1.  Fowler JF. The linearquadratic formula and progress in fractionated radiotherapy. Br J Radiol 1989;62:67994. [PUBMED] 
2.  Thames HD, Bentzen SM, Turesson I, Overgaard M, van den Bogaert W. Timedose factors in radiotherapy: a review of the human data. Radiother Oncol 1990;19:21935. [PUBMED] 
3.  Chougule A. Linear quadratic model compared with TDF and TSD concepts for malignancy of cervix uteri treated with twice daily fractionated schedules. Strahlenther Onkol 1994;170:305. [PUBMED] 
4.  Sood B, Garg M, Avadhani J, Gorla G, Malhotra H, Guha C, et al . Predictive value of linearquadratic model in the treatment of cervical cancer using highdoserate brachytherapy. Int J Radiat Oncol Biol Phys 2002;54:137787. [PUBMED] [FULLTEXT] 
5.  Joiner MC, van der Kogel AJ. The linearquadratic approach to fractionation and calculation of isoeffect relationships. In : Basic clinical radiobiology, 2nd edn. Steel GG, editor. London: Arnold; 1997. p. 10622. 
6.  Bentzen SM. Doseresponse relationships in radiotherapy. In : Basic clinical radiobiology, 2nd edn. Steel GG, editor. London: Arnold; 1997. p. 7886. 
7.  Hong JH, Chen MS, Lin FJ, Tang SG. Prognostic assessment of tumor regression after external irradiation for cervical cancer. Int J Radiat Oncol Biol Phys 1992;22:9137. [PUBMED] 
8.  Datta NR, Rajkumar A, Basu R. Tumor regression dynamics with external radiotherapy in cancer cervix and its implications. Ind J Cancer in press; 2004. 
9.  Chen SW, Liang JA, Yang SN, Ko HL, Lin FJ. The adverse effect of treatment prolongation in cervical cancer by highdoserate intracavitary brachytherapy. Radiother Oncol 2003;67:6976. [PUBMED] [FULLTEXT] 
10.  Dale RG, Fowler JF, Jones B. A new incompleterepair model based on a "reciprocaltime" pattern of sublethal damage repair. Acta Oncol 1999;38:91929. [PUBMED] [FULLTEXT] 
11.  Fowler JF. Clinical implications of hyperbolic (reciprocal straight line) repair of sublethal radiation damage. In : Progress in RadioOncology VI, Bolgna. Kogelink HD, Sedlmayer E, editors. Italy: Monduzzi Editore; 1998. p. 4719. 
12.  Fowler JF. Is repair of damage from ionizing radiation secondorder rather than first order? A simpler explanation of apparently multiexponential repair. Radiat Res 1999;152:12436. 
13.  Tarnawski R, Fowler J, Skladowski K, Swierniak A, Suwinski R, Maciejewski B, et al . How fast is repopulation of tumor cells during treatment gap? Int J Radiat Oncol Biol Phys 2002;54:22936. [PUBMED] [FULLTEXT] 
14.  Jones B, Dale RG. Mathematical models of tumor and normal tissue response. Acta Oncol 1999;38:8839. [PUBMED] [FULLTEXT] 
15.  Durand RE, AquinoParsons C. Predicting response to treatment in human cancers of the uterine cervix: sequential biopsies during external beam radiotherapy. Int J Radiat Oncol Biol Phys 2004;58:55560. [PUBMED] [FULLTEXT] 
16.  Travis EL, Tucker SL. Isoeffect models and fractionated radiation therapy. Int J Radiat Oncol Biol Phys 1987;13:2837. [PUBMED] 
17.  Zaider M, Wuu CS, Minerbo GN. The combined effects of sublethal damage repair, cellular repopulation and redistribution in the mitotic cycle. I. Survival probabilities after exposure to radiation. Radiat Res 1996;145:45766. 
18.  Hahnfeldt P, Hlatky L. Resensitization due to redistribution of cells in the phases of the cell cycle during arbitrary radiation protocols. Radiat Res 1996;145:13443. [PUBMED] 
19.  Wouters BG, Brown JM. Cells at intermediate oxygen levels can be more important than the "hypoxic fraction" in determining tumor response to fractionated radiotherapy. Radiat Res 1997;147:54150. [PUBMED] 
20.  Jones B, Dale RG. Cell loss factors and the linearquadratic model. Radiother Oncol 1995;37:1369. [PUBMED] [FULLTEXT] 
21.  Brenner DJ, Hlatky LR, Hahnfeldt PJ, Hall EJ, Sachs RK. A convenient extension of the linearquadratic model to include redistribution and reoxygenation. Int J Radiat Oncol Biol Phys 1995;32:37990. [PUBMED] [FULLTEXT] 
22.  Wein LM, Cohen JE, Wu JT. Dynamic optimization of a linearquadratic model with incomplete repair and volumedependent sensitivity and repopulation. Int J Radiat Oncol Biol Phys 2000;47:107383. [PUBMED] [FULLTEXT] 
23.  Bentzen SM, Ruifrok AC, Thames HD. Repair capacity and kinetics for human mucosa and epithelial tumors in head and neck: clinical data on the effect of changing the time interval between multiple fractions per day in radiotherapy. Radiother Oncol 1996;38:89101. [PUBMED] [FULLTEXT] 
24.  Lyng H, Sundfor K, Rofstad EK. Changes in tumor oxygen tension during radiotherapy of uterine cervical cancer: relationships to changes in vascular density, cell density and frequency of mitosis and apoptosis. Int J Radiat Oncol Biol Phys 2000;46:93546. [PUBMED] [FULLTEXT] 
25.  Coleman CN. International conference on translational research ICTR 2003 conference summary: marshalling resources in a complex time. Int J Radiat Oncol Biol Phys 2004;58:30719. [PUBMED] [FULLTEXT] 
Figures
[Figure  1], [Figure  2], [Figure  3] Tables
[Table  1], [Table  2], [Table  3]
This article has been cited by  1 
Optimal tumour control for earlystage nonsmallcell lung cancer: A radiobiological modelling perspective 

 Mohammed Alaswad,Christoph Kleefeld,Mark Foley   Physica Medica. 2019; 66: 55   [Pubmed]  [DOI]   2 
Overcoming radioresistance in WiDr cells with heavy ion irradiation and radiosensitization by 2deoxyglucose with photon irradiation 

 Felix Christian Hasse,Stefan Alexander Koerber,Elena Sophie Prigge,Jakob Liermann,Magnus von Knebel Doeberitz,Juergen Debus,Florian Sterzing   Clinical and Translational Radiation Oncology. 2019;   [Pubmed]  [DOI]   3 
The alfa and beta of tumours: a review of parameters of the linearquadratic model, derived from clinical radiotherapy studies 

 C. M. van Leeuwen,A. L. Oei,J. Crezee,A. Bel,N. A. P. Franken,L. J. A. Stalpers,H. P. Kok   Radiation Oncology. 2018; 13(1)   [Pubmed]  [DOI]   4 
Current Role of Magnetic Resonance Imaging in Evaluation and Radiotherapy in Locally Advanced Carcinoma Cervix 

 Saikat Das,Santanu Samanta,Manu Mathew,Reetu John,Abraham Peedicayil,Subhashini John   Indian Journal of Gynecologic Oncology. 2016; 14(2)   [Pubmed]  [DOI]  





