Purpose: Circulating tumor cells (CTC) might function as early markers for breast cancer metastasis or monitoring therapy efficacy. Enrichment and identification of CTCs are based on epithelial markers that might be modulated during epithelial–mesenchymal transition. Little is known about the expression of keratins in CTCs and whether all CTCs can be detected with antibodies directed against a limited panel of keratins.

Experimental Design: Protein expression of keratin 2, 4–10, 13–16, 18, and 19 were assessed by a cocktail of antibodies (C11, AE1, AE3, and K7) and keratin antibodies C11 and A45-B/B3 alone in 11 breast cancer cell lines and 50 primary breast carcinomas and their lymph node metastases. Furthermore, CTCs were assessed in blood of 70 metastatic breast cancer patients.

Results: Claudin-low cell lines did not show expression of normal breast epithelial keratins but were positive for K14 and K16, detected by the cocktail only. Primary breast carcinomas showed changes in keratin expression during metastatic progression to the lymph nodes. In 35 of 70 patients CTCs were identified, of which 83%, 40%, and 57% were identified by the cocktail, C11 and A45-B/B3, respectively. Identification of CTCs by the cocktail was associated with shorter survival (P < 0.01). In silico analyses revealed association between KRT16 expression and shorter relapse-free survival in metastatic breast cancer.

Conclusion: Breast cancer cells show a complex pattern of keratin expression with potential biologic relevance. Individual keratin antibodies recognizing only a limited set of keratins inherit the risk to miss biologically relevant CTCs in cancer patients, and antibody cocktails including these keratins are therefore recommended. Clin Cancer Res; 18(4); 993–1003. ©2012 AACR.

This article is featured in Highlights of This Issue, p. 915

Translational Relevance

Circulating tumor cells (CTC) found in blood might function as an early diagnostic marker for breast cancer metastasis and as a marker to monitor therapy efficacy. In this study, we show that keratins, the markers currently used for the identification of CTCs in breast cancer patients, can be downregulated during metastatic progression. The keratin antibodies used to date in CTC diagnostics (e.g., C11 in the U.S. Food and Drug Administration–approved CellSearch System) miss CTCs expressing certain keratin proteins (e.g., K16) associated with metastatic progression in breast cancer. The results presented here point to the use of broad-spectrum antibody cocktails including these biologically relevant keratins for CTC detection.

Tumor cells that are able to invade the body and enter the blood circulation may finally grow out to metastases in distant organs (1). Because cancer metastasis is the main cause of cancer-related death, early detection of (micro) metastases is of vital importance (2). Unfortunately, current imaging techniques are unable to detect micrometastases; therefore, different strategies have been explored to facilitate the detection of micrometastatic disease in early stage by identification of tumor cells in bone marrow or blood using immunocytochemistry (ICC; ref. 3). It has been shown that the identification of disseminated tumor cells (DTC) in the bone marrow of breast cancer patients covers prognostic information (4). In addition, identification of DTCs could help in decision making on the need of therapy and to predict the outcome of treatment (5). Alternatively to bone marrow samples, which are difficult to obtain in a routine fashion, blood has been shown to be a valuable source for the identification of dissociated tumor cells (6–8). Detection of these circulating tumor cells (CTC) might be used to monitor and predict therapy response in cancer treatment because repeated samples can be taken easily (9–11).

To be able to detect tumor cells in bone marrow or blood samples by monoclonal antibodies, different techniques and different markers are currently being used. The most commonly used marker to identify tumor cells from carcinomas are keratins, which have become the standard detection marker for DTCs and CTCs (1, 12–14). Twenty different keratins are described to be expressed in the cytoskeleton of epithelial cells. This family of intermediate filaments can be separated into neutral/basic (type I, K9-20) and acidic (type II, K1-8) keratins. Because hematopoietic cells rarely express keratin proteins detectable by immunostaining (15, 16), a specific detection can be achieved for cells from epithelial tumors, such as breast, prostate, colon, or lung carcinomas from blood cells. However, it is thought that tumor cells undergo epithelial–mesenchymal transition (EMT) to be able to migrate and invade the body; during this process the cytoskeleton is rearranged and epithelial markers such as E-cadherin, claudins, and keratins are downregulated (17–19). The question, whether current detection methods using keratin antibodies are able to detect all CTCs is still an open issue to be answered. If indeed not all CTCs are detected by the common detection methods, false-negative results will be the consequence and would hamper the clinical implementation of CTCs as diagnostic marker.

Because it has been reported that only a subset of keratins are downregulated during EMT, and others might be stably expressed or even upregulated, we hypothesized that an increased detection of tumor cells could be achieved by using more antibodies against a broader range of keratins than the current standards. In the study presented here, we have therefore investigated the keratin expression pattern of 11 breast cancer cell lines and 50 primary tumor specimens with their corresponding metastases. Moreover, we have studied the detection rate of CTCs in blood samples of 70 metastatic breast cancer patients, comparing the sensitivity of different keratin antibodies, and analyzed the clinical relevance of our findings by follow-up analyses of the outcome of cancer patients.

Antibodies

Protein expression of keratin 2, 4–10, 13–16, 18, and 19 was assessed by using a combination of pan-keratin antibodies C11 (catalog no. ab7753; Abcam), AE1 (catalog no. MAB1612; Millipore), AE3 (catalog no. MAB1611; Millipore), and keratin 7 antibody (catalog no. ab9021; Abcam; Table 1), further referred to as the keratin antibody cocktail. Furthermore, we have compared the keratin antibody cocktail with 2 of the most commonly used keratin antibodies for the detection of epithelial tumor cells: C11 and A45-B/B3 (catalog no. R002A; Micromet; refs. 13, 14, 20). The keratins detected by these antibodies are described in Table 1. Specificity of the antibodies has been confirmed by the manufacturers and in previous studies (21, 22).

Table 1.

Keratin antibodies and detection properties (K: Keratin)

KSize (kDa)C11A45-B/B3AE1AE3K7
2 65.5    •  
4 59 •   •  
5 58 •   •  
6 56 •   •  
7 54     • 
8 52.5 • •  •  
9 64   •   
10 56.5 •  •   
13 53 •     
14 50   •   
15 49   •   
16 48   •   
18 45 • •    
19 40  • •   
KSize (kDa)C11A45-B/B3AE1AE3K7
2 65.5    •  
4 59 •   •  
5 58 •   •  
6 56 •   •  
7 54     • 
8 52.5 • •  •  
9 64   •   
10 56.5 •  •   
13 53 •     
14 50   •   
15 49   •   
16 48   •   
18 45 • •    
19 40  • •   

Cell lines

The following breast cancer cell lines were acquired from American Type Culture Collection and cultured under the prescribed conditions: BT-20, BT-474, BT-549, GI-101, Hs 578T, MCF-7, MDA-MB-231, MDA-MB-435S, MDA-MB-468, SK-BR-3, and T-47D. Cells were grown in a 75-cm2 flask until confluency was reached. Cells were harvested by using trypsin/EDTA (catalog no. R001100; Gibco), washed, and resuspended in 1×PBS for either spiking experiments or cytospin preparation for direct staining.

Patients and blood sampling

Seventy metastatic breast cancer patients were included into the study (patient details are described in Supplementary Data S1). Patients were treated for metastatic breast cancer at the University Breast Center and received chemotherapy, endocrine treatment, or treatment with trastuzumab alone or in combination with chemotherapy based on national and international guidelines. After a median of 3.6 years after diagnosis of metastatic disease, patients were included into our study from 2009 to 2011; there were no exclusion criteria. All patients and healthy volunteers gave informed consent and the examination of blood samples was carried out blindly and was approved by the local ethics review board. Ten microliters of blood was drawn in an EDTA tube from all 70 metastatic breast cancer patients during their anticancer treatment at the hospital and from 10 healthy controls. Blood samples were stored at room temperature and processed within 24 hours after collection as follows: full blood was transferred to a 50-mL tube containing 40 mL HBSS (catalog no. L2045; Biochrom) and centrifuged at 400 ×g for 10 minutes at 4°C. Supernatant was removed by pipetting and pellet was resuspended in 30 mL 1×PBS. Cell suspension was added to 20 mL Ficoll (catalog no. 17-1440-03; GE Healthcare). Tube was spun at 400 ×g for 30 minutes at 4°C without acceleration and deceleration. The interface and supernatant, containing the mononuclear cells (i.e., leukocytes and tumor cells), were transferred to a new 50-mL tube. The tube was filled with 1×PBS and centrifuged at 400 ×g for 10 minutes at 4°C. Supernatant was discarded and cell pellet was resuspended in 1 mL 1×H-Lysis buffer (catalog no. WL1000; R&D Systems) and incubated for 1 to 3 minutes with gentle shaking at room temperature. Thirty microliters of PBS was added and sample was centrifuged again at 400 ×g for 10 minutes at 4°C. Supernatant was discarded and pellet was resuspended in 5 to 10 mL 1× PBS for cytospin preparation.

Cytospin preparation

For human blood and cell line cells, the following procedure was applied to prepare cytospins. The amount of viable cells was calculated by using a Neubauer cell counting chamber. Glass slides were placed into cytospin funnels of 220 or 120 mm2 in diameter and 5 × 105 or 5 × 104 cells were transferred for patient's blood or cell line experiments, respectively. Slides were spun at 170 ×g for 3 minutes at room temperature, and cytospins were allowed to air dry before staining was done.

Cytospin staining procedure

For staining of cytospins, primary antibodies A45-B/B3 (1:100), AE1 (1:100), AE3 (1:100), C11 (1:750), and CK7 (1:1000), or IgG1 isotype control clone MOPC21 (1:500; catalog no. M9269; Sigma) were used. Antibodies were diluted in 10% AB-Serum/PBS (catalog no. 805135; Biotest). The cell section was bordered using a Dako Pen (catalog no. S2002; Dako). The following solutions were added per slide, followed by incubation, removal of solution, and 3× washing with PBS for 3 minutes: (i) 4% PFA (pH 7.4), 10 minutes; (ii) 150 μL 10% AB-serum/PBS, 20 minutes, no washing; (iii) 125 μL primary antibody or MOPC21, 45 minutes; (iv) 125 μL secondary antibody (1:20; catalog no. Z0259; Dako), 30 minutes; (v) 125 μL APAAP complex (1:100; catalog no. D0651; Dako), 30 minutes; (vi) 150 μL [solution 2 (415 μL NaNO2, 166 μL 5% New Fuchsine, 3 minutes incubation in darkness) was added to solution 1 (78.8 mL H2O, 4.2 mL 1mol/L Tris pH 9.5, 200 μL Levamisole) and mixed by pipetting. Solution 3 (42 mg Naphthol AS-BI phosphate in 500 μL dimethylformamide) was added and mixed. The mixture was filtered with 0.2-μm filter], 20 minutes. Nuclei were counterstained with Mayer's hemalaun blue (1:5, catalog no. 1.09249; Merck) for 15 to 20 seconds. Slides were washed in H2O shortly, dipped 3× in 0.1% acetic acid and 3× in H2O, washed once in 0.1% NaHCO2 for 30 seconds, and finally washed 3× in H2O for 3 minutes.

Experimental setup

Four stainings were compared: (i) C11; (ii) A45-B/B3; (iii) Keratin antibody cocktail (AE1, AE3, C11, and CK7); and (iv) IgG1 isotype control (see Table 1 for antibody properties). For each cell line, 12 slides were prepared and each staining was carried out 3 times. Of each slide, 9 fields in the center of the cytospin were photographed; large cell clusters were avoided. Photographs were taken under 200× magnifications. Three individuals (S.J., J.S., and S.D.) each blindly reviewed 3 of the 9 pictures per slide and in total 36 photos per cell line. Positive and negative (as described above) cells were counted and the mean percentage of positive cells and SD per cell line were calculated.

Of the blood cytospins, a total of 10 slides were stained per blood sample. Nine slides were analyzed with A45-B/B3, C11, and keratin antibody cocktail (each 3 times as before) and one with IgG1 isotype control MOPC21. Slides were completely reviewed using the ChromaVision ACIS II automated microscope system. All positive hits were reviewed manually.

Western blot

To confirm immunocytochemical staining, Western blot was carried out as described elsewhere (23). A 5% SDS polyacrylamide gel for protein collection and a 10% SDS polyacrylamide gel for protein migration were used to separate 5 μg proteins, as determined with calibration series using Pierce BCA reagents (catalog no. 23225; Thermo Scientific). Proteins were transferred onto transfer membrane (Immobilon-psq, ISEQ00010), using transfer buffer containing 50 mmol/L Tris, 10 mmol/L Calycin, 0.4% SDS, and 20% Methanol. Next, blocking was carried out in 5% milk powder/TBST for 1 hour, followed by 3 × 5 minutes washing with TBST. Transfer membrane was incubated with primary antibodies over night as follows: A45-B/B3 (1:1,000) in 5% milk powder/TBST, AE1 (1:150,000) in 5% bovine serum albumin (BSA)/TBST, AE3 (1:150,000) in 5% BSA/TBST, C11 (1:1,000) in 5% BSA/TBST, and CK7 (1:1000) in 5% BSA/TBST. Next, 3 × 5 minutes washing was carried out and 1:1,500 secondary mouse antibody (Polyclonal Goat Anti Mouse Immunoglobulins/HRP Dako P0447) in 5% milk powder/TBST was incubated for 1 hour. After 3 × 5 minutes washing with TBST, Western blots were developed using 1:50 ECL Plus (Amersham).

Tissue microarrays

A commercially available tissue microarray (TMA) of 100 formalin-fixed, paraffin-embedded cores (1-mm diameter) was used in this study (catalog no. BR1005; US Biomax Inc.). The TMA contained 100 tissue cores from 50 breast cancer patients, unrelated to the earlier mentioned patients of this study. From each patient, one core originates from the primary breast carcinoma and one core originates from the corresponding carcinoma metastases in lymph nodes. Preceding the staining of the TMAs, slides were baked for 2 hours at 55°C and deparaffinized in Xylene for 15 minutes 2 times. Rehydration was done in Ethanol series (2× 100%, 2× 96%, 2× 80%, 2× 70%) 2 minutes each and finally 1×PBS. Antigen retrieval was done in boiling Sodium-citrate-buffer (10 mmol/L sodium citrate pH 6.0, 0.05% Tween 20) for 15 minutes. TMAs were washed in TBS twice and blocking was done with ready-to-use Protein block (catalog no. X0909; DakoCytomation) for 20 minutes at room temperature. Incubation of primary antibodies [A45-B/B3 (1:300), AE1 (1:300), AE3 (1:300), C11 (1:2250), and CK7 (1:3,000)] was done over night at 4°C. Following staining steps were carried out similarly as described for the cytospins. In total, 6 TMAs were stained with: the antibody cocktail, C11, and A45-B/B3, respectively, all in duplicate. Each slide was read independently by S.J. and J.H., blindly. Discrepancies were reevaluated and a consensus was found. Scoring was done as following: 0 (negative), no staining; 1+ (positive), a faint/barely perceptible staining; 2+ (positive), a weak to moderate staining; 3+ (positive), a strong staining was observed in more than 1% of the tumor cells.

Statistical analyses

Statistical significance between the fractions of positive cell line cells was calculated by ANOVA using each reviewed field as a measurement per staining (total fields per cell line per staining was 27) and reviewer as random variable. TMA results were categorized into positive and negative, and analyzed using the McNemar test. Survival analysis of patients tested for CTCs was done using the log-rank test after dividing the patient cohorts into CTC-positive and CTC-negative groups; HRs and 95% CI were calculated using Cox proportional hazards model. Statistical tools available from In-Silico Online, http://in-silico.net/statistics, 2012 (24).

In silico analysis

To examine a potential association of keratin gene expression with metastatic relapse in breast cancer, we examined the available gene expression data of 2 cohorts with a high rate of relapsed patients and information in metastasis-free survival (25, 26). Data were downloaded from GEO (GSE12276 and GSE2034), and all cases with relapse were used in our analysis. To retain the original series of the publications, there was an overlap of 12 cases in the 2 studies. For each individual set of keratin gene expression levels, data were sorted in decreasing order, along with the survival data of patients. The upper 35% of the data was labeled as “high keratin gene expression,” whereas the lower 35% of the data was labeled as “low keratin gene expression”; the remaining data (middle 30%) was discarded. The survival data of the low and high keratin gene expression groups were compared with each other using the log-rank test.

Keratin protein expression was assessed in breast cancer cell lines, primary breast tumors, lymph node metastases, and CTCs by keratin antibodies. The used antibodies were C11 (detects K4-6, 8, 10, 13, and 18), A45-B/B3 (detects K8, 18, and 19), and the keratin antibody cocktail (detects K2, 4–10, 13–16, 18, and 19), described in Table 1. In addition, keratin mRNA expression was determined by in silico analyses in metastatic breast cancer.

Keratin protein expression in breast cancer cell lines

Cytospins of 11 breast cancer cell lines were stained with C11, A45-B/B3, keratin antibody cocktail, and IgG1 isotype control (clone MOPC21) and reviewed under a 200× magnification. In total, 108 fields per cell line were reviewed and keratin-positive and keratin-negative cells were counted manually; a single field contained 80 cells on average [SD (s) = 25]. In 3 of 11 tested cell lines, a significant higher proportion of the cells was identified as positive with the keratin antibody cocktails as compared with C11 or A45-B/B3 alone (P < 0.01, ANOVA with Bonferroni–Holm correction; Table 2 and Supplementary Data S2). The biggest effect was seen in cell line MDA-MB-435S, in which 40% of the cells were stained using the keratin cocktail, whereas no cells were stained by C11 or A45-B/B3 alone. IgG1 isotype control was negative in all cases, showing the specificity of the stainings.

Table 2.

Mean percentage of positive cells and SD (s) of 11 breast cancer cell lines

Cell lineA45-B/B3 (s)C11 (s)KRT cocktail (s)P
BT-20 89.2 (5.9) 89.3 (6.8) 91.1 (5.6) >0.99 
BT-474 79.3 (10.6) 79.2 (7.0) 79.4 (9.0) >0.99 
BT-549 1.3 (1.9) <0.001 
GI-101 82.3 (8.9) 82.4 (6.1) 84.5 (5.9) >0.99 
Hs 578T 1.4 (1.7) <0.003 
MCF-7 86.0 (7.6) 84.9 (6.8) 84.6 (7.8) >0.99 
MDA-MB-231 96.7 (5.5) 99.1 (2.0) 98.4 (2.2) 0.38 
MDA-MB-435S 39.2 (11.1) <0.001 
MDA-MB-468 78.2 (10.2) 76.8 (10.1) 79.8 (7.8) >0.99 
SK-BR-3 92.3 (4.6) 91.9 (3.9) 93.4 (4.5) >0.99 
T-47D 94.5 (4.6) 94.0 (3.6) 96.1 (4.6) >0.99 
Cell lineA45-B/B3 (s)C11 (s)KRT cocktail (s)P
BT-20 89.2 (5.9) 89.3 (6.8) 91.1 (5.6) >0.99 
BT-474 79.3 (10.6) 79.2 (7.0) 79.4 (9.0) >0.99 
BT-549 1.3 (1.9) <0.001 
GI-101 82.3 (8.9) 82.4 (6.1) 84.5 (5.9) >0.99 
Hs 578T 1.4 (1.7) <0.003 
MCF-7 86.0 (7.6) 84.9 (6.8) 84.6 (7.8) >0.99 
MDA-MB-231 96.7 (5.5) 99.1 (2.0) 98.4 (2.2) 0.38 
MDA-MB-435S 39.2 (11.1) <0.001 
MDA-MB-468 78.2 (10.2) 76.8 (10.1) 79.8 (7.8) >0.99 
SK-BR-3 92.3 (4.6) 91.9 (3.9) 93.4 (4.5) >0.99 
T-47D 94.5 (4.6) 94.0 (3.6) 96.1 (4.6) >0.99 

NOTE: ANOVA was done between A45-B/B3, C11, and the KRT antibody cocktail, with the reviewer as random variable, to determine the significance in difference of positive cells.

P values are corrected for multiple testing (Bonferroni-Holm).

Next, Western blot was applied to investigate which specific keratins are expressed in cell lines BT-549, Hs 578T, and MDA-MB-435S. Protein extract from MDA-MB-231 was used as positive control, as this cell line had shown to be positive in all keratin subset mixes before (Table 2). Figure 1 shows the results of 5 Western blots, carried out with each of the individual pan-keratin antibodies and K7 antibody. As expected, MDA-MB-231 was positive for all applied antibodies. Cell line MDA-MB-435S, which was positive by ICC in our keratin cocktail, showed clear signals using AE1. The rare staining of cell lines BT-549 and Hs 578T using the keratin cocktail could be confirmed by the presence of a weak signal by AE1. These results showed that the AE1 pan-keratin antibody is important for the detection of keratins in breast cancer cells.

Figure 1.

Western blot done on cell lines BT-549, Hs 578T, MDA-MB-435S, and MDA-MB-231 (positive control). BT-549, Hs 578T, and MDA-MB-435S show signals at the height of K14, K16, and K19 as indicated by AE1.

Figure 1.

Western blot done on cell lines BT-549, Hs 578T, MDA-MB-435S, and MDA-MB-231 (positive control). BT-549, Hs 578T, and MDA-MB-435S show signals at the height of K14, K16, and K19 as indicated by AE1.

Close modal

Keratin gene expression in breast cancer cell lines (in silico analysis)

According to a study of Neve and colleagues, the 3 cell lines BT-549, Hs 578T, and MDA-MB-435S, in which no positive cells were detected using C11 or A45-B/B3, are all of the basal-like B (claudin-low) breast cancer subtype. MDA-MD-231 has also been classified as basal-like B but was positive in all stainings. All other cell lines from our study were classified as basal-like A or luminal subtype (27). This motivated us to analyze the keratin gene expression levels of the breast cancer cell lines studied by Neve and colleagues to investigate whether there are differences between keratin expression levels of claudin-low, basal-like, and luminal breast cancer cell lines. Complete data of 52 breast cancer cell lines were available, of which 25 were of the luminal subtype, 12 basal-like A, and 14 basal-like B (claudin-low).

Supplementary Data S3 summarizes the average gene expression levels of the different keratins per breast cancer subtype. As shown, KRT8 and KRT19 are significantly less expressed in basal-like B cell lines compared with basal-like A and luminal cell lines (P < 0.0005, Bonferonni). In addition, KRT18 is significantly lower expressed in all basal-like cell lines compared with luminal cell lines. KRT5, 6B, 14, and 16 are significantly higher expressed in basal-like A cell lines compared with luminal or basal-like B cell lines. Thus, the expression patterns of certain keratins are correlated to breast cancer subtypes.

Keratin protein expression in primary breast carcinomas and lymph node metastases

A TMA containing cores of 50 primary breast carcinomas and their corresponding lymph node metastases was investigated to assess the expression of keratins in tumor tissue and metastatic tumor cells. Concordance in staining intensity and positive cell percentage between the independent evaluations of the individual slides was approximately 90% per slide. Using C11 and A45-B/B3, 14 (28%) and 27 cases (54%) of the primary tumors were detected to be positive for keratin, respectively. When using the keratin antibody cocktail, all primary tumors (100%) were found to be positive (defined as a positive staining in >1% of the cells). For the corresponding metastasis to lymph nodes, 19 (38%) and 28 (56%) were positively stained using C11 and A45-B/B3, respectively; however, keratin status of the primary tumor differed frequently from the corresponding lymph node metastasis: 12 of 50 (24%) of the C11-stained cases and 15 of 50 (30%) of the A45-B/B3–stained cases showed deviation (Fig. 2). Forty-nine of 50 lymph node metastasis (98%) were stained positive using our combination of keratin antibodies (Fig. 2). These results suggested that the keratins detected by A45-B/B3 are expressed more often in breast cancer and breast cancer metastases in general, as the keratins detected by C11 or that by A45-B/B3 is more sensitive than C11. Moreover, by adding additional (pan-) keratin antibodies, we could stain significantly more tumors and lymph node metastases than with C11 or A45-B/B3 alone (P < 0.0001, McNemar test). Because keratin status in many of the lymph node metastases did not correspond to that of the primary tumor of same patient, detected by C11 or A45-B/B3, these results indicated that keratin expression changes during metastasis. This emphasizes the need of a broader spectrum of keratin antibodies to detect metastasized cancer cells.

Figure 2.

TMA containing core biopsies from 50 primary breast cancers (row A–E) and corresponding lymph node metastases (row F–J), stained by C11, A45-B/B3, and the keratin antibody cocktail. No staining (negative), white; a faint/barely perceptible staining (1+), light gray; a weak to moderate staining (2+), dark gray; a strong staining (3+), black. Percentage of positive cells as indicated or 100% (filled circles).

Figure 2.

TMA containing core biopsies from 50 primary breast cancers (row A–E) and corresponding lymph node metastases (row F–J), stained by C11, A45-B/B3, and the keratin antibody cocktail. No staining (negative), white; a faint/barely perceptible staining (1+), light gray; a weak to moderate staining (2+), dark gray; a strong staining (3+), black. Percentage of positive cells as indicated or 100% (filled circles).

Close modal

Keratin gene expression in metastasizing breast carcinomas (in silico analysis)

To investigate the relationship between the expression of keratins in metastatic breast cancer and relapse-free survival of patient, we carried out an in silico analysis on 2 previously published studies (25, 26). These 2 studies assessed mRNA gene expression by the HG-U133A GeneChip from Affymetrix in metastatic and nonmetastatic primary breast tumors. We extracted the expression levels of all keratins present in these studies from the cases with metastatic relapse only. Probes for all keratins investigated in our experimental setup were present in these data, except for KRT19. The study of Bos and colleagues contained primary tumors of 204 patients, who later developed metastatic disease, and the cohort of Wang and colleagues 107 metastatic cases.

After collecting the required data, survival analysis was done on those cases expressing keratins high or low (see Materials and Methods), for each keratin separately. In both studies of Bos and colleagues and Wang and colleagues, only expression of KRT16 showed a significant difference in survival (Fig. 3). For the cohort from the study of Bos and colleagues log-rank test resulted in P = 0.001 and cox regression in a HR of 1.69 (95% CI: 1.20–2.37). The median relapse-free survival time after primary diagnosis related to low and high expression of KRT16 was 23 and 14.8 months, respectively. For the cohort from the study of Wang and colleagues log-rank test resulted in P = 0.010 and cox regression in a HR of 1.79 (95% CI: 1.12–2.84). The median relapse-free survival related to low and high expression of KRT16 was 28 and 17.5 months. These results suggested that a high expression of KRT16 is associated with poor survival, whereas a low expression of KRT16 is correlated with a better survival in metastasizing breast cancer.

Figure 3.

Kaplan–Meier estimate of survival function of metastatic breast cancer patients separated on high expression of KRT16 (grey line) and low expression of KRT16 (black line). Data were taken from 2 published studies (25, 26). Longer survival correlated with low expression of KRT16 in both studies of Bos and colleagues (P = 0.001, n = 70) and Wang and colleagues (P = 0.010, n = 37) calculated with log-rank test.

Figure 3.

Kaplan–Meier estimate of survival function of metastatic breast cancer patients separated on high expression of KRT16 (grey line) and low expression of KRT16 (black line). Data were taken from 2 published studies (25, 26). Longer survival correlated with low expression of KRT16 in both studies of Bos and colleagues (P = 0.001, n = 70) and Wang and colleagues (P = 0.010, n = 37) calculated with log-rank test.

Close modal

Keratin protein expression in CTCs of breast cancer patients

Blood of an independent set of 70 metastatic breast cancer patients from our hospital was investigated for the presence of CTCs by capturing all mononuclear cells from 10 mL blood using Ficoll density gradient separation and subsequent staining with A45-B/B3, C11, and the keratin antibody cocktail. In total, CTCs were identified in 35 samples (50%), as defined by the detection of at least one positive cell; 57% of the positive patients was identified by A45-B/B2, 40% by C11, and 83% by the keratin antibody cocktail. For each individual patient, the fraction of positive cells detected with A45-B/B3, C11, and keratin antibody cocktail was calculated. Next, the relative fraction of positive cells for each staining over the whole cohort was calculated. On average, A45-B/B3 and C11 detected 26% and 14% of all the positive cells, respectively; the keratin antibody cocktail detected 61% (Supplementary Data S1). Figure 4A depicts the distribution of the number of cells detected per antibody in the metastatic breast cancer patients. These results suggested that more CTCs can be detected by using our keratin antibody cocktail as compared with A45-B/B3 or C11 alone (P < 0.01, 2-sided paired t test).

Figure 4.

CTCs in metastatic breast cancer patients (n = 70) identified by A45-B/B3, C11, and keratin antibody cocktail. A, total number of keratin-positive cells. By using the keratin antibody cocktail, more CTCs could be identified. B to D, Kaplan–Meier estimates of survival function of 70 metastatic breast cancer patients separated on CTC status as detected by A45-B/B3 [HR: 1.36 (95% CI: 0.53–3.46)], C11 [HR: 1.52 (95% CI: 0.55–4.20)] and keratin antibody cocktail [HR: 4.16 (95% CI: 1.49–11.59)]. Statistical significance determined by log-rank test. Between parentheses is the median survival in years of the corresponding patient group. CTC-negative groups are indicated by black lines, CTC-positive (≥1 cell) groups are indicated by grey lines, and censored patients are indicated by vertical bars (|).

Figure 4.

CTCs in metastatic breast cancer patients (n = 70) identified by A45-B/B3, C11, and keratin antibody cocktail. A, total number of keratin-positive cells. By using the keratin antibody cocktail, more CTCs could be identified. B to D, Kaplan–Meier estimates of survival function of 70 metastatic breast cancer patients separated on CTC status as detected by A45-B/B3 [HR: 1.36 (95% CI: 0.53–3.46)], C11 [HR: 1.52 (95% CI: 0.55–4.20)] and keratin antibody cocktail [HR: 4.16 (95% CI: 1.49–11.59)]. Statistical significance determined by log-rank test. Between parentheses is the median survival in years of the corresponding patient group. CTC-negative groups are indicated by black lines, CTC-positive (≥1 cell) groups are indicated by grey lines, and censored patients are indicated by vertical bars (|).

Close modal

To assess the rate of possible false positivity, blood of 10 apparently healthy volunteers was investigated. Of each person's blood, 20 cytospins containing each 500,000 cells were prepared and stained with keratin antibody cocktail. In all 200 slides, we have found no keratin-positive cells. This result confirmed that the specificity is as high as for C11 and A45-B/B3 that have already been shown in previous studies (13, 20).

Survival analysis of CTC-positive metastatic breast cancer patients

The 70 metastatic breast cancer patients were divided by their CTC status (positive or negative), detected by each of the tested keratin antibodies individually. Next, survival analyses starting from the day of blood sampling were carried out (median follow-up of 2.23 years), and statistical significance was calculated using the log-rank test; at time of analysis, 27% (19 of 70) of the women had died. No difference in survival was shown between the CTC-positive and CTC-negative patients as defined by A45-B/B3 [P = 0.8829, HR: 1.36 (95% CI: 0.53–3.46), Fig. 4B] or C11 [P = 0.8715, HR: 1.52 (95% CI: 0.55–4.20), Fig. 4C]. In contrast, a significant difference could be observed in the survival estimates of CTC-negative and CTC-positive patients as defined by the keratin antibody cocktail [P = 0.0035, HR: 4.16 (95% CI: 1.49–11.59), Fig. 4D]. Median survival time of CTC-positive patients was 1.78 years, whereas for the CTC-negative patients the median survival was 2.26 years. The results of this exploratory analysis indicated that detection of CTCs by the antibody cocktail is of prognostic relevance in metastatic breast cancer patients.

Identification of CTCs in the bloodstream of cancer patients can serve as a marker for early detection of micrometastasis, but it may also be used to monitor therapy response in early- and advance-stage patients as a “liquid biopsy” which may allow obtaining important information on the molecular make-up of occult and/or overt metastases, respectively. The detection of CTCs has the potential to become increasingly important in the clinical management of cancer patients. Currently, the most widely used markers for the identification of CTCs are keratins; however, little is known about the changes in the keratin expression pattern in tumor cells during metastatic spread and progression. It has been described that cells from metastasizing breast cancer may undergo EMT and that they may lose the expression of several keratins during this process (19, 28, 29). Usage of pan-keratin antibodies directed against a limited number of keratins may therefore lead to a failure in detecting CTCs and finally lead to false-negative results. The use of a cocktail of antibodies, directed against an even wider range of keratins that also cover the biologic behavior of keratins during the metastatic process and by this also keratins differentially regulated during EMT, may substantially improve the CTC detection rate. Indeed, it has been shown recently that an increased detection of DTCs in bone marrow can be achieved by applying a broader set of keratin antibodies (30). In this study, we evaluated keratin expression of primary tumor, lymph node metastasis, and CTCs. To further explore the dynamic changes in keratin expression during metastatic progression, future studies might also include keratin expression in distant metastases. Cells that leave the blood circulation and develop to metastases might have to undergo MET and thereby restoring keratin expression that was deregulated during EMT in CTCs (12, 18).

To test our hypothesis that a combination of available pan-keratin antibodies may improve the detection of CTCs, we first investigated 11 breast cancer cell lines with 2 of the most commonly used pan-keratin antibodies (C11 and A45-B/B3) and a keratin antibody cocktail, which was chosen on the basis of possible changes in keratin expression as a biologic feature during metastasis. Most cell lines expressed enough keratins to be clearly stained with either C11 or A45-B/B3 alone, except for the 3 cell lines BT-549, Hs 578T, and MDA-MB-435S. Using the keratin antibody cocktail, we can show that 1% to 2% of the BT-549 and HS 578T cells express keratin and approximately 40% of the MDA-MB-435S cells. Western blot analyses showed that these 3 cell lines (BT-549, Hs 578T, and MDA-MB-435S) could only be detected using AE1 and show expression of K16 and K19, although the latter could not be confirmed by A45-B/B3, possibly, due to its weak expression. In addition, K5 and K6 could be detected using AE3 but not with C11. This suggests that K5 and K6 might be below the detection level of C11.

According to the dataset of Neve and colleagues, the 3 cell lines BT-549, Hs 578T, and MDA-MB-435S are classified as breast cancer subtype basal-like B, which is also known as the claudin-low subtype (27). Claudin-low breast tumors originate from the myoepithelial cell layer that normally express K5, 14, and 17 only (31) and are characterized by low expression of epithelial markers such as keratins (32–34). Indeed, in silico analysis confirmed the lack of KRT5/14/17 in these 3 claudin-low cell lines; however, ICC analyses showed heterogeneity in keratin protein expression among the cells and K14 could be detected by Western blot. Because basal-like cell lines originate from luminal progenitor cells, these and luminal cell lines are expected to express K6-8/18/19 (35). The data from Neve and colleagues showed several differences between the luminal and basal-like cell lines in respect to these keratins. KRT6 was apparently downregulated in luminal cell lines, KRT18 was downregulated in basal-like cell lines, and unexpectedly, KRT16 was upregulated in basal-like breast cancer cell lines, which we could also detect by Western blot. Surprisingly, cell line MDA-MB-231 has also been classified as basal-like B but showed expression of many keratins, including those that would be expected in luminal breast cancer only.

Next, we investigated a TMA containing 100 core biopsies from primary tumors and corresponding lymph node metastases. We show that keratin expression is limited to only a set of keratins in primary tumors as well as in their lymph node metastases. In addition, keratin expression seems frequently to be altered between the primary tumor and the corresponding metastasis. However, in most of the investigated core biopsies, keratin expression could be identified by using the keratin antibody cocktail. These results indicate that keratin expression patterns are due to changes during invasion, or that keratin expression is heterogeneous within the primary tumor mass, similar to results obtained from our cell line experiments. On the basis of this heterogeneity in keratin expression and the results of others (19, 36, 37) it might be speculated that keratins can also be used in the future to identify the single-cell subpopulation that is able to initiate metastasis. It should be noted, however, that this investigation is limited by the fact that the staining procedure and tissue conditions can influence detection rates and that tissues identified as negative in this study might express low amounts of the respective keratins. Although we are confident that the staining procedure is optimal because the results were very reproducible, we could not control for the quality of the tissue. In addition, molecular subtype of the primary tumors is not known, a factor that is also correlated to expression of certain keratins.

As a next step to investigate whether broader keratin detection also leads to improved detection of CTCs, we investigated 70 blood samples of metastatic breast cancer patients for CTCs. To identify these rare cells, we first used a Ficoll density gradient to enrich for mononuclear cells. Although different methods exist to specifically enrich for tumor cells (38), enrichment by using cellular markers, such as EpCAM, might result in loss of cells of interest if they do not express these markers. Screening the samples of metastatic breast cancer patients for CTCs showed that an increased detection of tumor cells could be achieved with the keratin antibody cocktail. It should be noted that given the low number of detectable CTCs, it might be possible that CTCs have not been equally distributed among the investigated slides, but there is no possibility to correct for that. However, due to random staining and sufficient sample size, it seems not likely that there is any bias of the amount of CTCs among the stained slides. By carrying out the same screening procedure on healthy controls, with no findings of positive cells, we show the high specificity of the keratin antibodies used in this study and are convinced that we did not identify false positives in the cancer samples.

Most importantly, survival analyses on the patient cohort divided by positive and negative CTC status indicated that the CTC-negative patients (as defined by the keratin antibody cocktail) have a better survival. Apparently, more CTCs are detected by the keratin antibody cocktail as compared with C11 or A45-B/B3 alone that might be biologically relevant, resulting in a better separation between better and poor survival. A potential limitation of this investigation is that the patient cohort is not uniformly treated and that we have no untreated cohort to determine the prognostic impact of the examined factors. Therefore, it cannot be ruled out completely that also the time of blood sampling and the treatment influenced these results. Thus, to test the clinical significance of the cocktail, a larger cohort with uniform treatment and a longer follow-up will be required.

Investigating the possible association between survival and keratin expression in metastasizing breast cancer in silico, we found in published data sets from 2 earlier studies that patients who showed high expression of KRT16 in their primary tumor had a shorter relapse-free survival when compared with patients with KRT16 low expressing tumors (25, 26). In addition, we have seen that KRT16 upregulation is also a common phenomenon in basal-like breast cancer cell lines. These results could reflect the fact that keratin 16 expression is associated with a more aggressive course of breast cancer and would support our finding that CTCs detected by the antibody cocktail is associated with impaired survival in the metastatic setting. Our findings stress out the importance of the detection of K16, which is not recognized by either C11 or A45-B/B3. Expression of K16 in primary breast cancer and CTCs has to be further investigated for its potential clinical use.

In conclusion, our results indicate that breast tumors can alter expression of certain keratins during the process of metastatic development, which would make identification of CTCs of specifically these tumors by the currently used standard pan-keratin antibodies impossible. In this study, we show that by broadening up the spectrum of keratin detection, on the basis of the knowledge of the metastatic process compared with the use of single keratin antibodies, an increased detection of CTCs can be achieved, reducing the number of false-negative results. Our investigation also suggests that modulation of keratins due to EMT occurs frequently on CTCs in breast cancer patients and may be associated with an unfavorable outcome.

No potential conflicts of interest were disclosed.

The authors thank Stephanie Degwert and Susanne Hoppe for technical assistance.

This work has been supported by the German Ministry for Research and Education (BMBF), grant number 01ES0724.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Pantel
K
,
Brakenhoff
RH
. 
Dissecting the metastatic cascade
.
Nat Rev Cancer
2004
;
4
:
448
56
.
2.
Eccles
SA
,
Welch
DR
. 
Metastasis: recent discoveries and novel treatment strategies
.
Lancet
2007
;
369
:
1742
57
.
3.
Molloy
T
,
van't Veer
LJ
. 
Recent advances in metastasis research
.
Curr Opin Genet Dev
2008
;
18
:
35
41
.
4.
Braun
S
,
Vogl
FD
,
Naume
B
,
Janni
W
,
Osborne
MP
,
Coombes
RC
, et al
A pooled analysis of bone marrow micrometastasis in breast cancer
.
N Engl J Med
2005
;
353
:
793
802
.
5.
Braun
S
,
Auer
D
,
Marth
C
. 
The prognostic impact of bone marrow micrometastases in women with breast cancer
.
Cancer Invest
2009
;
27
:
598
603
.
6.
Muller
V
,
Stahmann
N
,
Riethdorf
S
,
Rau
T
,
Zabel
T
,
Goetz
A
, et al
Circulating tumor cells in breast cancer: correlation to bone marrow micrometastases, heterogeneous response to systemic therapy and low proliferative activity
.
Clin Cancer Res
2005
;
11
:
3678
85
.
7.
Nagaiah
G
,
Abraham
J
. 
Circulating tumor cells in the management of breast cancer
.
Clin Breast Cancer
2010
;
10
:
209
16
.
8.
Molloy
TJ
,
Bosma
AJ
,
Baumbusch
LO
,
Synnestvedt
M
,
Borgen
E
,
Russnes
HG
, et al
The prognostic significance of tumor cell detection in the peripheral blood versus the bone marrow in 733 early-stage breast cancer patients
.
Breast Cancer Res
2011
;
13
:
R61
.
9.
Nakamura
S
,
Yagata
H
,
Ohno
S
,
Yamaguchi
H
,
Iwata
H
,
Tsunoda
N
, et al
Multi-center study evaluating circulating tumor cells as a surrogate for response to treatment and overall survival in metastatic breast cancer
.
Breast Cancer
2010
;
17
:
199
204
.
10.
Gazzaniga
P
,
Naso
G
,
Gradilone
A
,
Cortesi
E
,
Gandini
O
,
Gianni
W
, et al
Chemosensitivity profile assay of circulating cancer cells: prognostic and predictive value in epithelial tumors
.
Int J Cancer
2010
;
126
:
2437
47
.
11.
Giuliano
M
,
Giordano
A
,
Jackson
S
,
Hess
KR
,
De
GU
,
Mego
M
, et al
Circulating tumor cells as prognostic and predictive markers in metastatic breast cancer patients receiving first-line systemic treatment
.
Breast Cancer Res
2011
;
13
:
R67
.
12.
Pantel
K
,
Brakenhoff
RH
,
Brandt
B
. 
Detection, clinical relevance and specific biological properties of disseminating tumour cells
.
Nat Rev Cancer
2008
;
8
:
329
40
.
13.
Braun
S
,
Pantel
K
,
Muller
P
,
Janni
W
,
Hepp
F
,
Kentenich
CR
, et al
Cytokeratin-positive cells in the bone marrow and survival of patients with stage I, II, or III breast cancer
.
N Engl J Med
2000
;
342
:
525
33
.
14.
Riethdorf
S
,
Fritsche
H
,
Muller
V
,
Rau
T
,
Schindlbeck
C
,
Rack
B
, et al
Detection of circulating tumor cells in peripheral blood of patients with metastatic breast cancer: a validation study of the CellSearch system
.
Clin Cancer Res
2007
;
13
:
920
8
.
15.
Borgen
E
,
Naume
B
,
Nesland
JM
,
Kvalheim
G
,
Beiske
K
,
Fodstad
O
, et al
Standardization of the immunocytochemical detection of cancer cells in BM and blood: I. establishment of objective criteria for the evaluation of immunostained cells
.
Cytotherapy
1999
;
1
:
377
88
.
16.
Jung
R
,
Petersen
K
,
Kruger
W
,
Wolf
M
,
Wagener
C
,
Zander
A
, et al
Detection of micrometastasis by cytokeratin 20 RT-PCR is limited due to stable background transcription in granulocytes
.
Br J Cancer
1999
;
81
:
870
3
.
17.
Radisky
DC
. 
Epithelial-mesenchymal transition
.
J Cell Sci
2005
;
118
:
4325
6
.
18.
Thiery
JP
,
Sleeman
JP
. 
Complex networks orchestrate epithelial-mesenchymal transitions
.
Nat Rev Mol Cell Biol
2006
;
7
:
131
42
.
19.
Willipinski-Stapelfeldt
B
,
Riethdorf
S
,
Assmann
V
,
Woelfle
U
,
Rau
T
,
Sauter
G
, et al
Changes in cytoskeletal protein composition indicative of an epithelial-mesenchymal transition in human micrometastatic and primary breast carcinoma cells
.
Clin Cancer Res
2005
;
11
:
8006
14
.
20.
Cristofanilli
M
,
Budd
GT
,
Ellis
MJ
,
Stopeck
A
,
Matera
J
,
Miller
MC
, et al
Circulating tumor cells, disease progression, and survival in metastatic breast cancer
.
N Engl J Med
2004
;
351
:
781
91
.
21.
Bartek
J
,
Bartkova
J
,
Schneider
J
,
Taylor-Papadimitriou
J
,
Kovarik
J
,
Rejthar
A
. 
Expression of monoclonal antibody-defined epitopes of keratin 19 in human tumours and cultured cells
.
Eur J Cancer Clin Oncol
1986
;
22
:
1441
52
.
22.
Conrad
K
,
Karsten
U
,
Frank
KH
,
Pasternak
G
. 
Immunoblotting of monoclonal anti-cytoskeletal antibodies with cytoskeletal extracts of HeLa and MCF-7 cells
.
Biomed Biochim Acta
1988
;
47
:
697
703
.
23.
Bartkowiak
K
,
Wieczorek
M
,
Buck
F
,
Harder
S
,
Moldenhauer
J
,
Effenberger
KE
, et al
Two-dimensional differential gel electrophoresis of a cell line derived from a breast cancer micrometastasis revealed a stem/progenitor cell protein profile
.
J Proteome Res
2009
;
8
:
2004
14
.
24.
in-silico.net [Internet]. In-Silico Online, statistical tools; [updated 2011 Dec; cited 2012 Jan]. Available from
: http://in-silico.net/statistics/.
25.
Bos
PD
,
Zhang
XH
,
Nadal
C
,
Shu
W
,
Gomis
RR
,
Nguyen
DX
, et al
Genes that mediate breast cancer metastasis to the brain
.
Nature
2009
;
459
:
1005
9
.
26.
Wang
Y
,
Klijn
JG
,
Zhang
Y
,
Sieuwerts
AM
,
Look
MP
,
Yang
F
, et al
Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer
.
Lancet
2005
;
365
:
671
9
.
27.
Neve
RM
,
Chin
K
,
Fridlyand
J
,
Yeh
J
,
Baehner
FL
,
Fevr
T
, et al
A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes
.
Cancer Cell
2006
;
10
:
515
27
.
28.
Raimondi
C
,
Gradilone
A
,
Naso
G
,
Vincenzi
B
,
Petracca
A
,
Nicolazzo
C
, et al
Epithelial-mesenchymal transition and stemness features in circulating tumor cells from breast cancer patients
.
Breast Cancer Res Treat
2011
;
130
:
449
55
.
29.
Gradilone
A
,
Raimondi
C
,
Nicolazzo
C
,
Petracca
A
,
Gandini
O
,
Vincenzi
B
, et al
Circulating tumor cells lacking cytokeratin in breast cancer: the importance of being mesenchymal
.
J Cell Mol Med
2011
;
15
:
1066
70
.
30.
Effenberger
KE
,
Borgen
E
,
Eulenburg
CZ
,
Bartkowiak
K
,
Grosser
A
,
Synnestvedt
M
, et al
Detection and clinical relevance of early disseminated breast cancer cells depend on their cytokeratin expression pattern
.
Breast Cancer Res Treat
2011
;
125
:
729
38
.
31.
Karantza
V
. 
Keratins in health and cancer: more than mere epithelial cell markers
.
Oncogene
2011
;
30
:
127
38
.
32.
Prat
A
,
Parker
JS
,
Karginova
O
,
Fan
C
,
Livasy
C
,
Herschkowitz
JI
, et al
Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer
.
Breast Cancer Res
2010
;
12
:
R68
.
33.
Prat
A
,
Perou
CM
. 
Mammary development meets cancer genomics
.
Nat Med
2009
;
15
:
842
4
.
34.
Sarrio
D
,
Rodriguez-Pinilla
SM
,
Hardisson
D
,
Cano
A
,
Moreno-Bueno
G
,
Palacios
J
. 
Epithelial-mesenchymal transition in breast cancer relates to the basal-like phenotype
.
Cancer Res
2008
;
68
:
989
97
.
35.
Groner
B
,
Vafaizadeh
V
,
Brill
B
,
Klemmt
P
. 
Stem cells of the breast and cancer therapy
.
Womens Health (Lond Engl)
2010
;
6
:
205
19
.
36.
Schaller
G
,
Fuchs
I
,
Pritze
W
,
Ebert
A
,
Herbst
H
,
Pantel
K
, et al
Elevated keratin 18 protein expression indicates a favorable prognosis in patients with breast cancer
.
Clin Cancer Res
1996
;
2
:
1879
85
.
37.
Alix-Panabieres
C
,
Vendrell
JP
,
Slijper
M
,
Pelle
O
,
Barbotte
E
,
Mercier
G
, et al
Full-length cytokeratin-19 is released by human tumor cells: a potential role in metastatic progression of breast cancer
.
Breast Cancer Res
2009
;
11
:
R39
.
38.
Talasaz
AH
,
Powell
AA
,
Huber
DE
,
Berbee
JG
,
Roh
KH
,
Yu
W
, et al
Isolating highly enriched populations of circulating epithelial cells and other rare cells from blood using a magnetic sweeper device
.
Proc Natl Acad Sci U S A
2009
;
106
:
3970
5
.