[Skip to Navigation]
Sign In
Figure 1.  Cohort Flowchart
Cohort Flowchart

aUsed to assess natural randomization.

Figure 2.  Overall Survival by Blood Donor Sex and Parity
Overall Survival by Blood Donor Sex and Parity

The shaded area of the graphs indicates 95% CIs.

Figure 3.  Subgroup Analyses by Patient Sex and Age
Subgroup Analyses by Patient Sex and Age
Table 1.  Pretransfusion Cohort Characteristics at Baselinea
Pretransfusion Cohort Characteristics at Baselinea
Table 2.  Patient Characteristics Before Each Transfusion During Follow-upa
Patient Characteristics Before Each Transfusion During Follow-upa
1.
Chapman  CE, Stainsby  D, Jones  H,  et al; Serious Hazards of Transfusion Steering Group.  Ten years of hemovigilance reports of transfusion-related acute lung injury in the United Kingdom and the impact of preferential use of male donor plasma.   Transfusion. 2009;49(3):440-452. doi:10.1111/j.1537-2995.2008.01948.x PubMedGoogle ScholarCrossref
2.
Tynell  E, Andersson  TM, Norda  R,  et al.  Should plasma from female donors be avoided? a population-based cohort study of plasma recipients in Sweden from 1990 through 2002.   Transfusion. 2010;50(6):1249-1256. doi:10.1111/j.1537-2995.2009.02568.x PubMedGoogle ScholarCrossref
3.
Middelburg  RA, Van Stein  D, Zupanska  B,  et al.  Female donors and transfusion-related acute lung injury: a case-referent study from the International TRALI Unisex Research Group.   Transfusion. 2010;50(11):2447-2454. doi:10.1111/j.1537-2995.2010.02715.x PubMedGoogle ScholarCrossref
4.
Triulzi  DJ, Kleinman  S, Kakaiya  RM,  et al.  The effect of previous pregnancy and transfusion on HLA alloimmunization in blood donors: implications for a transfusion-related acute lung injury risk reduction strategy.   Transfusion. 2009;49(9):1825-1835. doi:10.1111/j.1537-2995.2009.02206.x PubMedGoogle ScholarCrossref
5.
Vossoughi  S, Gorlin  J, Kessler  DA,  et al.  Ten years of TRALI mitigation: measuring our progress.   Transfusion. 2019;59(8):2567-2574. doi:10.1111/trf.15387 PubMedGoogle ScholarCrossref
6.
Chassé  M, Tinmouth  A, English  SW,  et al.  Association of blood donor age and sex with recipient survival after red blood cell transfusion.   JAMA Intern Med. 2016;176(9):1307-1314. doi:10.1001/jamainternmed.2016.3324 PubMedGoogle ScholarCrossref
7.
Middelburg  RA, Briët  E, van der Bom  JG.  Mortality after transfusions, relation to donor sex.   Vox Sang. 2011;101(3):221-229. doi:10.1111/j.1423-0410.2011.01487.x PubMedGoogle ScholarCrossref
8.
Caram-Deelder  C, Kreuger  AL, Evers  D,  et al.  Association of blood transfusion from female donors with and without a history of pregnancy with mortality among male and female transfusion recipients.   JAMA. 2017;318(15):1471-1478. doi:10.1001/jama.2017.14825 PubMedGoogle ScholarCrossref
9.
Zeller  MP, Rochwerg  B, Jamula  E,  et al.  Sex-mismatched red blood cell transfusions and mortality: a systematic review and meta-analysis.   Vox Sang. 2019;114(5):505-516. doi:10.1111/vox.12783 PubMedGoogle ScholarCrossref
10.
Heddle  NM, Cook  RJ, Liu  Y,  et al.  The association between blood donor sex and age and transfusion recipient mortality: an exploratory analysis.   Transfusion. 2019;59(2):482-491. doi:10.1111/trf.15011 PubMedGoogle ScholarCrossref
11.
Bjursten  H, Dardashti  A, Björk  J, Wierup  P, Algotsson  L, Ederoth  P.  Transfusion of sex-mismatched and non-leukocyte-depleted red blood cells in cardiac surgery increases mortality.   J Thorac Cardiovasc Surg. 2016;152(1):223-232. doi:10.1016/j.jtcvs.2015.12.022 PubMedGoogle ScholarCrossref
12.
Edgren  G, Murphy  EL, Brambilla  DJ,  et al; NHLBI Recipient Epidemiology and Donor Evaluation Study-III (REDS-III) Group.  Association of blood donor sex and prior pregnancy with mortality among red blood cell transfusion recipients.   JAMA. 2019;321(22):2183-2192. doi:10.1001/jama.2019.7084 PubMedGoogle ScholarCrossref
13.
Edgren  G, Ullum  H, Rostgaard  K,  et al.  Association of donor age and sex with survival of patients receiving transfusions.   JAMA Intern Med. 2017;177(6):854-860. doi:10.1001/jamainternmed.2017.0890 PubMedGoogle ScholarCrossref
14.
Desmarets  M, Bardiaux  L, Benzenine  E,  et al.  Effect of storage time and donor sex of transfused red blood cells on 1-year survival in patients undergoing cardiac surgery: an observational study.   Transfusion. 2016;56(5):1213-1222. doi:10.1111/trf.13537 PubMedGoogle ScholarCrossref
15.
Robins  JM, Hernán  MA, Brumback  B.  Marginal structural models and causal inference in epidemiology.   Epidemiology. 2000;11(5):550-560. doi:10.1097/00001648-200009000-00011 PubMedGoogle ScholarCrossref
16.
Karim  ME, Tremlett  H, Zhu  F, Petkau  J, Kingwell  E.  Dealing with treatment-confounder feedback and sparse follow-up in longitudinal studies: application of a marginal structural model in a multiple sclerosis cohort.   Am J Epidemiol. 2021;190(5):908-917. doi:10.1093/aje/kwaa243 PubMedGoogle ScholarCrossref
17.
Mansournia  MA, Etminan  M, Danaei  G, Kaufman  JS, Collins  G.  Handling time varying confounding in observational research.   BMJ. 2017;359:j4587. doi:10.1136/bmj.j4587 PubMedGoogle ScholarCrossref
18.
Tilling  K, Sterne  JA, Szklo  M.  Estimating the effect of cardiovascular risk factors on all-cause mortality and incidence of coronary heart disease using G-estimation: the atherosclerosis risk in communities study.   Am J Epidemiol. 2002;155(8):710-718. doi:10.1093/aje/155.8.710 PubMedGoogle ScholarCrossref
19.
Mansournia  MA, Altman  DG.  Inverse probability weighting.   BMJ. 2016;352:i189. doi:10.1136/bmj.i189 PubMedGoogle ScholarCrossref
20.
Craig  P, Katikireddi  SV, Leyland  A, Popham  F.  Natural experiments: an overview of methods, approaches, and contributions to public health intervention research.   Annu Rev Public Health. 2017;38(1):39-56. doi:10.1146/annurev-publhealth-031816-044327 PubMedGoogle ScholarCrossref
21.
Hernán  MA, Robins  JM.  Using big data to emulate a target trial when a randomized trial is not available.   Am J Epidemiol. 2016;183(8):758-764. doi:10.1093/aje/kwv254 PubMedGoogle ScholarCrossref
22.
Hernán  MA.  Methods of public health research—strengthening causal inference from observational data.   N Engl J Med. 2021;385(15):1345-1348. doi:10.1056/NEJMp2113319 PubMedGoogle ScholarCrossref
23.
Wikman  A, Larsson  S, Storry  J, Schött  U, Sato Folatre  JGS.  Blood components, special components and whole blood—what and to whom?.   Lakartidningen. 2021;118:20150.PubMedGoogle Scholar
24.
Ludvigsson  JF, Otterblad-Olausson  P, Pettersson  BU, Ekbom  A.  The Swedish personal identity number: possibilities and pitfalls in healthcare and medical research.   Eur J Epidemiol. 2009;24(11):659-667. doi:10.1007/s10654-009-9350-y PubMedGoogle ScholarCrossref
25.
Zhao  J, Rostgaard  K, Hjalgrim  H, Edgren  G.  The Swedish Scandinavian donations and transfusions database (SCANDAT3-S)—50 years of donor and recipient follow-up.   Transfusion. 2020;60(12):3019-3027. doi:10.1111/trf.16027 PubMedGoogle ScholarCrossref
26.
Elixhauser  A, Steiner  C, Harris  DR, Coffey  RM.  Comorbidity measures for use with administrative data.   Med Care. 1998;36(1):8-27. doi:10.1097/00005650-199801000-00004 PubMedGoogle ScholarCrossref
27.
Sharma  N, Schwendimann  R, Endrich  O, Ausserhofer  D, Simon  M.  Comparing Charlson and Elixhauser comorbidity indices with different weightings to predict in-hospital mortality: an analysis of national inpatient data.   BMC Health Serv Res. 2021;21(1):13. doi:10.1186/s12913-020-05999-5 PubMedGoogle ScholarCrossref
28.
Auvinen  MK, Zhao  J, Lassén  E,  et al.  Patterns of blood use in Sweden from 2008 to 2017: a nationwide cohort study.   Transfusion. 2020;60(11):2529-2536. doi:10.1111/trf.16092 PubMedGoogle ScholarCrossref
29.
Austin  PC, Stuart  EA.  Moving towards best practice when using inverse probability of treatment weighting (IPTW) using the propensity score to estimate causal treatment effects in observational studies.   Stat Med. 2015;34(28):3661-3679. doi:10.1002/sim.6607 PubMedGoogle ScholarCrossref
30.
Bandara  SN, Kennedy-Hendricks  A, Stuart  EA,  et al.  The effects of the Maryland Medicaid health home waiver on emergency department and inpatient utilization among individuals with serious mental illness.   Gen Hosp Psychiatry. 2020;64:99-104. doi:10.1016/j.genhosppsych.2019.12.004 PubMedGoogle ScholarCrossref
31.
Hernán  MA, Lanoy  E, Costagliola  D, Robins  JM.  Comparison of dynamic treatment regimes via inverse probability weighting.   Basic Clin Pharmacol Toxicol. 2006;98(3):237-242. doi:10.1111/j.1742-7843.2006.pto_329.x PubMedGoogle ScholarCrossref
32.
Cole  SR, Hernán  MA.  Constructing inverse probability weights for marginal structural models.   Am J Epidemiol. 2008;168(6):656-664. doi:10.1093/aje/kwn164 PubMedGoogle ScholarCrossref
33.
Williamson  T, Ravani  P.  Marginal structural models in clinical research: when and how to use them?   Nephrol Dial Transplant. 2017;32(suppl 2):ii84-ii90. doi:10.1093/ndt/gfw341 PubMedGoogle ScholarCrossref
34.
Shinozaki  T, Suzuki  E.  Understanding marginal structural models for time-varying exposures: pitfalls and tips.   J Epidemiol. 2020;30(9):377-389. doi:10.2188/jea.JE20200226 PubMedGoogle ScholarCrossref
35.
Dickerman  BA, García-Albéniz  X, Logan  RW, Denaxas  S, Hernán  MA.  Avoidable flaws in observational analyses: an application to statins and cancer.   Nat Med. 2019;25(10):1601-1606. doi:10.1038/s41591-019-0597-x PubMedGoogle ScholarCrossref
36.
Edgren  G, Hjalgrim  H.  Epidemiology of donors and recipients: lessons from the SCANDAT database.   Transfus Med. 2019;29(suppl 1):6-12. doi:10.1111/tme.12487 PubMedGoogle ScholarCrossref
37.
Pocock  SJ, Conthe  A.  Ever-pregnant female blood donors and mortality risk in male recipients.   JAMA. 2018;319(10):1048-1049. doi:10.1001/jama.2017.21868 PubMedGoogle ScholarCrossref
38.
Hernán  MA, Brumback  B, Robins  JM.  Marginal structural models to estimate the causal effect of zidovudine on the survival of HIV-positive men.   Epidemiology. 2000;11(5):561-570. doi:10.1097/00001648-200009000-00012 PubMedGoogle ScholarCrossref
39.
Hernán  MA, Brumback  BA, Robins  JM.  Estimating the causal effect of zidovudine on CD4 count with a marginal structural model for repeated measures.   Stat Med. 2002;21(12):1689-1709. doi:10.1002/sim.1144 PubMedGoogle ScholarCrossref
40.
Fergusson  DA, Chassé  M, Tinmouth  A,  et al.  Pragmatic, double-blind, randomised trial evaluating the impact of red blood cell donor sex on recipient mortality in an academic hospital population: the innovative Trial Assessing Donor Sex (iTADS) protocol.   BMJ Open. 2021;11(2):e049598.PubMedGoogle Scholar
41.
Chassé  M, Fergusson  DA, Acker  JP,  et al. An Innovative Trial Assessing Donor Sex on Recipient Survival (iTADS). Abstract presented at: American Association of Blood Banks Virtual Annual Meeting; October 18, 2021.
Original Investigation
June 13, 2022

Mortality Among Patients Undergoing Blood Transfusion in Relation to Donor Sex and Parity: A Natural Experiment

Author Affiliations
  • 1Clinical Epidemiology Division, Department of Medicine, Solna, Karolinska Institutet, Stockholm, Sweden
  • 2Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden
  • 3Department of Cardiology, Södersjukhuset, Stockholm, Sweden
JAMA Intern Med. 2022;182(7):747-756. doi:10.1001/jamainternmed.2022.2115
Key Points

Question  Does blood donor sex or parity affect mortality of patients undergoing transfusion?

Findings  In this nationwide cohort study involving a natural experiment of 368 778 patients undergoing red blood cell transfusion, the 2-year survival differences comparing transfusions from female and parous donors with male donors were both not statistically significant at −0.1% and 0.3%, respectively. Receiving blood from female donors was associated with an increased risk of additional transfusions owing to female donors having lower hemoglobin counts.

Meaning  Patients undergoing transfusion with blood from female or parous donors did not have a higher mortality compared with patients undergoing transfusion with blood from male donors; however, female blood units tended to have lower hemoglobin content.

Abstract

Importance  Prior evidence suggests that plasma and platelet transfusions from female and parous donors are associated with adverse clinical outcomes, which has led to the predominant use of male donors for plasma and platelets in many countries. However, studies on red blood cell transfusions have been conflicting.

Objective  To determine whether blood donor sex and parity affect mortality of patients undergoing transfusion with red blood cells.

Design, Setting, and Participants  This cohort study used nationwide blood bank and health register data in Sweden and involved a natural experiment created by donor sex and parity being concealed and randomly allocated. Patients were included if they were 18 to 90 years old, did not have a history of blood transfusion, and received a transfusion between January 1, 2010, and December 31, 2017. Patients were followed up from their first red blood cell transfusion until death, emigration, or end of study (June 30, 2018). Data analysis was performed between June 15 and December 15, 2021.

Exposures  (1) Female vs male donors and (2) parous or nonparous female vs male donors.

Main Outcomes and Measures  Overall survival up to 2 years estimated using inverse probability–weighted Kaplan-Meier estimates and relative risk for additional transfusions within 24 hours.

Results  Among the 368 778 included patients (mean [SD] age, 66.3 [17.7] years; 57.3% female), 2-year survival differences comparing red blood cell transfusions from female and parous donors to male donors were −0.1% (95% CI, −1.3% to 1.1%) and 0.3% (95% CI, −0.6% to 1.2%), respectively. No statistically significant survival differences were observed regardless of patient sex or age. Median (IQR) hemoglobin counts for female donors (13.5 [13.0-14.0] g/dL) were lower than for male donors (14.9 [14.4-15.5] g/dL). Red blood cell transfusions from female donors were associated with a relative risk of 1.12 (95% CI, 1.08-1.17) for additional transfusions within 24 hours but not after adjusting for donor hemoglobin counts (relative risk, 1.03; 95% CI, 0.98-1.08). Pretransfusion patient characteristics were naturally distributed as-if randomized.

Conclusions and Relevance  In this nationwide cohort study involving a natural experiment, after accounting for the lower hemoglobin values in blood from female donors, patients undergoing transfusion with blood from female or parous donors did not have higher 2-year mortality compared with recipients of blood from male donors.

Introduction

Female donor sex and parity are established risk factors for transfusion-related acute lung injury (TRALI) following plasma and platelet transfusions, which is a leading cause of transfusion-related mortality.1-4 Preferential use of male plasma and platelets in many countries has reduced the incidence for TRALI.1,5 However, recent observational studies have also shown increased mortality in patients receiving red blood cell transfusions from female donors,6,7 in male patients receiving transfusions of red blood cells from previously pregnant but not nonparous female donors,8 and in sex-mismatched transfusions.7,9-11 These findings are conflicting with other studies of red blood cell transfusions, which reported no such associations.12-14 However, it is not fully understood why results from previous studies have been conflicting, both between and within studies, despite being conducted in similar settings with similar methods.6,7,9-14 Limiting the use of red blood cell units from female donors may seriously jeopardize the blood supply.

Previous studies have not accounted for comparatively lower hemoglobin counts in female donors compared with male donors. Head-to-head comparisons are therefore likely inappropriate because units of blood from female donors essentially constitute lower “doses” of hemoglobin. The difference in hemoglobin content may create a type of time-dependent confounding known as treatment-confounder feedback, which can lead to unpredictable bias when using standard regression models.15 This is a source of bias that may explain the conflicting results both between and within previous studies.15-19

In this nationwide cohort study involving a natural experiment, we assessed if lower hemoglobin counts in female donors may result in treatment-confounder feedback and be a previously unrecognized source of bias. Furthermore, we assessed if blood donor parity and sex were naturally assigned as-if randomized and constituted a natural experiment, as blood allocation was blinded to donor sex or parity. Based on this, we emulated a nationwide randomized clinical trial to study the effect of donor sex and parity on mortality of patients undergoing transfusion of red blood cells, using appropriate analytical methods that mitigate bias owing to treatment-confounder feedback.

Methods
Study Design

In a nationwide study of transfusion-naive persons receiving red blood cell transfusions from January 1, 2010, through June 30, 2018, we separately assessed overall survival in 2 sets of red blood cell transfusion strategies: (1) female vs male donors and (2) nonparous or parous female vs male donors. No other outcomes or blood products were assessed.

To minimize confounding, we exploited a natural experiment in the blood allocation process with respect to donor sex and parity.20 The rationale for the natural randomization is that the choice of red blood cell units was solely based on blood group and otherwise followed a first-in, first-out policy. Donor sex and parity were also neither labeled on the blood product nor routinely communicated to the treating physician or patient. Therefore, donor sex and parity are expected to be as-if randomized with respect to important patient characteristics and are effectively double blinded. However, because both donor and patient characteristics may vary slightly with blood group, geographical region, and calendar year, these 3 factors also need to be accounted for. In analogy with a multicenter study, we collectively refer to blood group, region, and calendar year as center characteristics. After empirically confirming that the distribution of donor sex and parity were de facto independent of patient characteristics and that the balance between treatment arms further improved after adjusting for center characteristics, we used inverse probability weighting to emulate a randomized target trial.21,22 In the Supplement, see eMethods 1 for a graphical summary of the method, eMethods 2 for formal definitions, and eTable 1 for the study protocol of the target and emulated trial.

The conduct of this study and waiver for informed consent owing to use of deidentified data only was approved by the Regional Ethics Committee in Stockholm and the Swedish Ethics Review Authority. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines were used.

Data Sources and Study Population

The study population was defined using data collected from all blood banks and hospitals in Sweden and included all persons who received a first red blood cell transfusion at 18 to 90 years old between 2010 and 2017. Patients were followed up from the time of the first transfusion for a maximum of 2 years or until midnight on the date of death, emigration, end of study (June 30, 2018), or hour of the first transfusion of a unit from an unidentified donor or autologous unit, whichever came first. The TRALI-mitigation policies, such as preferential male plasma use or TRALI-related antibody screening, did not apply to red blood cell donors during the study period. We did not differentiate red blood cell units based on whether they were leukoreduced or irradiated. Universal leukoreduction was gradually introduced in Sweden starting in the mid-1990s.23

All blood donors, predonation hemoglobin counts, links between donors and patients, ABO and RhD blood groups, and the date and time of transfusions were extracted from blood bank databases. Unique personal identifiers for both donors and patients undergoing transfusion were used to ascertain comorbidities from the National Patient Register and Cancer Register. Data on migration, sex, and vital status were retrieved from national population registers.24 Donor parity was assessed through the Medical Birth Register, containing all births in Sweden with nationwide coverage starting 1973, and the Multi-generation Register starting 1932. Laboratory data were retrieved from all hospitals in Stockholm, Sweden, for all patients treated in this region. Reporting to all data sources used in this study is mandated by Swedish law, ensuring high completeness.

All data were assembled into a previously reported research database denoted the SCANDAT3-S database, which contains data on all electronically recorded donations and transfusions in Sweden from the late 1960s until July 1, 2018, with full nationwide coverage starting in the mid-1990s.25 Full coverage of relevant blood tests ordered for all patients in Stockholm (the Swedish capital and largest city), both inpatient and outpatient, was available during the entire follow-up. Laboratory data and diagnosis data were only used for evaluating the assumption of natural randomization.

Assessment of Natural Randomization

As a measure of whether blood donor sex and parity were naturally distributed as-if randomized, we assessed pretransfusion patient characteristics across treatment arms. We used descriptive statistics both at baseline and throughout follow-up for patient demographics; comorbidities, including cancer and 31 comorbidities as defined in the Elixhauser Comorbidity Index; results of the most recent laboratory tests (≤1 week before transfusion) such as complete blood cell count, C-reactive protein, and serum creatinine; and initial transfusion indication. The Elixhauser Comorbidity Index and transfusion indication was derived from discharge diagnoses using the eighth, ninth, and tenth revisions of the International Classification of Diseases.26-28 For inpatient episodes, the date of admission was used as date of diagnosis. To assess the influence of center characteristics, we compared absolute standardized mean differences of baseline patient characteristics across treatment arms, before and after weighting with the inverse probability of donor sex or parity accounting for center characteristics. This was calculated as the absolute value of the difference in means, divided by the variance of the compared groups.29 Weighting adjusts for the center characteristics by balancing them across treatment arms.15,30 Full details and additional diagnostics are available in eMethods 2 in the Supplement.

Treatment-Confounder Feedback Owing to Donor Hemoglobin Counts

Treatment-confounder feedback is the phenomenon by which a confounder is affected by previous treatment. Herein, treatment-confounder feedback between female donor sex (the treatment) and the number of transfusions (the confounder) is expected, because female red blood cell units tend to contain less hemoglobin, which may result in a need for additional transfusions. The number of transfusions is a critical confounder because it is directly related to being exposed and is an important predictor for comorbidity and death. Therefore, the number of transfusions is both a mediator for the relation between previous transfusion and death, but also a confounder for the present transfusion and death. This cannot be adjusted for using standard regression models such as Cox regression. While standard regression models would correctly remove confounding for the present transfusion and death, it would also attenuate any effect that prior transfusions may have on death, which creates bias toward the null. Additionally, such an adjustment creates a spurious (noncausal) statistical association between the previous transfusion and death through unmeasured confounders of the number of transfusions and death, which can produce a bias away from the null.15 These unmeasured confounders include frailty and comorbidity, which are common causes of both needing a transfusion and death. Bias owing to treatment-confounder feedback is therefore unpredictable and can result in a net bias effect both toward and away from the null.15

To directly demonstrate the presence of treatment-confounder feedback, we calculated the relative risk of receiving additional red blood cell units in the first 24 hours for patients receiving female-donor units compared with male-donor units only, with and without adjusting for donor hemoglobin count. To accommodate the high prevalence of double-unit transfusions during the study period,28 donor sex was assessed for a subcohort (62%), who received 2 transfusions within 6 hours, and additional transfusions were assessed in the ensuing 18 hours after the second transfusion. See eMethods 2 in the Supplement for additional details and a formal discussion of treatment-confounder feedback.

Statistical Analyses

We emulated a target trial by exploiting the natural randomization of donor sex and parity. To mitigate bias owing to treatment-confounder feedback, we used cloning, censoring, and inverse probability weighting. This approach has been used extensively within pharmacoepidemiology and is conceptually identical to a marginal structural model.15,22,31-34 For a nontechnical summary, as well as a detailed, formal discussion of the method and a directed acyclic graph, see eFigure 1 in the Supplement.

In brief, each patient was cloned and assigned every treatment strategy and subsequently censored when they no longer adhered to the treatment strategy. To account for selection bias owing to censoring of nonadherent patients, patients who continued to adhere to the treatment strategy were time-varyingly weighted upward. At every transfusion, patients who received the transfusion of the assigned donor characteristic with the probability p would henceforth be upweighted by a factor of 1/p. These probabilities were estimated nonparametrically from administrative database data on the calendar year, region, and the patient blood group. Because a transfusion could change the need for a future transfusion already by the next hour, follow-up was divided into hourly periods. Finally, the treatment strategies were compared using Kaplan-Meier curves estimated from the weighted populations. Under the verifiable assumption that the allocation of blood units is random after adjusting for center characteristics, the estimated survival difference is an unbiased estimate of the causal effect that would be observed in a randomized trial with full compliance.

To account for it being less likely to receive many consecutive units from donors with a less common characteristic, adherence was assessed up to the number of units that corresponded to the 99th percentile of the least common arm and to the 99.9th percentile in sensitivity analyses. Corresponding 95% CIs were estimated using bootstrapping with 1000 runs.35 We did not separately account for censoring owing to autologous red blood cell transfusions or units from nonidentified donors because it was rare (n = 564 [0.15%]). All data processing and statistical analyses were performed with SAS, version 9.4 (SAS Institute).

Results
Cohort Construction

From a total of 485 327 red blood cell recipients, 368 778 patients without a history of red blood cell transfusions were identified. In total, 1 936 757 red blood cell units from 402 190 donors were transfused during follow-up. Among the patients available for analysis, 189 375 received a red blood cell unit from a male donor as their first unit, and 125 358 received a female-donor unit. Of the latter, 32 698 received a unit from a nonparous female donor, and 84 970 received a unit from a parous female donor. The cohort construction is shown in Figure 1. Inverse probability–weighted analyses were conducted per hour of follow-up, totaling 1 852 976 482 person-hours for donor sex and 1 894 360 834 person-hours for donor parity. Descriptive statistics for transfusions and events are available in eTable 2 in the Supplement.

Assessment of Natural Randomization

Age, sex, proportion born outside of Sweden, year of first transfusion, blood group, transfusion indication, prevalence of comorbidities, and pretransfusion blood chemistry were similar across treatment arms, both at baseline and throughout follow-up (Tables 1 and 2). Absolute standardized mean differences of extended baseline covariates before and after weighting for center characteristics are presented in eFigures 2 and 3 in the Supplement. Of all 31 Elixhauser Comorbidity Index groups across all treatment arms, the absolute standardized difference compared with male donor strategy was at most 2.8% and decreased to 1.2% after weighting (additional diagnostics are available in eFigures 4 and 5 in the Supplement).

Assessment of Treatment-Confounder Feedback Owing to Donor Hemoglobin Counts

The median (IQR) hemoglobin count was higher among male donors (14.9 [14.4-15.5] g/dL) than female donors (13.5 [13.0-14.0] g/dL). Receiving red blood cell units from female donors was associated with a 12% increased risk of an additional red blood cell transfusion within 24 hours compared with male donors (relative risk, 1.12; 95% CI, 1.08-1.17). However, the difference was no longer statistically significant after adjusting for donor hemoglobin counts (relative risk, 1.03; 95% CI, 0.98-1.08).

Female vs Male Donors

There was no statistically significant difference in overall survival at 2 years between female and male donor strategies (survival difference, −0.1%; 95% CI, −1.3% to 1.1%; Figure 2A). Exposure assessment was truncated after the fifth red blood cell unit, which was the 99th percentile of the number of transfused units for patients receiving blood from female donors.

Parous Female and Nonparous Female vs Male Donors

Figure 2B shows the overall survival for parous and nonparous female vs male donor strategies. Exposure assessment was truncated after the third red blood cell unit, which was the 99th percentile in the nonparous female donor arm. At 2 years, there were no statistically significant differences in overall survival comparing parous female donors (survival difference, 0.3%; 95% CI, −0.6% to 1.2%) and nonparous female donors (survival difference, 0.1%; 95% CI, −3.8% to 3.2%) with male donor strategy.

Sex-Discordance, Subgroup, and Sensitivity Analyses

The survival difference after 2 years for the 2 sets of strategies in subgroups of patient sex and age was examined (Figure 3). There was no statistically significant difference in any subgroup. In sensitivity analyses, exposure was assessed to the 99.9th percentile, corresponding to the first 7 and 4 units for donor sex and parity, respectively. No statistically significant survival differences were found for any treatment strategy (eFigures 6 and 7 and eTable 3 in the Supplement).

Discussion

In this nationwide, natural experiment of adult patients who received red blood cell transfusions over an 8.5-year period, patient mortality did not differ in relation to the sex or parity of the blood donor. These findings were consistent in subgroups of patient sex and age. We demonstrated theoretical and empirical support for donor parity and sex being naturally assigned as-if randomized to patients undergoing transfusion, which we exploited to emulate a randomized trial. Furthermore, we demonstrated that red blood cell transfusions from female donors incurred a higher risk for additional transfusions owing to lower donor hemoglobin counts in female donors. We applied a suitable analytical framework that mitigated bias caused by the difference in red blood cell unit efficacy and were able to estimate the effect comparable to a randomized clinical trial. Given the use of nationwide blood bank data over an 8.5-year period with exposure and outcome ascertainment from nationwide health registers, these estimates were both precise and should be widely generalizable.

Previous observational studies have been conflicting. Chassé et al6 reported a survival benefit among recipients of red blood cell units from male donors compared with female donors, whereas Caram-Deelder et al8 reported lower survival for male patients undergoing transfusion with red blood cells from ever-pregnant donors but not never-pregnant female donors. Several studies and 1 meta-analysis have shown increased mortality for sex-discordant red blood cell transfusions.7,9-11 However, these studies have been criticized for residual confounding owing to model misspecification,6 large proportion of missing data,8 and introducing selection bias through informative censoring.36,37 These limitations were addressed by a trinational study that reported mostly null results.12 Still, results from the trinational study were partly inconsistent, and one of the cohorts reported increased mortality for sex-discordant transfusions in 1 subgroup.12

However, previous studies did not consider that female donors have lower hemoglobin counts and that their units contain a lower dose of hemoglobin. This creates treatment-confounder feedback, a type of time-dependent confounding that cannot be accounted for using standard regression models.15,22,31-34 Because previous observational studies have used standard regression models, they are to some extent inevitably biased toward or away from the null, with an unpredictable net effect, and are difficult to interpret causally.38 The present study used cloning, censoring, and weighting, a well-recognized method within pharmacoepidemiology, to remedy bias owing to treatment-confounder feedback.15,22,31-34,38,39 These results are also consistent with results from a yet-published randomized trial that did not show a mortality difference in relation to donor sex.40,41

Limitations

The study has limitations. We did not assess extreme cases or exposures that did not occur in a real-world setting and truncated exposure at the 99th percentile of the number of transfusions from the same donor sex or parity. However, sensitivity analyses with exposure assessment up to the 99.9th percentile also yielded null findings. We were not able to explicitly consider special cases of phenotype-matched transfusions, although this was rare and should not be contingent on donor sex or parity and is therefore unlikely to result in uncontrolled confounding. We did not evaluate other blood products and did not differentiate between red blood cell units with and without leukoreduction because the latter was uncommon. Lastly, we did not include pediatric patients or evaluate outcomes other than death.

Conclusions

In this nationwide cohort study involving a natural experiment, patients undergoing transfusion with blood from female or parous donors did not have a higher 2-year mortality compared with recipients of blood from male donors. Red blood cell units from female donors contained less hemoglobin compared with units from male donors, which is a hitherto unrecognized source of confounding that cannot be accounted for using standard regression models and, to our knowledge, has not been considered in previous studies.

Back to top
Article Information

Accepted for Publication: April 21, 2022.

Published Online: June 13, 2022. doi:10.1001/jamainternmed.2022.2115

Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2022 Zhao J et al. JAMA Internal Medicine.

Corresponding Author: Jingcheng Zhao, MD, PhD, Clinical Epidemiology Division, Department of Medicine, Solna, Karolinska Institutet, Karolinska University Hospital, 171 76 Stockholm, Sweden (jingcheng.zhao@ki.se).

Author Contributions: Drs Zhao and Edgren had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

Concept and design: All authors.

Acquisition, analysis, or interpretation of data: All authors.

Drafting of the manuscript: All authors.

Critical revision of the manuscript for important intellectual content: All authors.

Statistical analysis: All authors.

Obtained funding: Edgren.

Supervision: Sjölander, Edgren.

Conflict of Interest Disclosures: Dr Zhao reported grants from Karolinska Institutet, Region Stockholm, and Vetenskapsrådet during the conduct of the study. Dr Edgren reported grants from Vetenskapsrådet and Stockholm Region during the conduct of the study. No other disclosures were reported.

Funding/Support: The creation of the SCANDAT3-S database and the conduct of this study were made possible by a grant to Dr Edgren from the Swedish Research Council (2017-01954). Dr Edgren is supported by a clinical research appointment by Region Stockholm. Dr Zhao is supported by the Clinical Scientist Training Program and the Research Internship Program, both at Karolinska Institutet.

Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

References
1.
Chapman  CE, Stainsby  D, Jones  H,  et al; Serious Hazards of Transfusion Steering Group.  Ten years of hemovigilance reports of transfusion-related acute lung injury in the United Kingdom and the impact of preferential use of male donor plasma.   Transfusion. 2009;49(3):440-452. doi:10.1111/j.1537-2995.2008.01948.x PubMedGoogle ScholarCrossref
2.
Tynell  E, Andersson  TM, Norda  R,  et al.  Should plasma from female donors be avoided? a population-based cohort study of plasma recipients in Sweden from 1990 through 2002.   Transfusion. 2010;50(6):1249-1256. doi:10.1111/j.1537-2995.2009.02568.x PubMedGoogle ScholarCrossref
3.
Middelburg  RA, Van Stein  D, Zupanska  B,  et al.  Female donors and transfusion-related acute lung injury: a case-referent study from the International TRALI Unisex Research Group.   Transfusion. 2010;50(11):2447-2454. doi:10.1111/j.1537-2995.2010.02715.x PubMedGoogle ScholarCrossref
4.
Triulzi  DJ, Kleinman  S, Kakaiya  RM,  et al.  The effect of previous pregnancy and transfusion on HLA alloimmunization in blood donors: implications for a transfusion-related acute lung injury risk reduction strategy.   Transfusion. 2009;49(9):1825-1835. doi:10.1111/j.1537-2995.2009.02206.x PubMedGoogle ScholarCrossref
5.
Vossoughi  S, Gorlin  J, Kessler  DA,  et al.  Ten years of TRALI mitigation: measuring our progress.   Transfusion. 2019;59(8):2567-2574. doi:10.1111/trf.15387 PubMedGoogle ScholarCrossref
6.
Chassé  M, Tinmouth  A, English  SW,  et al.  Association of blood donor age and sex with recipient survival after red blood cell transfusion.   JAMA Intern Med. 2016;176(9):1307-1314. doi:10.1001/jamainternmed.2016.3324 PubMedGoogle ScholarCrossref
7.
Middelburg  RA, Briët  E, van der Bom  JG.  Mortality after transfusions, relation to donor sex.   Vox Sang. 2011;101(3):221-229. doi:10.1111/j.1423-0410.2011.01487.x PubMedGoogle ScholarCrossref
8.
Caram-Deelder  C, Kreuger  AL, Evers  D,  et al.  Association of blood transfusion from female donors with and without a history of pregnancy with mortality among male and female transfusion recipients.   JAMA. 2017;318(15):1471-1478. doi:10.1001/jama.2017.14825 PubMedGoogle ScholarCrossref
9.
Zeller  MP, Rochwerg  B, Jamula  E,  et al.  Sex-mismatched red blood cell transfusions and mortality: a systematic review and meta-analysis.   Vox Sang. 2019;114(5):505-516. doi:10.1111/vox.12783 PubMedGoogle ScholarCrossref
10.
Heddle  NM, Cook  RJ, Liu  Y,  et al.  The association between blood donor sex and age and transfusion recipient mortality: an exploratory analysis.   Transfusion. 2019;59(2):482-491. doi:10.1111/trf.15011 PubMedGoogle ScholarCrossref
11.
Bjursten  H, Dardashti  A, Björk  J, Wierup  P, Algotsson  L, Ederoth  P.  Transfusion of sex-mismatched and non-leukocyte-depleted red blood cells in cardiac surgery increases mortality.   J Thorac Cardiovasc Surg. 2016;152(1):223-232. doi:10.1016/j.jtcvs.2015.12.022 PubMedGoogle ScholarCrossref
12.
Edgren  G, Murphy  EL, Brambilla  DJ,  et al; NHLBI Recipient Epidemiology and Donor Evaluation Study-III (REDS-III) Group.  Association of blood donor sex and prior pregnancy with mortality among red blood cell transfusion recipients.   JAMA. 2019;321(22):2183-2192. doi:10.1001/jama.2019.7084 PubMedGoogle ScholarCrossref
13.
Edgren  G, Ullum  H, Rostgaard  K,  et al.  Association of donor age and sex with survival of patients receiving transfusions.   JAMA Intern Med. 2017;177(6):854-860. doi:10.1001/jamainternmed.2017.0890 PubMedGoogle ScholarCrossref
14.
Desmarets  M, Bardiaux  L, Benzenine  E,  et al.  Effect of storage time and donor sex of transfused red blood cells on 1-year survival in patients undergoing cardiac surgery: an observational study.   Transfusion. 2016;56(5):1213-1222. doi:10.1111/trf.13537 PubMedGoogle ScholarCrossref
15.
Robins  JM, Hernán  MA, Brumback  B.  Marginal structural models and causal inference in epidemiology.   Epidemiology. 2000;11(5):550-560. doi:10.1097/00001648-200009000-00011 PubMedGoogle ScholarCrossref
16.
Karim  ME, Tremlett  H, Zhu  F, Petkau  J, Kingwell  E.  Dealing with treatment-confounder feedback and sparse follow-up in longitudinal studies: application of a marginal structural model in a multiple sclerosis cohort.   Am J Epidemiol. 2021;190(5):908-917. doi:10.1093/aje/kwaa243 PubMedGoogle ScholarCrossref
17.
Mansournia  MA, Etminan  M, Danaei  G, Kaufman  JS, Collins  G.  Handling time varying confounding in observational research.   BMJ. 2017;359:j4587. doi:10.1136/bmj.j4587 PubMedGoogle ScholarCrossref
18.
Tilling  K, Sterne  JA, Szklo  M.  Estimating the effect of cardiovascular risk factors on all-cause mortality and incidence of coronary heart disease using G-estimation: the atherosclerosis risk in communities study.   Am J Epidemiol. 2002;155(8):710-718. doi:10.1093/aje/155.8.710 PubMedGoogle ScholarCrossref
19.
Mansournia  MA, Altman  DG.  Inverse probability weighting.   BMJ. 2016;352:i189. doi:10.1136/bmj.i189 PubMedGoogle ScholarCrossref
20.
Craig  P, Katikireddi  SV, Leyland  A, Popham  F.  Natural experiments: an overview of methods, approaches, and contributions to public health intervention research.   Annu Rev Public Health. 2017;38(1):39-56. doi:10.1146/annurev-publhealth-031816-044327 PubMedGoogle ScholarCrossref
21.
Hernán  MA, Robins  JM.  Using big data to emulate a target trial when a randomized trial is not available.   Am J Epidemiol. 2016;183(8):758-764. doi:10.1093/aje/kwv254 PubMedGoogle ScholarCrossref
22.
Hernán  MA.  Methods of public health research—strengthening causal inference from observational data.   N Engl J Med. 2021;385(15):1345-1348. doi:10.1056/NEJMp2113319 PubMedGoogle ScholarCrossref
23.
Wikman  A, Larsson  S, Storry  J, Schött  U, Sato Folatre  JGS.  Blood components, special components and whole blood—what and to whom?.   Lakartidningen. 2021;118:20150.PubMedGoogle Scholar
24.
Ludvigsson  JF, Otterblad-Olausson  P, Pettersson  BU, Ekbom  A.  The Swedish personal identity number: possibilities and pitfalls in healthcare and medical research.   Eur J Epidemiol. 2009;24(11):659-667. doi:10.1007/s10654-009-9350-y PubMedGoogle ScholarCrossref
25.
Zhao  J, Rostgaard  K, Hjalgrim  H, Edgren  G.  The Swedish Scandinavian donations and transfusions database (SCANDAT3-S)—50 years of donor and recipient follow-up.   Transfusion. 2020;60(12):3019-3027. doi:10.1111/trf.16027 PubMedGoogle ScholarCrossref
26.
Elixhauser  A, Steiner  C, Harris  DR, Coffey  RM.  Comorbidity measures for use with administrative data.   Med Care. 1998;36(1):8-27. doi:10.1097/00005650-199801000-00004 PubMedGoogle ScholarCrossref
27.
Sharma  N, Schwendimann  R, Endrich  O, Ausserhofer  D, Simon  M.  Comparing Charlson and Elixhauser comorbidity indices with different weightings to predict in-hospital mortality: an analysis of national inpatient data.   BMC Health Serv Res. 2021;21(1):13. doi:10.1186/s12913-020-05999-5 PubMedGoogle ScholarCrossref
28.
Auvinen  MK, Zhao  J, Lassén  E,  et al.  Patterns of blood use in Sweden from 2008 to 2017: a nationwide cohort study.   Transfusion. 2020;60(11):2529-2536. doi:10.1111/trf.16092 PubMedGoogle ScholarCrossref
29.
Austin  PC, Stuart  EA.  Moving towards best practice when using inverse probability of treatment weighting (IPTW) using the propensity score to estimate causal treatment effects in observational studies.   Stat Med. 2015;34(28):3661-3679. doi:10.1002/sim.6607 PubMedGoogle ScholarCrossref
30.
Bandara  SN, Kennedy-Hendricks  A, Stuart  EA,  et al.  The effects of the Maryland Medicaid health home waiver on emergency department and inpatient utilization among individuals with serious mental illness.   Gen Hosp Psychiatry. 2020;64:99-104. doi:10.1016/j.genhosppsych.2019.12.004 PubMedGoogle ScholarCrossref
31.
Hernán  MA, Lanoy  E, Costagliola  D, Robins  JM.  Comparison of dynamic treatment regimes via inverse probability weighting.   Basic Clin Pharmacol Toxicol. 2006;98(3):237-242. doi:10.1111/j.1742-7843.2006.pto_329.x PubMedGoogle ScholarCrossref
32.
Cole  SR, Hernán  MA.  Constructing inverse probability weights for marginal structural models.   Am J Epidemiol. 2008;168(6):656-664. doi:10.1093/aje/kwn164 PubMedGoogle ScholarCrossref
33.
Williamson  T, Ravani  P.  Marginal structural models in clinical research: when and how to use them?   Nephrol Dial Transplant. 2017;32(suppl 2):ii84-ii90. doi:10.1093/ndt/gfw341 PubMedGoogle ScholarCrossref
34.
Shinozaki  T, Suzuki  E.  Understanding marginal structural models for time-varying exposures: pitfalls and tips.   J Epidemiol. 2020;30(9):377-389. doi:10.2188/jea.JE20200226 PubMedGoogle ScholarCrossref
35.
Dickerman  BA, García-Albéniz  X, Logan  RW, Denaxas  S, Hernán  MA.  Avoidable flaws in observational analyses: an application to statins and cancer.   Nat Med. 2019;25(10):1601-1606. doi:10.1038/s41591-019-0597-x PubMedGoogle ScholarCrossref
36.
Edgren  G, Hjalgrim  H.  Epidemiology of donors and recipients: lessons from the SCANDAT database.   Transfus Med. 2019;29(suppl 1):6-12. doi:10.1111/tme.12487 PubMedGoogle ScholarCrossref
37.
Pocock  SJ, Conthe  A.  Ever-pregnant female blood donors and mortality risk in male recipients.   JAMA. 2018;319(10):1048-1049. doi:10.1001/jama.2017.21868 PubMedGoogle ScholarCrossref
38.
Hernán  MA, Brumback  B, Robins  JM.  Marginal structural models to estimate the causal effect of zidovudine on the survival of HIV-positive men.   Epidemiology. 2000;11(5):561-570. doi:10.1097/00001648-200009000-00012 PubMedGoogle ScholarCrossref
39.
Hernán  MA, Brumback  BA, Robins  JM.  Estimating the causal effect of zidovudine on CD4 count with a marginal structural model for repeated measures.   Stat Med. 2002;21(12):1689-1709. doi:10.1002/sim.1144 PubMedGoogle ScholarCrossref
40.
Fergusson  DA, Chassé  M, Tinmouth  A,  et al.  Pragmatic, double-blind, randomised trial evaluating the impact of red blood cell donor sex on recipient mortality in an academic hospital population: the innovative Trial Assessing Donor Sex (iTADS) protocol.   BMJ Open. 2021;11(2):e049598.PubMedGoogle Scholar
41.
Chassé  M, Fergusson  DA, Acker  JP,  et al. An Innovative Trial Assessing Donor Sex on Recipient Survival (iTADS). Abstract presented at: American Association of Blood Banks Virtual Annual Meeting; October 18, 2021.
×