- Open Access
Impact of tissue sampling on accuracy of Ki67 immunohistochemistry evaluation in breast cancer
© The Author(s). 2016
Received: 26 February 2016
Accepted: 31 July 2016
Published: 30 August 2016
Gene expression studies have identified molecular subtypes of breast cancer with implications to chemotherapy recommendations. For distinction of these types, a combination of immunohistochemistry (IHC) markers, including proliferative activity of tumor cells, estimated by Ki67 labeling index is used. Clinical studies are frequently based on IHC performed on tissue microarrays (TMA) with variable tissue sampling. This raises the need for evidence-based sampling criteria for individual IHC biomarker studies. We present a novel tissue sampling simulation model and demonstrate its application on Ki67 assessment in breast cancer tissue taking intratumoral heterogeneity into account.
Whole slide images (WSI) of 297 breast cancer sections, immunohistochemically stained for Ki67, were subjected to digital image analysis (DIA). Percentage of tumor cells stained for Ki67 was computed for hexagonal tiles super-imposed on the WSI. From this, intratumoral Ki67 heterogeneity indicators (Haralick’s entropy values) were extracted and used to dichotomize the tumors into homogeneous and heterogeneous subsets. Simulations with random selection of hexagons, equivalent to 0.75 mm circular diameter TMA cores, were performed. The tissue sampling requirements were investigated in relation to tumor heterogeneity using linear regression and extended error analysis.
The sampling requirements were dependent on the heterogeneity of the biomarker expression. To achieve a coefficient error of 10 %, 5–6 cores were needed for homogeneous cases, 11–12 cores for heterogeneous cases; in mixed tumor population 8 TMA cores were required. Similarly, to achieve the same accuracy, approximately 4,000 nuclei must be counted when the intratumor heterogeneity is mixed/unknown. Tumors of low proliferative activity would require larger sampling (10–12 TMA cores, or 6,250 nuclei) to achieve the same error measurement results as for highly proliferative tumors.
Our data show that optimal tissue sampling for IHC biomarker evaluation is dependent on the heterogeneity of the tissue under study and needs to be determined on a per use basis. We propose a method that can be applied to determine the sampling strategy for specific biomarkers, tissues and study targets. In addition, our findings highlight the benefit of high-capacity computer-based IHC measurement techniques to improve accuracy of the testing.
Gene expression studies have identified distinct molecular subtypes of breast cancer (Luminal A, Luminal B, HER2-enriched, basal-like and normal breast-like) with markedly different behavior and prognosis . Meanwhile, clinical practice of decision making largely relies on the definition of Luminal A-like and Luminal B-like disease, based on a combination of estrogen receptor (ER), progesterone receptor (PgR) and Ki67 immunohistochemistry (IHC) . Proliferative activity of tumor cells, estimated by Ki67 labeling index (Ki67 LI) is a key indicator to support this stratification and provides strong prognostic and predictive information on response to chemotherapy . Clinical utility of Ki67 LI is hampered by the lack of robust measurement methodologies and widely acknowledged issue of intratumor Ki67 heterogeneity expression. Consequently, it is hard to achieve consensus on cut-off values to stratify the patients for therapeutic decisions . Great effort has been made to standardize the techniques for manual and digital/automated Ki67 LI measurement, including criteria for tissue sampling, hotspot detection, and digital image analysis (DIA) tools [4–11].
Recently, Ki67 expression across distinct categories of breast cancer specimens including whole slide surgical specimens, needle core biopsies and tissue microarrays (TMA) was investigated by Knutsvik et al. . They found significant differences of Ki67 LI estimates across the different sample categories and suggested that specimen-specific cut-off values should be applied for practical use. While the recommendation is logical and may compensate for the inherent differences of the tissue sampling, its implementation requires better knowledge of measurement accuracy that can be achieved by the techniques, in general. Additionally, Going  has previously pointed out that the counting rules depend on level of mitotic activity in tumors. This dependency has not been investigated for tumors with varying Ki67 proliferation rates.
TMA has been often applied for discovery and clinical studies of IHC biomarkers. Initially proposed by Battifora , it enables multiple testing on numerous tissue samples in a standardized, tissue-sparing, and high-throughput manner by assembling small core biopsies from morphologically representative areas of tissues onto a single paraffin block . The approach was further refined into to a precise technique by Kononen . One inherent drawback of the TMA technique is related to the limited fraction of the original sample included, raising the need to achieve/be aware of adequate sampling requirements . Furthermore, TMA sampling requirements may vary depending on the target, lesion, tissue, and the goal of investigation. Therefore, it is important to determine the sampling parameters on a per-use basis. For instance, three cores of 0.6 mm diameter will have almost a similar area to one core of 1 mm diameter (0.85 mm2 versus 0.78 mm2), but provide different information about the specimen as they are likely to represent multiple areas [17, 18]. To address this issue, many studies have been performed to determine the impact of size and number of TMA cores [17, 19–28]. Most commonly, the recommended number of TMA cores varied from one to four with a diameter between 0.6 mm to 2 mm.
Determining optimal TMA sampling parameters by physical sampling of the cores, is not only time-consuming, but, more importantly, it limits the options of comprehensive statistical modeling and decomposes the original tissue sample to be used as the reference standard. To overcome these limitations, the concept of a virtual TMA was explored by utilizing digital whole slide images (WSI) to extract artificial TMA cores . The approach has been applied in several studies: Quintayo et al.  manually marked core positions on a low magnification image before acquiring images of the TMA cores at high magnification; they also matched core positions between H&E and IHC staining of the same tissue before the acquired cores were subsequently assembled to a virtual TMA of ductal carcinoma in situ. Pedersen et al.  reported a similar procedure, but used random sampling of six 1 mm diameter cores directly on 20x magnification images of both H&E and IHC slides before assembling the virtual TMA. The studies supported the principle that assembling a set of virtual TMAs by copying cores from digital images is a valuable approach in TMA-based tissue sampling modeling.
A methodology for comprehensive IHC evaluation with appraisal of intratumoral heterogeneity aspects in WSIs of Ki67-stained breast cancer tissue was recently proposed . It is based on systematic subsampling of DIA-generated data into a hexagonal tiling (HexT) arrays and enables computation of a comprehensive set of texture and distribution indicators for Ki67 intratumoral variability. While the primary aim of that study was to investigate intratumoral heterogeneity of Ki67 expression, in the current study we exploit the method for modeling tissue sampling precision in homogeneous and heterogeneous tumors dichotomized by spatial entropy of Ki67 expression: The hexagons in the HexT were chosen to simulate virtual TMA cores (or corresponding fields of view in conventional microscopy), with numbers of Ki67 positive and negative cells established by DIA. Using the spatial entropy extracted from the tiling as a spatial modeling of the Ki67 expression the impact tissue and cell sample size and tumor heterogeneity has on the accuracy of Ki67 LI measurement becomes possible to investigate. We present evidence that tumors with lower Ki67 LI as well as higher spatial heterogeneity of Ki67 expression require relatively larger sampling subsets to represent the global average of the biomarker expression in the tissue. Additionally, the results support the notion, that tumors at the low end of proliferation scale require higher cell counts .
Tissue and data
A data set consisting of primary breast cancer from 297 patients was used in this study. Details of the dataset are reported in . Briefly, 91 % of the tumors were invasive ductal carcinoma of the breast (270/297). Tissue samples were formalin-fixed, processed with standard paraffin embedding techniques. IHC for Ki67 was performed with antibody (clone MIB-1; DAKO, Glostrup, DK) and multimer technology-based detection system (ultraView Universal DAB, Ventana, Tucson, AZ, USA). Digital WSIs were recorded using a ScanScope XT Slide Scanner (Leica Aperio Technologies, Vista, CA, USA) under 20x objective magnification (0.5-μm resolution) and subsequently subjected to DIA by the Leica Aperio Genie Classifier v.1/Nuclear v.9 algorithm. This tool was previously calibrated based on tumor versus benign tissue recognition and positive versus negative cells detection. DIA algorithm was previously validated using a criterion standard achieved by stereological counting. The research was approved by the Vilnius Regional Biomedical Research Ethics committee (reference number NR.:40, date 2007-04-26). Additional informed consent was not required for the use of archived material.
TMA simulation using hexagonal tiling
The HexT methodology forming the basis of automated texture feature extraction is described in detail in . Briefly, the coordinates of positive and negative nuclei extracted by DIA were distributed into a dense HexT overlaid on each WSI. The HexT was randomly positioned within the invasive tumor area (Fig. 1, Middle). Hexagons containing no nuclear profiles by DIA were regarded as missing data; hexagons containing fewer than 100 nuclear profiles were regarded as insufficiently sampled. A minimum requirement of 30 informative hexagons per tumor was applied. Local Ki67 LI was calculated for each hexagon to construct co-occurrence matrix used to compute Haralick texture parameters.
The individual hexagons, with local Ki67 LI, were subsequently used as TMA cores for the random sampling simulations (Fig. 1, Right) and resembled approximately a TMA core of 0.75 mm circular diameter and 0.44 mm2 area. The tumors were dichotomized into homogeneous and heterogeneous groups based on the median entropy value obtained by the HexT methodology. The sampling simulations were carried out for all three tumor classes: all/mixed, homogeneous and heterogeneous.
In addition to giving insight about the minimum number of required TMA cores, the simulations can be used to infer error measurements according to how many nuclei are assessed. By dichotomizing the simulated cores by the number of nuclei contained, the error measurement can additionally be investigated as function of the nuclei count.
The experimental models and statistical methods
The impact caused by varying core number was investigated for a range of numbers feasible to punch out in practice. The chosen set of core numbers investigated is denoted HexN = (1, 2, …, 15).
Two different methods were used to simulate the impact of the number of hexagons/TMA cores on the precision of the sampling to represent the Ki67 LI reported by the DIA of the entire region of interest (ROI). First, the practice of “physical” TMA construction, in which a set of cores is sampled only once, was simulated by randomly sampling a subset of hexagons once. Single linear regression analysis was used to compare the data in a single random selection.
Both experiments were grouped by tumor heterogeneity and repeated for HexN = (1,2, …,15) with hexagons resembling a 0.75 mm diameter TMA core and the simulations were performed with 50,000 iterations.
From the simulations error measurements according to how many nuclei are assessed can be inferred as follows: for one tumor case 50,000 subsets of TMA cores are sampled of size HexN = (1, 2, …, 15). This yields a total of 750,000 subsets which are effectively grouped by HexN. By dichotomizing according to the number of nuclei sampled in each subset into bins of 250 nuclei (first bin [0;250), second bin [250;500) etc.), the error measurement can additionally be investigated as function of the nuclei count. To make it clear if CE is calculated according to hexagon area or nuclei number, the CE is denoted CE Area and CE Nuclei , respectively.
This is done for each choice of HexN, for a set of bins used for dichotomizing by nuclei count and for all three classes of heterogeneity (all/mixed, homogeneous and heterogeneous).
Statistical analysis was performed using R 3.1.2, GNU GCC 5.2.1, Open Office 4.1.2 and SAS 9.4 software.
Extensive dataset summary statistics of the Ki67 indicators, obtained by HexT methodology, are previously reported . Briefly, the global average of Ki67 LI values (in percentages) estimated by DIA of the WSIs was almost identical to the results obtained by HexT (mean: 32.5 ± 16.9 %, median 32.6 ± 17.4 % and sum 32.7 ± 17.3 %). Importantly, the HexT data provided a comprehensive set of intratissue variation indicators .
Single subsampling – linear regression analysis
Linear regression analysis results for hexagon size = 825 pixels (≈0.75 mm TMA core)
All tumor cases
To achieve R2 = 0.95 value in the regression models, random selection of at least four, three and twelve cores were required in the mixed, homogeneous and heterogeneous tumors, respectively.
Fit parameters for relative error CE_Area fitted to proliferation index for all three heterogeneity classes
Proliferation fit to relative error
Fit parameters for relative error CE_Nuclei fitted to proliferation index for all three heterogeneity classes
Proliferation fit to relative error
This study has exploited novel opportunities that digital microscopy images offers for virtual TMA modeling with incorporation of DIA results. Firstly, the virtual TMAs were modeled after the HexT methodology extracted both global texture information and local feature information from the WSI. Secondly, simulation of the TMA cores using the HexT dataset enabled multiple random sampling iterations bypassing the digital assembly of the virtual TMAs. This gave a much greater flexibility in investigating a wider range of sampling methodologies, parameters and error measurements. The added benefits do not impose any new limitations: if cores are needed for several stainings of the same tissue, cores can be sampled at the exact location in different images by applying mapping techniques similar to the ones reported by Quintayo et al. .
Previously, a similar approach was tested by Heus et al. , who utilized a dense grid of rectangular frames instead of hexagons. From each subsampled frame, a core was simulated by the largest circle contained within. This has a side-effect that tissue located at the corners of the frames will never be sampled; this effect is not independent of size of the simulated cores. The use of hexagons for virtual core simulation does not suffer from this: the dense HexT ensures that all parts of the tissue are considered with the same probability. Sampling without replacement further ensures that the same area is not represented by multiple cores in the subsets used in the simulations.
The analysis of core/cell sampling requirements in this study was made possible to group according to Haralick entropy texture feature extracted by the HexT methodology for each WSI. It must be noted that the Haralick entropy threshold value is not clearly defined. Therefore, the optimal method to split the dataset it into equal parts by median was chosen. In a similar study, the variance of the local Ki67 LI was used as entropy measurement, but without a complete error analysis for the entire dataset .
Combining Ki67 LI (or any other biomarker) from several cores is often needed in TMA studies. This introduces a risk of bias which involves assessing the number of positive and negative nuclei for the observed cores before recalculating the Ki67 LI. We evaluated this potential bias by comparing results from combined Ki67 LI from a set of simulated cores using the mean, the median and Ki67 LI calculated by using sum of nuclei in the sampled hexagons. We found that when larger set of cores were used, any bias with regard to the Ki67 LI calculation methods was negligible (Table 1), while calculation of the combined Ki67 LI, by assessing the core data first, is strongly advised where only a few TMA cores are used from heterogeneous tumors (Fig. 2, right).
The practical TMA construction, where cores were randomly chosen only once, was investigated using linear regression. This allowed comparison of the hexagonal simulation data to previous studies. For a tumor set with mixed heterogeneity, we found a number of cores to achieve R^2 = 0.95, to be four, in line with the previous reports [17, 25, 26, 31]. For homogeneous tumors, the optimal number of cores was three, depending whether the sum or mean calculations were used, respectively. A rather dramatic increase to the requirement of 12 cores was found in the heterogeneous tumors.
The single sampling brings some eventuality, because for each sampling, it is possible to obtain cores containing different tissue representation and thus biomarker expression level. The error is particularly important when considering tissue samples with varying degrees of heterogeneity, as it influences the representativeness of TMAs . A number of studies have shown that more cores will improve the agreement level and reduce the limitations due to the heterogeneity in various types of tumors and IHC biomarkers [25–27, 32, 33]. However, only the introduced method allows inference about the relative error caused by different TMA sampling parameters in combination with tumor heterogeneity. Our results show that to obtain a CE of 10 % the necessary number of cores in the dataset with mixed heterogeneity is 8; five cores hold sufficient information for Ki67 LI determination in homogeneous tumors, while heterogeneous tumors need at least 11–12 cores to be sampled.
In practical TMA applications, intratumoral heterogeneity of biomarker expression is usually unknown in advance; therefore, a more conservative approach would assume that all tumors in the study population are heterogeneous. On the other hand, Ki67 LI expression in breast cancer tissue is known for its spatial heterogeneity and may serve a reference standard for other biomarkers and tumors. In that sense, our study reveals that 11–12 random TMA cores of 0.75 mm diameter would sufficiently represent IHC biomarker expression in heterogeneous tumors. Our simulations also indicate that disagreements between different studies of TMA core numbers may in fact be due to unestablished differences in heterogeneity aspects. In general, our findings support the notion that heterogeneity information is crucial for optimizing TMA studies. Ideally, the presented method could be used in pilot studies to validate the optimal number of cores, or at least heterogeneity should be investigated from a larger set of cores, for instance by measuring a range of Ki67 LI between several TMA cores taken from the tissue.
Our study also provides evidence for minimum cell counting requirements to achieve robust Ki67 LI measurement, especially with regard to the limited capacity of manual counting procedures. Current clinical guidelines on the minimal number of cells to be counted are quite arbitrary, mostly set in the range of 500 and 2000 tumor cells . While small samples (e.g., needle core biopsies) may allow counting all the invasive tumor cells, it becomes impractical in larger samples. Therefore, to achieve adequate precision, it is recommended for the interpreting pathologist to score at least 1,000 cells, while 500 cells would be acceptable as the absolute minimum . Importantly, our findings reveal that to achieve 10 % CE approximately 4,000 nuclei must be counted when the intratumor heterogeneity is mixed/unknown (Fig. 4). These cell counts are rather large to accomplish in clinical practice for all breast carcinomas, but could be feasible for cases considered as “grey zone”, e.g. in the range of Ki67 LI 10-30 % . A visual scoring methodology proposed by Hida et al., might be used as method of choice for “low” (Ki67 LI <10 %) or “high” (Ki67 LI >30 %) proliferatively active cases, leaving behind “grey zone” cases, which requires more precise methodologies .
The inverse relationship between relative estimation error and mitotic activity previously highlighted by Going  was confirmed to also exist between each of the two error estimates (CE Area , CE Nuclei ) and the Ki67 LI proliferation activity indicator (Figs. 5 and 6). This dependency of CE on Ki67 LI shows that tumor cases with low proliferation rate contribute most of the CE in Fig. 4 which is a set of “mixed” proliferation rate. Consequently, when scoring a single case with unknown Ki67 LI one may need to evaluate a higher cell count or larger TMA sample to ensure a 10 % CE at a specific grey zone. Specifically, the tumors at the lower scale of proliferative activity (Ki67 LI < 20 %, Fig. 5, left) will for a mixed/unknown heterogeneity case require larger sampling (at least 10–11 TMA cores) to achieve the same error measurement (10 % CE Area ) results as for highly proliferative tumors (4–6 TMA cores). Similarly, for cases with Ki67 LI < 20 % (Fig. 5, right) at least 6,250 nuclei are necessary (for 10 % CE Nuclei ) As such, Figs. 5 and 6 may aid determining practical sampling requirements of individual cases for acceptable CE at specific grey zones.
In general, the results of our study suggest that adequate accuracy levels of Ki67 LI measurement can hardly be achieved by manual counts and argue in favor of DIA-based techniques to benefit from the high-capacity methods. In addition, automated hotspot detection with standard definitions by DIA, which was out of scope in the present study, would provide another advantage compared to the visual evaluation by conventional microscopy or inspection of WSI.
Several aspects raised in this study relate to the evaluation of Ki67 immunohistochemistry in breast cancer in clinical research and practice. Firstly, obtaining an optimal number of TMA cores/cell number needed for biomarker research studies depends on the tissue, especially its intratissue heterogeneity and level of expression. For Ki67 LI in breast cancer, we found 5–6 cores sufficient for homogeneous expression in the tissue, 8 cores for tumors with mixed heterogeneity and at least 11 cores for heterogeneous tumors. Secondly, our findings reveal that to achieve low error estimates when evaluating by cell counting, approximately 4,000 nuclei must be evaluated when the intratumor heterogeneity is mixed/unknown. In breast cancer cases of the lower proliferative activity (Ki67 LI < 20 %) larger sampling is required to achieve the same error measurement results as for highly proliferative tumors. The presented data may aid in defining practical sampling requirements of individual cases and specific grey zones.
The wide range of the number of cores/nuclei needed supports the notion that optimal sampling requirements must be determined on a peruse basis and that heterogeneity information must be assessed in the study. The method presented can be applied for individual pilot study measurements. In addition, our findings highlight the importance of high-capacity computer-based IHC measurement techniques to improve accuracy of the testing.
CE, coefficient of error; CE Area , coefficient of error calculated according to hexagon area; CE Nuclei :coefficient of error calculated according to nuclei number; DIA, digital image analysis; ER, estrogen receptor; H&E, hematoxylin and eozin staining; HexN, set of hexagonal cores; HexT, hexagonal tiling; IHC, immunohistochemistry; Ki67LI, Ki67 labeling index; PgR, progesterone receptor; TMA, tissue microarray; TMAs, tissue microarrays; WSI, whole slide images
This research is fully funded by the European Social Fund under the Global Grant measure, Grant #VP1-3.1-SMM-07-K-03-051.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article and its Additional file 2.
BP, AR constructed and performed TMA simulation procedures and statistical analysis systems. AL, AR, BP, RA, and PH performed statistical analyses. JB, in collaboration with AR and AL, drafted essential parts of the manuscript. AG, AiL, AR, IE reviewed and edited the manuscript. All authors participated in conception and design of the study, reviewing the analysis results, critically revised and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The research was approved by the Vilnius Regional Biomedical Research Ethics committee (reference number NR.:40, date 2007-04-26). Additional informed consent was not required for the use of archived material.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Knutsvik G, Stefansson IM, Aziz S, Arnes J, Eide J, Collett K, Akslen LA. Evaluation of Ki67 expression across distinct categories of breast cancer specimens: a population-based study of matched surgical specimens, core needle biopsies and tissue microarrays. PLoS One. 2014;9(11), e112121. doi:10.1371/journal.pone.0112121.View ArticlePubMedPubMed CentralGoogle Scholar
- Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thurlimann B, Senn HJ. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Ann Oncol. 2013;24(9):2206–23. doi:10.1093/annonc/mdt303.View ArticlePubMedPubMed CentralGoogle Scholar
- Untch M, Harbeck N, Huober J, von Minckwitz G, Gerber B, Kreipe HH, Liedtke C, Marschner N, Mobus V, Scheithauer H, Schneeweiss A, Thomssen C, Jackisch C, Beckmann MW, Blohmer JU, Costa SD, Decker T, Diel I, Fasching PA, Fehm T, Janni W, Luck HJ, Maass N, Scharl A, Loibl S. Primary Therapy of Patients with Early Breast Cancer: Evidence, Controversies, Consensus: Opinions of German Specialists to the 14th St. Gallen International Breast Cancer Conference 2015 (Vienna 2015). Geburtshilfe Frauenheilkd. 2015;75(6):556–65. doi:10.1055/s-0035-1546120.View ArticlePubMedPubMed CentralGoogle Scholar
- Romero Q, Bendahl PO, Ferno M, Grabau D, Borgquist S. A novel model for Ki67 assessment in breast cancer. Diagn Pathol. 2014;9:118. doi:10.1186/1746-1596-9-118.View ArticlePubMedPubMed CentralGoogle Scholar
- Lu H, Papathomas TG, van Zessen D, Palli I, de Krijger RR, van der Spek PJ, Dinjens W, Stubbs AP. Automated Selection of Hotspots (ASH): enhanced automated segmentation and adaptive step finding for Ki67 hotspot detection in adrenal cortical cancer. Diagn Pathol. 2014;9(1):216. doi:10.1186/s13000-014-0216-6.View ArticlePubMedPubMed CentralGoogle Scholar
- Potts SJ, Krueger JS, Landis ND, Eberhard DA, Young GD, Schmechel SC, Lange H. Evaluating tumor heterogeneity in immunohistochemistry-stained breast cancer tissue. Lab Invest. 2012;92(9):1342–57. doi:10.1038/labinvest.2012.91.View ArticlePubMedGoogle Scholar
- Haroske G, Dimmer V, Steindorf D, Schilling U, Theissig F, Kunze KD. Cellular sociology of proliferating tumor cells in invasive ductal breast cancer. Anal Quant Cytol Histol. 1996;18(3):191–8.PubMedGoogle Scholar
- Gudlaugsson E, Skaland I, Janssen EA, Smaaland R, Shao Z, Malpica A, Voorhorst F, Baak JP. Comparison of the effect of different techniques for measurement of Ki67 proliferation on reproducibility and prognosis prediction accuracy in breast cancer. Histopathology. 2012;61(6):1134–44. doi:10.1111/j.1365-2559.2012.04329.x.View ArticlePubMedGoogle Scholar
- Dowsett M, Nielsen TO, A'Hern R, Bartlett J, Coombes RC, Cuzick J, Ellis M, Henry NL, Hugh JC, Lively T, McShane L, Paik S, Penault-Llorca F, Prudkin L, Regan M, Salter J, Sotiriou C, Smith IE, Viale G, Zujewski JA, Hayes DF, International Ki-67 in Breast Cancer Working G. Assessment of Ki67 in breast cancer: recommendations from the International Ki67 in Breast Cancer working group. J Natl Cancer Inst. 2011;103(22):1656–64. doi:10.1093/jnci/djr393.View ArticlePubMedPubMed CentralGoogle Scholar
- Laurinavicius A, Plancoulaine B, Laurinaviciene A, Herlin P, Meskauskas R, Baltrusaityte I, Besusparis J, Dasevicius D, Elie N, Iqbal Y, Bor C. A methodology to ensure and improve accuracy of Ki67 labelling index estimation by automated digital image analysis in breast cancer tissue. Breast Cancer Res. 2014;16(2):R35. doi:10.1186/bcr3639.View ArticlePubMedPubMed CentralGoogle Scholar
- Laurinavicius A, Plancoulaine B, Laurinaviciene A, Herlin P, Meskauskas R, Baltrusaityte I, Besusparis J, Elie N, Belhomme P, Iqbal Y, Bor-Angelier C (2013) A methodology to ensure and improve accuracy of Ki67 digital immunohistochemistry analysis in breast cancer tissue. Molecular Cancer Research 11. doi:10.1158/1557-3125.ADVBC-B116
- Going JJ. Techniques of mitosis counting. Hum Pathol. 1993;24(1):113–4.View ArticlePubMedGoogle Scholar
- Battifora H. The multitumor (sausage) tissue block: novel method for immunohistochemical antibody testing. Lab Invest. 1986;55(2):244–8.PubMedGoogle Scholar
- Permuth-Wey J, Boulware D, Valkov N, Livingston S, Nicosia S, Lee JH, Sutphen R, Schildkraut J, Narod S, Parker A, Coppola D, Sellers T, Pal T. Sampling Strategies for Tissue Microarrays to evaluate biomarkers in Ovarian Cancer. Cancer Epidemiol Biomarkers Prev. 2009;18(1):28–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Kononen J, Bubendorf L, Kallioniemi A, Barlund M, Schraml P, Leighton S, Torhorst J, Mihatsch MJ, Sauter G, Kallioniemi OP. Tissue microarrays for high-throughput molecular profiling of tumor specimens. Nat Med. 1998;4(7):844–7.View ArticlePubMedGoogle Scholar
- Kallioniemi O-P, Wagner U, Kononen J, Sauter G. Tissue microarray technology for high-throughput molecular profiling of cancer. Hum Mol Genet. 2001;10(7):657–62. doi:10.1093/hmg/10.7.657.View ArticlePubMedGoogle Scholar
- Ilyas M, Grabsch H, Ellis IO, Womack C, Brown R, Berney D, Fennell D, Salto-Tellez M, Jenkins M, Landberg G, Byers R, Treanor D, Harrison D, Green AR, Ball G, Hamilton P. Guidelines and considerations for conducting experiments using tissue microarrays. Histopathology. 2013;62(6):827–39. doi:10.1111/his.12118.View ArticlePubMedGoogle Scholar
- Anagnostou VK, Lowery FJ, Syrigos KN, Cagle PT, Rimm DL. Quantitative evaluation of protein expression as a function of tissue microarray core diameter: is a large (1.5 mm) core better than a small (0.6 mm) core? Arch Pathol Lab Med. 2010;134(4):613–9. doi:10.1043/1543-2165-134.4.613.PubMedGoogle Scholar
- Quintayo MA, Starczynski J, Yan FJ, Wedad H, Nofech-Mozes S, Rakovitch E, Bartlett JM. Virtual tissue microarrays: a novel and viable approach to optimizing tissue microarrays for biomarker research applied to ductal carcinoma in situ. Histopathology. 2013. doi:10.1111/his.12336.Google Scholar
- Schmidt LH, Biesterfeld S, Kummel A, Faldum A, Sebastian M, Taube C, Buhll R, Wiewrodt R. Tissue microarrays are reliable tools for the clinicopathological characterization of lung cancer tissue. Anticancer Res. 2009;29(1):201–9.PubMedGoogle Scholar
- Zhang D, Salto-Tellez M, Putti TC, Do E, Koay ES. Reliability of tissue microarrays in detecting protein expression and gene amplification in breast cancer. Mod Pathol. 2003;16(1):79–84. doi:10.1097/01.MP.0000047307.96344.93.View ArticlePubMedGoogle Scholar
- Torhorst J, Bucher C, Kononen J, Haas P, Zuber M, Kochli OR, Mross F, Dieterich H, Moch H, Mihatsch M, Kallioniemi OP, Sauter G. Tissue microarrays for rapid linking of molecular changes to clinical endpoints. Am J Pathol. 2001;159(6):2249–56. doi:10.1016/S0002-9440(10)63075-1.View ArticlePubMedPubMed CentralGoogle Scholar
- Mucci NR, Akdas G, Manely S, Rubin MA. Neuroendocrine expression in metastatic prostate cancer: evaluation of high throughput tissue microarrays to detect heterogeneous protein expression. Hum Pathol. 2000;31(4):406–14. doi:10.1053/hp.2000.7295.View ArticlePubMedGoogle Scholar
- Camp RL, Charette LA, Rimm DL. Validation of tissue microarray technology in breast carcinoma. Lab Invest. 2000;80(12):1943–9.View ArticlePubMedGoogle Scholar
- Goethals L, Perneel C, Debucquoy A, De Schutter H, Borghys D, Ectors N, Geboes K, McBride WH, Haustermans KM. A new approach to the validation of tissue microarrays. J Pathol. 2006;208(5):607–14. doi:10.1002/path.1934.View ArticlePubMedGoogle Scholar
- Gulmann C, Butler D, Kay E, Grace A, Leader M. Biopsy of a biopsy: validation of immunoprofiling in gastric cancer biopsy tissue microarrays. Histopathology. 2003;42(1):70–6.View ArticlePubMedGoogle Scholar
- Alkushi A. Validation of tissue microarray biomarker expression of breast carcinomas in Saudi women. Hematol Oncol Stem Cell Ther. 2009;2(3):394–8.View ArticlePubMedGoogle Scholar
- Pedersen MB, Riber-Hansen R, Nielsen PS, Bendix K, Hamilton-Dutoit SJ, D'Amore F, Steiniche T. Digital pathology for the validation of tissue microarrays in peripheral T-cell lymphomas. Appl Immunohistochem Mol Morphol. 2014;22(8):577–84. doi:10.1097/PAI.0b013e3182a7d16d.View ArticlePubMedGoogle Scholar
- Plancoulaine B, Laurinaviciene A, Herlin P, Besusparis J, Meskauskas R, Baltrusaityte I, Iqbal Y, Laurinavicius A. A methodology for comprehensive breast cancer Ki67 labeling index with intra-tumor heterogeneity appraisal based on hexagonal tiling of digital image analysis data. Virchows Arch. 2015;467(6):711–22. doi:10.1007/s00428-015-1865-x.View ArticleGoogle Scholar
- Heus R (2009) Approches virtuelles dédiées à la technologie des puces à tissus "Tissue MicroArrays " TMA : Application à l'étude de la transformation tumorale du tissu colorectal. PhD thesis, Life Sciences, Université Joseph-Fourier - Grenoble I French https://tel.archives-ouvertes.fr/tel-00429056
- Rubin MA, Dunn R, Strawderman M, Pienta KJ. Tissue microarray sampling strategy for prostate cancer biomarker analysis. Am J Surg Pathol. 2002;26(3):312–9.View ArticlePubMedGoogle Scholar
- Lax SF, Pizer ES, Ronnett BM, Kurman RJ. Comparison of estrogen and progesterone receptor, Ki-67, and p53 immunoreactivity in uterine endometrioid carcinoma and endometrioid carcinoma with squamous, mucinous, secretory, and ciliated cell differentiation. Hum Pathol. 1998;29(9):924–31.View ArticlePubMedGoogle Scholar
- Ruiz C, Seibt S, Al KK, Siraj AK, Mirlacher M, Schraml P, Maurer R, Spichtin H, Torhorst J, Popovska S, Simon R, Sauter G. Tissue microarrays for comparing molecular features with proliferation activity in breast cancer. IntJCancer. 2006;118(9):2190–4.Google Scholar
- Hida AI, Oshiro Y, Inoue H, Kawaguchi H, Yamashita N, Moriya T. Visual assessment of Ki67 at a glance is an easy method to exclude many luminal-type breast cancers from counting 1000 cells. Breast Cancer. 2015;22(2):129–34. doi:10.1007/s12282-013-0460-8.View ArticlePubMedGoogle Scholar