This retrospective study was approved by the institutional ethics committee 2024YX070 and conducted in accordance with the Declaration of Helsinki. Given the non-interventional, de-identified nature of the dataset and minimal risk, informed consent was waived per institutional policy. Data were stored on secure, access-controlled servers, and analyses used an anonymized, locked dataset with an auditable trail and a pre-specified data-freeze date.
Study design and setting
We performed a single-center retrospective cohort study including patients with breast cancer (BC) who developed a second primary lung cancer (pLC) (BC-pLC) and contemporaneous patients with BC only. All patients underwent curative-intent surgery between January 2012 and January 2017, with follow-up through June 30, 2023. The BC-pLC cohort required pathologic confirmation of both primary BC and primary lung cancer; the BC cohort included patients with primary BC only. Propensity score matching was not performed and is acknowledged as a limitation; to mitigate selection bias, covariates were pre-specified for multivariable adjustment and are reported transparently. Synchronous BC-pLC was defined as a diagnostic interval of ≤6 months between primaries and metachronous BC-pLC as >6 months. A flow diagram summarizing screening, exclusions, inclusion, and analyses is presented in Figure 1.
Eligibility criteria and case adjudication
Eligible participants were aged ≥18 years with pathologically confirmed early-stage primary BC, available core clinicopathologic variables, imaging for review, and blocks/slides for central pathology review when indicated. Patients were excluded if they had a concurrent malignancy at another site, radiologic or pathologic evidence of pulmonary metastasis from BC, metastatic or recurrent BC at the index diagnosis, or incomplete follow-up data. Classification of lung lesions as primary lung cancer rather than metastatic BC was based on histomorphology and, when available, an immunophenotypic profile consistent with lung origin (e.g., TTF-1/Napsin A for adenocarcinoma and p40 for squamous differentiation). Discordant or indeterminate cases were resolved by multidisciplinary consensus (breast surgery, thoracic surgery, radiology, pathology).
Imaging protocol and surveillance
During BC care, chest imaging followed institutional pathways that included baseline staging CT at diagnosis and surveillance CT typically every 6-12 months during the first 2-3 years, and annually thereafter, or earlier if symptoms or abnormal findings prompted imaging. Acquisition parameters were standardized per radiology protocols. Because surveillance intensity can affect detection and stage distribution, potential detection/lead-time bias was addressed in sensitivity analyses by considering imaging intensity and calendar period.
Pathology workflow and immunohistochemistry
Tumor specimens were fixed in neutral buffered formalin, 4 %, typically 6-24 min/mm of tissue thickness (conventionally totaling 6-24 h for routine specimens), dehydrated, and paraffin-embedded before sectioning. Sections of 3-4 µm were cut and mounted on charged slides. Heat-induced epitope retrieval was performed using either citrate buffer (pH 6.0) or EDTA buffer (pH 9.0) at 95-100 °C for approximately 20 min, followed by cooling at room temperature for 10 min. Primary antibodies for ER, PR, HER2, and Ki-67 were applied according to validated laboratory protocols with typical incubation of 30-60 min at room temperature, polymer-based secondary detection, and DAB development for 5-10 min; slides were counterstained with hematoxylin for 30-60 s, dehydrated, cleared, and mounted. ER and PR positivity were defined as ≥1% tumor-cell nuclear staining. HER2 followed ASCO/CAP 2018 criteria (3+ by IHC or amplification by in situ hybridization when 2+ and equivocal). Ki-67 was dichotomized at a pre-specified cut-off of 14%. Two board-certified pathologists independently reviewed all slides, and discrepancies were resolved by consensus to standardize scoring.
Biospecimen storage conditions
Formalin-fixed, paraffin-embedded (FFPE) blocks were stored at controlled room temperature (18-25 °C), dry and dark, with ambient humidity maintained below 60%. Blocks were retained for ≥10 years per institutional policy, and re-cutting from the original block was performed whenever section age exceeded predefined stability windows. Unstained sections (3-4 µm) were stored in sealed slide boxes with desiccant, protected from light. For immunohistochemistry (ER, PR, Ki-67), unstained slides were refrigerated at 2-8 °C and stained within 8 weeks of sectioning; slides older than this threshold were discarded and re-cut. For in situ hybridization/FISH confirmation of equivocal HER2 results, unstained slides were used within 2-4 weeks of sectioning at 2-8 °C storage, or re-cut if older. Storage temperature and sectioning dates were logged in the laboratory information system, and any deviations triggered re-cutting prior to staining.
Data collection and management
Demographic and clinicopathologic variables were abstracted from the electronic medical record per a predefined data dictionary. Variables included age at diagnosis, menopausal status, hypertension, family history, primary tumor size, histologic subtype and grade, nodal status and number of examined nodes, lymphovascular invasion, and biomarker status for ER, PR, HER2, and Ki-67. Histologic subtypes and TNM staging were summarized descriptively in Table 1 and Table 2, with clarifying footnotes for abbreviations and subgroup definitions. Data quality procedures comprised range and logic checks and source verification by a second reviewer. Analyses were reproduced from a locked dataset with a documented codebook and manifest.
Outcomes and follow-up
Overall survival (OS) was measured from the date of cancer diagnosis to either death from any cause or the most recent follow-up. Invasive disease-free survival (iDFS) was calculated from the end of curative treatment to the occurrence of an invasive recurrence, the development of a new primary malignancy, or death from unrelated causes. Patients alive without events were censored at last contact. Follow-up was conducted via telephone contacts, clinic visits, and electronic medical record review through June 30, 2023. Cause-specific mortality was not adjudicated; competing-risks analyses were therefore not performed, and this limitation is noted in the Discussion.
Statistical analysis
Continuous variables were evaluated for normality using the Shapiro-Wilk test and summarized as median (IQR) or mean (SD) as appropriate, and categorical variables as n (%). Group comparisons used the t test or Wilcoxon rank-sum test and the chi-square test or Fisher's exact test where appropriate. For correlates of BC-pLC, univariable logistic regression was followed by multivariable logistic regression using a priori variables (tumor size category, ER, PR, Ki-67), with events-per-variable considerations, variance inflation factors for collinearity, and standard goodness-of-fit diagnostics. Survival was assessed with Kaplan-Meier estimates and log-rank tests and with Cox proportional-hazards models prespecifying age, tumor size, nodal status, ER, PR, HER2, Ki-67, and cohort (BC-pLC vs BC). Proportional-hazards assumptions were evaluated using Schoenfeld residuals. Missing data were handled with complete-case analysis, and the extent of missingness is reflected in denominators. Given multiple biomarker comparisons across ER, PR, HER2, and Ki-67, family-wise error control used Bonferroni adjustment with an adjusted ɑ of 0.0125, and nominal two-sided P values are also reported. All P values in the text and tables are paired with effect sizes and 95% confidence intervals.
Sensitivity and robustness analyses
Robustness was examined through interval-stratified survival analyses that contrasted synchronous cases with metachronous cases categorized as ≤12 months, 13-36 months, and >36 months, with adjusted hazard ratios reported alongside 95% confidence intervals and nominal P values. A 12 month landmark approach was used to mitigate guarantee-time bias where postoperative or adjuvant therapies might otherwise distort time-to-event estimates. Additional models adjusted for treatment classes where available, including lung surgery type, radiotherapy exposure, and systemic therapy categories (chemotherapy, targeted therapy, immunotherapy). Sensitivity analyses further considered imaging surveillance intensity and calendar period to address potential detection and secular-trend biases.
Reproducibility checkpoints
Key checkpoints influencing reproducibility included explicit adjudication criteria distinguishing primary lung cancer from metastatic BC, prespecified thresholds for ER, PR, HER2, and Ki-67, dual independent pathology reads with consensus resolution, transparent reporting of imaging surveillance intensity, and a priori covariate specification with multiplicity control. Any deviations from protocol were documented, including the rationale and an assessment of their potential impact on results.