Scientists Map Hidden Lymphatic Damage in Breast Cancer Survivors Using AI

Breast cancer saves lives, but for roughly one in five survivors, the treatment that cures the tumor leaves behind a quieter, chronic burden: breast cancer-related lymphedema, or BCRL. This condition arises when surgery or radiotherapy disrupts the lymphatic drainage of the arm, causing fluid to accumulate, fat tissue to proliferate, and lymph to reflux backward into the skin, a phenomenon clinicians call dermal backflow. Now, a research team spanning the University of Auckland in New Zealand and Macquarie University in Australia has transformed hundreds of clinician-drawn diagnostic diagrams into something unprecedented: a spatial atlas of where lymphatic dysfunction strikes the upper limb, paired with machine learning models that can stage disease severity objectively. The study, published in Breast Cancer Research and Treatment, draws on one of the largest imaging cohorts ever assembled for this purpose and could reshape how a debilitating complication is detected and monitored.

The data came from the Australian Lymphedema Education, Research and Treatment (ALERT) Centre Databank at Macquarie University in Sydney. The researchers extracted records from 282 female patients with confirmed BCRL, covering 289 affected limbs, including seven patients with lymphedema on both sides. Each patient had undergone indocyanine green (ICG) lymphography, an imaging technique in which a fluorescent tracer is injected into the skin and visualized with a near-infrared camera. During these sessions, four trained clinicians drew the observed lymphatic drainage pathways directly onto paper templates of the patient’s arm, marking regions of dermal backflow in red once the signal stabilized after roughly an hour. These annotated summary diagrams, digitized as high-resolution images, became the raw material for the atlas. Alongside the images, the databank provided rich clinical information: age, duration of lymphedema, cancer treatment history, limb volume difference, bioimpedance spectroscopy scores known as L-Dex, body mass index, and the clinician-assigned MD Anderson Cancer Center (MDACC) stage of lymphedema severity.

Turning hand-drawn diagrams from hundreds of patients into a single coherent map posed a genuine computational challenge. The team chose a left-arm template as the reference space and horizontally flipped diagrams from right-sided cases, an approach justified by the known symmetry of normal upper-limb lymphatic drainage when assessed by nodal basin. Each diagram was then aligned to the template using feature-based image co-registration: the ORB algorithm, which combines the FAST keypoint detector with the BRIEF descriptor, identified matching landmarks, while the RANSAC algorithm estimated the geometric transformation needed to overlay every image onto the common frame. Segmentation of the red dermal backflow regions was performed semi-automatically by thresholding in the CIELAB perceptual color space, exploiting the fact that the A-channel captures the green-red axis of an image. Morphological operations cleaned up small artifacts, and the researchers manually corrected masks where necessary using image editing software before quality checks.

Once every diagram was registered and segmented, the binary masks, marking background as zero and dermal backflow as one, were stacked into a three-dimensional volume with one slice per patient study. A pixel-wise mean across the stack produced a frequency map: for every anatomical location on the arm, the proportion of patients exhibiting dermal backflow at that spot. The result is, in effect, a probabilistic heat map of lymphatic failure painted onto the anatomy of the upper limb. The team also generated subgroup-specific heat maps stratified by MDACC stage, L-Dex ratio, limb volume difference, lymphatic drainage region, and the presence of backflow in the hand, allowing them to see how the geography of disease shifts as severity increases.

The atlas revealed striking and reproducible spatial patterns. In the anterior view, dermal backflow reached its highest prevalence, up to roughly 80 percent, in the medial wrist and forearm, then faded progressively toward the elbow and upper arm. In the posterior view, prevalence climbed even higher, approaching 90 percent in the olecranon region and forearm. By contrast, the fingers, the shoulder, and, notably, the region near the axilla showed the lowest frequencies, below 10 percent. The thenar eminence at the base of the thumb was the most commonly affected part of the hand. This distal-to-proximal gradient suggests that lymphatic dysfunction in BCRL is not a diffuse, uniform process but a regionally patterned one, with certain anatomical zones acting as consistent hotspots of failure across hundreds of individuals.

The subgroup analyses added clinical depth. As MDACC stage increased from 1 to 4, dermal backflow became more extensive and spread more proximally up the limb. Higher L-Dex ratios and greater limb volume differences were likewise associated with broader regions of involvement. Patients whose lymph drained through alternative routes beyond the predominant ipsilateral axilla showed distinct patterns, and those with hand involvement displayed more distally concentrated backflow with less upper-arm involvement. Perhaps most clinically provocative, dermal backflow was frequently observed in limbs with minimal or even negative limb volume difference, hinting that lymphatic dysfunction can precede measurable swelling. Because negative volume differences can arise from arm dominance and baseline variability, the finding underscores a key limitation of conventional volume- and impedance-based assessments: they may underestimate early or spatially localized disease that a spatial heat map can capture as a continuous probability signal.

Why the forearm and wrist should be so vulnerable remains an open question, but the authors point to the underlying organization of upper-extremity lymphatic drainage and the lymphosome concept, which divides the body into discrete lymphatic territories. The antecubital region is associated with drainage toward the epitrochlear nodes en route to the axilla, while the low prevalence near the axilla itself may simply reflect the normal convergence of superficial lymphatics into deeper nodal structures as they approach their destination. The patterns also echo earlier physiological work: a study by Mayrovitz and colleagues reported regional variability in superficial tissue water content in women treated for breast cancer even before overt lymphedema appeared. The new atlas suggests that this regional susceptibility is preserved, and possibly amplified, in established disease, with distal regions showing roughly 70 percent prevalence compared with about 30 percent proximally.

The second half of the study translated these imaging insights into prediction. Using Spearman’s rank correlation, the team narrowed 28 clinical variables down to 13 statistically significant and clinically relevant features. Dermal backflow percentage, calculated as the segmented backflow area divided by the total upper-limb area across anterior and posterior views, showed the strongest correlation with MDACC stage at rho equals 0.90, far exceeding L-Dex at 0.67 and limb volume difference at 0.65. Five machine learning models were then trained to classify stage: multinomial logistic regression, random forest, multi-layer perceptron, k-nearest neighbors, and support vector machine, evaluated with five-fold stratified cross-validation under a one-versus-rest framework.

The results were clear. When dermal backflow percentage was included, random forest achieved the highest area under the curve at 0.94 with precision of 0.768, while multinomial logistic regression posted the best F1 score at 0.926 and sensitivity of 0.735, both with accuracy around 0.737 and specificity above 0.92. The neural network and support vector machine performed at an intermediate level, and k-nearest neighbors lagged throughout. Critically, removing dermal backflow from the feature set degraded every model, with random forest’s AUC falling to 0.832. The authors note that the simpler, interpretable models outperforming deep architectures is good news for clinical translation, where understanding why a model assigns a stage matters as much as the assignment itself. Interestingly, the non-linear models showed higher sensitivity for advanced stage 4 disease, though by that point lymphedema is typically clinically obvious anyway.

The study has limitations the authors acknowledge candidly. The atlas reflects clinician interpretation of ICG studies rather than raw imaging data, and although standardized training reduced variability among the four annotators, differences in interpretation may have influenced the maps. A U-Net deep learning segmentation model was trialed to automate the pipeline but degraded on the full dataset, particularly in the small and variable hand regions, so manual correction remained necessary. The two-dimensional diagrams also approximate, rather than perfectly capture, the curved three-dimensional skin surface. Still, the framework opens concrete paths forward: external validation in multi-centre cohorts, fully automated registration and segmentation, three-dimensional volumetric analysis, and investigation of whether sentinel versus axillary lymph node dissection produce different backflow patterns. For now, the atlas offers clinicians a reference for what typical lymphatic failure looks like across a large cohort, and a quantitative, reproducible measure of severity that could complement, rather than replace, existing assessment, potentially catching the disease before the tape measure can.

Subject of Research: Spatial mapping of dermal backflow and machine learning-based staging of breast cancer-related lymphedema

Article Title: A dermal backflow atlas of lymphatic dysfunction in breast cancer-related lymphedema with machine learning-based stage prediction

Article References: Sutejitsiri, L., Suami, H., Nielsen, P. M. F., Phillips, A. R. J., & Reynolds, H. M. (2026). A dermal backflow atlas of lymphatic dysfunction in breast cancer-related lymphedema with machine learning-based stage prediction. Breast Cancer Research and Treatment, 219(3), Article 28. https://doi.org/10.1007/s10549-026-08093-2

Image Credits: AI Generated

DOI: 10.1007/s10549-026-08093-2

Keywords: breast cancer-related lymphedema, dermal backflow, ICG lymphography, machine learning, spatial atlas, MDACC staging, lymphatic drainage, random forest, bioimpedance spectroscopy, upper limb, image segmentation, early detection