We thank Giovannetti and colleagues for their thoughtful comments about our recent article (1) entitled “Associations between Single-Nucleotide Polymorphisms in the PI3K-PTEN-AKT-mTOR Pathway and Increased Risk of Brain Metastasis in Patients with Non–Small Cell Lung Cancer.” In their letter, the authors express that the sample size should be determined according to the “rule of tens.” This is a simple and effective method used to evaluate prior events. In fact, for prospective studies, the method used commonly is based on the sample size estimation formula of survival analysis [|$(z_\beta + z_{1 - \alpha } )^2 /(P_A P_B \log \,_e^2 \Delta )$|; refs. 2, 3]. Here, P is the proportion of each genotype in our study. According to this formula, the sample size is associated with the gene mutation frequency of the SNP. On the basis of our study and public data, there should be approximately 80 events for analysis of rs2498804 and rs2494732, and we have 99 patients with brain metastases. The rule of tens is more about the ability to control for other variables in a multivariate analysis. We agree that the data may have been too sparse for the control variables we included. However, the unadjusted hazard ratio (HR) is similar (Table 1).
Polymorphisms and genotypes . | No. of patients (all) . | No. of events (%) . | HR (95% CI) . | P . | No. of patients without BM at diagnosis . | No. of events (%) . | HR (95% CI) . | P . |
---|---|---|---|---|---|---|---|---|
AKT1: rs2498804 | ||||||||
TT | 125 | 29 (23) | 1.000 | 122 | 26 (21) | 1.000 | ||
GT+GG | 192 | 70 (37) | 1.720 (1.115–2.651) | 0.014 | 181 | 59 (33) | 1.640 (1.034–2.602) | 0.036 |
AKT1: rs2494732 | ||||||||
CC | 156 | 37 (24) | 1.000 | 151 | 32 (21) | 1.000 | ||
CT+TT | 159 | 62 (39) | 1.783 (1.187–2.680) | 0.005 | 150 | 53 (35) | 1.786 (1.152–2.770) | 0.010 |
PIK3CA: rs2699887 | ||||||||
GG | 272 | 79 (29) | 1.000 | 263 | 70 (27) | 1.000 | ||
AG+AA | 45 | 20 (44) | 1.701 (1.041–2.778) | 0.034 | 40 | 15 (37) | 1.466 (0.839–2.560) | 0.179 |
Polymorphisms and genotypes . | No. of patients (all) . | No. of events (%) . | HR (95% CI) . | P . | No. of patients without BM at diagnosis . | No. of events (%) . | HR (95% CI) . | P . |
---|---|---|---|---|---|---|---|---|
AKT1: rs2498804 | ||||||||
TT | 125 | 29 (23) | 1.000 | 122 | 26 (21) | 1.000 | ||
GT+GG | 192 | 70 (37) | 1.720 (1.115–2.651) | 0.014 | 181 | 59 (33) | 1.640 (1.034–2.602) | 0.036 |
AKT1: rs2494732 | ||||||||
CC | 156 | 37 (24) | 1.000 | 151 | 32 (21) | 1.000 | ||
CT+TT | 159 | 62 (39) | 1.783 (1.187–2.680) | 0.005 | 150 | 53 (35) | 1.786 (1.152–2.770) | 0.010 |
PIK3CA: rs2699887 | ||||||||
GG | 272 | 79 (29) | 1.000 | 263 | 70 (27) | 1.000 | ||
AG+AA | 45 | 20 (44) | 1.701 (1.041–2.778) | 0.034 | 40 | 15 (37) | 1.466 (0.839–2.560) | 0.179 |
Abbreviations: BM, brain metastases; CI, confidence interval.
Another issue is whether the Bonferroni correction should be used. It is important to correct the P value of each gene when performing a statistical test on a group of genes. GeneSpring offers four types of multiple testing corrections, and the Bonferroni correction is the most stringent. The more stringent a multiple testing correction, the less false-positive genes are allowed. The trade-off of a stringent multiple testing correction is that the rate of false negatives is very high. To reduce false negatives, P values were not corrected. This may be debatable but is done in practice (3, 4). Furthermore, the association remained significant after the Bonferroni correction for AKT1:rs2494732 (i.e., <0.05/16 = 0.003) in our study. Otherwise, AKT1:rs2498804 and rs2494732 were in strong linkage disequilibrium (D = 0.95), and we believe that AKT1: rs2498804 is not a false-positive gene (1). In addition, we are pleased to find that the association of AKT1:rs2498804 was confirmed in your independent cohorts. We think that this is a very good method to avoid false positives.
Giovannetti and colleagues also raised concerns about the treatment details in our study. We agree with them on this issue. In fact, we have described that the treatment details included surgery, chemotherapy, and radiotherapy. Unfortunately, we have no information on EGFR mutation status or ALK translocations, so we are unable to address this important comment. However, we can report early findings from an ongoing study of 70 patients receiving oral erlotinib in a program supported by the China Charity Federation. The results were the same as those in our study (Table 2). The data have been provided in our response to the reviewers.
Characteristic . | No. of patients (all) . | No. of events (%) . | Univariate analysis HR (95% CI) . | P . | Multivariate analysis HR (95% CI) . | P . | Pa . |
---|---|---|---|---|---|---|---|
AKT1: rs2494804 | |||||||
TT | 27 | 2 (8) | 1.000 | 1.000 | |||
GT+GG | 43 | 15 (35) | 5.228 (1.195–22.873) | 0.028 | 5.721 (1.264–25.896) | 0.024 | 0.010 |
AKT1: rs2494732 | |||||||
CC | 33 | 2 (6) | 1.000 | 1.000 | |||
CT+TT | 35 | 15 (43) | 8.345 (1.907–36.515) | 0.005 | 12.261 (2.649–56.753) | 0.001 | 0.001 |
PIK3CA: rs2699887 | |||||||
GG | 58 | 14 (24) | 1.000 | 1.000 | |||
AG+AA | 12 | 3 (25) | 1.089 (0.313–3.792) | 0.893 | 2.492 (0.655–9.482) | 0.181 | 1.000 |
Characteristic . | No. of patients (all) . | No. of events (%) . | Univariate analysis HR (95% CI) . | P . | Multivariate analysis HR (95% CI) . | P . | Pa . |
---|---|---|---|---|---|---|---|
AKT1: rs2494804 | |||||||
TT | 27 | 2 (8) | 1.000 | 1.000 | |||
GT+GG | 43 | 15 (35) | 5.228 (1.195–22.873) | 0.028 | 5.721 (1.264–25.896) | 0.024 | 0.010 |
AKT1: rs2494732 | |||||||
CC | 33 | 2 (6) | 1.000 | 1.000 | |||
CT+TT | 35 | 15 (43) | 8.345 (1.907–36.515) | 0.005 | 12.261 (2.649–56.753) | 0.001 | 0.001 |
PIK3CA: rs2699887 | |||||||
GG | 58 | 14 (24) | 1.000 | 1.000 | |||
AG+AA | 12 | 3 (25) | 1.089 (0.313–3.792) | 0.893 | 2.492 (0.655–9.482) | 0.181 | 1.000 |
NOTE. Multivariate analyses in this table were adjusted for sex, patient age, disease stage, tumor histology, KPS, and smoking status.
Abbreviations: BM, brain metastases; CI, confidence interval.
aP values were calculated by the Fisher exact test.
Finally, we admit that there are limitations with our study and agree with Giovannetti and colleagues that further research is needed to unravel the functional significance of these polymorphisms. We are currently using immunohistochemistry to detect AKT protein expression in NSCLC primary tissues, and additional functional validation is ongoing.
See the original Letter to the Editor, p. 3623
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Grant Support
This study was funded by three grants from the National Natural Science Foundation of China (grants 81071832, 81272492, and 81101691).