Full field digital mammography and breast density: comparison of calibrated and noncalibrated measurements

Reproductive Factors Related to Childbearing and a Novel Automated Mammographic Measure, V

BACKGROUND

We investigated the associations between several reproductive factors related to childbearing and the variation (V) measure (a novel, objective, single summary measure of breast image intensity) by menopausal status.

METHODS

Our study included 3,814 cancer-free women within the Nurses' Health Study (NHS) and NHSII cohorts. The data on reproductive variables and covariates were obtained from biennial questionnaires closest to the mammogram date. V-measures were obtained from mammographic images using a previously developed algorithm capturing the standard deviation of pixel values. We used multivariate linear regression to examine the associations of parity, age at first birth, time between menarche and first birth, time since last pregnancy, and lifetime breastfeeding duration with V-measure, adjusting for breast cancer risk factors, including the percentage of mammographic density (PMD). We further examined whether these associations were statistically accounted for (mediated) by PMD.

RESULTS

Among premenopausal women, none of the reproductive factors were associated with V. Among postmenopausal women, inverse associations of parity and positive associations of age at first birth with V were mediated by PMD (percent mediated: nulliparity: 66.7%, P < 0.0001; parity: 50.5%, P < 0.01; age at first birth 76.1%, P < 0.001) and were no longer significant in PMD-adjusted models. Lifetime duration of breastfeeding was positively associated with V [>36 vs. 0 ≤1 months β = 0.29; 95% confidence interval (CI) 0.07; 0.52, Ptrend < 0.01], independent of PMD.

CONCLUSIONS

Parity, age at first birth, and breastfeeding were associated with postmenopausal V.

IMPACT

This study highlights associations of reproductive factors with mammographic image intensity.

The Role of the Biosafety Cabinet in Preventing Infection in the Clinical Laboratory

Clinical laboratories are essential in healthcare to better diagnose, treat, and track medical diseases. However, handling infectious organisms and possibly infectious materials in these laboratories puts the safety of laboratory workers and the general public at risk. By controlling the distribution of infectious substances and stopping the spread of diseases, biosafety cabinets (BSCs) have become crucial tools in guaranteeing laboratory safety. The prevention of infections is most important in medical and laboratory settings. In clinical laboratories, biological and infectious agents are handled, posing threats to healthcare workers and the general public. To avoid infections, proper training of the BSC is essential. Laboratory employees are instructed in aseptic procedures, proper hand posture, and efficient personal protection when working in the cabinet. These instructions decrease the chance of contaminating the surrounding area. Additionally, user ergonomics are taken into account while designing BSC, reducing operator fatigue, and guaranteeing that staff can execute tasks precisely for extended periods. This review highlights the importance of biosafety cabinets in maintaining a secure laboratory environment and explains their crucial function in infection control.

Breast Density Analysis Using Digital Breast Tomosynthesis

We evaluated an automated percentage of breast density (BD) technique (PDa) with digital breast tomosynthesis (DBT) data. The approach is based on the wavelet expansion followed by analyzing signal dependent noise. Several measures were investigated as risk factors: normalized volumetric; total dense volume; average of the DBT slices (slice-mean); a two-dimensional (2D) metric applied to the synthetic images; and the mean and standard deviations of the pixel values. Volumetric measures were derived theoretically, and PDa was modeled as a function of compressed breast thickness. An alternative method for constructing synthetic 2D mammograms was investigated using the volume results. A matched case-control study (n = 426 pairs) was analyzed. Conditional logistic regression modeling, controlling body mass index and ethnicity, was used to estimate odds ratios (ORs) for each measure with 95% confidence intervals provided parenthetically. There were several significant findings: volumetric measure [OR = 1.43 (1.18, 1.72)], which produced an identical OR as the slice-mean measure as predicted; [OR =1.44 (1.18, 1.75)] when applied to the synthetic images; and mean of the pixel values (volume or 2D synthetic) [ORs ~ 1.31 (1.09, 1.57)]. PDa was modeled as 2nd degree polynomial (concave-down): its maximum value occurred at 0.41×(compressed breast thickness), which was similar across case-control groups, and was significant from this position [OR = 1.47 (1.21, 1.78)]. A standardized 2D synthetic image was produced, where each pixel value represents the percentage of BD above its location. The significant findings indicate the validity of the technique. Derivations supported by empirical analyses produced a new synthetic 2D standardized image technique. Ancillary to the objectives, the results provide evidence for understanding the percentage of BD measure applied to 2D mammograms. Notwithstanding the findings, the study design provides a template for investigating other measures such as texture.

Associations of Oral Contraceptives with Mammographic Breast Density in Premenopausal Women

BACKGROUND

We investigated the associations of oral contraceptives (OC) with percent breast density (PD), absolute dense area (DA), nondense area (NDA), and a novel image intensity variation (V) measure in premenopausal women.

METHODS

This study included 1,233 controls from a nested case-control study within Nurses' Health Study II cohort. Information on OCs was collected in 1989 and updated biennially. OC use was defined from the questionnaire closest to the mammogram date. PD, DA, and NDA were measured from digitized film mammograms using a computer-assisted thresholding technique; the V measure was obtained with a previously developed algorithm measuring the SD of pixel values in the eroded breast region. Generalized linear regression was used to assess associations between OCs and density measures (square root-transformed PD, DA, and NDA, and -untransformed V).

RESULTS

OC use was not associated with PD [current vs. never: β = -0.06; 95% confidence interval (CI), -0.37-0.24; past vs. never: β = 0.10; 95% CI, -0.09-0.29], DA (current vs. never: β = -0.20; 95% CI -0.59-0.18; past vs. never: β = 0.13; 95% CI, -0.12-0.39), and NDA (current vs. never: β = -0.19; 95% CI, -0.56-0.18; past vs. never: β = -0.01; 95% CI, -0.28-0.25). Women with younger age at initiation had significantly greater V-measure (<20 years vs. never: β = 26.88; 95% CI, 3.18-50.58; 20-24 years vs. never: β = 20.23; 95% CI, -4.24-44.71; 25-29 years vs. never: β = 2.61; 95% CI -29.00-34.23; ≥30 years vs. never: β = 0.28; 95% CI, -34.16-34.72, P _trend = 0.03).

CONCLUSIONS

Our findings suggest that an earlier age at first OC use was associated with significantly greater V.

IMPACT

These findings could guide decisions about the age for OC initiation.

Mammographic Variation Measures, Breast Density, and Breast Cancer Risk

OBJECTIVE. Our previous work showed that variation measures, which represent breast architecture derived from mammograms, were significantly associated with breast cancer. For replication purposes, we examined the association of three variation measures (variation [V], which is measured in the image domain, and P₁ and p₁ [a normalized version of P₁], which are derived from restricted regions in the Fourier domain) with breast cancer risk in an independent population. We also compared these measures to volumetric density measures (volumetric percent density [VPD] and dense volume [DV]) from a commercial product. MATERIALS AND METHODS. We examined 514 patients with breast cancer and 1377 control patients from a screening practice who were matched for age, date of examination, mammography unit, facility, and state of residence. Spearman rank-order correlation was used to evaluate the monotonic association between measures. Breast cancer associations were estimated using conditional logistic regression, after adjustment for age and body mass index. Odds ratios were calculated per SD increment in mammographic measure. RESULTS. These variation measures were strongly correlated with VPD (correlation, 0.68-0.80) but not with DV (correlation, 0.31-0.48). Similar to previous findings, all variation measures were significantly associated with breast cancer (odds ratio per SD: 1.30 [95% CI, 1.16-1.46] for V, 1.55 [95% CI, 1.35-1.77] for P₁, and 1.51 [95% CI, 1.33-1.72] for p₁). Associations of volumetric density measures with breast cancer were similar (odds ratio per SD: 1.54 [95% CI, 1.33-1.78] for VPD and 1.34 [95% CI, 1.20-1.50] for DV). When DV was included with each variation measure in the same model, all measures retained significance. CONCLUSION. Variation measures were significantly associated with breast cancer risk (comparable to the volumetric density measures) but were independent of the DV.

Empirically-derived synthetic populations to mitigate small sample sizes

Limited sample sizes can lead to spurious modeling findings in biomedical research. The objective of this work is to present a new method to generate synthetic populations (SPs) from limited samples using matched case-control data (n = 180 pairs), considered as two separate limited samples. SPs were generated with multivariate kernel density estimations (KDEs) with unconstrained bandwidth matrices. We included four continuous variables and one categorical variable for each individual. Bandwidth matrices were determined with Differential Evolution (DE) optimization by covariance comparisons. Four synthetic samples (n = 180) were derived from their respective SPs. Similarity between observed samples with synthetic samples was compared assuming their empirical probability density functions (EPDFs) were similar. EPDFs were compared with the maximum mean discrepancy (MMD) test statistic based on the Kernel Two-Sample Test. To evaluate similarity within a modeling context, EPDFs derived from the Principal Component Analysis (PCA) scores and residuals were summarized with the distance to the model in X-space (DModX) as additional comparisons. Four SPs were generated from each sample. The probability of selecting a replicate when randomly constructing synthetic samples (n = 180) was infinitesimally small. MMD tests indicated that the observed sample EPDFs were similar to the respective synthetic EPDFs. For the samples, PCA scores and residuals did not deviate significantly when compared with their respective synthetic samples. The feasibility of this approach was demonstrated by producing synthetic data at the individual level, statistically similar to the observed samples. The methodology coupled KDE with DE optimization and deployed novel similarity metrics derived from PCA. This approach could be used to generate larger-sized synthetic samples. To develop this approach into a research tool for data exploration purposes, additional evaluation with increased dimensionality is required. Moreover, given a fully specified population, the degree to which individuals can be discarded while synthesizing the respective population accurately will be investigated. When these objectives are addressed, comparisons with other techniques such as bootstrapping will be required for a complete evaluation.

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.

Calibrated Breast Density Measurements

RATIONALE AND OBJECTIVES

Mammographic density is an important risk factor for breast cancer, but translation to the clinic requires assurance that prior work based on mammography is applicable to current technologies. The purpose of this work is to evaluate whether a calibration methodology developed previously produces breast density metrics predictive of breast cancer risk when applied to a case-control study.

MATERIALS AND METHODS

A matched case control study (n = 319 pairs) was used to evaluate two calibrated measures of breast density. Two-dimensional mammograms were acquired from six Hologic mammography units: three conventional Selenia two-dimensional full-field digital mammography systems and three Dimensions digital breast tomosynthesis systems. We evaluated the capability of two calibrated breast density measures to quantify breast cancer risk: the mean (PGm) and standard deviation (PGsd) of the calibrated pixels. Matching variables included age, hormone replacement therapy usage/duration, screening history, and mammography unit. Calibrated measures were compared to the percentage of breast density (PD) determined with the operator-assisted Cumulus method. Conditional logistic regression was used to generate odds ratios (ORs) from continuous and quartile (Q) models with 95% confidence intervals. The area under the receiver operating characteristic curve (Az) was also used as a comparison metric. Both univariate models and models adjusted for body mass index and ethnicity were evaluated.

RESULTS

In adjusted models, both PGsd and PD were statistically significantly associated with breast cancer with similar Az of 0.61-0.62. The corresponding ORs and confidence intervals were also similar. For PGsd, the OR was 1.34 (1.09, 1.66) for the continuous measure and 1.83 (1.11, 3.02), 2.19 (1.28, 3.73), and 2.20 (1.26, 3.85) for Q2-Q4. For PD, the OR was 1.43 (1.16, 1.76) for the continuous measure and 0.84 (0.52, 1.38), 1.96 (1.19, 3.23), and 2.27 (1.29, 4.00) for Q2-Q4. The results for PGm were slightly attenuated and not statistically significant. The OR was 1.22 (0.99, 1.51) with Az = 0.60 for the continuous measure and 1.24 (0.78, 1.97), 0.98 (0.60, 1.61), and 1.26, (0.77, 2.07) for Q2-Q4 with Az = 0.60.

CONCLUSION

The calibrated PGsd measure provided significant associations with breast cancer comparable to those given by PD. The calibrated PGm performed slightly worse. These findings indicate that the calibration approach developed previously replicates under more general conditions.

Spatial Correlation and Breast Cancer Risk

We present a novel method for evaluating local spatial correlation structure in two-dimensional (2D) mammograms and evaluate its capability for risk prediction as one possible application. Two matched case-control studies were analyzed. Study 1 included women (N = 588 pairs) with mammograms acquired with either Hologic Selenia full field digital mammography (FFDM) units or Hologic Dimensions digital breast tomosynthesis units. Study 2 included women (N =180 pairs) with mammograms acquired with a General Electric Senographe 2000D FFDM unit. Matching variables included age, HRT usage/duration, screening history, and mammography unit. Local autocorrelation functions were determined with Fourier analysis and compared with a template defined as a 2D double-sided exponential function with one spatial extent parameter: n = 4, 12, 24, 50, 74, 100, and 124, where (n+1)×(n+1) is the area of the local spatial extent measured in pixels. The difference between the local correlation and template was gauged within an adjustable parameter kernel and summarized, producing two measures: the mean (m_n+1)_, and standard deviation (s_n+1). Both adjustable parameters were varied in Study 1. Select measures that produced significant associations with breast cancer were translated to Study 2. Breast cancer associations were evaluated with conditional logistic regression, adjusted for body mass index and ethnicity. Odds ratios (ORs) were estimated as per standard deviation increment with 95% confidence intervals (CIs). Two measures were selected for breast cancer association analysis in Study 1: m₇₅ and s₂₅. Both measures revealed significant associations with breast cancer: OR = 1.45 (1.23, 1.66) for m₇₅ and OR = 1.30 (1.14, 1.49) for s_25. When translating to Study 2, these measures also revealed significant associations: OR = 1.49 (1.12, 1.96) for m₇₅ and OR = 1.34 (1.06, 1.69) for s₂₅. Novel correlation metrics presented in this work produced significant associations with breast cancer risk. This approach is general and may have applications beyond mammography.

Technical challenges in generalizing calibration techniques for breast density measurements

PURPOSE

We are developing a calibration methodology for full-field digital mammography (FFDM). Calibration compensates for image acquisition technique influences on the pixel representation, ideally producing improved inter-image breast density estimates. This approach relies on establishing references with rigid breast tissue-equivalent phantoms (BTEs) and requires an accurate estimate of the compressed breast thickness because the system readout is nominal. There is also an attenuation mismatch between adipose breast tissue and the adipose BTE that was noted in our previous work. It is referred to as the "attenuation anomaly" and addressed in this report. The objectives are to evaluate methods to correct for the compressed breast thickness and compensate for the attenuation anomaly.

METHODS

Thickness correction surfaces were established with a deformable phantom (DP) using both image and physical measurements for three direct x-ray conversion FFDM units. The Cumulative Sum serial quality control procedure was established to ensure the thickness correction measurements were stable over time by imaging and calibrating DPs biweekly in lieu of physical measurements. The attenuation anomaly was addressed by evaluating adipose image regions coupled with an optimization technique to adjust the adipose calibration data. We compared calibration consistency across matched left and right cranial caudal (CC) mammographic views (n = 199) with and without corrections using Bland-Altman plots. These plots were complemented by comparing the right and left breast calibrated average (μ_a ) and population distribution mean (m_a ) with 95% confidence intervals and difference distribution variances with the F-test for uncorrected and corrected data.

RESULTS

Thickness correction surfaces were well approximated as tilted planes and were dependent upon compression force. A correction was developed for the attenuation anomaly. All paddles (large and small paddles for all units) exhibited similar tilt as a function of force. Without correction, m_a  = 0.92 (-1.77, 3.62) was not significantly different from zero with many negative μ_a samples. The thickness correction produced a significant shift in the μ_a distribution in the positive direction with m_a  = 13.99 (11.17, 16.80) and reduced the difference distribution variance significantly (P < 0.0001). Applying both corrections in tandem gave m_a  = 22.83 (20.32, 25.34), representing another significant positive shift in comparison with the thickness correction in isolation. Thickness corrections were stable over approximately a 2-year timeframe for all units.

CONCLUSION

These correction techniques are valid approaches for addressing technical problems with calibration that relies on reference phantoms. The efficacy of the calibration methodology will require validation with clinical endpoints in future studies.

Generalized breast density metrics

Mammograms represent data that can inform future risk of breast cancer. Data from two case-control study populations were analyzed. Population 1 included women (N  =  180 age matched case-control pairs) with mammograms acquired with one indirect x-ray conversion mammography unit. Population 2 included women (N  =  319 age matched case-control pairs) with mammograms acquired from 6 direct x-ray conversion units. The Fourier domain was decomposed into n concentric rings (radial spatial frequency bands). The power in each ring was summarized giving a set of measures. We investigated images in raw, for presentation (processed) and calibrated representations and made comparison with the percentage of breast density (BD) determined with the operator assisted Cumulus method. Breast cancer associations were evaluated with conditional logistic regression, adjusted for body mass index and ethnicity. Odds ratios (ORs), per standard deviation increase derived from the respective breast density distributions and 95% confidence intervals (CIs) were estimated. A measure from a lower radial frequency ring, corresponding 0.083-0.166 cycles mm^-1 and BD had significant associations with risk in both populations. In Population 1, the Fourier measure produced significant associations in each representation: OR  =  1.76 (1.33, 2.32) for raw; OR  =  1.43 (1.09, 1.87) for processed; and OR  =  1.68 (1.26, 2.25) for calibrated. BD also provided significant associations in Population 1: OR  =  1.72 (1.27, 2.33). In Population 2, the Fourier measure produced significant associations for each representation as well: OR  =  1.47 (1.19, 1.80) for raw; OR  =  1.38 (1.15, 1.67) for processed; and OR  =  1.42 (1.15, 1.75) for calibrated. BD provided significant associations in Population 2: OR  =  1.43 (1.17, 1.76). Other coincident spectral regions were also predictive of case-control status. In sum, generalized breast density measures were significantly associated with breast cancer in both FFDM technologies.

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.

Segmentation of breast masses on dedicated breast computed tomography and three-dimensional breast ultrasound images

We present and evaluate a method for the three-dimensional (3-D) segmentation of breast masses on dedicated breast computed tomography (bCT) and automated 3-D breast ultrasound images. The segmentation method, refined from our previous segmentation method for masses on contrast-enhanced bCT, includes two steps: (1) initial contour estimation and (2) active contour-based segmentation to further evolve and refine the initial contour by adding a local energy term to the level-set equation. Segmentation performance was assessed in terms of Dice coefficients (DICE) for 129 lesions on noncontrast bCT, 38 lesions on contrast-enhanced bCT, and 98 lesions on 3-D breast ultrasound (US) images. For bCT, DICE values of 0.82 and 0.80 were obtained on contrast-enhanced and noncontrast images, respectively. The improvement in segmentation performance with respect to that of our previous method was statistically significant ( p = 0.002 ). Moreover, segmentation appeared robust with respect to the presence of glandular tissue. For 3-D breast US, the DICE value was 0.71. Hence, our method obtained promising results for both 3-D imaging modalities, laying a solid foundation for further quantitative image analysis and potential future expansion to other 3-D imaging modalities.

A comparison of calibration data from full field digital mammography units for breast density measurements

Background

Breast density is a significant breast cancer risk factor measured from mammograms. The most appropriate method for measuring breast density for risk applications is still under investigation. Calibration standardizes mammograms to account for acquisition technique differences prior to making breast density measurements. We evaluated whether a calibration methodology developed for an indirect x-ray conversion full field digital mammography (FFDM) technology applies to direct x-ray conversion FFDM systems.

Methods

Breast tissue equivalent (BTE) phantom images were used to establish calibration datasets for three similar direct x-ray conversion FFDM systems. The calibration dataset for each unit is a function of the target/filter combination, x-ray tube voltage, current × time (mAs), phantom height, and two detector fields of view (FOVs). Methods were investigated to reduce the amount of calibration data by restricting the height, mAs, and FOV sampling. Calibration accuracy was evaluated with mixture phantoms. We also compared both intra- and inter-system calibration characteristics and accuracy.

Results

Calibration methods developed previously apply to direct x-ray conversion systems with modification. Calibration accuracy was largely within the acceptable range of ± 4 standardized units from the ideal value over the entire acquisition parameter space for the direct conversion units. Acceptable calibration accuracy was maintained with a cubic-spline height interpolation, representing a modification to previous work. Calibration data is unit specific, can be acquired with the large FOV, and requires a minimum of one reference mAs sample. The mAs sampling, calibration accuracy, and the necessity for machine specific calibration data are common characteristics and in agreement with our previous work.

Conclusion

The generality of our calibration approach was established under ideal conditions. Evaluation with patient data using breast cancer status as the endpoint is required to demonstrate that the approach produces a breast density measure associated with breast cancer.

Beyond mammography: new frontiers in breast cancer screening

Breast cancer screening remains a subject of intense and, at times, passionate debate. Mammography has long been the mainstay of breast cancer detection and is the only screening test proven to reduce mortality. Although it remains the gold standard of breast cancer screening, there is increasing awareness of subpopulations of women for whom mammography has reduced sensitivity. Mammography also has undergone increased scrutiny for false positives and excessive biopsies, which increase radiation dose, cost, and patient anxiety. In response to these challenges, new technologies for breast cancer screening have been developed, including low-dose mammography, contrast-enhanced mammography, tomosynthesis, automated whole breast ultrasound, molecular imaging, and magnetic resonance imaging. Here we examine some of the current controversies and promising new technologies that may improve detection of breast cancer both in the general population and in high-risk groups, such as women with dense breasts. We propose that optimal breast cancer screening will ultimately require a personalized approach based on metrics of cancer risk with selective application of specific screening technologies best suited to the individual's age, risk, and breast density.

Reported mammographic density: film-screen versus digital acquisition

PURPOSE

To test the hypothesis that American College of Radiology Breast Imaging Reporting and Data System (BI-RADS) categories for breast density reported by radiologists are lower when digital mammography is used than those reported when film-screen (FS) mammography is used.

MATERIALS AND METHODS

This study was institutional review board approved and HIPAA compliant. Demographic data, risk factors, and BI-RADS breast density categories were collected from five mammography registries that were part of the Breast Cancer Surveillance Consortium. Active, passive, or waiver of consent was obtained for all participants. Women aged 40 years and older who underwent at least two screening mammographic examinations less than 36 months apart between January 1, 2000, and December 31, 2009, were included. Women with prior breast cancer, augmentation, or use of agents known to affect density were excluded. The main sample included 89 639 women with both FS and digital mammograms. The comparison group included 259 046 women with two FS mammograms and 87 066 women with two digital mammograms. BI-RADS density was cross-tabulated according to the order in which the two types of mammogram were acquired and by the first versus second interpretation.

RESULTS

Regardless of acquisition method, the percentage of women with a change in density from one reading to the next was similar. Breast density was lower in 19.8% of the women who underwent FS before digital mammography and 17.1% of those who underwent digital before FS mammography. Similarly, lower density classifications were reported on the basis of the second mammographic examination regardless of acquisition method (15.8%-19.8%). The percentage of agreement between density readings was similar regardless of mammographic types paired (67.3%-71.0%).

CONCLUSION

The study results showed no difference in reported BI-RADS breast density categories according to acquisition method. Reported BI-RADS density categories may be useful in the development of breast cancer risk models in which FS, digital, or both acquisition methods are used.

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.

Estimation of breast percent density in raw and processed full field digital mammography images via adaptive fuzzy c-means clustering and support vector machine segmentation

PURPOSE

The amount of fibroglandular tissue content in the breast as estimated mammographically, commonly referred to as breast percent density (PD%), is one of the most significant risk factors for developing breast cancer. Approaches to quantify breast density commonly focus on either semiautomated methods or visual assessment, both of which are highly subjective. Furthermore, most studies published to date investigating computer-aided assessment of breast PD% have been performed using digitized screen-film mammograms, while digital mammography is increasingly replacing screen-film mammography in breast cancer screening protocols. Digital mammography imaging generates two types of images for analysis, raw (i.e., "FOR PROCESSING") and vendor postprocessed (i.e., "FOR PRESENTATION"), of which postprocessed images are commonly used in clinical practice. Development of an algorithm which effectively estimates breast PD% in both raw and postprocessed digital mammography images would be beneficial in terms of direct clinical application and retrospective analysis.

METHODS

This work proposes a new algorithm for fully automated quantification of breast PD% based on adaptive multiclass fuzzy c-means (FCM) clustering and support vector machine (SVM) classification, optimized for the imaging characteristics of both raw and processed digital mammography images as well as for individual patient and image characteristics. Our algorithm first delineates the breast region within the mammogram via an automated thresholding scheme to identify background air followed by a straight line Hough transform to extract the pectoral muscle region. The algorithm then applies adaptive FCM clustering based on an optimal number of clusters derived from image properties of the specific mammogram to subdivide the breast into regions of similar gray-level intensity. Finally, a SVM classifier is trained to identify which clusters within the breast tissue are likely fibroglandular, which are then aggregated into a final dense tissue segmentation that is used to compute breast PD%. Our method is validated on a group of 81 women for whom bilateral, mediolateral oblique, raw and processed screening digital mammograms were available, and agreement is assessed with both continuous and categorical density estimates made by a trained breast-imaging radiologist.

RESULTS

Strong association between algorithm-estimated and radiologist-provided breast PD% was detected for both raw (r = 0.82, p < 0.001) and processed (r = 0.85, p < 0.001) digital mammograms on a per-breast basis. Stronger agreement was found when overall breast density was assessed on a per-woman basis for both raw (r = 0.85, p < 0.001) and processed (0.89, p < 0.001) mammograms. Strong agreement between categorical density estimates was also seen (weighted Cohen's κ ≥ 0.79). Repeated measures analysis of variance demonstrated no statistically significant differences between the PD% estimates (p > 0.1) due to either presentation of the image (raw vs processed) or method of PD% assessment (radiologist vs algorithm).

CONCLUSIONS

The proposed fully automated algorithm was successful in estimating breast percent density from both raw and processed digital mammographic images. Accurate assessment of a woman's breast density is critical in order for the estimate to be incorporated into risk assessment models. These results show promise for the clinical application of the algorithm in quantifying breast density in a repeatable manner, both at time of imaging as well as in retrospective studies.

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.