Acessibilidade / Reportar erro

Nuclear entropy, angular second moment, variance and texture correlation of thymus cortical and medullar lymphocytes: Grey level co-occurrence matrix analysis

Abstracts

Grey level co-occurrence matrix analysis (GLCM) is a well-known mathematical method for quantification of cell and tissue textural properties, such as homogeneity, complexity and level of disorder. Recently, it was demonstrated that this method is capable of evaluating fine structural changes in nuclear structure that otherwise are undetectable during standard microscopy analysis. In this article, we present the results indicating that entropy, angular second moment, variance, and texture correlation of lymphocyte nuclear structure determined by GLCM method are different in thymus cortex when compared to medulla. A total of 300 thymus lymphocyte nuclei from 10 one-month-old mice were analyzed: 150 nuclei from cortex and 150 nuclei from medullar regions of thymus. Nuclear GLCM analysis was carried out using National Institutes of Health ImageJ software. For each nucleus, entropy, angular second moment, variance and texture correlation were determined. Cortical lymphocytes had significantly higher chromatin angular second moment (p < 0.001) and texture correlation (p < 0.05) compared to medullar lymphocytes. Nuclear GLCM entropy and variance of cortical lymphocytes were on the other hand significantly lower than in medullar lymphocytes (p < 0.001). These results suggest that GLCM as a method might have a certain potential in detecting discrete changes in nuclear structure associated with lymphocyte migration and maturation in thymus.

texture analysis; GLCM matrix; CD4; CD8; gene


Análise da matriz de co-ocorrência de níveis de cinza (GLCM) é um método matemático bem conhecido para a quantificação das propriedades estruturais das células e tecidos, incluindo homogeneidade, complexidade e nível de desordem. Recentemente foi demonstrado que esse método é capaz de avaliar as mudanças estruturais finas do núcleo que, de outra forma não são detectáveis pela análise microscópica habitual. Nesse artigo, apresentamos resultados indicando que a entropia, segundo momento angular, variância, e correlação de textura da estrutura nuclear de linfócitos, determinada pelo método de GLCM, diferem na córtex do timo em comparação a medula. Um total de 300 núcleos de linfócitos do timo de 10 camundongos de um mês de idade foram analisados: 150 núcleos da córtex e 150 núcleos da medula do timo. Análise de GLCM nuclear foi realizada usando o software ImageJ do National Institutes of Health (NIH). Para cada núcleo, entropia, segundo momento angular, variância e correlação de textura foram determinadas. Linfócitos corticais continham significativamente maiores os parâmetros segundo momento angular da cromatina (p<0.001) e correlação de textura (p<0.05) comparados a linfócitos medulares. A entropia GLCM nuclear e variância de linfócitos corticais eram por outro lado significativamente menores que nos linfócitos medulares (p<0.001). Esses resultados sugerem que GLCM como método tenha um potencial de detector mudanças discretas na estrutura nuclear associadas a migração e maturação de linfócitos no timo.

análise de textura; matriz GLCM; CD4; CD8; gene


INTRODUCTION

During the past several decades, there have been numerous research efforts to design and implement a simple, accurate, and affordable computational biology method for quantification and evaluation of nuclear architecture. Many of the developed techniques include image analysis of cell cultures and stained tissue micrographs, and they combine standard morphometric measurements with second and higher order statistical parameters. One of the recently applied techniques that include second order statistical texture analysis is the so-called grey level co-occurrence matrix method (GLCM). GLCM is a rapid mathematical method for assessing image structural properties, such as homogeneity, complexity and level of disorder. It was introduced in the work of Haralick et al. (1973)Haralick RM, Shanmugam K and Dinstein I. 1973. Textural Features for Image Classification. IEEE Transactions on Systems Man Cybernetics 3: 610-621., who have established a total of 14 parameters based on pixel pair grey level occurrences within a grey scale image. In time, GLCM was successfully applied in nuclear magnetic resonance imaging (Li et al. 2009Li X, Pai A, Blumenkrantz G, Carballido-Gamio J, Link T, Ma B, Ries M and Majumdar S. 2009. Spatial distribution and relationship of T1 rho and T2 relaxation times in knee cartilage with osteoarthritis. Magn Reson Med 61: 1310-1318.), computed tomography (Huber et al. 2011Huber MB, Nagarajan MB, Leinsinger G, Eibel R, Ray LA and WismÜller A. 2011. Performance of topological texture features to classify fibrotic interstitial lung disease patterns. Med Phys 38: 2035-2044) and other clinical research areas. In fundamental medical and biology research, Shamir et al. demonstrated that tissue GLCM entropy is an indicator of tissue age-related structural degradation (Shamir et al. 2009Shamir L, Wolkow CA and Goldberg IG. 2009. Quantitative measurement of aging using image texture entropy. Bioinformatics 25: 3060-3063.). Losa et al. suggested that in cell cultures, GLCM parameters of nuclear structure may be important for evaluation of chromatin structural changes during apoptosis (Losa and Castelli 2005Losa GA and Castelli C. 2005. Nuclear patterns of human breast cancer cells during apoptosis: characterisation by fractal dimension and co-occurrence matrix statistics. Cell Tissue Res 322: 257-267.).

Today, of all GLCM parameters, entropy, angular second moment (ASM), variance and texture correlation are perhaps most commonly used in experimental and clinical medicine (Shamir et al. 2009Shamir L, Wolkow CA and Goldberg IG. 2009. Quantitative measurement of aging using image texture entropy. Bioinformatics 25: 3060-3063., Pantic and Pantic 2012Pantic I and Pantic S. 2012. Germinal center texture entropy as possible indicator of humoral immune response: immunophysiology viewpoint. Mol Imaging Biol 14: 534-540., Pantic et al. 2012aPantic I, Pantic S and Basta-Jovanovic G. 2012a. Grey level co-occurrence matrix (GLCM) texture analysis of germinal center light zone lymphocyte nuclei: physiology viewpoint with focus on apoptosis. Microsc Microanal 18: 470-475., Losa and Castelli 2005Losa GA and Castelli C. 2005. Nuclear patterns of human breast cancer cells during apoptosis: characterisation by fractal dimension and co-occurrence matrix statistics. Cell Tissue Res 322: 257-267., Li et al. 2009Li X, Pai A, Blumenkrantz G, Carballido-Gamio J, Link T, Ma B, Ries M and Majumdar S. 2009. Spatial distribution and relationship of T1 rho and T2 relaxation times in knee cartilage with osteoarthritis. Magn Reson Med 61: 1310-1318., Alvarenga et al. 2010Alvarenga AV, Teixeira CA, Ruano MG and Pereira WC. 2010. Influence of temperature variations on the entropy and correlation of the Grey-Level Co-occurrence Matrix from B-Mode images. Ultrasonics 50: 290-293.). In recent research carried out by our laboratory, we have suggested that GLCM entropy and angular second moment can be successfully used to evaluate structural properties of a secondary lymphoid organ. In immunized guinea-pigs, entropy of spleen germinal centers was negatively correlated with humoral immune response, indicating that entropy may be a good predictor of structural deterioration and degradation that occurs in spleen germinal center tissue (Pantic and Pantic 2012Pantic I and Pantic S. 2012. Germinal center texture entropy as possible indicator of humoral immune response: immunophysiology viewpoint. Mol Imaging Biol 14: 534-540.). Also, in another study, we found that in germinal center light zones (where the process of B-lymphocyte selection takes place), entropy and angular second moment of lymphocyte nuclear structures are related to the nuclear size (Pantic et al. 2012aPantic I, Pantic S and Basta-Jovanovic G. 2012a. Grey level co-occurrence matrix (GLCM) texture analysis of germinal center light zone lymphocyte nuclei: physiology viewpoint with focus on apoptosis. Microsc Microanal 18: 470-475.). These results suggested that nuclear GLCM parameters may be good indicators of structural changes that follow the shrinkage of nuclei during early stages of programmed cell death (apoptosis) in spleen white pulp.

Unlike secondary lymphoid organs, there are virtually no data regarding primary lymphoid organ nuclear GLCM properties. In this article, we present the results indicating that lymphocytes located in thymus cortex have significantly different nuclear entropy, ASM and correlation, when compared to the lymphocytes located in thymus medulla. These findings suggest that GLCM as a method might be capable of detecting fine structural changes in cell nuclei associated with migration and maturation of thymus lymphocytes.

MATERIALS AND METHODS

Thymus tissue was obtained from 10 male one-month-old Swiss white albino outbred mice. All animals were previously bread and kept under same laboratory and environmental conditions (standard laboratory diet, ambient temperature of 22°C, natural light cycle). The research followed the rules presented in the Guide for the Care and Use of Laboratory Animals, published by the U.S. National Institute of Health (NIH Publication No. 85 – 23, revised 1985). The principal investigator (IP) had an authorization from the Ethical Commission of the University of Belgrade, School of Medicine, and the Serbian Ministry of agriculture for the experimental work with the laboratory animals (Approval No. 323-07-03985/2012-05/1).

Thymus tissue was fixated in Carnoy solution, mounted on glass slides (5 micrometer thickness), and stained with hematoxylin/eosin (H/E). An Olympus C-5060 Wide Zoom digital camera instrument with an oil immersion objective (×1,000 magnification) of Olympus BX41 microscope (1 pixel equaled 0.048 µm) was used for micrograph acquisition similarly to our previous study (Pantic et al. 2012aPantic I, Pantic S and Basta-Jovanovic G. 2012a. Grey level co-occurrence matrix (GLCM) texture analysis of germinal center light zone lymphocyte nuclei: physiology viewpoint with focus on apoptosis. Microsc Microanal 18: 470-475.). A total of 300 thymus lymphocyte nuclei (30 randomly selected non-overlapping lymphocyte nuclei from each animal) were analyzed: 150 from cortex (15 per animal) and 150 from medullar regions (15 per animal) of thymus (figure 1).

Figure 1
Digital micrographs (1,000x magnification) of thymus cortex (A) and medulla (B).

Nuclear image acquisition was done after creating an automatic threshold of the nuclear boundaries in 8-bit format using National Institutes of Health ImageJ (version 1.44p) software. The following protocol was used for each nucleus. The digital micrographs of thymus cortex/medulla in tiff format were converted to 8-bit (ImageJ>Type>8-bit). The converted 8-bit micrographs were then processed using the function ImageJ>Adjust>Threshold. During the visualization of the nuclei in ImageJ software, two independent histology experts determined the values of the upper and lower threshold level settings, which in their opinion, provided the most accurate quality of the segmentation. Upper threshold level was set to 110 while lower threshold level was 0. Although this visual evaluation of automatic image processing is subjective in nature, we assume that it has not significantly influenced the final results, as the values of GLCM parameters depend much more on textural patterns of the interior of the nuclei, than on precision of nuclear boundary determination. The sample of non-overlapping lymphocyte nuclei was selected by [Wand tracing tool] (Mode: Legacy, Tolerance:0.0) and cropped (ImageJ>Image>Crop). After the cropping process, individual nuclei images were analyzed by GLCM method.

GLCM method was first introduced in 1973, by Haralick et al. who developed a method able to quantify spatial relationship between grey tones (pixel brightness values) in an image. This method uses [second order] statistics, which means that it estimates the relationship between groups of two resolution cells (pixels) separated by a defined distance d. GLCM analysis in our study was carried out by ImageJ and its texture analysis plugins developed by Julio E. Cabrera, and Toby C. Cornish.

Grey level co-occurrence matrix entropy, angular second moment, variance and correlation for each nucleus were calculated based on the following formulas:

where i and j are coordinates of the co-occurrence matrix, whereas σ and µ represent means and standard deviations of px and py for a selected distance d and angle θ (Haralick et al. 1973Haralick RM, Shanmugam K and Dinstein I. 1973. Textural Features for Image Classification. IEEE Transactions on Systems Man Cybernetics 3: 610-621.).

GLCM parameters were determined separately for direction angles θ equal to 0°, 45°, 90° and 135°, as well as distances d = 1, d = 2 and d = 3, after which the average values were calculated similarly to the work of Losa and Castelli (2005)Losa GA and Castelli C. 2005. Nuclear patterns of human breast cancer cells during apoptosis: characterisation by fractal dimension and co-occurrence matrix statistics. Cell Tissue Res 322: 257-267.. This way, a total of 3,600 individual GLCM measurements were calculated. For details regarding the mathematical protocol of the study, the reader is referred to our previous articles on this subject (Pantic and Pantic 2012Pantic I and Pantic S. 2012. Germinal center texture entropy as possible indicator of humoral immune response: immunophysiology viewpoint. Mol Imaging Biol 14: 534-540., Pantic et al. 2012aPantic I, Pantic S and Basta-Jovanovic G. 2012a. Grey level co-occurrence matrix (GLCM) texture analysis of germinal center light zone lymphocyte nuclei: physiology viewpoint with focus on apoptosis. Microsc Microanal 18: 470-475.), the work of Haralick et al. (1973)Haralick RM, Shanmugam K and Dinstein I. 1973. Textural Features for Image Classification. IEEE Transactions on Systems Man Cybernetics 3: 610-621., as well as technical data of National Institutes of Health ImageJ software.

Statistical analysis was done using GraphPad Prism software (La Jolla, CA, USA) and SPSS statistical package (v. 10; SPSS Inc., Chicago, IL, USA). Student's t-test was used to determine possible statistically significant difference between cortical and medullar lymphocytes.

RESULTS

Average values of GLCM entropy, angular second moment and correlation for different GLCM distances d are presented in table I. All parameters presented statistically significant difference.

TABLE I
Average values of chromatin entropy, angular second moment, variance and correlation for different GLCM distances (d).

For distance d = 1, mean entropy was 6.389 ± 0.253 for cortical lymphocytes and 6.510 ± 0.280 for medullar lymphocytes. The difference between the two lymphocyte populations was statistically highly significant (p < 0.001, figure 2). Average angular second moment for cortical lymphocytes was 0.0021 ± 0.0006, and it was significantly higher than the ASM of medullar lymphocytes (0.0019 ± 0.0007, p < 0.001, figure 2). Mean GLCM variance in cortical and medullar lymphocytes were 174.3 ± 102.8 and 216.4 ± 108.5, respectively (p < 0.001). Mean texture correlation of cortical lymphocytes was significantly higher than correlation of medullar lymphocytes (0.0065 ± 0.0030 versus 0.0055 ± 0.0038, p < 0.05, figure 2).

Figure 2
GLCM entropy (a), angular second moment (b), correlation (c) and variance (d) of thymus cortical and medullar lymphocyte nuclei.

For distance d = 2, mean entropy was 6.537 ± 0.264 for cortical lymphocytes and 6.675 ± 0.294 for medullar lymphocytes. The entropy of medullar lymphocytes was significantly higher (p < 0.001, figure 2). Similar results were found for textural variance (p < 0.001; 170.8 ± 98.1 versus 212.1 ± 104.3). Cortical lymphocyte population had higher ASM (0.0018 ± 0.0005 versus 0.0016 ± 0.0006, p < 0.001, figure 2) and texture correlation (0.0058 ± 0.0025 versus 0.0050 ± 0.0033, p < 0.05, figure 2).

For distance d = 3, mean ASM was 0.0016 ± 0.0005 for cortical lymphocytes and 0.0014 ± 0.0006 for medullar lymphocytes. The difference between the two lymphocyte populations was statistically highly significant (p < 0.001, figure 2). Cortical lymphocyte population had lower entropy (6.611 ± 0.270 versus 6.756 ± 0.301, p < 0.001, figure 2) and variance (170.8 ± 96.4 versus 211.4 ± 102.9) as well as higher texture correlation (0.0048 ± 0.0019 versus 0.0043 ± 0.0027, p < 0.05, figure 2).

Of all determined parameters, entropy in both thymus regions showed the highest degree of homogeneity with inter and intra-individual variability lower than 10%. There was no significant difference in GLCM parameters between individual animals (ANOVA, p > 0.05).

For each GLCM parameter a receiver operating characteristic (ROC) curve was constructed to illustrate the performance of the parameter as a binary classifier of lymphocyte nuclei in terms of their location in thymus. As presented in figure 3, the areas under the ROC curves were 0.652 for entropy, 0.649 for angular second moment, 0.626 for correlation and 0.637 for variance.

Figure 3
Receiver operating characteristic (ROC) curves for GLCM parameters.

DISCUSSION

The results of our study suggest that lymphocytes of thymus medulla might have higher nuclear GLCM entropy and variance than lymphocytes in thymus cortex. Angular second moment and texture correlation are on the other hand lower in medullar lymphocytes. In biology research, these GLCM parameters are generally seen as indicators of cell/tissue structural complexity and heterogeneity. High entropy and variance usually implies elevated level of disorder and disorganization, whereas angular second moment is a measure of homogeneity (Shamir et al. 2009Shamir L, Wolkow CA and Goldberg IG. 2009. Quantitative measurement of aging using image texture entropy. Bioinformatics 25: 3060-3063., Haralick et al. 1973Haralick RM, Shanmugam K and Dinstein I. 1973. Textural Features for Image Classification. IEEE Transactions on Systems Man Cybernetics 3: 610-621.). Texture correlation is a measure of grey level linear dependencies in a micrograph (Haralick et al. 1973Haralick RM, Shanmugam K and Dinstein I. 1973. Textural Features for Image Classification. IEEE Transactions on Systems Man Cybernetics 3: 610-621.) and can have values from 0 (uncorrelated micrograph) to +1 and -1 (perfectly correlated). GLCM parameters are today commonly used for quantification of cell/tissue texture, and for detection of fine cell/tissue structural changes that are otherwise undetectable during standard microscopy analysis.

Thymus is a primary lymphoid organ with various biological functions regarding immune system development and its ability to respond to infections. Cortex and medulla are two major thymus histological compartments. Both compartments consist primarily of T-lymphocytes (T-cells), cells which are essential for adaptive immune responses. One of the main events that take place in thymus environment is selection and maturation of T-cells. In cortex, these cells undergo the so-called [T-cell receptor gene rearrangement], a process involving random joining of gene segments in nucleus in order to form a genetic sequence responsible for creation of cell surface T-cell receptor (TCR), which is an essential property of a functional T-lymphocyte antigen recognition ability (Sleckman 2005Sleckman BP. 2005. Lymphocyte antigen receptor gene assembly: multiple layers of regulation. Immunol Res 32: 253-258., Nitta et al. 2008Nitta T, Murata S, Ueno T, Tanaka K and Takahama Y. 2008. Thymic microenvironments for T-cell repertoire formation. Adv Immunol 99: 59-94., Krangel 2009Krangel MS. 2009. Mechanics of T cell receptor gene rearrangement. Curr Opin Immunol 21: 133-139., Seitan et al. 2011). After this TCR gene rearrangement, the T-cell expresses two surface TCR co-receptor CD molecules: CD4 and CD8 (double positive CD4+CD8+ immature T-cell).

Another process that takes place in thymus cortex is positive selection of immature T-cells (Teng et al. 2011, Takahama et al. 2010Takahama Y, Nitta T, Mat Ripen A, Nitta S, Murata S and Tanaka K. 2010. Role of thymic cortex-specific self-peptides in positive selection of T cells. Semin Immunol 22: 287-293., Klein et al. 2009Klein L, Hinterberger M, Wirnsberger G and Kyewski B. 2009. Antigen presentation in the thymus for positive selection and central tolerance induction. Nat Rev Immunol 9: 833-844., Gommeaux et al. 2009Gommeaux J, Grégoire C, Nguessan P, Richelme M, Malissen M, Guerder S, Malissen B and Carrier A. 2009. Thymus-specific serine protease regulates positive selection of a subset of CD4+ thymocytes. Eur J Immunol 39: 956-964.). Lymphocytes that have T-cell receptor capable of binding to a major histocompatibility complex (MHC) molecule, receive pro-survival signal. It is thought that during positive selection approximately 90% of immature T-cells (double positive CD4+CD8+) fail to receive this signal as their TCR is unable to interact with MHC molecule (Clarke et al. 2009Clarke RL, Thiemann S, Refaeli Y, Werlen G and Potter TA. 2009. A new function for LAT and CD8 during CD8- mediated apoptosis that is independent of TCR signal transduction. Eur J Immunol 39: 1619-1631.). As a result, these cells go through programmed cell death (apoptosis). The remaining, surviving T-cells become either CD4+CD8- or CD4-CD8+ (single positive) through the process of downregulation and upregulation. This lineage commitment is followed by migration of single positive T-cells to thymus medulla (Teng et al. 2011, Kurobe et al. 2006, Choi et al. 2008Choi YI, Duke-Cohan JS, Ahmed WB, Handley MA, Mann F, Epstein JA, Clayton LK and Reinherz EL. 2008. PlexinD1 glycoprotein controls migration of positively selected thymocytes into the medulla. Immunity 29: 888-898.).

In thymus medulla, the single positive T-cells that show affinity to self-antigens presented on MHC molecules receive pro-apoptotic signal (Le Borgne et al. 2009, Qiu et al. 2010Qiu Q, Ravens I, Seth S, Rathinasamy A, Maier MK, Davalos-Misslitz A, Forster R and Bernhardt G. 2010. CD155 is involved in negative selection and is required to retain terminally maturing CD8 T cells in thymus. J Immunol 184: 1681-1689., Atibalentja et al. 2009Atibalentja DF, Byersdorfer CA and Unanue ER. 2009. Thymus-blood protein interactions are highly effective in negative selection and regulatory T cell induction. J Immunol 183: 7909-7918., McCaughtry et al. 2008McCaughtry TM, Baldwin TA, Wilken MS and Hogquist KA. 2008. Clonal deletion of thymocytes can occur in the cortex with no involvement of the medulla. J Exp Med 205: 2575-2584.). This elimination of auto-reactive T-cells is called negative selection, and its main purpose is to prevent the development of autoimmunity. Although negative selection can also occur in double positive T-cells in cortex, it is much more pronounced and efficient in medulla.

Therefore, cortical and medullar lymphocytes, although all being T-cells, have some major differences in their functional and maturational status. First, these two populations have different expression of TCR co-receptor CD molecules: cortical lymphocytes are predominantly double negative (CD4-CD8-) or double positive (CD4+CD8+), while medullar, more mature lymphocytes are mostly single positive (CD4+CD8- or CD4-CD8+). Second, the microenvironment of thymus cortex is substantially different than in medulla; this refers to different interleukin levels, presence or absence of certain growth factors, as well as differences in epithelial reticular cell functions (Nitta et al. 2008Nitta T, Murata S, Ueno T, Tanaka K and Takahama Y. 2008. Thymic microenvironments for T-cell repertoire formation. Adv Immunol 99: 59-94., Benz et al. 2004Benz C, Heinzel K and Bleul CC. 2004. Homing of immature thymocytes to the subcapsular microenvironment within the thymus is not an absolute requirement for T cell development. Eur J Immunol 34: 3652-3563., Ma et al. 2013Ma D, Wei Y and Liu F. 2013. Regulatory mechanisms of thymus and T cell development. Dev Comp Immunol 39: 91-102.). All these factors can have an impact on nuclear structural organization of the two lymphocyte populations.

Our results indicate that GLCM as a method may be capable of detecting discrete nuclear structural changes that are associated with thymus lymphocyte migration and maturation. This is in accordance with the study of Losa and Castelli (2005)Losa GA and Castelli C. 2005. Nuclear patterns of human breast cancer cells during apoptosis: characterisation by fractal dimension and co-occurrence matrix statistics. Cell Tissue Res 322: 257-267. who suggested that some GLCM parameters (i.e. textural correlation) may be very efficient in identifying chromatin organizational changes in cell cultures, in some cases with greater sensitivity than conventional molecular biology methods such as TdT-mediated dUTP nick end labeling assay, FACScan Sub G0/G1 peak, or Annexin-V labelling. In our study, of analyzed GLCM parameters, changes in entropy, variance and angular second moment were more pronounced than in correlation. In our study, in addition to the standard statistical difference tests, the ability of the GLCM parameters to discriminate the two lymphocyte populations was tested by ROC analysis. In clinical sciences, ROC analysis is essential for evaluating potential value of the studied diagnostic test. ROC curve is constructed based on the fraction of true positives out of the positives (sensitivity) and the fraction of false positives out of the negatives (1-specificity), taking into account different thresholds. In our work, in SPSS software [state variable]=1 was assigned to the cortical lymphocytes, while [state variable]=0 was assigned to the medullar lymphocytes. Therefore, cortical location would be equivalent to [positive] diagnosis in a clinical diagnostic test.

Potential discriminatory performance of the test is evaluated by analyzing the area under the ROC curve. A test with no diagnostic value would have the area under the ROC curve approximately higher than 0.5 and lower than 0.6. Area between 0.6 and 0.7 indicates [poor] performance, area between 0.7 and 0.8 [fair] performance, and area between 0.8 and 0.9 [good] performance (Zweig and Campbell 1993Zweig MH and Campbell G. 1993. Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine. Clin Chem 39: 561-577., Sandelowsky et al. 2011Sandelowsky H, Ställberg B, Nager A and Hasselström J. 2011. The prevalence of undiagnosed chronic obstructive pulmonary disease in a primary care population with respiratory tract infections - a case finding study. BMC Fam Pract 12, 122 p.). Tests that belong to the category [excellent] usually have the ROC areas higher than 0.9.

In this study, all four tested GLCM parameters had ROC areas between 0.6 and 0.7, with entropy having the highest accuracy and correlation the lowest. Although these values show relatively low discriminatory performance of GLCM parameters, it should be noted that at present there is no known image analysis method that would represent a gold standard for nuclear structure classification in thymus. Also, we cannot exclude the possibility that in some future study, a scoring system combining various GLCM parameters could be designed and would have much higher performance than individual GLCM parameters taken separately. Alternatively, one could combine GLCM method with other known image complexity analysis techniques, such as fractal analysis, particularly having in mind that some computational biology programs (i.e. ImageJ) can be easily upgraded to apply both techniques on same digital micrographs (Pantic et al. 2012bPantic I, Harhaji-Trajkovic L, Pantovic A, Milosevic NT and Trajkovic V. 2012b. Changes in fractal dimension and lacunarity as early markers of UV-induced apoptosis. J Theor Biol 303: 87-92.).

The presented findings regarding entropy and angular second moment are also in accordance with the results of our previous study concerning GLCM analysis of B-lymphocytes (B-cells) in spleen germinal center (Pantic et al. 2012aPantic I, Pantic S and Basta-Jovanovic G. 2012a. Grey level co-occurrence matrix (GLCM) texture analysis of germinal center light zone lymphocyte nuclei: physiology viewpoint with focus on apoptosis. Microsc Microanal 18: 470-475.). We suggest that nuclear structural entropy and ASM of spleen germinal center light zone B-cell nuclei are correlated with nuclear structural reorganization during apoptosis. This reorganization may have included chromatin condensation, marginalization and fluctuations in euchromatin / heterochromatin ratio. The chromatin changes that occur during lymphocyte maturation in thymus (and migration from cortex to medulla) are probably more discrete than the ones taking place in germinal center, which makes GLCM as a method potentially even more sensitive than previously thought.

Our study had several potential limitations. First, all 300 analyzed lymphocyte nuclei originated from thymus tissues of one-month-old mice. It is well known that with age, thymus as a primary lymphoid organ begins to atrophy, followed by its decreased function (Calder et al. 2011Calder AE, Hince MN, Dudakov JA, Chidgey AP and Boyd RL. 2011. Thymic involution: where endocrinology meets immunology. Neuroimmunomodulation 18: 281-289., Gruver et al. 2007Gruver AL, Hudson LL and Sempowski GD. 2007. Immuno-senescence of ageing. J Pathol 211: 144-156.). It remains unclear whether GLCM method would detect the same differences between cortical and medullar lymphocyte chromatin in older tissues. Perhaps, in the future, it would be of interest to conduct a similar, but larger study which would investigate possible effects of aging on T-cell chromatin GLCM parameters. Second, if we consider the total number of T-cells that are present in thymus at a given time, the sample of 300 nuclei is somewhat limited. Nevertheless, we should have in mind that GLCM as a mathematical and computational biology method is relatively exact and when combined with NIH ImageJ software, it significantly lowers possible measurement variations that would otherwise be the result of subjective interpretation, researcher bias, or small sample size.

In conclusion, our findings suggest that grey level co-occurrence matrix entropy, angular second moment, variance and correlation might be useful in detection of discrete structural changes of lymphocyte nuclei in thymus. Although discriminatory performance of the calculated GLCM parameters is relatively poor, we assume that in the future, a more sensitive scoring system based on a combination of individual GLCM features could be designed. Along with results of other studies, it would represent a basis for further research regarding applicability of GLCM method in medical and biological sciences.

The authors are grateful to The Ministry of Education and Science, Republic of Serbia, Research Projects oi-175059 and iii-41027, as well as the staff at The Institute of Physiology Dragisa Nikitovic and Dragana Pavlovic for their help in obtaining the laboratory animals.

The authors report no conflict of interests regarding this manuscript.

REFERENCES

  • Alvarenga AV, Teixeira CA, Ruano MG and Pereira WC. 2010. Influence of temperature variations on the entropy and correlation of the Grey-Level Co-occurrence Matrix from B-Mode images. Ultrasonics 50: 290-293.
  • Atibalentja DF, Byersdorfer CA and Unanue ER. 2009. Thymus-blood protein interactions are highly effective in negative selection and regulatory T cell induction. J Immunol 183: 7909-7918.
  • Benz C, Heinzel K and Bleul CC. 2004. Homing of immature thymocytes to the subcapsular microenvironment within the thymus is not an absolute requirement for T cell development. Eur J Immunol 34: 3652-3563.
  • Calder AE, Hince MN, Dudakov JA, Chidgey AP and Boyd RL. 2011. Thymic involution: where endocrinology meets immunology. Neuroimmunomodulation 18: 281-289.
  • Choi YI, Duke-Cohan JS, Ahmed WB, Handley MA, Mann F, Epstein JA, Clayton LK and Reinherz EL. 2008. PlexinD1 glycoprotein controls migration of positively selected thymocytes into the medulla. Immunity 29: 888-898.
  • Clarke RL, Thiemann S, Refaeli Y, Werlen G and Potter TA. 2009. A new function for LAT and CD8 during CD8- mediated apoptosis that is independent of TCR signal transduction. Eur J Immunol 39: 1619-1631.
  • Gommeaux J, Grégoire C, Nguessan P, Richelme M, Malissen M, Guerder S, Malissen B and Carrier A. 2009. Thymus-specific serine protease regulates positive selection of a subset of CD4+ thymocytes. Eur J Immunol 39: 956-964.
  • Gruver AL, Hudson LL and Sempowski GD. 2007. Immuno-senescence of ageing. J Pathol 211: 144-156.
  • Haralick RM, Shanmugam K and Dinstein I. 1973. Textural Features for Image Classification. IEEE Transactions on Systems Man Cybernetics 3: 610-621.
  • Huber MB, Nagarajan MB, Leinsinger G, Eibel R, Ray LA and WismÜller A. 2011. Performance of topological texture features to classify fibrotic interstitial lung disease patterns. Med Phys 38: 2035-2044
  • Image J. 2011. National Institutes of Health (NIH). http://rsbweb.nih.gov/ij Accessed 26 December 2011.
    » http://rsbweb.nih.gov/ij
  • Klein L, Hinterberger M, Wirnsberger G and Kyewski B. 2009. Antigen presentation in the thymus for positive selection and central tolerance induction. Nat Rev Immunol 9: 833-844.
  • Krangel MS. 2009. Mechanics of T cell receptor gene rearrangement. Curr Opin Immunol 21: 133-139.
  • Kurobe H et al. 2006. CCR7-dependent cortex-to-medulla migration of positively selected thymocytes is essential for establishing central tolerance. Immunity 24: 165-177.
  • Le Borgne M, Ladi E, Dzhagalov I, Herzmark P, Liao YF, Chakraborty AK and Robey EA. 2009. The impact of negative selection on thymocyte migration in the medulla. Nat Immunol 10: 823-830.
  • Li X, Pai A, Blumenkrantz G, Carballido-Gamio J, Link T, Ma B, Ries M and Majumdar S. 2009. Spatial distribution and relationship of T1 rho and T2 relaxation times in knee cartilage with osteoarthritis. Magn Reson Med 61: 1310-1318.
  • Losa GA and Castelli C. 2005. Nuclear patterns of human breast cancer cells during apoptosis: characterisation by fractal dimension and co-occurrence matrix statistics. Cell Tissue Res 322: 257-267.
  • Ma D, Wei Y and Liu F. 2013. Regulatory mechanisms of thymus and T cell development. Dev Comp Immunol 39: 91-102.
  • McCaughtry TM, Baldwin TA, Wilken MS and Hogquist KA. 2008. Clonal deletion of thymocytes can occur in the cortex with no involvement of the medulla. J Exp Med 205: 2575-2584.
  • Nitta T, Murata S, Ueno T, Tanaka K and Takahama Y. 2008. Thymic microenvironments for T-cell repertoire formation. Adv Immunol 99: 59-94.
  • Pantic I, Harhaji-Trajkovic L, Pantovic A, Milosevic NT and Trajkovic V. 2012b. Changes in fractal dimension and lacunarity as early markers of UV-induced apoptosis. J Theor Biol 303: 87-92.
  • Pantic I and Pantic S. 2012. Germinal center texture entropy as possible indicator of humoral immune response: immunophysiology viewpoint. Mol Imaging Biol 14: 534-540.
  • Pantic I, Pantic S and Basta-Jovanovic G. 2012a. Grey level co-occurrence matrix (GLCM) texture analysis of germinal center light zone lymphocyte nuclei: physiology viewpoint with focus on apoptosis. Microsc Microanal 18: 470-475.
  • Qiu Q, Ravens I, Seth S, Rathinasamy A, Maier MK, Davalos-Misslitz A, Forster R and Bernhardt G. 2010. CD155 is involved in negative selection and is required to retain terminally maturing CD8 T cells in thymus. J Immunol 184: 1681-1689.
  • Sandelowsky H, Ställberg B, Nager A and Hasselström J. 2011. The prevalence of undiagnosed chronic obstructive pulmonary disease in a primary care population with respiratory tract infections - a case finding study. BMC Fam Pract 12, 122 p.
  • Seitan VC et al. 2011. A role for cohesin in T-cell-receptor rearrangement and thymocyte differentiation. Nature 476: 467-471.
  • Shamir L, Wolkow CA and Goldberg IG. 2009. Quantitative measurement of aging using image texture entropy. Bioinformatics 25: 3060-3063.
  • Sleckman BP. 2005. Lymphocyte antigen receptor gene assembly: multiple layers of regulation. Immunol Res 32: 253-258.
  • Takahama Y, Nitta T, Mat Ripen A, Nitta S, Murata S and Tanaka K. 2010. Role of thymic cortex-specific self-peptides in positive selection of T cells. Semin Immunol 22: 287-293.
  • Teng F et al. 2011. The molecular signature underlying the thymic migration and maturation of TCRαβ+ CD4+ CD8 thymocytes. PLoS One 6: e25567.
  • Zweig MH and Campbell G. 1993. Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine. Clin Chem 39: 561-577.

Publication Dates

  • Publication in this collection
    13 Aug 2013
  • Date of issue
    Sept 2013

History

  • Received
    3 Apr 2012
  • Accepted
    11 Oct 2012
Academia Brasileira de Ciências Rua Anfilófio de Carvalho, 29, 3º andar, 20030-060 Rio de Janeiro RJ Brasil, Tel: +55 21 3907-8100 - Rio de Janeiro - RJ - Brazil
E-mail: aabc@abc.org.br