Calibrated measures for breast density estimation

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.

A genome-wide association study of mammographic texture variation

Background

Breast parenchymal texture features, including grayscale variation (V), capture the patterns of texture variation on a mammogram and are associated with breast cancer risk, independent of mammographic density (MD). However, our knowledge on the genetic basis of these texture features is limited.

Methods

We conducted a genome-wide association study of V in 7040 European-ancestry women. V assessments were generated from digitized film mammograms. We used linear regression to test the single-nucleotide polymorphism (SNP)-phenotype associations adjusting for age, body mass index (BMI), MD phenotypes, and the top four genetic principal components. We further calculated genetic correlations and performed SNP-set tests of V with MD, breast cancer risk, and other breast cancer risk factors.

Results

We identified three genome-wide significant loci associated with V: rs138141444 (6q24.1) in ECT2L, rs79670367 (8q24.22) in LINC01591, and rs113174754 (12q22) near PGAM1P5. 6q24.1 and 8q24.22 have not previously been associated with MD phenotypes or breast cancer risk, while 12q22 is a known locus for both MD and breast cancer risk. Among known MD and breast cancer risk SNPs, we identified four variants that were associated with V at the Bonferroni-corrected thresholds accounting for the number of SNPs tested: rs335189 (5q23.2) in PRDM6, rs13256025 (8p21.2) in EBF2, rs11836164 (12p12.1) near SSPN, and rs17817449 (16q12.2) in FTO. We observed significant genetic correlations between V and mammographic dense area (r_g = 0.79, P = 5.91 × 10⁻⁵), percent density (r_g = 0.73, P = 1.00 × 10⁻⁴), and adult BMI (r_g =  − 0.36, P = 3.88 × 10⁻⁷). Additional significant relationships were observed for non-dense area (z =  − 4.14, P = 3.42 × 10⁻⁵), estrogen receptor-positive breast cancer (z = 3.41, P = 6.41 × 10⁻⁴), and childhood body fatness (z =  − 4.91, P = 9.05 × 10⁻⁷) from the SNP-set tests.

Conclusions

These findings provide new insights into the genetic basis of mammographic texture variation and their associations with MD, breast cancer risk, and other breast cancer risk factors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-022-01570-8.

Automated percent mammographic density, mammographic texture variation, and risk of breast cancer: a nested case-control study

Percent mammographic density (PMD) is a strong breast cancer risk factor, however, other mammographic features, such as V, the standard deviation (SD) of pixel intensity, may be associated with risk. We assessed whether PMD, automated PMD (APD), and V, yielded independent associations with breast cancer risk. We included 1900 breast cancer cases and 3921 matched controls from the Nurses' Health Study (NHS) and the NHSII. Using digitized film mammograms, we estimated PMD using a computer-assisted thresholding technique. APD and V were determined using an automated computer algorithm. We used logistic regression to generate odds ratios (ORs) and 95% confidence intervals (CIs). Median time from mammogram to diagnosis was 4.1 years (interquartile range: 1.6-6.8 years). PMD (OR per SD:1.52, 95% CI: 1.42, 1.63), APD (OR per SD:1.32, 95% CI: 1.24, 1.41), and V (OR per SD:1.32, 95% CI: 1.24, 1.40) were positively associated with breast cancer risk. Associations for APD were attenuated but remained statistically significant after mutual adjustment for PMD or V. Women in the highest quartile of both APD and V (OR vs Q1/Q1: 2.49, 95% CI: 2.02, 3.06), or PMD and V (OR vs Q1/Q1: 3.57, 95% CI: 2.79, 4.58) had increased breast cancer risk. An automated method of PMD assessment is feasible and yields similar, but somewhat weaker, estimates to a manual measure. PMD, APD and V are each independently, positively associated with breast cancer risk. Women with dense breasts and greater texture variation are at the highest relative risk of breast cancer.

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.

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.

Volumetric breast density measurement for personalized screening: accuracy, reproducibility, consistency, and agreement with visual assessment

Assessment of breast density at the point of mammographic examination could lead to optimized breast cancer screening pathways. The onsite breast density information may offer guidance of when to recommend supplemental imaging for women in a screening program. A software application (Insight BD, Siemens Healthcare GmbH) for fast onsite quantification of volumetric breast density is evaluated. The accuracy of the method is assessed using breast tissue equivalent phantom experiments resulting in a mean absolute error of 3.84%. Reproducibility of measurement results is analyzed using 8427 exams in total, comparing for each exam (if available) the densities determined from left and right views, from cranio-caudal and medio-lateral oblique views, from full-field digital mammograms (FFDM) and digital breast tomosynthesis (DBT) data and from two subsequent exams of the same breast. Pearson correlation coefficients of 0.937, 0.926, 0.950, and 0.995 are obtained. Consistency of the results is demonstrated by evaluating the dependency of the breast density on women's age. Furthermore, the agreement between breast density categories computed by the software with those determined visually by 32 radiologists is shown by an overall percentage agreement of 69.5% for FFDM and by 64.6% for DBT data. These results demonstrate that the software delivers accurate, reproducible, and consistent measurements that agree well with the visual assessment of breast density by radiologists.

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.

Mammographic Breast Density in Chinese Women: Spatial Distribution and Autocorrelation Patterns

Mammographic breast density (MBD) is a strong risk factor for breast cancer. The spatial distribution of MBD in the breast is variable and dependent on physiological, genetic, environmental and pathological factors. This pilot study aims to define the spatial distribution and autocorrelation patterns of MBD in Chinese women aged 40–60. By analyzing their digital mammographic images using a public domain Java image processing program for segmentation and quantification of MBD, we found their left and right breasts were symmetric to each other in regard to their breast size (Total Breast Area), the amount of BMD (overall PD) and Moran's I values. Their MBD was also spatially autocorrelated together in the anterior part of the breast in those with a smaller breast size, while those with a larger breast size tend to have their MBD clustered near the posterior part of the breast. Finally, we observed that the autocorrelation pattern of MBD was dispersed after a 3-year observation period.

Preliminary evaluation of the publicly available Laboratory for Breast Radiodensity Assessment (LIBRA) software tool: comparison of fully automated area and volumetric density measures in a case-control study with digital mammography

INTRODUCTION

Breast density, commonly quantified as the percentage of mammographically dense tissue area, is a strong breast cancer risk factor. We investigated associations between breast cancer and fully automated measures of breast density made by a new publicly available software tool, the Laboratory for Individualized Breast Radiodensity Assessment (LIBRA).

METHODS

Digital mammograms from 106 invasive breast cancer cases and 318 age-matched controls were retrospectively analyzed. Density estimates acquired by LIBRA were compared with commercially available software and standard Breast Imaging-Reporting and Data System (BI-RADS) density estimates. Associations between the different density measures and breast cancer were evaluated by using logistic regression after adjustment for Gail risk factors and body mass index (BMI). Area under the curve (AUC) of the receiver operating characteristic (ROC) was used to assess discriminatory capacity, and odds ratios (ORs) for each density measure are provided.

RESULTS

All automated density measures had a significant association with breast cancer (OR = 1.47-2.23, AUC = 0.59-0.71, P < 0.01) which was strengthened after adjustment for Gail risk factors and BMI (OR = 1.96-2.64, AUC = 0.82-0.85, P < 0.001). In multivariable analysis, absolute dense area (OR = 1.84, P < 0.001) and absolute dense volume (OR = 1.67, P = 0.003) were jointly associated with breast cancer (AUC = 0.77, P < 0.01), having a larger discriminatory capacity than models considering the Gail risk factors alone (AUC = 0.64, P < 0.001) or the Gail risk factors plus standard area percent density (AUC = 0.68, P = 0.01). After BMI was further adjusted for, absolute dense area retained significance (OR = 2.18, P < 0.001) and volume percent density approached significance (OR = 1.47, P = 0.06). This combined area-volume density model also had a significantly (P < 0.001) improved discriminatory capacity (AUC = 0.86) relative to a model considering the Gail risk factors plus BMI (AUC = 0.80).

CONCLUSIONS

Our study suggests that new automated density measures may ultimately augment the current standard breast cancer risk factors. In addition, the ability to fully automate density estimation with digital mammography, particularly through the use of publically available breast density estimation software, could accelerate the translation of density reporting in routine breast cancer screening and surveillance protocols and facilitate broader research into the use of breast density as a risk factor for breast cancer.

Associations between breast density and a panel of single nucleotide polymorphisms linked to breast cancer risk: a cohort study with digital mammography

Background

Breast density and single-nucleotide polymorphisms (SNPs) have both been associated with breast cancer risk. To determine the extent to which these two breast cancer risk factors are associated, we investigate the association between a panel of validated SNPs related to breast cancer and quantitative measures of mammographic density in a cohort of Caucasian and African-American women.

Methods

In this IRB-approved, HIPAA-compliant study, we analyzed a screening population of 639 women (250 African American and 389 Caucasian) who were tested with a validated panel assay of 12 SNPs previously associated to breast cancer risk. Each woman underwent digital mammography as part of routine screening and all were interpreted as negative. Both absolute and percent estimates of area and volumetric density were quantified on a per-woman basis using validated software. Associations between the number of risk alleles in each SNP and the density measures were assessed through a race-stratified linear regression analysis, adjusted for age, BMI, and Gail lifetime risk.

Results

The majority of SNPs were not found to be associated with any measure of breast density. SNP rs3817198 (in LSP1) was significantly associated with both absolute area (p = 0.004) and volumetric (p = 0.019) breast density in Caucasian women. In African-American women, SNPs rs3803662 (in TNRC9/TOX3) and rs4973768 (in NEK10) were significantly associated with absolute (p = 0.042) and percent (p = 0.028) volume density respectively.

Conclusions

The majority of SNPs investigated in our study were not found to be significantly associated with breast density, even when accounting for age, BMI, and Gail risk, suggesting that these two different risk factors contain potentially independent information regarding a woman’s risk to develop breast cancer. Additionally, the few statistically significant associations between breast density and SNPs were different for Caucasian versus African American women. Larger prospective studies are warranted to validate our findings and determine potential implications for breast cancer risk assessment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1159-3) contains supplementary material, which is available to authorized users.

Measurement of breast density with digital breast tomosynthesis--a systematic review

Digital breast tomosynthesis (DBT) has gained acceptance as an adjunct to digital mammography in screening. Now that breast density reporting is mandated in several states in the USA, it is increasingly important that the methods of breast density measurement be robust, reliable and consistent. Breast density assessment with DBT needs some consideration since quantitative methods are modelled for two-dimensional (2D) mammography. A review of methods used for breast density assessment with DBT was performed. Existing evidence shows Cumulus has better reproducibility than that of the breast imaging reporting and data system (BI-RADS®) but still suffers from subjective variability; MedDensity is limited by image noise, whilst Volpara and Quantra are robust and consistent. The reported BI-RADs inter-reader breast density agreement (k) ranged from 0.65 to 0.91, with inter-reader correlation (r) ranging from 0.70 to 0.93. The correlation (r) between BI-RADS and Cumulus ranged from 0.54-0.94, whilst that of BI-RADs and MedDensity ranged from 0.48-0.78. The reported agreement (k) between BI-RADs and Volpara is 0.953. Breast density correlation between DBT and 2D mammography ranged from 0.73 to 0.97, with agreement (k) ranging from 0.56 to 0.96. To avoid variability and provide more reliable breast density information for clinicians, automated volumetric methods are preferred.

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.

Automated Volumetric Breast Density derived by Shape and Appearance Modeling

The image shape and texture (appearance) estimation designed for facial recognition is a novel and promising approach for application in breast imaging. The purpose of this study was to apply a shape and appearance model to automatically estimate percent breast fibroglandular volume (%FGV) using digital mammograms. We built a shape and appearance model using 2000 full-field digital mammograms from the San Francisco Mammography Registry with known %FGV measured by single energy absorptiometry method. An affine transformation was used to remove rotation, translation and scale. Principal Component Analysis (PCA) was applied to extract significant and uncorrelated components of %FGV. To build an appearance model, we transformed the breast images into the mean texture image by piecewise linear image transformation. Using PCA the image pixels grey-scale values were converted into a reduced set of the shape and texture features. The stepwise regression with forward selection and backward elimination was used to estimate the outcome %FGV with shape and appearance features and other system parameters. The shape and appearance scores were found to correlate moderately to breast %FGV, dense tissue volume and actual breast volume, body mass index (BMI) and age. The highest Pearson correlation coefficient was equal 0.77 for the first shape PCA component and actual breast volume. The stepwise regression method with ten-fold cross-validation to predict %FGV from shape and appearance variables and other system outcome parameters generated a model with a correlation of r² = 0.8. In conclusion, a shape and appearance model demonstrated excellent feasibility to extract variables useful for automatic %FGV estimation. Further exploring and testing of this approach is warranted.

Agreement of mammographic measures of volumetric breast density to MRI

Background

Clinical scores of mammographic breast density are highly subjective. Automated technologies for mammography exist to quantify breast density objectively, but the technique that most accurately measures the quantity of breast fibroglandular tissue is not known.

Purpose

To compare the agreement of three automated mammographic techniques for measuring volumetric breast density with a quantitative volumetric MRI-based technique in a screening population.

Materials and Methods

Women were selected from the UCSF Medical Center screening population that had received both a screening MRI and digital mammogram within one year of each other, had Breast Imaging Reporting and Data System (BI-RADS) assessments of normal or benign finding, and no history of breast cancer or surgery. Agreement was assessed of three mammographic techniques (Single-energy X-ray Absorptiometry [SXA], Quantra, and Volpara) with MRI for percent fibroglandular tissue volume, absolute fibroglandular tissue volume, and total breast volume.

Results

Among 99 women, the automated mammographic density techniques were correlated with MRI measures with R² values ranging from 0.40 (log fibroglandular volume) to 0.91 (total breast volume). Substantial agreement measured by kappa statistic was found between all percent fibroglandular tissue measures (0.72 to 0.63), but only moderate agreement for log fibroglandular volumes. The kappa statistics for all percent density measures were highest in the comparisons of the SXA and MRI results. The largest error source between MRI and the mammography techniques was found to be differences in measures of total breast volume.

Conclusion

Automated volumetric fibroglandular tissue measures from screening digital mammograms were in substantial agreement with MRI and if associated with breast cancer could be used in clinical practice to enhance risk assessment and prevention.

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.

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.

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

RATIONALE AND OBJECTIVES

Mammographic breast density is an important and widely accepted risk factor for breast cancer. A statement about breast density in the mammographic report is becoming a requirement in many States. However, there is significant inter-observer variation between radiologists in their interpretation of breast density. A properly designed automated system could provide benefits in maintaining consistency and reproducibility. We have developed a new automated and calibrated measure of breast density using full field digital mammography (FFDM). This new measure assesses spatial variation within a mammogram and produced significant associations with breast cancer in a small study. The costs of this automation are delays from advanced image and data analyses before the study can be processed. We evaluated this new calibrated variation measure using a larger dataset than previously. We also explored the possibility of developing an automated measure from unprocessed (raw data) mammograms as an approximation for this calibrated breast density measure.

MATERIALS AND METHODS

A case-control study comprised of 160 cases and 160 controls matched by age, screening history, and hormone replacement therapy was used to compare the calibrated variation measure of breast density with three variants of a noncalibrated measure of spatial variation. The operator-assisted percentage of breast density measure (PD) was used as a standard reference for comparison. Odds ratio (OR) quartile analysis was used to compare these measures. Linear regression analysis was applied to assess the calibration's impact on the raw pixel distribution.

RESULTS

All breast density measures showed significant breast cancer associations. The calibrated spatial variation measure produced the strongest associations (OR: 1.0 [ref.], 4.6, 4.3, 7.4). The associations for PD were diminished in comparison (OR: 1.0 [ref.], 2.7, 2.9, 5.2). Two additional non-calibrated measures restricted in region size also showed significant associations (OR: 1.0 [ref.], 2.9, 4.4, 5.4), and (OR: 1.0 [ref.], 3.5, 3.1, 4.9). Regression analyses indicated the raw image mean is influenced by the calibration more so than its standard deviation.

CONCLUSION

Breast density measures can be automated. The associated calibration produced risk information not retrievable from the raw data representation. Although the calibrated measure produced the stronger association, the non-calibrated measures may offer an alternative to PD and other operator based methods after further evaluation, because they can be implemented automatically with a simple processing algorithm.

Mammographic density and breast cancer risk: current understanding and future prospects

Variations in percent mammographic density (PMD) reflect variations in the amounts of collagen and number of epithelial and non-epithelial cells in the breast. Extensive PMD is associated with a markedly increased risk of invasive breast cancer. The PMD phenotype is important in the context of breast cancer prevention because extensive PMD is common in the population, is strongly associated with risk of the disease, and, unlike most breast cancer risk factors, can be changed. Work now in progress makes it likely that measurement of PMD will be improved in the near future and that understanding of the genetics and biological basis of the association of PMD with breast cancer risk will also improve. Future prospects for the application of PMD include mammographic screening, risk prediction in individuals, breast cancer prevention research, and clinical decision making.