Reproductive risk factors for breast cancer and association with novel breast density measurements among Hispanic, Black, and White women

PURPOSE

There are differences in the distributions of breast cancer incidence and risk factors by race and ethnicity. Given the strong association between breast density and breast cancer, it is of interest describe racial and ethnic variation in the determinants of breast density.

METHODS

We characterized racial and ethnic variation in reproductive history and several measures of breast density for Hispanic (n = 286), non-Hispanic Black (n = 255), and non-Hispanic White (n = 1694) women imaged at a single hospital. We quantified associations between reproductive factors and percent volumetric density (PVD), dense volume (DV), non-dense volume (NDV), and a novel measure of pixel intensity variation (V) using multivariable-adjusted linear regression, and tested for statistical heterogeneity by race and ethnicity.

RESULTS

Reproductive factors most strongly associated with breast density were age at menarche, parity, and oral contraceptive use. Variation by race and ethnicity was most evident for the associations between reproductive factors and NDV (minimum p-heterogeneity:0.008) and V (minimum p-heterogeneity:0.004) and least evident for PVD (minimum p-heterogeneity:0.042) and DV (minimum p-heterogeneity:0.041).

CONCLUSION

Reproductive choices, particularly those related to childbearing and oral contraceptive use, may contribute to racial and ethnic variation in breast density.

Association of Daily Alcohol Intake, Volumetric Breast Density, and Breast Cancer Risk

High alcohol intake and breast density increase breast cancer (BC) risk, but their interrelationship is unknown. We examined whether volumetric density modifies and/or mediates the alcohol-BC association. BC cases (n = 2233) diagnosed from 2006 to 2013 in the San Francisco Bay area had screening mammograms 6 or more months before diagnosis; controls (n = 4562) were matched on age, mammogram date, race or ethnicity, facility, and mammography machine. Logistic regression was used to estimate alcohol-BC associations adjusted for age, body mass index, and menopause; interaction terms assessed modification. Percent mediation was quantified as the ratio of log (odds ratios [ORs]) from models with and without density measures. Alcohol consumption was associated with increased BC risk (2-sided P_trend = .004), as were volumetric percent density (OR = 1.45 per SD, 95% confidence interval [CI] = 1.36 to 1.56) and dense volume (OR = 1.30, 95% CI = 1.24 to 1.37). Breast density did not modify the alcohol-BC association (2-sided P > .10 for all). Dense volume mediated 25.0% (95% CI = 5.5% to 44.4%) of the alcohol-BC association (2-sided P = .01), suggesting alcohol may partially increase BC risk by increasing fibroglandular tissue.

Abstract P5-08-02: Association of daily alcohol intake, volumetric density and breast cancer risk

Background: Alcohol intake and breast density are established risk factors for breast cancer (BC). A few studies suggest breast density is on the causal pathway between alcohol and BC but the literature is not consistent. We examine the interrelationship between alcohol, volumetric breast density and BC risk, specifically whether there is a stronger alcohol and BC association among women with dense breasts and/or whether breast density mediates the alcohol and BC association.

Methods: We conducted a nested case-control study within the San Francisco Mammography Registry consisting of four breast screening centers. From 2006-2013, through linkages with the California cancer registry, we identified 2572 women diagnosed with BC who had screening mammograms performed at least six months prior to diagnosis. Controls (5119) were matched to cases on age, date of earliest mammogram, race/ethnicity, facility, and mammography machine. We obtained the raw format of digital mammograms on average 3.1 (standard deviation (SD) =1.7) years prior to diagnosis or corresponding date for controls. We ascertained usual daily alcohol intake and other risk factors from a clinical questionnaire at time of mammography. We obtained BI-RADS density from clinical records and used VolparaTM to assess volumetric percent density, dense volume and non-dense volume from the mammogram. We examined the associations of daily alcohol intake (none, one or less, two or more drinks per day) and volumetric density phenotypes (per 1 SD) with BC risk using logistic regression (odds ratios, OR; 95% confidence intervals, CI; and trend tests). We examined deviation from multiplicative interaction using chi-squared tests. We evaluated mediation of the alcohol and BC association by volumetric density measures using logistic regression to estimate the association between alcohol use and BC with and without adjustment for density measures. Percent mediation was estimated using the differences in the log OR estimates from the two models. All models were adjusted for age, 1/BMI and menopause and matching factors. Analyses were also stratified by menopausal status.

Results: Alcohol intake was available on 2233 cases and 4562 controls, 88% of those eligible. BC cases and controls had similar age (57.2 years (SD=11.5) vs. 57.1 years (SD=11.5)), BMI (25.3 kg/m2 (SD=5.3) vs. 24.9 kg/ m2 (SD=5.1)) and race (69.0% vs. 68.2% Caucasian). Cases were more likely to drink alcohol daily than controls (52.1% vs. 49.0%), in particular two or more drinks per day (14.8% vs. 13.2%). Alcohol was associated with increased BC risk (OR=1.14, 95% CI:1.02-1.27, for one or less drinks per day and OR=1.22, 95% CI:1.05-1.42 for 2 or more drinks per day) compared to non-drinkers (p-trend=0.004). Percent volumetric density (OR=1.45 per SD, 95%CI: 1.36-1.56,) and dense volume (OR=1.30, 95% CI: 1.24-1.37) were also positively associated with BC risk; non-dense volume was inversely associated (OR=0.93, 95%CI: 0.86-1.01). Associations were similar by menopausal subgroup. There was no evidence for a differential association of alcohol and breast cancer risk by dense breasts assessed using any of the density phenotypes examined (all P’s>0.1). However, the association between alcohol and overall risk of BC was partially mediated by dense volume among all women (percent mediated=25%, P=0.01) and postmenopausal women (percent mediated=19%, P=0.03).

Conclusions: The association of daily alcohol intake and breast cancer risk was similar among women with dense and non-dense breasts. However, dense volume partially mediated the association between alcohol and risk of breast cancer, particularly among postmenopausal women, suggesting that alcohol partially influences breast cancer risk through changes in breast tissue composition.

Citation Format: Celine Marie Vachon, Christopher G. Scott, Stacey J Winham, John A Shepherd, Kathleen R Brandt, Matthew R Jensen, Carrie B Hruska, John J Heine, V. Shane Pankratz, Karla Kerlikowske. Association of daily alcohol intake, volumetric density and breast cancer risk [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P5-08-02.

Early-Life and Adult Anthropometrics in Relation to Mammographic Image Intensity Variation in the Nurses' Health Studies

BACKGROUND

The V measure captures grayscale intensity variation on a mammogram and is positively associated with breast cancer risk, independent of percent mammographic density (PMD), an established marker of breast cancer risk. We examined whether anthropometrics are associated with V, independent of PMD.

METHODS

The analysis included 1,700 premenopausal and 1,947 postmenopausal women without breast cancer within the Nurses' Health Study (NHS) and NHSII. Participants recalled their body fatness at ages 5, 10, and 20 years using a 9-level pictogram (level 1: most lean) and reported weight at age 18 years, current adult weight, and adult height. V was estimated by calculating standard deviation of pixels on screening mammograms. Linear mixed models were used to estimate beta coefficients (ß) and 95% confidence intervals (CI) for the relationships between anthropometric measures and V, adjusting for confounders and PMD.

RESULTS

V and PMD were positively correlated (Spearman r = 0.60). Higher average body fatness at ages 5 to 10 years (level ≥ 4.5 vs. 1) was significantly associated with lower V in premenopausal (ß = -0.32; 95% CI, -0.48 to -0.16) and postmenopausal (ß = -0.24; 95% CI, -0.37 to -0.10) women, independent of current body mass index (BMI) and PMD. Similar inverse associations were observed with average body fatness at ages 10 to 20 years and BMI at age 18 years. Current BMI was inversely associated with V, but the associations were largely attenuated after adjustment for PMD. Height was not associated with V.

CONCLUSIONS

Our data suggest that early-life body fatness may reflect lifelong impact on breast tissue architecture beyond breast density. However, further studies are needed to confirm the results.

IMPACT

This study highlights strong inverse associations of early-life adiposity with mammographic image intensity variation.

Radiologic Features of Small Pulmonary Nodules and Lung Cancer Risk in the National Lung Screening Trial: A Nested Case-Control Study

Purpose To extract radiologic features from small pulmonary nodules (SPNs) that did not meet the original criteria for a positive screening test and identify features associated with lung cancer risk by using data and images from the National Lung Screening Trial (NLST). Materials and Methods Radiologic features in SPNs in baseline low-dose computed tomography (CT) screening studies that did not meet NLST criteria to be considered a positive screening examination were extracted. SPNs were identified for 73 incident case patients who were given a diagnosis of lung cancer at either the first or second follow-up screening study and for 157 control subjects who had undergone three consecutive negative screening studies. Multivariable logistic regression was used to assess the association between radiologic features and lung cancer risk. All statistical tests were two sided. Results Nine features were significantly different between case patients and control subjects. Backward elimination followed by bootstrap resampling identified a reduced model of highly informative radiologic features with an area under the receiver operating characteristic curve of 0.932 (95% confidence interval [CI]: 0.88, 0.96), a specificity of 92.38% (95% CI: 52.22%, 84.91%), and a sensitivity of 76.55% (95% CI: 87.50%, 95.35%) that included total emphysema score (odds ratio [OR] = 1.71; 95% CI: 1.39, 2.01), attachment to vessel (OR = 2.41; 95% CI: 0.99, 5.81), nodule location (OR = 3.25; 95% CI: 1.09, 8.55), border definition (OR = 7.56; 95% CI: 1.89, 30.8), and concavity (OR = 2.58; 95% CI: 0.89, 5.64). Conclusion A set of clinically relevant radiologic features were identified that that can be easily scored in the clinical setting and may be of use to determine lung cancer risk among participants with SPNs. ^© RSNA, 2017 Online supplemental material is available for this article.

Longitudinal Changes in Volumetric Breast Density with Tamoxifen and Aromatase Inhibitors

Background: Reductions in breast density with tamoxifen and aromatase inhibitors may be an intermediate marker of treatment response. We compare changes in volumetric breast density among breast cancer cases using tamoxifen or aromatase inhibitors (AI) to untreated women without breast cancer.Methods: Breast cancer cases with a digital mammogram prior to diagnosis and after initiation of tamoxifen (n = 366) or AI (n = 403) and a sample of controls (n = 2170) were identified from the Mayo Clinic Mammography Practice and San Francisco Mammography Registry. Volumetric percent density (VPD) and dense breast volume (DV) were measured using Volpara (Matakina Technology) and Quantra (Hologic) software. Linear regression estimated the effect of treatment on annualized changes in density.Results: Premenopausal women using tamoxifen experienced annualized declines in VPD of 1.17% to 1.70% compared with 0.30% to 0.56% for controls and declines in DV of 7.43 to 15.13 cm³ compared with 0.28 to 0.63 cm³ in controls, for Volpara and Quantra, respectively. The greatest reductions were observed among women with ≥10% baseline density. Postmenopausal AI users had greater declines in VPD than controls (Volpara P = 0.02; Quantra P = 0.03), and reductions were greatest among women with ≥10% baseline density. Declines in VPD among postmenopausal women using tamoxifen were only statistically greater than controls when measured with Quantra.Conclusions: Automated software can detect volumetric breast density changes among women on tamoxifen and AI.Impact: If declines in volumetric density predict breast cancer outcomes, these measures may be used as interim prognostic indicators. Cancer Epidemiol Biomarkers Prev; 26(6); 930-7. ©2017 AACR.

Abstract 3424: Longitudinal changes in volumetric breast density with adjuvant endocrine therapy among women with breast cancer

Breast density represents the stromal and epithelial tissue in the breast and is a strong risk factor for breast cancer. Reductions in breast density attributable to tamoxifen (TAM) and aromatase inhibitors (AI) may be associated with reduced risk of first primary and subsequent breast cancer. Studies assessing breast density changes have principally used two-dimensional measures. We assess changes in breast density following initiation of TAM and AI using two automated volumetric density measures that have shown strong associations to breast cancer risk.

Breast cancer cases with a full field digital mammogram (FFDM) prior to diagnosis (index mammogram) and after initiation of TAM (n = 379) or AI (n = 425) were identified from the Mayo Clinic Mammography Practice and the San Francisco Mammography Registry. Volumetric percent density (VPD) and dense volume (DV) were measured on 4-view FFDM using VolparaTM (Matakina Technology) and QuantraTM (Hologic) automated software. We used linear regression to assess the effect of treatment on mean annualized change in VPD and DV (change from index to last mammogram divided by time interval) for each software type, adjusting for age, body mass index (BMI) and density at index mammogram.

The median time between index mammogram and cancer diagnosis was 0.6 months (IQR: 0.2, 2.2) and median time between index and last mammogram was 3 years (IQR: 2.0, 3.9). Women on TAM were younger, had lower BMI and higher baseline VPD and DV relative to women on AI (Table). An annual decrease in VPD and DV was observed with both TAM and AI. Both Volpara and Quantra estimated a similar magnitude of change in VPD in women on TAM and AI, and a greater change in DV with TAM.

Our findings suggest that both Volpara and Quantra can assess volumetric changes in breast density among women on hormone therapy. If declines in volumetric density correlate with a reduction in breast cancer risk, these automated measures could be used in clinical practice to assess response to therapy. Annualized changes in volumetric breast density estimated by linear regression.Tamoxifen (n = 379)Aromatase Inhibitors (n = 425)Baseline Median (IQR)Annualized Change (95% CI)*Baseline Median (IQR)Annualized Change (95% CI)*Age at Diagnosis50.0 (45.0, 60.0)–63.0 (58.0, 71.0)–Body Mass Index (BMI)23.6 (21.5, 26.8)–25.7 (22.7, 29.9)–Time Interval¥3.0 (2.1, 3.9)–3.0 (2.1, 3.9)–VolparaPercent Density (VPD,%)11.6 (6.8, 18.8)-0.17 (-0.27, -0.10)7.2 (5.0, 11.0)-0.19 (-0.29, -0.12)Dense Volume (DV, cm3)64.7 (45.4, 90.9)-0.90 (-1.45, -0.48)51.9 (38.9, 69.9)-0.52 (-0.93, -0.23)QuantraPercent Density (VPD,%)14.5 (9.2, 20.2)-0.42 (-0.59, -0.28)9.9 (7.1, 14.5)-0.38 (-0.54, -0.25)Dense Volume (DV, cm3)94.0 (58.0, 144.0)-2.20 (-3.52, -1.19)80.0 (49.0, 128.0)-0.95 (-1.85, -0.35)IQR = Interquartile range ¥ Median number of years between index mammogram and last mammogram post-initiation of therapy. *Annualized change estimated as change from index to last mammogram divided by time interval and adjusted for study site, age at diagnosis, BMI and density at index mammogram.

Citation Format: Natalie J. Engmann, Celine M. Vachon, Christopher G. Scott, Matthew R. Jensen, Lin Ma, Kathleen R. Brandt, Amir P. Mahmoudzadeh, Serghei Malkov, Dana H. Whaley, Carrie B. Hruska, Fang F. Wu, Stacey J. Winham, Diana L. Miglioretti, Aaron D. Norman, John J. Heine, John Shepherd, V Shane Pankratz, Karla Kerlikowske. Longitudinal changes in volumetric breast density with adjuvant endocrine therapy among women with breast cancer. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3424.

Association Between Computed Tomographic Features and Kirsten Rat Sarcoma Viral Oncogene Mutations in Patients With Stage I Lung Adenocarcinoma and Their Prognostic Value

BACKGROUND

We investigated the association between computed tomographic (CT) features and Kirsten rat sarcoma viral oncogene (KRAS) mutations in patients with stage I lung adenocarcinoma and their prognostic value.

PATIENTS AND METHODS

A total of 79 patients with pathologic stage I lung adenocarcinoma, available KRAS mutational status, preoperative CT images available, and survival data were included in the present study. Seven CT features, including spiculation, concavity, ground-glass opacity, bubble-like lucency, air bronchogram, pleural retraction, and pleural attachment, were evaluated. The association among the clinical characteristics, CT features, and mutational status was analyzed using Student's t test, the χ(2) test or Fisher's exact test, and logistic regression. The association among CT features, mutational status, and overall survival was analyzed using Kaplan-Meier survival curves with the log-rank test and Cox proportional hazard regression.

RESULTS

The prevalence of KRAS mutations was 41.77%. Spiculation was significantly associated with the presence of KRAS mutations (odds ratio, 2.99; 95% confidence interval [CI], 1.16-7.68). Although KRAS mutational status was not significantly associated with overall survival, the presence of pleural attachment was associated with an increased risk of death (hazard ratio, 2.46; 95% CI, 1.09-5.53). When analyzing KRAS mutational status and pleural attachment combined, patients with wild-type KRAS and no pleural attachment had significantly better survival than did those with wild-type KRAS and pleural attachment (P = .014).

CONCLUSION

These data suggest that spiculation is associated with KRAS mutations and pleural attachment is associated with overall survival in patients with stage I lung adenocarcinoma. Combining the analysis of KRAS mutational status and CT features could better predict survival.

Calibrated breast density methods for full field digital mammography: a system for serial quality control and inter-system generalization

PURPOSE

The authors are developing a system for calibrated breast density measurements using full field digital mammography (FFDM). Breast tissue equivalent (BTE) phantom images are used to establish baseline (BL) calibration curves at time zero. For a given FFDM unit, the full BL dataset is comprised of approximately 160 phantom images, acquired prior to calibrating prospective patient mammograms. BL curves are monitored serially to ensure they produce accurate calibration and require updating when calibration accuracy degrades beyond an acceptable tolerance, rather than acquiring full BL datasets repeatedly. BL updating is a special case of generalizing calibration datasets across FFDM units, referred to as cross-calibration. Serial monitoring, BL updating, and cross-calibration techniques were developed and evaluated.

METHODS

BL curves were established for three Hologic Selenia FFDM units at time zero. In addition, one set of serial phantom images, comprised of equal proportions of adipose and fibroglandular BTE materials (50/50 compositions) of a fixed height, was acquired biweekly and monitored with the cumulative sum (Cusum) technique. These 50/50 composition images were used to update the BL curves when the calibration accuracy degraded beyond a preset tolerance of ±4 standardized units. A second set of serial images, comprised of a wide-range of BTE compositions, was acquired biweekly to evaluate serial monitoring, BL updating, and cross-calibration techniques.

RESULTS

Calibration accuracy can degrade serially and is a function of acquisition technique and phantom height. The authors demonstrated that all heights could be monitored simultaneously while acquiring images of a 50/50 phantom with a fixed height for each acquisition technique biweekly, translating into approximately 16 image acquisitions biweekly per FFDM unit. The same serial images are sufficient for serial monitoring, BL updating, and cross-calibration. Serial calibration accuracy was maintained within ±4 standardized unit variation from the ideal when applying BL updating. BL updating is a special case of cross-calibration; the BL dataset of unit 1 can be converted to the BL dataset for another similar unit (i.e., unit 2) at any given time point using the 16 serial monitoring 50/50 phantom images of unit 2 (or vice versa) acquired near this time point while maintaining the ±4 standardized unit tolerance.

CONCLUSIONS

A methodology for monitoring and maintaining serial calibration accuracy for breast density measurements was evaluated. Calibration datasets for a given unit can be translated forward in time with minimal phantom imaging effort. Similarly, cross-calibration is a method for generalizing calibration datasets across similar units without additional phantom imaging. This methodology will require further evaluation with mammograms for complete validation.

Breast Imaging Reporting and Data System (BI-RADS) breast composition descriptors: automated measurement development for full field digital mammography

PURPOSE

The Breast Imaging Reporting and Data System (BI-RADS) breast composition descriptors are used for standardized mammographic reporting and are assessed visually. This reporting is clinically relevant because breast composition can impact mammographic sensitivity and is a breast cancer risk factor. New techniques are presented and evaluated for generating automated BI-RADS breast composition descriptors using both raw and calibrated full field digital mammography (FFDM) image data.

METHODS

A matched case-control dataset with FFDM images was used to develop three automated measures for the BI-RADS breast composition descriptors. Histograms of each calibrated mammogram in the percent glandular (pg) representation were processed to create the new BR(pg) measure. Two previously validated measures of breast density derived from calibrated and raw mammograms were converted to the new BR(vc) and BR(vr) measures, respectively. These three measures were compared with the radiologist-reported BI-RADS compositions assessments from the patient records. The authors used two optimization strategies with differential evolution to create these measures: method-1 used breast cancer status; and method-2 matched the reported BI-RADS descriptors. Weighted kappa (κ) analysis was used to assess the agreement between the new measures and the reported measures. Each measure's association with breast cancer was evaluated with odds ratios (ORs) adjusted for body mass index, breast area, and menopausal status. ORs were estimated as per unit increase with 95% confidence intervals.

RESULTS

The three BI-RADS measures generated by method-1 had κ between 0.25-0.34. These measures were significantly associated with breast cancer status in the adjusted models: (a) OR = 1.87 (1.34, 2.59) for BR(pg); (b) OR = 1.93 (1.36, 2.74) for BR(vc); and (c) OR = 1.37 (1.05, 1.80) for BR(vr). The measures generated by method-2 had κ between 0.42-0.45. Two of these measures were significantly associated with breast cancer status in the adjusted models: (a) OR = 1.95 (1.24, 3.09) for BR(pg); (b) OR = 1.42 (0.87, 2.32) for BR(vc); and (c) OR = 2.13 (1.22, 3.72) for BR(vr). The radiologist-reported measures from the patient records showed a similar association, OR = 1.49 (0.99, 2.24), although only borderline statistically significant.

CONCLUSIONS

A general framework was developed and validated for converting calibrated mammograms and continuous measures of breast density to fully automated approximations for the BI-RADS breast composition descriptors. The techniques are general and suitable for a broad range of clinical and research applications.

Automated Percentage of Breast Density Measurements for Full-field Digital Mammography Applications

RATIONALE AND OBJECTIVES

Increased mammographic breast density is a significant risk factor for breast cancer. A reproducible, accurate, and automated breast density measurement is required for full-field digital mammography (FFDM) to support clinical applications. We evaluated a novel automated percentage of breast density measure (PDa) and made comparisons with the standard operator-assisted measure (PD) using FFDM data.

METHODS

We used a nested breast cancer case-control study matched on age, year of mammogram and diagnosis with images acquired from a specific direct x-ray conversion FFDM technology. PDa was applied to the raw and clinical display (or processed) representation images. We evaluated the transformation (pixel mapping) of the raw image, giving a third representation (raw-transformed), to improve the PDa performance using differential evolution optimization. We applied PD to the raw and clinical display images as a standard for measurement comparison. Conditional logistic regression was used to estimate the odd ratios (ORs) for breast cancer with 95% confidence intervals (CI) for all measurements; analyses were adjusted for body mass index. PDa operates by evaluating signal-dependent noise (SDN), captured as local signal variation. Therefore, we characterized the SDN relationship to understand the PDa performance as a function of data representation and investigated a variation analysis of the transformation.

RESULTS

The associations of the quartiles of operator-assisted PD with breast cancer were similar for the raw (OR: 1.00 [ref.]; 1.59 [95% CI, 0.93-2.70]; 1.70 [95% CI, 0.95-3.04]; 2.04 [95% CI, 1.13-3.67]) and clinical display (OR: 1.00 [ref.]; 1.31 [95% CI, 0.79-2.18]; 1.14 [95% CI, 0.65-1.98]; 1.95 [95% CI, 1.09-3.47]) images. PDa could not be assessed on the raw images without preprocessing. However, PDa had similar associations with breast cancer when assessed on 1) raw-transformed (OR: 1.00 [ref.]; 1.27 [95% CI, 0.74-2.19]; 1.86 [95% CI, 1.05-3.28]; 3.00 [95% CI, 1.67-5.38]) and 2) clinical display (OR: 1.00 [ref.]; 1.79 [95% CI, 1.04-3.11]; 1.61 [95% CI, 0.90-2.88]; 2.94 [95% CI, 1.66-5.19]) images. The SDN analysis showed that a nonlinear relationship between the mammographic signal and its variation (ie, the biomarker for the breast density) is required for PDa. Although variability in the transform influenced the respective PDa distribution, it did not affect the measurement's association with breast cancer.

CONCLUSIONS

PDa assessed on either raw-transformed or clinical display images is a valid automated breast density measurement for a specific FFDM technology and compares well against PD. Further work is required for measurement generalization.

Comparison of percent density from raw and processed full-field digital mammography data

Introduction

Mammographic density has been established as a strong risk factor for breast cancer, primarily using digitized film mammograms. Full-field digital mammography (FFDM) is replacing film mammography, has different properties than film, and provides both raw and processed clinical display representation images. We evaluated and compared FFDM raw and processed breast density measures and their associations with breast cancer.

Methods

A case-control study of 180 cases and 180 controls matched by age, postmenopausal hormone use, and screening history was conducted. Mammograms were acquired from a General Electric Senographe 2000D FFDM unit. Percent density (PD) was assessed for each FFDM representation using the operator-assisted Cumulus method. Reproducibility within image type (n = 80) was assessed using Lin's concordance correlation coefficient (r_c). Correlation of PD between image representations (n = 360) was evaluated using Pearson's correlation coefficient (r) on the continuous measures and the weighted kappa statistic (κ) for quartiles. Conditional logistic regression was used to estimate odds ratios (ORs) for the PD and breast cancer associations for both image representations with 95% confidence intervals. The area under the receiver operating characteristic curve (AUC) was used to assess the discriminatory accuracy.

Results

Percent density from the two representations provided similar intra-reader reproducibility (r_c= 0.92 for raw and r_c= 0.87 for processed images) and was correlated (r = 0.82 and κ = 0.64). When controlling for body mass index, the associations of quartiles of PD with breast cancer and discriminatory accuracy were similar for the raw (OR: 1.0 (ref.), 2.6 (1.2 to 5.4), 3.1 (1.4 to 6.8), 4.7 (2.1 to 10.6); AUC = 0.63) and processed representations (OR: 1.0 (ref.), 2.2 (1.1 to 4.1), 2.2 (1.1 to 4.4), 3.1 (1.5 to 6.6); AUC = 0.64).

Conclusions

Percent density measured with an operator-assisted method from raw and processed FFDM images is reproducible and correlated. Both percent density measures provide similar associations with breast cancer.

The influence of mammogram acquisition on the mammographic density and breast cancer association in the Mayo Mammography Health Study cohort

Introduction

Mammographic density is a strong risk factor for breast cancer. Image acquisition technique varies across mammograms to limit radiation and produce a clinically useful image. We examined whether acquisition technique parameters at the time of mammography were associated with mammographic density and whether the acquisition parameters confounded the density and breast cancer association.

Methods

We examined this question within the Mayo Mammography Health Study (MMHS) cohort, comprised of 19,924 women (51.2% of eligible) seen in the Mayo Clinic mammography screening practice from 2003 to 2006. A case-cohort design, comprising 318 incident breast cancers diagnosed through December 2009 and a random subcohort of 2,259, was used to examine potential confounding of mammogram acquisition technique parameters (x-ray tube voltage peak (kVp), milliampere-seconds (mAs), thickness and compression force) on the density and breast cancer association. The Breast Imaging Reporting and Data System four-category tissue composition measure (BI-RADS) and percent density (PD) (Cumulus program) were estimated from screen-film mammograms at time of enrollment. Spearman correlation coefficients (r) and means (standard deviations) were used to examine the relationship of density measures with acquisition parameters. Hazard ratios (HR) and C-statistics were estimated using Cox proportional hazards regression, adjusting for age, menopausal status, body mass index and postmenopausal hormones. A change in the HR of at least 15% indicated confounding.

Results

Adjusted PD and BI-RADS density were associated with breast cancer (p-trends < 0.001), with a 3 to 4-fold increased risk in the extremely dense vs. fatty BI-RADS categories (HR: 3.0, 95% CI, 1.7 - 5.1) and the ≥ 25% vs. ≤ 5% PD categories (HR: 3.8, 95% CI, 2.5 - 5.9). Of the acquisition parameters, kVp was not correlated with PD (r = 0.04, p = 0.07). Although thickness (r = -0.27, p < 0.001), compression force (r = -0.16, p < 0.001), and mAs (r = -0.06, p = 0.008) were inversely correlated with PD, they did not confound the PD or BI-RADS associations with breast cancer and their inclusion did not improve discriminatory accuracy. Results were similar for associations of dense and non-dense area with breast cancer.

Conclusions

We confirmed a strong association between mammographic density and breast cancer risk that was not confounded by mammogram acquisition technique.

A novel automated mammographic density measure and breast cancer risk

BACKGROUND

Mammographic breast density is a strong breast cancer risk factor but is not used in the clinical setting, partly because of a lack of standardization and automation. We developed an automated and objective measurement of the grayscale value variation within a mammogram, evaluated its association with breast cancer, and compared its performance with that of percent density (PD).

METHODS

Three clinic-based studies were included: a case-cohort study of 217 breast cancer case subjects and 2094 non-case subjects and two case-control studies comprising 928 case subjects and 1039 control subjects and 246 case subjects and 516 control subjects, respectively. Percent density was estimated from digitized mammograms using the computer-assisted Cumulus thresholding program, and variation was estimated from an automated algorithm. We estimated hazards ratios (HRs), odds ratios (ORs), the area under the receiver operating characteristic curve (AUC), and 95% confidence intervals (CIs) using Cox proportional hazards models for the cohort and logistic regression for case-control studies, with adjustment for age and body mass index. We performed a meta-analysis using random study effects to obtain pooled estimates of the associations between the two mammographic measures and breast cancer. All statistical tests were two-sided.

RESULTS

The variation measure was statistically significantly associated with the risk of breast cancer in all three studies (highest vs lowest quartile: HR = 2.0 [95% CI = 1.3 to 3.1]; OR = 2.7 [95% CI = 2.1 to 3.6]; OR = 2.4 [95% CI = 1.4 to 3.9]; [corrected] all P (trend) < .001). [corrected]. The risk estimates and AUCs for the variation measure were similar to [corrected] those for percent density (AUCs for variation = 0.60-0.62 and [corrected] AUCs for percent density = 0.61-0.65). [corrected]. A meta-analysis of the three studies demonstrated similar associations [corrected] between variation and breast cancer (highest vs lowest quartile: RR = 1.8, 95% CI = 1.4 to 2.3) and [corrected] percent density and breast cancer (highest vs lowest quartile: RR = 2.3, 95% CI = 1.9 to 2.9).

CONCLUSION

The association between the automated variation measure and the risk of breast cancer is at least as strong as that for percent density. Efforts to further evaluate and translate the variation measure to the clinical setting are warranted.

Perturbation and candidate analysis to combat overfitting of gene expression microarray data

Analysis of gene expression microarray datasets presents the high risk of over-fitting (spurious patterns) because of their feature-rich but case-poor nature. This paper describes our ongoing efforts to develop a method to combat over-fitting and determine the strongest signal in the dataset. A GA-SVM hybrid along with Gaussian noise (manual noise gain) is used to discover feature sets of minimal size that accurately classifies the cases under cross-validation. Initial results on a colorectal cancer dataset shows that the strongest signal (modest number of candidates) can be found by a binary search.

Statistical learning methods as a preprocessing step for survival analysis: evaluation of concept using lung cancer data

BACKGROUND

Statistical learning (SL) techniques can address non-linear relationships and small datasets but do not provide an output that has an epidemiologic interpretation.

METHODS

A small set of clinical variables (CVs) for stage-1 non-small cell lung cancer patients was used to evaluate an approach for using SL methods as a preprocessing step for survival analysis. A stochastic method of training a probabilistic neural network (PNN) was used with differential evolution (DE) optimization. Survival scores were derived stochastically by combining CVs with the PNN. Patients (n = 151) were dichotomized into favorable (n = 92) and unfavorable (n = 59) survival outcome groups. These PNN derived scores were used with logistic regression (LR) modeling to predict favorable survival outcome and were integrated into the survival analysis (i.e. Kaplan-Meier analysis and Cox regression). The hybrid modeling was compared with the respective modeling using raw CVs. The area under the receiver operating characteristic curve (Az) was used to compare model predictive capability. Odds ratios (ORs) and hazard ratios (HRs) were used to compare disease associations with 95% confidence intervals (CIs).

RESULTS

The LR model with the best predictive capability gave Az = 0.703. While controlling for gender and tumor grade, the OR = 0.63 (CI: 0.43, 0.91) per standard deviation (SD) increase in age indicates increasing age confers unfavorable outcome. The hybrid LR model gave Az = 0.778 by combining age and tumor grade with the PNN and controlling for gender. The PNN score and age translate inversely with respect to risk. The OR = 0.27 (CI: 0.14, 0.53) per SD increase in PNN score indicates those patients with decreased score confer unfavorable outcome. The tumor grade adjusted hazard for patients above the median age compared with those below the median was HR = 1.78 (CI: 1.06, 3.02), whereas the hazard for those patients below the median PNN score compared to those above the median was HR = 4.0 (CI: 2.13, 7.14).

CONCLUSION

We have provided preliminary evidence showing that the SL preprocessing may provide benefits in comparison with accepted approaches. The work will require further evaluation with varying datasets to confirm these findings.

Calibrated measures for breast density estimation

RATIONALE AND OBJECTIVES

Breast density is a significant breast cancer risk factor measured from mammograms. Evidence suggests that the spatial variation in mammograms may also be associated with risk. We investigated the variation in calibrated mammograms as a breast cancer risk factor and explored its relationship with other measures of breast density using full field digital mammography (FFDM).

MATERIALS AND METHODS

A matched case-control analysis was used to assess a spatial variation breast density measure in calibrated FFDM images, normalized for the image acquisition technique variation. Three measures of breast density were compared between cases and controls: (a) the calibrated average measure, (b) the calibrated variation measure, and (c) the standard percentage of breast density (PD) measure derived from operator-assisted labeling. Linear correlation and statistical relationships between these three breast density measures were also investigated.

RESULTS

Risk estimates associated with the lowest to highest quartiles for the calibrated variation measure were greater in magnitude (odds ratios: 1.0 [ref.], 3.5, 6.3, and 11.3) than the corresponding risk estimates for quartiles of the standard PD measure (odds ratios: 1.0 [ref.], 2.3, 5.6, and 6.5) and the calibrated average measure (odds ratios: 1.0 [ref.], 2.4, 2.3, and 4.4). The three breast density measures were highly correlated, showed an inverse relationship with breast area, and related by a mixed distribution relationship.

CONCLUSION

The three measures of breast density capture different attributes of the same data field. These preliminary findings indicate the variation measure is a viable automated method for assessing breast density. Insights gained by this work may be used to develop a standard for measuring breast density.

A quantitative description of the percentage of breast density measurement using full-field digital mammography

RATIONALE AND OBJECTIVES

Breast density is a significant breast cancer risk factor that is measured from mammograms. However, uncertainty remains in both understanding its underlying physical properties as it relates to the breast and determining the optimal method for its measurement. A quantitative description of the information captured by the standard operator-assisted percentage of breast density (PD) measure was developed using full-field digital mammography (FFDM) images that were calibrated to adjust for interimage acquisition technique differences.

MATERIALS AND METHODS

The information captured by the standard PD measure was quantified by developing a similar measure of breast density (PD(c)) from calibrated mammograms automatically by applying a static threshold to each image. The specific threshold was estimated by first sampling the probability distributions for breast tissue in calibrated mammograms. A percent glandular (PG) measure of breast density was also derived from calibrated mammograms. The PD, PD(c), and PG breast density measures were compared using both linear correlation (R) and quartile odds ratio measures derived from a matched case-control study.

RESULTS

The standard PD measure is an estimate of the number of pixel values above a fixed idealized x-ray attenuation fraction. There was significant correlation (P < .0001) between the PD(c)-PD (r = 0.78), PD(c)-PG (r = 0.87), and PD-PG (r = 0.71) measures of breast density. Risk estimates associated with the lowest to highest quartiles for the PD(c) measure (odds ratio [OR]: 1.0 ref., 3.4, 3.6, and 5.6), and the standard PD measure (OR 1.0 ref., 2.9, 4.8, and 5.1) were similar and greater than that of the calibrated PG measure (OR 1.0 ref., 2.0, 2.4, and 2.4).

CONCLUSIONS

The information captured by the standard PD measure was quantified as it relates to calibrated mammograms and used to develop an automated method for measuring breast density. These findings represent an initial step for developing an automated measure built on an established calibration platform. A fully developed automated measure may be useful for both research- and clinical-based risk applications.

Abstract 3716: The Mayo Mammography Health Study (MMHS): A prospective cohort study on mammographic breast density and breast cancer

INTRODUCTION: Mammographic breast density is a strong risk factor for breast cancer (BC). We established the Mayo Mammography Health Study (MMHS) cohort study at the Mayo Clinic in Rochester, Minnesota (MN) to examine the association of breast density with BC. We evaluated the influence of the image acquisition technique on the density and BC association.

METHODS: From October 2003 to September 2006, all women scheduled for screening mammography at the Mayo Clinic were invited to participate. Eligible women were residents of MN, Iowa (IA) or Wisconsin (WI); age 35+; and had no personal history of BC. A risk factor questionnaire, consent form and permission to link to tumor registries were obtained. Incident BC was identified through 2009 by linkage to the Mayo and state cancer registries. A case-cohort of all incident BCs and 2300 randomly selected women (the subcohort) were used to examine the association of breast density and BC using digitized film mammograms at the time of enrollment while controlling for the influence of the acquisition parameters (peak kilovoltage, milliampere-second, and compressed breast thickness). Two density measures were considered: a quantitative percent density (PD) measure estimated using the computer-assisted thresholding program, Cumulus (University of Toronto), and the 4-category clinical BI-RADS measure. Proportional hazards regression was used to calculate hazards ratios (HR) and 95% confidence intervals (CI) for BC associated with quartiles of PD (0-6.2%,6.3-14.9%, 15.0-25.7% and 25.8%+) and BI-RADS categories 1-4 (almost entirely fatty to extremely dense). All models included age, postmenopausal hormone use (PMH), BMI, and menopausal status. The influence of acquisition parameters was evaluated by examining models with and without their inclusion.

RESULTS: A total of 20,982 (50%) women participated in MMHS; 1058 (5%) with a prior history of BC were excluded, for a total cohort of 19,924. Compared to nonresponders, responders were younger (57.5 vs. 58.4 yrs), more likely to have ever used PMH (45% vs. 33%), and more likely to have a BC family history (19% vs 16%). The case-cohort consisted of 317 incident cases and 2300 women in the subcohort. Women were excluded from the current analysis if PD could not be estimated or if acquisition parameters were not available, leaving 249 cases and 1937 in the subcohort. As expected, PD was associated with BC [HR (95% CI): 1.0 (REF), 2.1 (1.4-3.1), 3.0 (2.0-4.5), and 4.6 (3.0-7.0) for quartiles; p-trend&lt;0.001]. Controlling for acquisition parameters attenuated the association [HR (95% CI): 1.0 (REF), 2.3 (1.5-3.4), 2.4 (1.6-3.7), and 3.0 (1.8-5.0) for quartiles; p-trend&lt;0.001]. Results for BI-RADS density were similar to those for PD.

CONCLUSION: This study confirms that breast density is a significant risk factor for BC and demonstrates that the acquisition technique confounds the density and BC risk association.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3716. doi:10.1158/1538-7445.AM2011-3716

Statistical learning techniques applied to epidemiology: a simulated case-control comparison study with logistic regression

Background

When investigating covariate interactions and group associations with standard regression analyses, the relationship between the response variable and exposure may be difficult to characterize. When the relationship is nonlinear, linear modeling techniques do not capture the nonlinear information content. Statistical learning (SL) techniques with kernels are capable of addressing nonlinear problems without making parametric assumptions. However, these techniques do not produce findings relevant for epidemiologic interpretations. A simulated case-control study was used to contrast the information embedding characteristics and separation boundaries produced by a specific SL technique with logistic regression (LR) modeling representing a parametric approach. The SL technique was comprised of a kernel mapping in combination with a perceptron neural network. Because the LR model has an important epidemiologic interpretation, the SL method was modified to produce the analogous interpretation and generate odds ratios for comparison.

Results

The SL approach is capable of generating odds ratios for main effects and risk factor interactions that better capture nonlinear relationships between exposure variables and outcome in comparison with LR.

Conclusions

The integration of SL methods in epidemiology may improve both the understanding and interpretation of complex exposure/disease relationships.

Effective radiation attenuation calibration for breast density: compression thickness influences and correction

Background

Calibrating mammograms to produce a standardized breast density measurement for breast cancer risk analysis requires an accurate spatial measure of the compressed breast thickness. Thickness inaccuracies due to the nominal system readout value and compression paddle orientation induce unacceptable errors in the calibration.

Method

A thickness correction was developed and evaluated using a fully specified two-component surrogate breast model. A previously developed calibration approach based on effective radiation attenuation coefficient measurements was used in the analysis. Water and oil were used to construct phantoms to replicate the deformable properties of the breast. Phantoms consisting of measured proportions of water and oil were used to estimate calibration errors without correction, evaluate the thickness correction, and investigate the reproducibility of the various calibration representations under compression thickness variations.

Results

The average thickness uncertainty due to compression paddle warp was characterized to within 0.5 mm. The relative calibration error was reduced to 7% from 48-68% with the correction. The normalized effective radiation attenuation coefficient (planar) representation was reproducible under intra-sample compression thickness variations compared with calibrated volume measures.

Conclusion

Incorporating this thickness correction into the rigid breast tissue equivalent calibration method should improve the calibration accuracy of mammograms for risk assessments using the reproducible planar calibration measure.

Cumulative sum quality control for calibrated breast density measurements

PURPOSE

Breast density is a significant breast cancer risk factor. Although various methods are used to estimate breast density, there is no standard measurement for this important factor. The authors are developing a breast density standardization method for use in full field digital mammography (FFDM). The approach calibrates for interpatient acquisition technique differences. The calibration produces a normalized breast density pixel value scale. The method relies on first generating a baseline (BL) calibration dataset, which required extensive phantom imaging. Standardizing prospective mammograms with calibration data generated in the past could introduce unanticipated error in the standardized output if the calibration dataset is no longer valid.

METHODS

Sample points from the BL calibration dataset were imaged approximately biweekly over an extended timeframe. These serial samples were used to evaluate the BL dataset reproducibility and quantify the serial calibration accuracy. The cumulative sum (Cusum) quality control method was used to evaluate the serial sampling.

RESULTS

There is considerable drift in the serial sample points from the BL calibration dataset that is x-ray beam dependent. Systematic deviation from the BL dataset caused significant calibration errors. This system drift was not captured with routine system quality control measures. Cusum analysis indicated that the drift is a sign of system wear and eventual x-ray tube failure.

CONCLUSIONS

The BL calibration dataset must be monitored and periodically updated, when necessary, to account for sustained system variations to maintain the calibration accuracy.

Texture features from mammographic images and risk of breast cancer

Mammographic percent density (PD) is a strong risk factor for breast cancer, but there has been relatively little systematic evaluation of other features in mammographic images that might additionally predict breast cancer risk. We evaluated the association of a large number of image texture features with risk of breast cancer using a clinic-based case-control study of digitized film mammograms, all with screening mammograms before breast cancer diagnosis. The sample was split into training (123 cases and 258 controls) and validation (123 cases and 264 controls) data sets. Age-adjusted and body mass index (BMI)-adjusted odds ratios (OR) per SD change in the feature, 95% confidence intervals, and the area under the receiver operator characteristic curve (AUC) were obtained using logistic regression. A bootstrap approach was used to identify the strongest features in the training data set, and results for features that validated in the second half of the sample were reported using the full data set. The mean age at mammography was 64.0+/-10.2 years, and the mean time from mammography to breast cancer was 3.7+/-1.0 (range, 2.0-5.9 years). PD was associated with breast cancer risk (OR, 1.49; 95% confidence interval, 1.25-1.78). The strongest features that validated from each of several classes (Markovian, run length, Laws, wavelet, and Fourier) showed similar ORs as PD and predicted breast cancer at a similar magnitude (AUC=0.58-0.60) as PD (AUC=0.58). All of these features were automatically calculated (unlike PD) and measure texture at a coarse scale. These features were moderately correlated with PD (r=0.39-0.76), and after adjustment for PD, each of the features attenuated only slightly and retained statistical significance. However, simultaneous inclusion of these features in a model with PD did not significantly improve the ability to predict breast cancer.

An automated approach for estimation of breast density

Breast density is a strong risk factor for breast cancer; however, no standard assessment method exists. An automated breast density method was modified and compared with a semi-automated, user-assisted thresholding method (Cumulus method) and the Breast Imaging Reporting and Data System four-category tissue composition measure for their ability to predict future breast cancer risk. The three estimation methods were evaluated in a matched breast cancer case-control (n = 372 and n = 713, respectively) study at the Mayo Clinic using digitized film mammograms. Mammograms from the craniocaudal view of the noncancerous breast were acquired on average 7 years before diagnosis. Two controls with no previous history of breast cancer from the screening practice were matched to each case on age, number of previous screening mammograms, final screening exam date, menopausal status at this date, interval between earliest and latest available mammograms, and residence. Both Pearson linear correlation (R) and Spearman rank correlation (r) coefficients were used for comparing the three methods as appropriate. Conditional logistic regression was used to estimate the risk for breast cancer (odds ratios and 95% confidence intervals) associated with the quartiles of percent breast density (automated breast density method, Cumulus method) or Breast Imaging Reporting and Data System categories. The area under the receiver operator characteristic curve was estimated and used to compare the discriminatory capabilities of each approach. The continuous measures (automated breast density method and Cumulus method) were highly correlated with each other (R = 0.70) but less with Breast Imaging Reporting and Data System (r = 0.49 for automated breast density method and r = 0.57 for Cumulus method). Risk estimates associated with the lowest to highest quartiles of automated breast density method were greater in magnitude [odds ratios: 1.0 (reference), 2.3, 3.0, 5.2; P trend < 0.001] than the corresponding quartiles for the Cumulus method [odds ratios: 1.0 (reference), 1.7, 2.1, and 3.8; P trend < 0.001] and Breast Imaging Reporting and Data System [odds ratios: 1.0 (reference), 1.6, 1.5, 2.6; P trend < 0.001] method. However, all methods similarly discriminated between case and control status; areas under the receiver operator characteristic curve were 0.64, 0.63, and 0.61 for automated breast density method, Cumulus method, and Breast Imaging Reporting and Data System, respectively. The automated breast density method is a viable option for quantitatively assessing breast density from digitized film mammograms.

Effective x-ray attenuation coefficient measurements from two full field digital mammography systems for data calibration applications

Background

Breast density is a significant breast cancer risk factor. Currently, there is no standard method for measuring this important factor. Work presented here represents an essential component of an ongoing project that seeks to determine the appropriate method for calibrating (standardizing) mammography image data to account for the x-ray image acquisition influences. Longer term goals of this project are to make accurate breast density measurements in support of risk studies.

Methods

Logarithmic response calibration curves and effective x-ray attenuation coefficients were measured from two full field digital mammography (FFDM) systems with breast tissue equivalent phantom imaging and compared. Normalization methods were studied to assess the possibility of reducing the amount of calibration data collection. The percent glandular calibration map functional form was investigated. Spatial variations in the calibration data were used to assess the uncertainty in the calibration application by applying error propagation analyses.

Results

Logarithmic response curves are well approximated as linear. Measured effective x-ray attenuation coefficients are characteristic quantities independent of the imaging system and are in agreement with those predicted numerically. Calibration data collection can be reduced by applying a simple normalization technique. The calibration map is well approximated as linear. Intrasystem calibration variation was on the order of four percent, which was approximately half of the intersystem variation.

Conclusion

FFDM systems provide a quantitative output, and the calibration quantities presented here may be used for data acquired on similar FFDM systems.

New statistical learning theory paradigms adapted to breast cancer diagnosis/classification using image and non-image clinical data

The automated decision paradigms presented in this work address the false positive (FP) biopsy occurrence in diagnostic mammography. An EP/ES stochastic hybrid and two kernelized Partial Least Squares (K-PLS) paradigms were investigated with following studies: methodology performance comparisonsautomated diagnostic accuracy assessments with two data sets. The findings showed: the new hybrid produced comparable results more rapidlythe new K-PLS paradigms train and operate Essentially in real time for the data sets studied. Both advancements are essential components for eventually achieving the FP reduction goal, while maintaining acceptable diagnostic sensitivities.

Effective x-ray attenuation measurements with full field digital mammography

This work shows that effective x-ray attenuation coefficients may be estimated by applying Beer's Law to phantom image data acquired with the General Electric Senographe 2000D full field digital mammography system. Theoretical developments are provided indicating that an approximate form of the Beer's relation holds for polychromatic x-ray beams. The theoretical values were compared with experimentally determined measured values, which were estimated at various detector locations. The measured effective attenuation coefficients are in agreement with those estimated with theoretical developments and numerical integration. The work shows that the measured quantities show little spatial variation. The main ideas are demonstrated with polymethylmethacrylate and breast tissue equivalent phantom imaging experiments. The work suggests that the effective attenuation coefficients may be used as known values for radiometric standardization applications that compensate for the image acquisition influences. The work indicates that it is possible to make quantitative attenuation coefficient measurements from a system designed for clinical purposes.

Aspects of signal-dependent noise characterization

An automated method is presented for analyzing signal-dependent noise. Signal-dependent noise is present in many types of data-acquisition processes and has been investigated by other researchers with various methods. Regardless of the noise analysis methods, often the starting point is based on a particular signal-dependent noise model, which also forms the basis for our study. The approach determines whether the estimated noise variance is dependent on the signal by approximating the functional relation within the constraints of the assumed signal-noise model. The method relies on the Fourier attributes of the signal and noise and uses the wavelet expansion for separating these components. The technique does not rely on the underlying noise and signal probability distributions. Two-dimensional simulations as well as mammography data are used to illustrate the merits of the approach.

Improved Interpretation of Digitized Mammography with Wavelet Processing:A Localization Response Operating Characteristic Study

OBJECTIVE

Our objective was the implementation and evaluation of a novel enhancement technique for improved interpretation of high-resolution digitized mammograms from computer monitors.

MATERIALS AND METHODS

A wavelet algorithm was designed to attenuate the image spectral characteristics responsible for the long-range image correlation that often interferes with digital display. The algorithm was evaluated with a localization response operating characteristic (LROC) experiment with 500 negative, benign, and cancer cases with masses and calcification clusters. Three observers reviewed the original and wavelet-enhanced images on a 5-Mpixel monitor using a custom-made workstation user interface.

RESULTS

Performance indexes were estimated for four different case combinations, each observer, and each interpretation mode. Wavelet enhancement improved the performance of all observers in all case combinations. Detection accuracy ranged from 0.678 to 0.827 for the unprocessed original data and 0.709-0.871 for the enhanced cases. Localization accuracy ranged from 0.547 to 0.785 for the original images and 0.568-0.847 for the enhanced cases, yielding increases of 5-15%. The difference between enhanced and original performances was statistically significant at the 0.10 level and in a few combinations at the 0.05 level.

CONCLUSION

Soft-copy digitized mammography could replace standard film mammography under appropriate display parameters and conditions. The optimization of the soft-copy quality is expected to require more advanced processing techniques than standard gray-scale adjustments. Wavelet-based algorithms, such as the one proposed here, offer better soft-copy quality than the originals and a better starting point for additional manual gray-scale adjustments or automated postprocessing.

Spectral analysis of full field digital mammography data

The spectral content of mammograms acquired from using a full field digital mammography (FFDM) system are analyzed. Fourier methods are used to show that the FFDM image power spectra obey an inverse power law; in an average sense, the images may be considered as 1/f fields. Two data representations are analyzed and compared (1) the raw data, and (2) the logarithm of the raw data. Two methods are employed to analyze the power spectra (1) a technique based on integrating the Fourier plane with octave ring sectioning developed previously, and (2) an approach based on integrating the Fourier plane using rings of constant width developed for this work. Both methods allow theoretical modeling. Numerical analysis indicates that the effects due to the transformation influence the power spectra measurements in a statistically significant manner in the high frequency range. However, this effect has little influence on the inverse power law estimation for a given image regardless of the data representation or the theoretical analysis approach. The analysis is presented from two points of view (1) each image is treated independently with the results presented as distributions, and (2) for a given representation, the entire image collection is treated as an ensemble with the results presented as expected values. In general, the constant ring width analysis forms the foundation for a spectral comparison method for finding spectral differences, from an image distribution sense, after applying a nonlinear transformation to the data. The work also shows that power law estimation may be influenced due to the presence of noise in the higher frequency range, which is consistent with the known attributes of the detector efficiency. The spectral modeling and inverse power law determinations obtained here are in agreement with that obtained from the analysis of digitized film-screen images presented previously. The form of the power spectrum for a given image is approximately l/f2beta with beta approximately 1.4-1.5.

Mammographic tissue, breast cancer risk, serial image analysis, and digital mammography. Part 1. Tissue and related risk factors

This work is presented as a sequence of two parts. In this leading section, a review of the breast tissue-risk research is provided. Although controversy remains, there is substantial evidence indicating that dense mammographic tissue (a) is a breast cancer risk factor that is at least similar, if not greater, in magnitude with the other known breast cancer risk factors and (b) may be a partial biomarker for some of the other risk factors. Understanding these influences may provide a mechanism for measuring the dynamics of breast cancer risk. The totality of this work is to provide support for an automated serial mammography study under way at the authors' institution, where digital mammographic images are acquired with a full-field digital mammography system. This is a filmless imaging system, where the image is acquired in digital format. This electronic imaging acquisition system provides a prime opportunity to easily couple and manipulate the image data with patient information such as risk probability analysis or other pertinent personal history data for improved automated decision making. In this leading section, the main focus is on understanding elements that will assist in fusing risk probability analysis with automated computer-aided diagnosis. The evidence indicates that there are many factors that influence breast tissue at any given time and thus have the ability to alter the associated radiographic image appearance over time. At the initiation of the serial study it was clear that the authors did not fully understand the nature of the problem: automatically comparing similar mammographic scenes acquired at different times. In the second part of this sequence, the more time-related tissue influences are reviewed.

Mammographic tissue, breast cancer risk, serial image analysis, and digital mammography. Part 2. Serial breast tissue change and related temporal influences

The work presented herein is the second part of a review of breast tissue-related cancer-risk research. Briefly, in part 1, the tissue-risk research is discussed. In this part, factors that influence temporal breast tissue change are reviewed. Since breast composition is correlated with some of the known risk factors, understanding these influences may provide a mechanism for measuring the dynamics of breast cancer risk. The purpose of this work is to provide support for an automated serial mammography study under way at the authors' institution, where the digital mammographic images are acquired with a full-field digital mammography imaging system. At the initiation of the serial study, it was clear that the authors did not fully understand the nature of the problem: automatically comparing similar mammographic scenes acquired at different times. The evidence indicates that there are many factors that influence breast tissue at any given time and thus have the ability to alter the associated radiographic appearance over time. In general, the topics considered herein include aging; involution; breast development; exogenous and endogenous hormonal interactions such as hormone replacement therapy, oral contraceptive use, and menstrual timing; screening sensitivity issues and interval cancers; tumor growth rates; sojourn times; and lifestyle factors such as diet and exercise. Throughout this work, commentaries and suggestive strategies for automated serial image analysis are provided.

A statistical methodology for mammographic density detection

A statistical methodology is presented based on a chi-square probability analysis that allows the automated discrimination of radiolucent tissue (fat) from radiographic densities (fibroglandular tissue) in digitized mammograms. The method is based on earlier work developed at this facility that shows mammograms may be considered as evolving from a linear filtering operation where a random input field is passed through a 1/f filtering process. The filtering process is reversible which allows the solution of the input field with knowledge obtained from the raw image (the output). The input field solution is analogous to a prewhitening technique or deconvolution. This field contains all the information of the raw image in a much simplified format that can be approximated and analyzed with parametric methods. In the work presented here evidence indicates that there are two random events occurring in the input field with differing variances: (1) one relating to fat tissue with the smaller variance, and (2) the second relating to all other tissue with the larger variance. A statistical comparison of the variances is made by scanning the image with a small search window. A relaxation method allows for making a reliable estimate of the smaller variance which is considered as the global reference. If a local variance deviates significantly from the reference variance, based on chi-square analysis, it is labeled as nonfat; otherwise it is labeled as fat. This statistical test procedure results in a region by region continuous labeling of fat and nonfat tissue across the image. In the work presented here, the emphasis is on the methodology development with supporting preliminary results that are very encouraging. It is widely accepted that mammographic density is a breast cancer risk factor. An important application of this work is to incorporate density-based risk analysis into the ongoing statistical-based detection work developed at this facility. Additional applications include risk analysis dependent on either percentages or total amounts of fat or dense tissue. This work may be considered as the initial step in introducing many of the known breast cancer risk factors into the actual image data analysis.

On the statistical nature of mammograms

We show that digitized mammograms can be considered as evolving from a simple process. A given image results from passing a random input field through a linear filtering operation, where the filter transfer function has a self-similar characteristic. By estimating the functional form of the filter and solving the corresponding filtering equation, the analysis shows that the input field gray value distribution and spectral content can be approximated with parametric methods. The work gives a simple explanation for the variegated image appearance and multimodal character of the gray value distribution common to mammograms. Using the image analysis as a guide, a simulated mammogram is generated that has many statistical characteristics of real mammograms. Additional benefits may follow from understanding the functional form of the filter in conjunction with the input field characteristics that include the approximate parametric description of mammograms, showing the distinction between homogeneously dense and nondense images, and the development of mass analysis methods.

Review and evaluation of MRI nonuniformity corrections for brain tumor response measurements

Current MRI nonuniformity correction techniques are reviewed and investigated. Many approaches are used to remedy this artifact, but it is not clear which method is the most appropriate in a given situation, as the applications have been with different MRI coils and different clinical applications. In this work four widely used nonuniformity correction techniques are investigated in order to assess the effect on tumor response measurements (change in tumor volume over time): a phantom correction method, an image smoothing technique, homomorphic filtering, and surface fitting approach. Six brain tumor cases with baseline and follow-up MRIs after treatment with varying degrees of difficulty of segmentation were analyzed without and with each of the nonuniformity corrections. Different methods give significantly different correction images, indicating that rf nonuniformity correction is not yet well understood. No improvement in tumor segmentation or in tumor growth/shrinkage assessment was achieved using any of the evaluated corrections.

Multiresolution statistical analysis of high-resolution digital mammograms

A multiresolution statistical method for identifying clinically normal tissue in digitized mammograms is used to construct an algorithm for separating normal regions from potentially abnormal regions; that is, small regions that may contain isolated calcifications. This is the initial phase of the development of a general method for the automatic recognition of normal mammograms. The first step is to decompose the image with a wavelet expansion that yields a sum of independent images, each containing different levels of image detail. When calcifications are present, there is strong empirical evidence that only some of the image components are necessary for the purpose of detecting a deviation from normal. The underlying statistic for each of the selected expansion components can be modeled with a simple parametric probability distribution function. This function serves as an instrument for the development of a statistical test that allows for the recognition of normal tissue regions. The distribution function depends on only one parameter, and this parameter itself has an underlying statistical distribution. The values of this parameter define a summary statistic that can be used to set detection error rates. Once the summary statistic is determined, spatial filters that are matched to resolution are applied independently to each selected expansion image. Regions of the image that correlate with the normal statistical model are discarded and regions in disagreement (suspicious areas) are flagged. These results are combined to produce a detection output image consisting only of suspicious areas. This type of detection output is amenable to further processing that may ultimately lead to a fully automated algorithm for the identification of normal mammograms.

MRI segmentation: methods and applications

The current literature on MRI segmentation methods is reviewed. Particular emphasis is placed on the relative merits of single image versus multispectral segmentation, and supervised versus unsupervised segmentation methods. Image pre-processing and registration are discussed, as well as methods of validation. The application of MRI segmentation for tumor volume measurements during the course of therapy is presented here as an example, illustrating problems associated with inter- and intra-observer variations inherent to supervised methods.