Skip to main content

BRIEF RESEARCH REPORT article

Front. Physiol., 05 May 2021
Sec. Fractal Physiology

Loss of Mammographic Tissue Homeostasis in Invasive Lobular and Ductal Breast Carcinomas vs. Benign Lesions

  • 1Laboratory of Physical Foundations of Strength, Institute of Continuous Media Mechanics UB RAS, Perm, Russia
  • 2CompuMAINE Laboratory, University of Maine, Orono, ME, United States
  • 3Department of Pathology, Perm State Medical University Named After Academician E. A. Wagner, Perm, Russia
  • 4Laboratoire Ondes et Matière d’Aquitaine, Universite de Bordeaux, Bordeaux, France
  • 5Department of Chemical and Biomedical Engineering, University of Maine, Orono, ME, United States

The 2D wavelet transform modulus maxima (WTMM) method is used to perform a comparison of the spatial fluctuations of mammographic breast tissue from patients with invasive lobular carcinoma, those with invasive ductal carcinoma, and those with benign lesions. We follow a procedure developed and validated in a previous study, in which a sliding window protocol is used to analyze thousands of small subregions in a given mammogram. These subregions are categorized according to their Hurst exponent values (H): fatty tissue (H ≤ 0.45), dense tissue (H ≥ 0.55), and disrupted tissue potentially linked with tumor-associated loss of homeostasis (0.45 < H < 0.55). Following this categorization scheme, we compare the mammographic tissue composition of the breasts. First, we show that cancerous breasts are significantly different than breasts with a benign lesion (p-value ∼ 0.002). Second, the asymmetry between a patient’s cancerous breast and its contralateral counterpart, when compared to the asymmetry from patients with benign lesions, is also statistically significant (p-value ∼ 0.006). And finally, we show that lobular and ductal cancerous breasts show similar levels of disruption and similar levels of asymmetry. This study demonstrates reproducibility of the WTMM sliding-window approach to help detect and characterize tumor-associated breast tissue disruption from standard mammography. It also shows promise to help with the detection lobular lesions that typically go undetected via standard screening mammography at a much higher rate than ductal lesions. Here both types are assessed similarly.

Introduction

Breast cancer is the second-most occurring cancer type, and is ranked as the fifth in terms of mortality (Siegel et al., 2015; Bray et al., 2018). The “gold standard” for assessing the state of the breast is X-ray screening mammography (Lannin and Wang, 2017). The primary and basic radiographic signs of breast cancer are masses and microcalcifications. Microcalcifications are indicative of the presence of calcium oxalate and calcium phosphate within the breast tissue, and have a small-cell character (1 mm or less in size), resembling grains of sand (Henrot et al., 2014; Wilkinson et al., 2017). Often, microcalcifications are the only radiologic manifestation of early breast cancer (Nalawade, 2009). However, this sign is not pathognomonic because in some histological forms of breast cancer, for example, lobular carcinoma, microcalcifications rarely occur (Fischmann, 2008). On the other hand, microcalcifications can occur in such benign processes as sclerosing adenosis or some fibroadenomas (Fischmann, 2008; Henrot et al., 2014). These findings are driving us to “think outside of the tumor” (Gerasimova-Chechkina et al., 2016) and to develop a computational approach to study and quantitatively characterize tissue microenvironment throughout the whole breast (Marin et al., 2017). Indeed, the breast tumor microenvironment plays a key role in early tumorigenesis. When the microenvironment is structurally sound, the tumor, if any, is under control (sleeping tumor; Bissell and Hines, 2011). If tissue structure is altered through cell cycle disruptions, the microenvironment may in fact promote tumor growth by selectively favoring the survival of cancer stem cells and protecting them from therapy or treatment (Maguer-Satta, 2011). A 2017 study presented “preliminary evidence that tissue disruption and loss of homeostasis in breast tissue microenvironment and breast bilateral asymmetry could be quantitatively and objectively assessed from mammography via a localized, wavelet-based multifractal analysis of the whole breast” (Marin et al., 2017). Tissue homeostatic balance emerges from the integration of multiple subcellular, intercellular, extracellular, chemical, and physical signals, and constraints (Rejniak, 2012). Tissue disruption is a term used in this paper and in Marin et al. (2017) to characterize what we infer as a larger-scale tissue architecture alteration caused by loss of tissue homeostasis. Contralateral asymmetry refers to imbalanced proportions of tissue organization as measured with the metrics used in this article.

The wavelet transform modulus maxima (WTMM) method is a multifractal formalism used to analyze complex 1D signals (Ivanov et al., 1999; Gerasimova et al., 2013, 2014), 2D images (Arneodo et al., 2000, 2003; Decoster et al., 2000; Roux et al., 2000; Kestener et al., 2001; Khalil et al., 2006; Batchelder et al., 2014; Marin et al., 2017), 3D images (Kestener and Arneodo, 2003), and vector fields (Kestener and Arneodo, 2004). The strategy previously developed to characterize loss of tissue homeostasis and tissue disruption (Marin et al., 2017) consists in deploying the 2D WTMM method to calculate the so-called (monofractal) Hurst exponent, H, by applying a sliding window approach to each mammogram. Then each subregion is color-coded based on its corresponding H value. Subregions where H ≤ 0.45 (fatty tissue) are colored blue, 0.45 < H < 0.55 (disrupted tissue) are yellow, and H ≥ 0.55 (dense tissue) are red. That pilot study demonstrated that “disrupted regions associated with loss of tissue homeostasis, as quantified by H∼1/2, and loss of breast symmetry, were found significantly more in tumorous cases when compared to normal cases” (Marin et al., 2017). Physical signatures associated with density fluctuations that are uncorrelated (H∼1/2) include randomness and free diffusion, which underpins the hypothesis that mammographic tumor-associated microenvironment displaying these uncorrelated fluctuations would be linked to a disruptive nature of the tissue and a perturbed homeostasis (Marin et al., 2017). Note that this methodology contrasts with most existing computer-aided detection/diagnostic methods (Karahaliou et al., 2008; Ayer et al., 2010; Tsai et al., 2011; Haberle et al., 2012; Meselhy Eltoukhy et al., 2012) that are predisposed for texture analysis or feature extraction.

A primary goal of this study is to demonstrate the reproducibility of the approach on a different set of mammograms than those analyzed in Marin et al. (2017). A further goal is to explore whether the tissue disruption associated with loss of homeostasis are similarly or differently characterized on patients with invasive lobular carcinomas (ILC) vs. those with invasive ductal carcinomas (IDC) [also known as invasive mammary carcinoma of no special type (Sinn and Kreipe, 2013)]. We perform a computational analysis on the mediolateral oblique (MLO) mammographic views from 81 patients with a malignant tumor (43 ILC and 38 IDC) and 23 patients with a benign tumor (12 fibroadenoma and 11 fibrocystic mastopathy).

Materials and Methods

The Data

Screening digital mammograms were obtained from Perm Regional Oncological Dispensary in Russia, following approval by the institution’s ethics committee. The mammographic procedure involved taking two X-ray images of the two compressed breasts for each patient, the standard MLO and craniocaudal (CC) views using an Alpha ST Mammograph (GE Healthcare). Only the MLO view was analyzed in this study. The MLO view is well-known to radiologists to include more breast tissue than any other single view (Bassett and Conner, 2003). The spatial resolution of the images is 50 μm per pixel. The data consist of mammograms with a pathology-proven diagnostic from 104 women, aged 47 to 72 (average 63.6 years old). See Table 1.

TABLE 1
www.frontiersin.org

Table 1. Study design and population.

WTMM Multifractal Analysis of Mammograms

The reader is referred to 2D WTMM references (Arneodo et al., 2000, 2003; Khalil et al., 2006, 2009; McAteer et al., 2010; Batchelder et al., 2014; Marin et al., 2017), specifically Marin et al. (2017), where for the first time an automated selection of critical 2D WTMM algorithmic parameters was presented. This same automated selection method was used here [see section 2.D. in Marin et al. (2017)]. A very brief overview of the 2D WTMM method follows. Using the Gaussian function θ(r) = exp(−r2/2), where r2 = x2 + y2, and x,y, the wavelet transform (WT) of a function (i.e., an image) f(x,y) ∈ L2() with respect to the two wavelets ψ1(x,y) = θ(x,y)/∂⁡x and ψ2(x,y) = θ(x,y)/∂⁡y is the vector

T ψ [ f ] ( b , a ) = { T θ [ f ] ( b , a ) } ,

where Tθ[f](b,a) = a−2∫∫f(r)θ(a−1(rb))d2r. At any given scale a > 0, the WTMM are the positions b in the image where the WT modulus Mψ[f](b,a) is locally maximum. For images such as those analyzed here (i.e., everywhere continuous but nowhere differentiable rough surfaces) these WTMM form connected chains. On these chains, further maxima are found, which are themselves connected through scales to form a collection of maxima lines that are part of the so-called WT space-scale skeleton, L(a). The partition functions

Z ( q , a ) = l L ( a ) ( sup ( x , a ) L , a a M ψ [ f ] ( b , a ) ) q ,

where q are statistical order moments, allow us to obtain the scaling exponents τ(q) via the relation Z(q,a)∼aτ(q), from which we can then obtain the singularity spectrum D(h)=minq(qh-τ(q). Here D is the fractal dimension of the set of points in the domain of f where the local Holder roughness exponent is h, i.e., the set of all points r where f(r+ϵ)−f(r)∼ϵh(r),|ϵ|→0. When the Holder exponent is identical throughout the domain of f, then the global Hurst roughness exponent, H is used to characterize the monofractal scaling properties. As discussed in Marin et al. (2017), the mammogram subregions almost always display monofractal scaling properties (those that do not are discarded from the analysis), and therefore, a single Hurst exponent value is associated to each subregion.

Image Sliding Window Analysis

As was done in Marin et al. (2017), each mammogram was divided into several thousands of overlapping subregions of size 360 × 360-pixel, of which only the central 256 × 256 wavelet-transformed portion was kept for analysis, by sliding a window by 32-pixel increments, going from the top left portion of the mammogram, down to the bottom right. Then the 2D WTMM method, for which the scaling parameters were set automatically via the objective method presented in Marin et al. (2017), was used compute the Hurst exponents for each subregion, except when the automated method failed to identify scaling parameters satisfying user-inputted minimal threshold requirements, as discussed in Marin et al. (2017). These latter cases were discarded from the analysis and shown in gray in Figure 1. Otherwise, the algorithm associated a color to the value of H: subregions where H ≤ 0.45 correspond to anti-correlated density fluctuations, which were found to be associated with fatty tissue, were colored blue; subregions where H ≥ 0.55 correspond to long-range correlated density fluctuation, which were found to be associated with dense tissue, were red; and subregions where 0.45 < H < 0.55 correspond to uncorrelated density fluctuation, which we refer to as disrupted tissue, were yellow (Marin et al., 2017). Sample cases are presented in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Color-coded MLO views of the tumorous and opposite breast from three sample patients. Each pixel represents “a 360 × 360-pixel mammogram subregion colored according to its H value. Subregions where H ≤ 0.45 (fatty) are colored blue, 0.45 < H < 0.55 (disrupted) are yellow, and H ≥ 0.55 (dense) are red. Gray pixels correspond to rejected subregions” [see (Marin et al., 2017)]. TOP: Patient with an invasive lobular carcinoma. MIDDLE: Patient with an invasive ductal carcinoma. BOTTOM: Patient with a benign lesion (fibrocystic mastopathy). For the two cancer patients, the two MLO panels show evident differences in terms of yellow (disrupted tissue, 0.45 < H < 0.55) subregions. However, for the benign patient, both breasts do not display evident visual differences in tissue disruption.

Metrics

The following metrics were calculated for each mammogram: the percentage of blue subregions per mammogram (%B), the percentage of yellow subregions per mammogram (%Y), the percentage of red subregions per mammogram (%R), and a combination score: %Y/%B. The latter metric was implemented based on empirical evidence that tumorous breasts had not only higher yellow subregions concentrations, but also lower blue subregions concentrations (Figure 2). And finally, to measure similarity between the tumorous breast and its contralateral counterpart, another metric, coined “YB Factor,” calculates the ratio of the %Y/%B scores from the tumorous breast to the contralateral opposite:

FIGURE 2
www.frontiersin.org

Figure 2. Box plots representing the distributions of %B (in blue), %Y (in yellow), %R (in red), and %Y/%B (in green). For the top row, all cases are included for both cancer and benign patients. In the middle row, the IDC population is compared to the ILC population. In the bottom row, the two benign subgroups, fib_a and fib_m, are compared. The p-values shown at the top of each plot were obtained by running a non-parametric Wilcoxon Rank Sum test.

YB Factor = % Y tumor / % B tumor % Y opposite / % B opposite .

The introduction of the YB Factor brings two advantages. It yields a single score per patient, which allows us to perform population statistics on patients (in addition to statistics on populations of breasts). And since the YB Factor is a ratio, its value is meaningful: YB Factor ∼ 1 is representative of similarity between the two breasts and any value away from 1 represents and quantifies asymmetry. Radiologists typically consider contralateral asymmetry as a feature of tumor development because breast asymmetry was shown to be greater in healthy women who later developed breast cancer than in women who did not (Scutt et al., 2006).

Statistical Tests

Using the Shapiro–Wilk’s test of normality, we found that some distributions were normal and others were not (Royston, 1982). All statistical distribution analyses and Wilcoxon Rank Sum tests (yielding the unadjusted p-values presented in this article) were performed using the R language, version 3.6.2 (R Core-Team, 2019).

Results

The four metrics were calculated for each mammogram (Figure 2), and the YB Factors were calculated for each patient (Figure 3). Statistical significance tests were performed using the non-parametric Wilcoxon Rank Sum test and box plots were created to display the differences and similarities between different types of breast and patient populations. The box plots and p-values are presented in Figures 2, 3, where for the latter one star represents a p-value < 0.05 and two stars represents p-value < 0.01. In these two figures, blue represents %B, yellow represents %Y, red represents %R, green represents %Y/%B, and purple represents the YB Factor. In Figure 2 the top row represents all cancer breasts vs. all benign breasts, the middle row represents the invasive ductal carcinoma (IDC1) breasts vs. the ILC breasts, and the bottom row represents the two benign breast categories: fibroadenoma (fib_a) vs. fibrocystic mastopathy (fib_m). In these figures, “Cancer” refers to the combined data from both IDC and ILC and “Benign” refers to the combined data from fib_a and fib_m. In Figure 3, the IDC data are broken down into their pathological grades (IDC1, IDC2, IDC3 referring to grades 1, 2, 3, respectively, – see Table 1).

FIGURE 3
www.frontiersin.org

Figure 3. Box plots representing the distribution of YB Factor for all subgroups of patient data. A horizontal dashed line is included at YB Factor = 1, an indicator of similarity between the tumorous and the opposite breast. One star represents p-value < 0.05 and two stars represent p-value < 0.01.

Overall, there is a statistically significant difference between all tumorous breasts with cancer (IDC and ILC combined) when compared to all breasts with a benign lesion (fib_a and fib_m combined). Cancer breasts have higher levels of disrupted tissue than benign breasts using both the %Y metric (p-value ∼ 0.003 – Figure 2, top yellow) and the disrupted tissue relative to fatty tissue metric, i.e., %Y/%B, (p-value ∼ 0.002 – Figure 2, top green).

On a patient level, when comparing breast asymmetry via the YB Factor, again, there is a statistically significant difference between the YB Factors from cancer patients vs. the YB Factors from benign patients (p-value ∼ 0.0062). In addition to measuring significant differences, the YB Factors yield a direct assessment of asymmetry: in terms of proportions, approximately 75% of the cancer patients had a YB Factor > 1 vs. only ∼50% for the benign cases (Figure 3).

Individually, the IDC patients and the ILC patients are significantly different than the benign patients when considering the YB Factor (Figure 3). The p-value for IDC vs. benign is ∼0.003 and the p-value for ILC vs. benign is ∼0.035. However, and quite notably, there is no significant difference between IDC breasts vs. ILC breasts.

And finally, we further note that the breasts fibrocystic mastopathy (fib_m) have higher levels of disrupted tissue than breasts with fibroadenoma (fib_a). This is true based on the %Y metric (Figure 2, bottom yellow, p-value ∼ 0.0426) as well as the %Y/%B metric (Figure 2, bottom green, p-value ∼ 0.0503). However, it is critically important to note that there is no significant asymmetry with the benign cases, as assessed by the YB Factor (Figure 3).

Discussion

Using a previously published computational methodology to assess tumor-associated loss of homeostasis via mammography, referred to as “tissue disruption,” we performed an analysis of tumorous breasts (benign vs. cancer) as well as an analysis of patients (tumorous breast vs. the contralateral opposite). Mammograms of cancerous breasts were found to have significantly higher levels of tissue disruption compared to benign breasts. Additionally, on a patient population level, the breast disruption asymmetry, as measured here via the newly introduced YB Factor, showed that while benign patients have relatively similar disruption levels in both breasts (YB Factor ∼ 1), that ratio is significantly higher in the tumorous breast for cancer patients. This means that when comparing the tumorous breasts to the healthy opposites, only the malignant breasts are significantly different from their opposites, not the benign breasts. This not only adds obvious diagnostic value to the tumorous vs. opposite breast approach, but, through further study, may lead to a better biophysical and physiological understanding of difference in loss of tissue homeostasis in malignant vs. benign microenvironments.

Patients with IDC and those with the much rarer (but more difficult to detect) ILC were individually significantly different than the benign population. Our investigation into potential differences between IDC and ILC showed that no significant differences existed between the two populations. This leads us to hypothesize that our methodology may be particularly well-suited to characterize tissue disruption associated with ILC and not just IDC. ILC is known for being difficult to detect radiologically (Krecke and Gisvold, 1993; Brem et al., 2009; Lopez and Bassett, 2009). It is typically detected at a later stage than IDC (Chen et al., 2017) and the mammographic appearance of ILC is often dangerously subtle (Johnson et al., 2015). Although ILC only accounts for 5–10% of breast cancers (Pestalozzi et al., 2008), it has an inherently invasive growth pattern and it disproportionately has greater metastasis compared to IDC (Arpino et al., 2004; Ferlicot et al., 2004; Chen et al., 2017; Nayeem et al., 2018). Moreover, ILC is more likely to lead to mastectomy and a lower long-term (5+ years) survival rate than IDC (Pestalozzi et al., 2008; Chen et al., 2017). Therefore, further validation of our methodology on a larger cohort could eventually translate to a major advancement in earlier detection of lobular breast cancer.

Important distinctions separate the analysis performed in previous work (Marin et al., 2017) vs. what is reported here. In the former, the analysis was performed on digitized mammogram films (scans) obtained from the Digital Database for Screening Mammography (DDSM; Heath et al., 1998, 2001), while here digitally acquired mammograms were analyzed. Although the pixel sizes are similar in both studies (∼50 microns/pixel), the dynamic range is different: lossless LJPEG 12-bit data (pixel values from 0 to 4095) for the DDSM data and uncompressed 8-bit BMP data (0–255) for the current study. Nonetheless, disrupted tissue regions (yellow squares, where the Hurst exponent, H, is between 0.45 and 0.55) are significantly found in larger numbers in tumorous breasts in both studies. This is important for two reasons: it demonstrates reproducibility of the approach (different patient populations were used in each study), and the same method offers a similar diagnostic potential in scans and digitally acquired mammograms.

A more refined study will be undertaken to investigate more patients, with additional metrics [also see Marin et al. (2017)], and to explore whether or not (and to what extent) different histopathological types of tumors could be discriminated with our methodology. In particular, future investigations and explorations of combinations of metrics are recommended for larger cohorts, ideally with longitudinal datasets.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Local Ethics Committee of Perm State Medical University named after E.A. Wagner (protocol No 10 from 28/11/2018). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

EG-C, GF, and IA: data acquisition. BW: image pre-processing. BT and AK: WTMM analysis. EG-C, KB, and AK: statistical analyses and figure preparation. EG-C, KB, and AK: manuscript writing. All authors have read and approved of the manuscript.

Funding

The authors thank the Ministry of Education and Science of Russia (grant No 14.W01.17.2674-MK). We acknowledge the support of NVIDIA Corporation with the donation of the Titan X Pascal GPU. Research reported in this manuscript was partially supported by National Cancer Institute of the National Institutes of Health under award number R15CA246335. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Conflict of Interest

We declare the United States Patent 10,467,755 Methods of Cancer Detection 11/05/2019. Inventors: AK and KB.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank roengenologist Marina Alexandrenok from the Department of Roentgenology, Perm (Russia) Regional Oncological Dispensary, for her help in selecting mammograms and her expert conclusions. We also thank Amy Harrow and Christine Houser for helpful discussions.

Footnotes

  1. ^ Also referred to as invasive mammary carcinoma of no special type (Sinn and Kreipe, 2013).

References

Siegel, R. L., Miller, K. D., and Jemal, A. (2015). Cancer statistics, 2015. CA Cancer J Clin 65, 5–29. doi: 10.3322/caac.21254

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: Cancer J. Clin 68, 394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

Lannin, D. R., and Wang, S. (2017). Are small breast cancers good because they are small or small because they are good? N Engl J Med 376, 2286–2291. doi: 10.1056/nejmsr1613680

PubMed Abstract | CrossRef Full Text | Google Scholar

Henrot, P., Leroux, A., Barlier, C., and Génin, P. (2014). Breast microcalcifications: The lesions in anatomical pathology. Diagn Interv Radiol 95, 141–152. doi: 10.1016/j.diii.2013.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkinson, L., Thomas, V., and Sharma, N. (2017). Microcalcification on mammography: approaches to interpretation and biopsy. Br J Radiol 90, 20160594. doi: 10.1259/bjr.20160594

PubMed Abstract | CrossRef Full Text | Google Scholar

Nalawade, Y. V. (2009). Evaluation of breast calcifications. Indian J Radiol Imaging 19, 282–286.

Google Scholar

Fischmann, A. (2008). “Microcalcification in breast lesions: radiography and histopathology,” in Cancer Imaging, ed. M. A. Hayat (San Diego., CA: Academic Press), 383–392. doi: 10.1016/b978-012374212-4.50042-0

CrossRef Full Text | Google Scholar

Gerasimova-Chechkina, E., Toner, B., Marin, Z., Audit, B., Roux, S., Argoul, F., et al. (2016). Combining multifractal analyses of digital mammograms and infrared thermograms to assist in early breast cancer diagnosis. AIP Conference Proceedings 1760, 020018.

Google Scholar

Marin, Z., Batchelder, K. A., Toner, B. C., Guimond, L., Gerasimova-Chechkina, E., Harrow, A. R., et al. (2017). Mammographic evidence of microenvironment changes in tumorous breasts. Med Phys 44, 1324–1336. doi: 10.1002/mp.12120

PubMed Abstract | CrossRef Full Text | Google Scholar

Bissell, M. J., and Hines, W. C. (2011). Why don’t we get more cancer? A proposed role of the microenvironment in restraining cancer progression. Nat Med 17, 320–329. doi: 10.1038/nm.2328

PubMed Abstract | CrossRef Full Text | Google Scholar

Maguer-Satta, V. (2011). “The Stem Cell Niche: The Black Master of Cancer,” in Cancer Stem Cells Theories and Practice, ed. S. Shostak (Rijeka: InTech), 215–240.

Google Scholar

Rejniak, K. A. (2012). Homeostatic Imbalance in Epithelial Ducts and Its Role in Carcinogenesis. Scientifica 2012, 132978.

Google Scholar

Gerasimova, E., Audit, B., Roux, S. G., Khalil, A., Argoul, F., Naimark, O., et al. (2013). Multifractal analysis of dynamic infrared imaging of breast cancer. EPL 104, 68001. doi: 10.1209/0295-5075/104/68001

CrossRef Full Text | Google Scholar

Gerasimova, E., Audit, B., Roux, S. G., Khalil, A., Gileva, O., Argoul, F., et al. (2014). Wavelet-based multifractal analysis of dynamic infrared thermograms to assist in early breast cancer diagnosis. Front. Physiol. 5:176. doi: 10.3389/fphys.2014.00176

PubMed Abstract | CrossRef Full Text | Google Scholar

Ivanov, P. C., Amaral, L. A., Goldberger, A. L., Havlin, S., Rosenblum, M. G., Struzik, Z. R., et al. (1999). Multifractality in human heartbeat dynamics. Nature 399, 461–465. doi: 10.1038/20924

PubMed Abstract | CrossRef Full Text | Google Scholar

Arneodo, A., Decoster, N., and Roux, S. (2000). A wavelet-based method for multifractal image analysis. I. Methodology and test applications on isotropic and anisotropic random rough surfaces. Eur Phys J B 15, 567–600. doi: 10.1007/s100510051161

CrossRef Full Text | Google Scholar

Decoster, N., Arneodo, A., and Roux, S. (2000). A wavelet-based method for multifractal image analysis. II. Applications to synthetic multifractal rough surfaces. Eur Phys J B 15, 739–764. doi: 10.1007/s100510051179

CrossRef Full Text | Google Scholar

Roux, S., Decoster, N., and Arneodo, A. (2000). A wavelet-based method for multifractal image analysis. III. Applications to high-resolution satellite images of cloud structure. Eur Phys J B 15, 765–786. doi: 10.1007/s100510051180

CrossRef Full Text | Google Scholar

Kestener, P., Lina, J. M., St-Jean, P., and Arneodo, A. (2001). Wavelet-based multifractal formalism to assist in diagnosis in digitized mammograms. Image Anal Stereol 20, 169–174. doi: 10.5566/ias.v20.p169-174

CrossRef Full Text | Google Scholar

Arneodo, A., Decoster, N., Kestener, P., and Roux, S. G. (2003). A wavelet-based method for multifractal image analysis: From theoretical concepts to experimental applications. Adv Imaging Electr Phys 126, 1–92. doi: 10.1016/s1076-5670(03)80014-9

CrossRef Full Text | Google Scholar

Khalil, A., Joncas, G., Nekka, F., Kestener, P., and Arneodo, A. (2006). Morphological analysis of HI features. II. Wavelet-based multifractal formalism. Astrophys J Suppl S 165, 512–550. doi: 10.1086/505144

CrossRef Full Text | Google Scholar

Batchelder, K. A., Tanenbaum, A. B., Albert, S., Guimond, L., Kestener, P., Arneodo, A., et al. (2014). Wavelet-based 3D reconstruction of microcalcification clusters from two mammographic views: New evidence that fractal tumors are malignant and Euclidean tumors are benign. PLoS One 9:e107580. doi: 10.1371/journal.pone.0107580

PubMed Abstract | CrossRef Full Text | Google Scholar

Kestener, P., and Arneodo, A. (2003). Three-dimensional wavelet-based multifractal method: The need for revisiting the multifractal description of turbulence dissipation data. Phys Rev Lett 91, 194501.

Google Scholar

Kestener, P., and Arneodo, A. (2004). Generalizing the wavelet-based multifractal formalism to vector-valued random fields: application to turbulent velocity and vorticity 3D numerical data. Phys Rev Lett 93, 044501.

Google Scholar

Karahaliou, A. N., Boniatis, I. S., Skiadopoulos, S. G., Sakellaropoulos, F. N., Arikidis, N. S., Likaki, E. A., et al. (2008). Breast cancer diagnosis: analyzing texture of tissue surrounding microcalcifications. IEEE Trans Inf Technol Biomed 12, 731–738. doi: 10.1109/titb.2008.920634

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayer, T., Ayvaci, M. U., Liu, Z. X., Alagoz, O., and Burnside, E. S. (2010). Computer-aided diagnostic models in breast cancer screening. Imaging Med 2, 313–323. doi: 10.2217/iim.10.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, N. C., Chen, H. W., and Hsu, S. L. (2011). Computer-aided diagnosis for early-stage breast cancer by using Wavelet Transform. Comput Med Imaging Graph 35, 1–8. doi: 10.1016/j.compmedimag.2010.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Haberle, L., Wagner, F., Fasching, P. A., Jud, S. M., Heusinger, K., Loehberg, C. R., et al. (2012). Characterizing mammographic images by using generic texture features. Breast Cancer Res 14, R59.

Google Scholar

Meselhy Eltoukhy, M., Faye, I., and Belhaouari Samir, B. (2012). A statistical based feature extraction method for breast cancer diagnosis in digital mammogram using multiresolution representation. Comput Biol Med 42, 123–128. doi: 10.1016/j.compbiomed.2011.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Bassett, L., and Conner, K. (2003). “Types of Mammography,” in Holland-Frei Cancer Medicine, 6th Edn, eds D. W. Kufe, R. E. Pollock, and R. R. Weichselbaum (Hamilton ON: BC Decker).

Google Scholar

McAteer, R. T. J., Kestener, P., Arneodo, A., and Khalil, A. (2010). Automated detection of coronal loops using a wavelet transform modulus maxima method. Sol Phys 262, 387–397. doi: 10.1007/s11207-010-9530-7

CrossRef Full Text | Google Scholar

Khalil, A., Aponte, C., Zhang, R., Davisson, T., Dickey, I., Engelman, D., et al. (2009). Image analysis of soft-tissue in-growth and attachment into highly porous alumina ceramic foam metals. Med Eng Phys 31, 775–783. doi: 10.1016/j.medengphy.2009.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Scutt, D., Lancaster, G. A., and Manning, J. T. (2006). Breast asymmetry and predisposition to breast cancer. Breast Cancer Research 8, R14.

Google Scholar

Royston, P. (1982). An extension of Shapiro and Wilk’s W test for normality to large samples. Applied Statistics 31, 115–124. doi: 10.2307/2347973

CrossRef Full Text | Google Scholar

R Core-Team (2019). R: A language and environment for statistical computing. Available online at: https://www.R-project.org/ (accessed October 22, 2020)

Google Scholar

Krecke, K. N., and Gisvold, J. J. (1993). Invasive lobular carcinoma of the breast: mammographic findings and extent of disease at diagnosis in 184 patients. AJR Am J Roentgenol 161, 957–960. doi: 10.2214/ajr.161.5.8273634

PubMed Abstract | CrossRef Full Text | Google Scholar

Brem, R. F., Ioffe, M., Rapelyea, J. A., Yost, K. G., Weigert, J. M., Bertrand, M. L., et al. (2009). Invasive lobular carcinoma: detection with mammography, sonography, MRI, and breast-specific gamma imaging. AJR Am J Roentgenol 192, 379–383. doi: 10.2214/ajr.07.3827

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, J. K., and Bassett, L. W. (2009). Invasive Lobular Carcinoma of the Breast: Spectrum of Mammographic, US, and MR Imaging Findings. RadioGraphics 29, 165–176. doi: 10.1148/rg.291085100

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Yang, J., Li, S., Lv, M., Shen, Y., Wang, B., et al. (2017). Invasive lobular carcinoma of the breast: A special histological type compared with invasive ductal carcinoma. PloS one 12:e0182397. doi: 10.1371/journal.pone.0182397

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, K., Sarma, D., and Hwang, E. S. (2015). Lobular breast cancer series: imaging. Breast Cancer Research 17, 94.

Google Scholar

Pestalozzi, B. C., Zahrieh, D., Mallon, E., Gusterson, B. A., Price, K. N., Gelber, R. D., et al. (2008). Distinct clinical and prognostic features of infiltrating lobular carcinoma of the breast: combined results of 15 International Breast Cancer Study Group clinical trials. J Clin Oncol 26, 3006–3014. doi: 10.1200/jco.2007.14.9336

PubMed Abstract | CrossRef Full Text | Google Scholar

Nayeem, M., Dang, M.-T., Saleh, H., and Abu-Rashed, A. (2018). The Breast in the End: Metastatic Lobular Breast Cancer Found in the Rectum: 1544. American Journal of Gastroenterology 113, S888.

Google Scholar

Ferlicot, S., Vincent-Salomon, A., Médioni, J., Genin, P., Rosty, C., Sigal-Zafrani, B., et al. (2004). Wide metastatic spreading in infiltrating lobular carcinoma of the breast. European Journal of Cancer 40, 336–341. doi: 10.1016/j.ejca.2003.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Arpino, G., Bardou, V. J., Clark, G. M., and Elledge, R. M. (2004). Infiltrating lobular carcinoma of the breast: tumor characteristics and clinical outcome. Breast cancer research : BCR 6, R149–R156.

Google Scholar

Heath, M., Bowyer, K., Kopans, D., Kegelmeyer, W. P., Moore, R., Chang, K., et al. (1998). “Current status of the Digital Database for Screening Mammography,” in Proceedings of the 4th International Workshop on Digital Mammography, (Nijmegen).

Google Scholar

Heath, M., Bowyer, K., Kopans, D., Moore, R., and Kegelmeyer, W. P. (2001). “The Digital Database for Screening Mammography,” in Proceedings of the 5th International Workshop on Digital Mammography, (Madison, WI).

Google Scholar

Sinn, H.-P., and Kreipe, H. (2013). A Brief Overview of the WHO Classification of Breast Tumors, 4th Edition, Focusing on Issues and Updates from the 3rd Edition. Breast care (Basel) 8, 149–154. doi: 10.1159/000350774

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Radiomics, tissue homeostasis, mammography, wavelets, breast density, fractals, multifractals, Hurst exponent (H)

Citation: Gerasimova-Chechkina E, Toner BC, Batchelder KA, White B, Freynd G, Antipev I, Arneodo A and Khalil A (2021) Loss of Mammographic Tissue Homeostasis in Invasive Lobular and Ductal Breast Carcinomas vs. Benign Lesions. Front. Physiol. 12:660883. doi: 10.3389/fphys.2021.660883

Received: 29 January 2021; Accepted: 09 April 2021;
Published: 05 May 2021.

Edited by:

Paul Bogdan, University of Southern California, Los Angeles, United States

Reviewed by:

Yih-Kuen Jan, University of Illinois at Urbana-Champaign, United States
Corina Stefania Drapaca, Pennsylvania State University, United States

Copyright © 2021 Gerasimova-Chechkina, Toner, Batchelder, White, Freynd, Antipev, Arneodo and Khalil. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Andre Khalil, andre.khalil@maine.edu

These authors have contributed equally to this work

Deceased

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.