The death of a child is arguably one of the most stressful events that a parent can experience. Grief is the natural emotional response to such loss and has been described across human populations (Lindemann, 1944; Bowlby, 1980; Parkes, 2014; Boyden et al., 2014; Einarsdóttir, 2004). Intense or prolonged grief can induce behaviours and physiological states that are detrimental to health and survival. Thus, child loss has been associated with elevated risks of psychiatric symptoms (Kreicbergs et al., 2004), psychiatric hospitalizations (Li et al., 2005) and sickness absence (Wilcox et al., 2015), cardiovascular disease (Li et al., 2002a; Li et al., 2003a), and some cancers (Fang et al., 2011; Huang et al., 2013; Li et al., 2002b). Studies also suggest that parents, particularly mothers, who lose a child face increased risk of death (Li et al., 2002a; Li et al., 2003a; Levav et al., 1988; Li et al., 2003b; Qin and Mortensen, 2003; Agerbo, 2005; Rostila et al., 2012; Rostila et al., 2015; Espinosa and Evans, 2013; Harper et al., 2011; Chen et al., 2012; Schorr et al., 2016), both from suicide and natural causes. The aforementioned studies were all conducted in affluent developed populations, where infant mortality rates are low and life expectancy is high. In contemporary developing nations, infant mortality rates (number of deaths during first year per 1000 live births) >50 are common, whereas in developed nations they are usually lower than eight (United Nations, 2015). Parent-child attachment in pre-industrial and developing populations has been a matter of debate, with some historians arguing for parental reluctance to invest emotions in children where mortality rates are high, while others argue that parents' attachment and grief after child loss is similar across human populations, past and present (Magnússon, 2010; Woods, 2003). One way to shed light on this issue is to study the impact of child loss on parent mortality in the same population before and after a demographic transition that dramatically altered infant mortality rates.

Here, we examine the association between child loss and parent mortality in Iceland from 1817 until 2015, a period that spans a drastic transition from a poor agricultural society to an affluent industrial society. One marker of this drastic change is that infant mortality rates in Iceland (per 1000 births) dropped from one of the highest known values in history, 238 in 1881–1885 (Jónsson and Magnússon, 1997), to one of the lowest, 2.07 in 2005–2010 (Statistics Iceland, 2019). We used a comprehensive genealogical database to assess whether the relative rate of parent mortality after child loss changed during this transition, and to determine whether this rate was further affected by factors such as parent’s sex and age at loss, the number of children as well as the age of the lost child. In line with some previous studies (Li et al., 2003b), we used parent mortality as a measure of the cumulative health response to the emotional and physiological strain of child loss. While mortality obviously does not capture the entire range of such reactions, it has the advantage of being reliably recorded at different times in a genealogical database.

We followed all Icelandic parents, born 1800 or later, from the birth of their first child (earliest from 1817) until death or to the end of 2015 (see further in Materials and methods). From 1817 to 2015, 64 044 parents lost at least one child (Supplementary file 1 - Table 1). The proportion of parents who experienced child loss decreased from 61.1% among parents born 1800–1880 to 5.2% among parents born 1931–1996. To adjust for potential confounding by socioeconomic status or other familial factors that might affect the risk of mortality, our primary analyses were focused on the 47,711 parents who lost a child and had at least one sibling who did not lose a child (N = 126,342). Compared to their siblings who did not lose a child, parents who lost a child were more likely to be women, have more children, be younger at first childbirth and older at cohort entry, and have a longer follow-up time (p<0.05; Table 1). Similar patterns were observed for mothers and fathers who lost a child, respectively (Supplementary file 1 - Table 2).

Parent mortality rate after child loss across two centuries

We used a conditional Cox proportional hazards model to contrast the mortality rates of parents who lost at least one child to their respective siblings who did not (Figure 1). Throughout the entire study period, we observed an elevated hazard ratio (HR) of overall mortality rate among parents who lost a child (1.04; 95% confidence interval [CI] 1.02–1.06), which was mostly driven by maternal (HR 1.07; 95% CI 1.04–1.10) rather than paternal mortality rates (HR 1.02; 95% CI 0.99–1.05; p-for-difference=0.0096). Figure 1 shows that there are no obvious temporal trends in the hazard ratios, except for an increase in mortality rate after child loss for parents born after 1930 (p-for-difference = 0.002 when compared to parents born 1800–1930). Similar results were obtained from the population-based matched cohort and a subset of this cohort with identifiable siblings, separately, (Figure 1—figure supplement 1), albeit with greater HRs revealing the importance of the sibling-controlled design to adjust for unmeasured familial confounding factors. The HRs changed minimally after exclusion of parent-child pairs with the same date of death (Figure 1—figure supplement 2).

Figure 1 with 2 supplements see all

Download asset Open asset

Child loss and premature maternal mortality rate

As the average life expectancy of Icelanders increased considerably from 1817 to the present, we performed separate analyses of parent mortality before and after the age of 50 years (Figure 2). Before the age of 50, mothers who lost a child experienced increased mortality (HR 1.36; 95% CI 1.27–1.45) during the whole study period with the HRs consistently elevated across all seven birth cohorts (Figure 2a). For mothers, the relative mortality rate increase ranged from 35% (born 1800–1930) to 64% (born after 1930), yet was not statistically significantly different between these two cohorts (p=0.0641). In contrast, no such excess mortality was observed for fathers who lost a child in any of the birth cohorts (HR for entire observation period 1.04; 95% CI 0.97–1.11). Interestingly, the excess maternal mortality rates before the age of 50 after child loss seem limited to loss of a child younger than 18 years, while the elevation in paternal mortality rates are more notable after loss of an adult child (Table 2). The difference between mothers and fathers in excess mortality rate before the age of 50 is particularly marked after loss of infants, representing 59% of child losses during the study period. These results were robust to the exclusion of mothers who died within a week after giving birth (Figure 2—figure supplement 1), making it unlikely that the higher mortality rate in mothers is explained by delivery or postpartum complications that resulted in the deaths of both newborn and mother (Nour, 2008).

Figure 2 with 6 supplements see all

Download asset Open asset

The study period spans two centuries of drastic demographic transition in Iceland, from a poor developing nation with a particularly high infant mortality rate, to a modern affluent nation with one of the lowest rates of infant mortality in history. Although we observed a more modest elevation in overall parent mortality rates after child loss before 1930 compared to thereafter, the most striking finding is the consistent increase in premature maternal deaths in response to child loss throughout the 199 years of observation.

This study leverages a unique nationwide data source of genealogies with virtually complete follow-up over a 199 year period spanning an extreme change in indicators of societal development, for example infant and child mortality rate. Although based on an unusually large set of data, our study has several limitations that should be noted. First, we have no information about the causes of death for children or parents. Such data could be used to help verify the causal relationship between deaths of children and parents and to uncover patterns in the causes of death among parents (e.g. suicide or cardiovascular deaths) after child loss. Some impact of reverse causation, for example the ill-health of parents increasing the mortality risk of their children, cannot be ruled out. However, the health and welfare of children was even more dependent on their parents, particularly fathers, in pre-demographic transition (pre-welfare state times) societies. The relatively stable effect estimates of parent mortality rates after child loss across calendar periods and the consistently lower impact of child loss on paternal mortality argues against a strong influence of reverse causality. Second, we did not have direct information about the socioeconomic status of parents. However, as socioeconomic status tends to cluster in families, our sibling-controlled analysis is likely to adjust for much of the confounding due to this source. Third, it is possible that the fragmented registration of births in earlier times (particularly in the case of infant deaths) resulted in misclassification of some parents who had truly lost a child and, thus, somewhat conservative estimates of excess parent mortality rates in earlier times. Last, in the sibling cohort, parents who lost a child differed somewhat from those that did not lose a child with respect to sex, number of children, age at first birth, age at entry to the cohort, length of follow-up time and age at death. We therefore made every effort to adjust for these covariates in our analyses. We adjusted for age at follow-up and sex in all analyses and presented separate results for the different birth cohorts. We further performed subgroup analyses by age of parents at loss (or matching), parent age at first childbirth, follow-up time and number of children alive at loss, and demonstrated that the association between child loss and parent mortality is observed across the range of these covariates. Indeed, parents who lost a child may have to live longer to experience such an event, however, such survivor bias (if any) due to a selection of survivors to the exposed group would rather yield conservative estimates of parent mortality after a child loss. On the other hand, the siblings used for comparison that were not yet parents at the age as when the index parents lost a child may have a survival advantage, which could yield exaggerated estimates of parent mortality after child loss. However, we observed similar results in a sensitivity analysis that was restricted to unexposed siblings who were already parents at the same age as when the index parents lost a child.

The large increase in the rate of premature maternal mortality after child loss in 19^th century Iceland suggests that mothers formed strong emotional bonds with their children – in spite of very harsh conditions and a strong expectation of child loss. This is consistent with accounts of intense grieving among mothers who lost a child in a contemporary developing society (Einarsdóttir, 2004). In contrast, the paternal mortality response to child loss is much weaker throughout the observation period and only statistically significant for older fathers (i.e. parent mortality after 50 years of age) born after 1930.

What could account for this stable sex difference in parent mortality rate after child loss? One possibility is a difference in the attachment of mothers and fathers to their children, which could result in a greater emotional response of mothers and concomitant psychological and physiological consequences. Throughout the 199-year study period, mothers may have had greater opportunity than fathers to form early emotional attachment to their infants through gestation and breastfeeding. Although historical accounts suggest that breastfeeding may have been rarer in 19th to early 20th century Iceland (Garðarsdóttir, 2005) than in the present, it does not seem to have prevented mothers from forming strong emotional bonds with their children as illustrated by the elevated risks of maternal death after child loss. It is possible that females are born with a tendency to form stronger emotional bonds with children than males. Although our results cannot disentangle the impact of nature and nurture in this context, it is interesting to note that the mortality response of fathers to child loss before the age of 50 has not significantly increased despite the greater contribution of fathers to the postnatal care of children in recent decades (Huerta et al., 2014).

Another possible explanation for our findings is a sex difference in the emotional response to traumatic events in general. Women have been reported to have a greater risk of mental disorders after exposure to emotional trauma than men (Tolin and Foa, 2006). The intensity of such responses may result in a multitude of other mental and physical ailments which in turn can affect the risk of mortality through, for example suicide (Wilcox et al., 2009; Panagioti et al., 2012; Gradus et al., 2010) and cardiovascular disease (Song et al., 2019; Edmondson et al., 2013; Gradus et al., 2015). Previous studies in contemporary developed populations have reported evidence for greater risk of psychiatric hospitalizations (Li et al., 2005) and mortality (Li et al., 2003b) after child loss in mothers than in fathers. Interestingly, an opposite sex difference is seen for mortality rate after spousal loss, where mortality rate elevation is more pronounced among widowers than widows (Moon et al., 2011; Shor et al., 2012). However, this difference may in part be due to a greater dependence of older males on social and emotional support provided by spouses. A similar mechanism is indicated in our results, namely the increased mortality of older parents after child loss in more recent times.

Our results show that even in the poor agricultural population of 19th century Iceland, with one of the highest infant mortality rates measured in humans, there was a considerable impact of child loss on premature maternal mortality rates. This contradicts claims that mothers in poor cultures adapt to high infant mortality rates with reduced emotional investment (as reviewed in Woods, 2003). Indeed, one implication of our findings is that child loss is likely to constitute a major threat to the survival of mothers in societies with high infant mortality rates. Taken together these findings highlight the importance of the mother-child bond and the extensive health threats that arise when it is broken.

The data used in this study are compiled in the Genealogy database (Íslendingabók) at deCODE Genetics and are available for review on site, on request. The use of these data is in accordance with Icelandic law and with permission of the Icelandic Bioethics Committee. Therefore, the authors cannot make the dataset publicly available. On the other hand, interested researchers can obtain data access upon approvals by the Icelandic Bioethics Committee (https://www.vsn.is/en) and by contacting the corresponding authors.

1. 1. Agerbo E

   (2005) Midlife suicide risk, partner's psychiatric illness, spouse and child bereavement by suicide or other modes of death: a gender specific study

   Journal of Epidemiology & Community Health 59:407–412.

   https://doi.org/10.1136/jech.2004.024950

   • Google Scholar
2. 1. Bowlby J

   (1980)

   Attachment and Loss

   Loss, Sadness and Depression, Attachment and Loss, 3, New York, Basic Books.

   • Google Scholar
3. 1. Boyden JY
   2. Kavanaugh K
   3. Issel LM
   4. Eldeirawi K
   5. Meert KL

   (2014) Experiences of african american parents following perinatal or pediatric death: a literature review

   Death Studies 38:374–380.

   https://doi.org/10.1080/07481187.2013.766656

   • PubMed
   • Google Scholar
4. 1. Chen CC
   2. Kuo CC
   3. Wu TN
   4. Yang CY

   (2012) Death of a son is associated with risk of suicide among parous women in Taiwan: a nested case-control study

   Journal of Epidemiology 22:532–536.

   https://doi.org/10.2188/jea.JE20120060

   • PubMed
   • Google Scholar
5. 1. Edmondson D
   2. Kronish IM
   3. Shaffer JA
   4. Falzon L
   5. Burg MM

   (2013) Posttraumatic stress disorder and risk for coronary heart disease: a meta-analytic review

   American Heart Journal 166:806–814.

   https://doi.org/10.1016/j.ahj.2013.07.031

   • PubMed
   • Google Scholar
6. 1. Einarsdóttir J

   (2004)

   Mother Love, Child Death and Poverty in Guinea-Bissau

   Tired of Weeping, Mother Love, Child Death and Poverty in Guinea-Bissau, 2nd edition, Madison, Wisconsin, University of Wisconsin Press.

   • Google Scholar
7. 1. Espinosa J
   2. Evans WN

   (2013) Maternal bereavement: the heightened mortality of mothers after the death of a child

   Economics & Human Biology 11:371–381.

   https://doi.org/10.1016/j.ehb.2012.06.002

   • PubMed
   • Google Scholar
8. 1. Fang F
   2. Fall K
   3. Sparén P
   4. Adami HO
   5. Valdimarsdóttir HB
   6. Lambe M
   7. Valdimarsdóttir U

   (2011) Risk of infection-related cancers after the loss of a child: a follow-up study in Sweden

   Cancer Research 71:116–122.

   https://doi.org/10.1158/0008-5472.CAN-10-0470

   • PubMed
   • Google Scholar
9. 1. Garðarsdóttir Ó

   (2005)

   Brjóstagjöf og áhrif hennar á lífslíkur ungbarna á Íslandi 1910-1925

   Ljósmæðrablaðið 83:19–27.

   • Google Scholar
10. 1. Gradus JL
    2. Qin P
    3. Lincoln AK
    4. Miller M
    5. Lawler E
    6. Sørensen HT
    7. Lash TL

    (2010) Posttraumatic stress disorder and completed suicide

    American Journal of Epidemiology 171:721–727.

    https://doi.org/10.1093/aje/kwp456

    • PubMed
    • Google Scholar
11. 1. Gradus JL
    2. Farkas DK
    3. Svensson E
    4. Ehrenstein V
    5. Lash TL
    6. Milstein A
    7. Adler N
    8. Sørensen HT

    (2015) Associations between stress disorders and cardiovascular disease events in the danish population

    BMJ Open 5:e009334.

    https://doi.org/10.1136/bmjopen-2015-009334

    • PubMed
    • Google Scholar
12. 1. Gudbjartsson DF
    2. Helgason H
    3. Gudjonsson SA
    4. Zink F
    5. Oddson A
    6. Gylfason A
    7. Besenbacher S
    8. Magnusson G
    9. Halldorsson BV
    10. Hjartarson E
    11. Sigurdsson GT
    12. Stacey SN
    13. Frigge ML
    14. Holm H
    15. Saemundsdottir J
    16. Helgadottir HT
    17. Johannsdottir H
    18. Sigfusson G
    19. Thorgeirsson G
    20. Sverrisson JT
    21. Gretarsdottir S
    22. Walters GB
    23. Rafnar T
    24. Thjodleifsson B
    25. Bjornsson ES
    26. Olafsson S
    27. Thorarinsdottir H
    28. Steingrimsdottir T
    29. Gudmundsdottir TS
    30. Theodors A
    31. Jonasson JG
    32. Sigurdsson A
    33. Bjornsdottir G
    34. Jonsson JJ
    35. Thorarensen O
    36. Ludvigsson P
    37. Gudbjartsson H
    38. Eyjolfsson GI
    39. Sigurdardottir O
    40. Olafsson I
    41. Arnar DO
    42. Magnusson OT
    43. Kong A
    44. Masson G
    45. Thorsteinsdottir U
    46. Helgason A
    47. Sulem P
    48. Stefansson K

    (2015) Large-scale whole-genome sequencing of the icelandic population

    Nature Genetics 47:435–444.

    https://doi.org/10.1038/ng.3247

    • PubMed
    • Google Scholar
13. 1. Halldorsson BV
    2. Palsson G
    3. Stefansson OA
    4. Jonsson H
    5. Hardarson MT
    6. Eggertsson HP
    7. Gunnarsson B
    8. Oddsson A
    9. Halldorsson GH
    10. Zink F
    11. Gudjonsson SA
    12. Frigge ML
    13. Thorleifsson G
    14. Sigurdsson A
    15. Stacey SN
    16. Sulem P
    17. Masson G
    18. Helgason A
    19. Gudbjartsson DF
    20. Thorsteinsdottir U
    21. Stefansson K

    (2019) Characterizing mutagenic effects of recombination through a sequence-level genetic map

    Science 363:eaau1043.

    https://doi.org/10.1126/science.aau1043

    • PubMed
    • Google Scholar
14. 1. Harper M
    2. O'Connor RC
    3. O'Carroll RE

    (2011) Increased mortality in parents bereaved in the first year of their child's life

    BMJ Supportive & Palliative Care 1:306–309.

    https://doi.org/10.1136/bmjspcare-2011-000025

    • PubMed
    • Google Scholar
15. 1. Helgason A
    2. Einarsson AW
    3. Guðmundsdóttir VB
    4. Sigurðsson Á
    5. Gunnarsdóttir ED
    6. Jagadeesan A
    7. Ebenesersdóttir SS
    8. Kong A
    9. Stefánsson K

    (2015) The Y-chromosome point mutation rate in humans

    Nature Genetics 47:453–457.

    https://doi.org/10.1038/ng.3171

    • PubMed
    • Google Scholar
16. 1. Huang J
    2. Valdimarsdóttir U
    3. Fall K
    4. Ye W
    5. Fang F

    (2013) Pancreatic Cancer risk after loss of a child: a register-based study in Sweden during 1991-2009

    American Journal of Epidemiology 178:582–589.

    https://doi.org/10.1093/aje/kwt045

    • PubMed
    • Google Scholar
17. 1. Huerta MC
    2. Adema W
    3. Baxter J
    4. Han WJ
    5. Lausten M
    6. Lee R
    7. Waldfogel J

    (2014) Fathers' Leave and fathers' Involvement: evidence from four OECD countries

    European Journal of Social Security 16:308–346.

    https://doi.org/10.1177/138826271401600403

    • PubMed
    • Google Scholar
18. Book

    1. Jónsson G
    2. Magnússon MS

    (1997)

    Hagskinna: Icelandic Historical Statistics

    Reykjavik, Iceland: Hagstofa Íslands.

    • Google Scholar
19. 1. Kreicbergs U
    2. Valdimarsdóttir U
    3. Onelöv E
    4. Henter JI
    5. Steineck G

    (2004) Anxiety and depression in parents 4-9 years after the loss of a child owing to a malignancy: a population-based follow-up

    Psychological Medicine 34:1431–1441.

    https://doi.org/10.1017/S0033291704002740

    • PubMed
    • Google Scholar
20. 1. Levav I
    2. Friedlander Y
    3. Kark JD
    4. Peritz E

    (1988) An epidemiologic study of mortality among bereaved parents

    New England Journal of Medicine 319:457–461.

    https://doi.org/10.1056/NEJM198808253190801

    • PubMed
    • Google Scholar
21. 1. Li J
    2. Hansen D
    3. Mortensen PB
    4. Olsen J

    (2002a) Myocardial infarction in parents who lost a child: a nationwide prospective cohort study in Denmark

    Circulation 106:1634–1639.

    https://doi.org/10.1161/01.cir.0000031569.45667.58

    • PubMed
    • Google Scholar
22. 1. Li J
    2. Johansen C
    3. Hansen D
    4. Olsen J

    (2002b) Cancer incidence in parents who lost a child: a nationwide study in Denmark

    Cancer 95:2237–2242.

    https://doi.org/10.1002/cncr.10943

    • PubMed
    • Google Scholar
23. 1. Li J
    2. Johnsen SP
    3. Olsen J

    (2003a) Stroke in parents who lost a child: a nationwide follow-up study in Denmark

    Neuroepidemiology 22:211–216.

    https://doi.org/10.1159/000069898

    • PubMed
    • Google Scholar
24. 1. Li J
    2. Precht DH
    3. Mortensen PB
    4. Olsen J

    (2003b) Mortality in parents after death of a child in Denmark: a nationwide follow-up study

    The Lancet 361:363–367.

    https://doi.org/10.1016/S0140-6736(03)12387-2

    • Google Scholar
25. 1. Li J
    2. Laursen TM
    3. Precht DH
    4. Olsen J
    5. Mortensen PB

    (2005) Hospitalization for mental illness among parents after the death of a child

    New England Journal of Medicine 352:1190–1196.

    https://doi.org/10.1056/NEJMoa033160

    • PubMed
    • Google Scholar
26. 1. Lindemann E

    (1944) Symptomatology and management of acute grief

    American Journal of Psychiatry 101:141–148.

    https://doi.org/10.1176/ajp.101.2.141

    • Google Scholar
27. 1. Magnússon SG

    (2010)

    A Social History of Iceland

    99–111,  Wasteland with Words, A Social History of Iceland, London, Reaktion Books.

    • Google Scholar
28. 1. Moon JR
    2. Kondo N
    3. Glymour MM
    4. Subramanian SV

    (2011) Widowhood and mortality: a meta-analysis

    PLOS ONE 6:e23465.

    https://doi.org/10.1371/journal.pone.0023465

    • PubMed
    • Google Scholar
29. 1. Nour NM

    (2008)

    An introduction to maternal mortality

    Reviews in Obstetrics & Gynecology 1:77–81.

    • PubMed
    • Google Scholar
30. 1. Panagioti M
    2. Gooding PA
    3. Tarrier N

    (2012) A meta-analysis of the association between posttraumatic stress disorder and suicidality: the role of comorbid depression

    Comprehensive Psychiatry 53:915–930.

    https://doi.org/10.1016/j.comppsych.2012.02.009

    • PubMed
    • Google Scholar
31. Book

    1. Parkes C

    (2014)

    The Price of Love: The Selected Works of Colin Murray Parkes

    London: Routledge.

    • Google Scholar
32. 1. Qin P
    2. Mortensen PB

    (2003) The impact of parental status on the risk of completed suicide

    Archives of General Psychiatry 60:797–802.

    https://doi.org/10.1001/archpsyc.60.8.797

    • PubMed
    • Google Scholar
33. 1. Rostila M
    2. Saarela J
    3. Kawachi I

    (2012) Mortality in parents following the death of a child: a nationwide follow-up study from Sweden

    Journal of Epidemiology and Community Health 66:927–933.

    https://doi.org/10.1136/jech-2011-200339

    • Google Scholar
34. 1. Rostila M
    2. Saarela J
    3. Kawachi I
    4. Hjern A

    (2015) Testing the anniversary reaction: causal effects of bereavement in a nationwide follow-up study from Sweden

    European Journal of Epidemiology 30:239–247.

    https://doi.org/10.1007/s10654-015-9989-5

    • PubMed
    • Google Scholar
35. 1. Schorr L
    2. Burger A
    3. Hochner H
    4. Calderon R
    5. Manor O
    6. Friedlander Y
    7. Lawrence GM
    8. Paltiel O

    (2016) Mortality, Cancer incidence, and survival in parents after bereavement

    Annals of Epidemiology 26:115–121.

    https://doi.org/10.1016/j.annepidem.2015.12.008

    • PubMed
    • Google Scholar
36. 1. Shor E
    2. Roelfs DJ
    3. Curreli M
    4. Clemow L
    5. Burg MM
    6. Schwartz JE

    (2012) Widowhood and mortality: a meta-analysis and meta-regression

    Demography 49:575–606.

    https://doi.org/10.1007/s13524-012-0096-x

    • PubMed
    • Google Scholar
37. 1. Song H
    2. Fang F
    3. Arnberg FK
    4. Mataix-Cols D
    5. Fernández de la Cruz L
    6. Almqvist C
    7. Fall K
    8. Lichtenstein P
    9. Thorgeirsson G
    10. Valdimarsdóttir UA

    (2019) Stress related disorders and risk of cardiovascular disease: population based, sibling controlled cohort study

    BMJ 365:l1255.

    https://doi.org/10.1136/bmj.l1255

    • PubMed
    • Google Scholar
38. Website

    1. Statistics Iceland

    (2019) Infant mortality and late fetal deaths 1951-2018

    Accessed October 30, 2019.

    https://px.hagstofa.is/pxen/pxweb/en/Ibuar/Ibuar__Faeddirdanir__danir__danir/MAN05321.px
39. 1. Sun JX
    2. Helgason A
    3. Masson G
    4. Ebenesersdóttir SS
    5. Li H
    6. Mallick S
    7. Gnerre S
    8. Patterson N
    9. Kong A
    10. Reich D
    11. Stefansson K

    (2012) A direct characterization of human mutation based on microsatellites

    Nature Genetics 44:1161–1165.

    https://doi.org/10.1038/ng.2398

    • PubMed
    • Google Scholar
40. 1. Tolin DF
    2. Foa EB

    (2006) Sex differences in trauma and posttraumatic stress disorder: a quantitative review of 25 years of research

    Psychological Bulletin 132:959–992.

    https://doi.org/10.1037/0033-2909.132.6.959

    • PubMed
    • Google Scholar
41. Book

    1. United Nations

    (2015)

    Statistical Yearbook 2012, Fifty-Seventh Issue

    New York, United States: United Nations.

    • Google Scholar
42. 1. Wilcox HC
    2. Storr CL
    3. Breslau N

    (2009) Posttraumatic stress disorder and suicide attempts in a community sample of urban american young adults

    Archives of General Psychiatry 66:305–311.

    https://doi.org/10.1001/archgenpsychiatry.2008.557

    • PubMed
    • Google Scholar
43. 1. Wilcox HC
    2. Mittendorfer-Rutz E
    3. Kjeldgård L
    4. Alexanderson K
    5. Runeson B

    (2015) Functional impairment due to bereavement after the death of adolescent or young adult offspring in a national population study of 1,051,515 parents

    Social Psychiatry and Psychiatric Epidemiology 50:1249–1256.

    https://doi.org/10.1007/s00127-014-0997-7

    • PubMed
    • Google Scholar
44. 1. Woods R

    (2003) Did montaigne love his children? demography and the hypothesis of parental indifference

    Journal of Interdisciplinary History 33:421–442.

    https://doi.org/10.1162/002219502320815181

    • Google Scholar

Acceptance summary:

This is an interesting and valuable analysis that uses the historic registry data in Iceland to ask whether losing a child increases the risk of death for the parents. As the authors state, several papers have shown loss of a child is associated – in observational studies – with a higher risk of parental death, particularly for mothers. However, the novelty in this study is using a sibling comparison design over a 200-year time frame with a large sample size. Hence they are able to answer the question "is the phenomenon of interest – increased risk of parental death following the death of a child – particular to industrialized human populations in which child death is quite rare." The answer from this paper is clearly 'no' – even in a context in which child death affects more than half of parents, the authors detect a raised risk of death among mothers and (in a more specific context) fathers.

Decision letter after peer review:

Thank you for submitting your article "The mother's risk of premature death after child loss across two centuries" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by M. Dawn Teare as the Reviewing Editor and Eduardo Franco as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Susan C Alberts (Reviewer #1).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Essential revisions:

1) There are some puzzling features of this study that it would be helpful to understand better. Parents who lost a child and had more children, were younger at first birth, were more likely to be women and had longer follow-up that those who did not lose a child (Table 1). Was the difference in% women for parents with and without loss significant? If so, please add to the text. More generally, please explain these differences between the parents with and without loss, and consider whether the differences are an artefact of the design/implementation that might bias the estimates.

2) The difference in age between parents with loss and their siblings could be because follow-up time started when the child died for the parents with loss, and started for the sibling when the sibling became a parent. As such, for the sibling the time between their nephew/niece dying and the birth of their first child might be "immortal" time (1) that gives the sibling a survival advantage. For example, if the parent and their childless sibling both died a week after the child died, then the sibling would not be eligible for the study, making the surviving eligible siblings more strongly selected survivors than the parents with loss. Please consider whether such an issue could affect the estimates and amend the method or add a limitation explaining this point.

3) Please explain the adjustment for sex in the sibling design in more detail. Was the sibling design adjusted for sex composition of the sibling pairs, i.e., both female, both male, etc., because men and women have children at different ages, giving a different exposure time till age 50 years?

4) Given the above concerns about the handling of exposure time, would it be helpful to add a twin comparison possibly with exposure time handled in a different way?

5) Number of children could be a confounder (common cause of exposure and outcome) or a mediator, because a higher number of children affects parental survival. Obviously, the more children in a family the greater risk of a child dying. Alternatively, parents might have a 'replacement' child after a death. Would it be possible to distinguish between these possibilities and conduct the analyses accordingly? Would a marginal structural model help? Alternatively, it would be helpful to see estimates adjusted for number of children in the family.

6) The authors allude to an important point in the Discussion when they say "…we have no information about the causes of death for children or parents. Such data could be used to help verify the causal relationship between deaths of children and parents and to uncover patterns in the causes of death among parents (e.g. suicide or cardiovascular deaths) after child loss." This seems a highly salient point: Indeed, a causal relationship in the other direction is quite possible, and mutual causal arrows are possible. That is, child death may be greatly increased in cases where parents are unwell or in poor condition – implying a causal arrow from parent health to child death. Conversely, the risk of parental death after child death may be particularly high for parents in poor health. The authors need to lay out these alternative causal pathways more explicitly. Is there any data that could provide insight into these questions?

7) It would be good to know a bit more about the completeness (or lack thereof) of the dataset and whether the authors used any specific techniques to deal with this incompleteness. For instance, The authors make the statement "it is possible that the fragmented registration of births in earlier times (particularly in the case of infant deaths) resulted in misclassification of some parents who had truly lost a child and, thus, somewhat conservative estimates of excess parent mortality rates in earlier times." And "The database contains an almost complete record of the ancestors of contemporary Icelanders back to around 1650 and, after 1880, is based on death records from a geographically comprehensive set of parish records." It is difficult to get a full picture of dataset completeness from this relatively scanty information. What proportion of records are incomplete? A table with information about various types of missing data would be advisable.

8) It would be helpful if the discussion focused more on the strengths and weaknesses of the study.

[Editors' note: further revisions were requested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "The mother's risk of premature death after child loss across two centuries" for further consideration at eLife. Your revised article has been favorably evaluated by a Senior Editor and a Reviewing Editor, and two reviewers.

The manuscript has been improved but there is one remaining issue that needs to be addressed before acceptance. Reviewer #2 still has concerns about survivor bias. Comparing people who survived at least one child with those who did not may result in selection bias. Adjusting for potential confounding variables will not protect against this. Essentially the people who have child loss may have to live longer, which may obscure the effect of losing a child. The authors need to elaborate on this issue. This could be done by explaining why this is not a problem, performing an additional comparison or highlighting the issues under limitations. We transcribe below the reviewers' comments.

Reviewer #1:

The authors have thoroughly and clearly addressed the concerns and completed the essential revisions. I believe the paper is ready for publication and will make a strong contribution to the literature.

Reviewer #2:

This is an interesting study looking at the relation of offspring death with survival over 200 years in a unique data set.

My previous concern with this study was that the people who had lost a child are different from the comparison group who did not lose a child, which might cause bias, due to some form of survivor bias, i.e., selection bias. Specifically, the parents with loss are more likely to be women, younger at first birth, to have more children and to have longer follow-up. The authors have responded by saying that they addressed these issues by adjustment or stratification, which are means of dealing with confounding not selection bias.

Essentially, the study appears to be comparing people who survived at least one of their children with people who did not survive any of their children. It is easier to survive at least one of your children if you are a woman (because women live longer than men), are younger when the children are born, if you have more children and are followed up for longer which explains why the people with child loss are more likely to be women, younger at first birth, with more children and to have longer follow-up.

I remain concerned as to whether comparing people who survived at least one child with those who did not gives the effect of child loss on survival. Essentially the people who have child loss may have to live longer, which may obscure the effect of losing a child. Of course, it is possible that these issues do not matter, but it would be helpful if the authors could explain why not, or amend accordingly.

Essential revisions:

1) There are some puzzling features of this study that it would be helpful to understand better. Parents who lost a child had more children, were younger at first birth, were more likely to be women and had longer follow-up that those who did not lose a child (Table 1). Was the difference in% women for parents with and without loss significant? If so, please add to the text. More generally, please explain these differences between the parents with and without loss, and consider whether the differences are an artefact of the design/implementation that might bias the estimates.

Given the large sample size, any differences between parents who did vs. did not lose a child is likely to be statistically significant even when the absolute difference is small. We have now added P for difference for all comparisons in Table 1 and revised the text accordingly. The P for difference across parental groups was <0.05 for all covariates, including: parents’ sex, number of children, age at first birth, age at entry to the cohort, length of follow-up time and age at death.

Although we note a statistically significant difference between sibling groups on all covariates we have directly addressed each of these factors in our analyses and, therefore, are confident that they do not account for our results. We adjusted for age at follow-up (as the underlying timescale) and sex in all analyses and present results for each of the studied birth cohorts throughout. In Table 2, we further present the association between child loss and premature parent mortality rate stratified by parental age at loss/matching and follow-up time, and sub-grouped by the number of children alive at loss. These results are not consistent with the notion that our findings are driven by differences between the groups on the abovementioned factors. Although varying somewhat, the association between loss and survival is evident across the range of these covariates. The considerable differences in average follow-up time between parents who did vs. did not lose a child (31.79 vs. 25.45 years, respectively) is likely explained by the parents who did not lose a child at the time of matching but later encountered a child loss and were then censored as end of follow-up in that matching set. This is relatively common in the earlier birth cohorts due to the higher mortality rates for children in that period. This concern is partly alleviated by Figure 2, where the follow-up of parents is restricted to 50 years across all birth cohorts. Moreover, the results in Table 2 do not suggest a considerable variation in parent mortality across extended periods of follow-up time from child loss/matching. This further argues against the possibility that the difference of average follow-up length between two groups may account for our findings.

Modified manuscript text:

Results section: “Compared to their siblings who did not lose a child, parents who lost a child were more likely to be women, have more children, be younger at first birth and older at cohort entry, and have a longer follow-up time (P<0.05; Table 1).”

Discussion section: “Last, in the sibling cohort, parents who lost a child differed somewhat from those that did not lose a child with respect to sex, number of children, age at first birth, age at entry to the cohort, length of follow-up time and age at death. We therefore made every effort to adjust for these covariates in our analyses. We adjusted for age at follow-up and sex in all analyses and presented separate results for the different birth cohorts. We further stratified analyses by age of parents at loss (or matching), follow-up time and number of children alive at loss, and demonstrated that the association between child loss and parent mortality is observed across the range of these covariates.”

Materials and methods: “In this analysis, we used a chi-square test for the categorical and a t-test for the continuous covariates.”

See also Table 1.

2) The difference in age between parents with loss and their siblings could be because follow-up time started when the child died for the parents with loss, and started for the sibling when the sibling became a parent. As such, for the sibling the time between their nephew/niece dying and the birth of their first child might be "immortal" time (1) that gives the sibling a survival advantage. For example, if the parent and their childless sibling both died a week after the child died, then the sibling would not be eligible for the study, making the surviving eligible siblings more strongly selected survivors than the parents with loss. Please consider whether such an issue could affect the estimates and amend the method or add a limitation explaining this point.

We agree with the reviewer that the sibling design needs further clarification in the manuscript. First, in order to be at risk of losing a child, one needs to have one. Thus, our analysis of index persons who lost a child and their siblings only considers follow-up time where both of these groups were parents. We started to follow the siblings who already had a child from the same age as their index-sibling when he/she lost a child. If the siblings did not have any children at the age when the index-persons lost a child, these were not included in the analysis until they subsequently had one, and then only from the age when they became parents. Indeed, we agree with the reviewer that the time between the sibling’s nephew/niece dying and the birth of their first child is “immortal” and was therefore not included in our analysis and should not bias our estimates. Indeed, both groups (individuals with or without child loss) need to survive long enough to become parents to be considered for analysis. The reviewer notes that the parents who lost a child are a bit further into his/her life history than their siblings who are still childless at the time of the loss. However, the population cohort – where even higher hazard ratios are observed – does not have that problem. Please also refer to the results of our sensitivity analysis below (Point #4), restricting the sibling cohort to unexposed siblings born within 5 years of the parents who lost a child.

We have now clarified this important issue in the Materials and methods section.

Modified manuscript text:

“Specifically, siblings who did not lose a child started to contribute to the follow-up at the same age as the parent who lost a child or at a later age when they became parents.”

3) Please explain the adjustment for sex in the sibling design in more detail. Was the sibling design adjusted for sex composition of the sibling pairs, i.e., both female, both male, etc., because men and women have children at different ages, giving a different exposure time till age 50 years?

We did not restrict the sibling sets to all females or all males since that would inevitably have further raised the bar for selection of siblings. Instead, we matched a mother who lost a child to all her eligible siblings with children, including both her sisters and brothers, who had not lost a child at the same age. We adjusted for sex in the sibling-based analysis which accounts for the effect of sex on mortality rate across sibling sets. Moreover, we used attained age as the underlying time scale and therefore the age difference was also controlled for. We have now further clarified the sibling-based analysis in the Materials and methods.

Modified manuscript text:

“Of note, we did not match on parental sex but adjusted for it in all models and thereby accounting for the effect of sex on mortality rate across the sibling groups.”

4) Given the above concerns about the handling of exposure time, would it be helpful to add a twin comparison possibly with exposure time handled in a different way?

Unfortunately, a twin-based comparison would be underpowered due to the low number of pairs. However, we have now conducted a sensitivity analysis restricting to 96,041 (55.2%) siblings born five years before or after the index parents. The results remain largely unchanged, as presented in Figure 2—figure supplement 5.

Modified manuscript text:

Results section: “In a sensitivity analysis, restricting parent-sibling pairs born within five years from the parents who lost a child or adjusting for the number of children alive at time of matching, we observed similar estimates of parental mortality rates after child loss (Figure 2—figure supplements 5-6). The only notable difference was that, in the birth cohort 1931-1996, the impact of child loss on maternal mortality after 50 years of age became statistically significant.”

Materials and methods: “In a sensitivity analyses, we restricted to parent-sibling pairs born within five years from the parents who lost a child (N=96,041, 55.2%) to reduce the potential differences between siblings.”

5) Number of children could be a confounder (common cause of exposure and outcome) or a mediator, because a higher number of children affects parental survival. Obviously, the more children in a family the greater risk of a child dying. Alternatively, parents might have a 'replacement' child after a death. Would it be possible to distinguish between these possibilities and conduct the analyses accordingly? Would a marginal structural model help? Alternatively, it would be helpful to see estimates adjusted for number of children in the family.

We agree with the reviewer that the number of children could be a confounder or a mediator in the association between child loss and parent mortality. This is why we did not adjust for it in our primary analysis. Instead, we performed a subgroup analysis by the number of children alive at the time of loss among parents who lost a child (Table 2), where we did not observe varying effect sizes across subgroups. In accordance with the reviewer’s note, we now present results from a sensitivity analysis, adjusting for the number of children alive at loss/matching. The results are consistent with those of the main analyses, except that maternal mortality after 50 years of age, in the birth cohort 1931-1996, becomes statistically significant (Figure 2—figure supplement 6). Adjusting for or matching on children born after loss (“replacement children”) is challenging, since the parents need to survive long enough to have them. This would introduce a survival bias. We have therefore opted not to adjust for the number of any post-loss/matching children in our analysis.

Modified manuscript text:

Results section:

In a sensitivity analysis, restricting to parent-sibling pairs born within five years from the parents who lost a child or adjusting for the number of children alive at time of matching, we observed similar estimates of parental mortality rates after child loss (Figure 2—figure supplements 5-6). The only notable difference was that, in the birth cohort 1931-1996, the impact of child loss on maternal mortality after 50 years of age became statistically significant.

Materials and methods: To further address the potential effect of number of children, in another sensitivity analysis, we additionally adjusted for number of alive children at the time of loss/matching.

6) The authors allude to an important point in the Discussion when they say "…we have no information about the causes of death for children or parents. Such data could be used to help verify the causal relationship between deaths of children and parents and to uncover patterns in the causes of death among parents (e.g. suicide or cardiovascular deaths) after child loss." This seems a highly salient point: Indeed, a causal relationship in the other direction is quite possible, and mutual causal arrows are possible. That is, child death may be greatly increased in cases where parents are unwell or in poor condition – implying a causal arrow from parent health to child death. Conversely, the risk of parental death after child death may be particularly high for parents in poor health. The authors need to lay out these alternative causal pathways more explicitly. Is there any data that could provide insight into these questions?

This is a valid point that requires further clarification. We have no data on causes of death while previous studies (e.g. Li et al., 2003) suggest that child loss in modern times (after 1980) increases the risk of maternal death, particularly of unnatural causes but also natural causes. We agree with the reviewer that some risk of reverse causation (parent’s ill health increasing risk of child’s death and subsequent death of parent) cannot be ruled out. However, if deteriorating parental health was strongly associated with the risk of child death, we might expect the association between child loss and parent mortality to be particularly evident during the first years of follow-up. However, the relatively stable association between child loss and parent mortality across follow-up times (as shown in Table 2) does not indicate a strong influence of reverse causality in these data. Indeed, families (mothers and children) were extremely financially dependent on fathers, such that if fathers died or were unable to work, then families were often separated and children fostered – which usually meant very hard work and diminished living conditions for these children. Thus, if parents' ill-health was a driving factor, we would expect a stronger association between child loss and father's mortality. However, it is indeed maternal mortality rate, and not paternal mortality rate after child loss, which is consistently elevated across all follow-up periods.

Modified manuscript text:

Discussion: “Some impact of reverse causation, e.g. the ill-health of parents increasing the mortality risk of their children, cannot be ruled out. However, the health and welfare of children was even more dependent on their parents, particularly fathers, in pre-demographic transition (pre-welfare state times) societies. The relatively stable effect estimates of parent mortality rates after child loss across calendar periods and the consistently lower impact of child loss on paternal mortality argues against a strong influence of reverse causality.”

7) It would be good to know a bit more about the completeness (or lack thereof) of the dataset and whether the authors used any specific techniques to deal with this incompleteness. For instance, The authors make the statement "it is possible that the fragmented registration of births in earlier times (particularly in the case of infant deaths) resulted in misclassification of some parents who had truly lost a child and, thus, somewhat conservative estimates of excess parent mortality rates in earlier times." And "The database contains an almost complete record of the ancestors of contemporary Icelanders back to around 1650 and, after 1880, is based on death records from a geographically comprehensive set of parish records." It is difficult to get a full picture of dataset completeness from this relatively scanty information. What proportion of records are incomplete? A table with information about various types of missing data would be advisable.

In terms of linkages between parents and offspring, the completeness of the genealogical database across the entire observation period is 99%. The proportion of individuals without identifiable one or both parents in the database is: 2.4% in individuals born 1800-1880, 0.6% in individuals born 1881-1930, and 1.5% in individuals 1931-1996. Moreover, the database has been extensively evaluated for the accuracy of such linkages using genetic data – (Helgason et al., 2015, Sun et al., 2012, Halldorsson et al., 2019)

The individuals missing from the genealogical database are primarily children born prior to 1900, who died during the first few weeks of life and were not always recorded in censuses or parish records. This would lead to a misclassification of some parents, which in our study are considered as unexposed (parents that have not lost a child), but in fact should have been considered as exposed (parents who lost a child). We note that the impact of such misclassification on our results is in the form of a conservative bias. Thus, if there is a difference in survival between exposed and unexposed parents, such misclassification would only lead us to underestimate this difference and is therefore not a major concern for this study. We now provide this information in the Materials and methods section.

Modified manuscript text:

“In terms of linkages between parents and offspring, the completeness of the genealogical database across the entire observation period is 99%. The database has been extensively evaluated for the accuracy of such linkage using genetic data (40-42).”

8) It would be helpful if the discussion focused more on the strengths and weaknesses of the study.

Based on the points raised above we have now extended the discussion of limitations and moved that part to the beginning of the Discussion section.

[Editors' note: further revisions were requested prior to acceptance, as described below.]

The manuscript has been improved but there is one remaining issue that needs to be addressed before acceptance. Reviewer #2 still has concerns about survivor bias. Comparing people who survived at least one child with those who did not may result in selection bias. Adjusting for potential confounding variables will not protect against this. Essentially the people who have child loss may have to live longer, which may obscure the effect of losing a child. The authors need to elaborate on this issue. This could be done by explaining why this is not a problem, performing an additional comparison or highlighting the issues under limitations. We transcribe below the reviewers' comments.

Thank you for the opportunity for a second resubmission of our paper. Please find below our responses to the issues raised along with the corresponding modifications to the manuscript. We believe that we have now addressed reviewer #2’s main issue, i.e. survivor bias, and we hope you now find our manuscript acceptable for publication in eLife. Please do not hesitate to contact us with any remaining issues.

Reviewer #1:

The authors have thoroughly and clearly addressed the concerns and completed the essential revisions. I believe the paper is ready for publication and will make a strong contribution to the literature.

Thank you for the positive comment on our work and our paper.

Reviewer #2:

This is an interesting study looking at the relation of offspring death with survival over a 200 years in a unique data set.

My previous concern with this study was that the people who had lost a child are different from the comparison group who did not lose a child, which might cause bias, due to some form of survivor bias, i.e., selection bias. Specifically, the parents with loss are more likely to be women, younger at first birth, to have more children and to have longer follow-up. The authors have responded by saying that they addressed these issues by adjustment or stratification, which are means of dealing with confounding not selection bias.

Essentially, the study appears to be comparing people who survived at least one of their children with people who did not survive any of their children. It is easier to survive at least one of your children if you are a woman (because women live longer than men), are younger when the children are born, if you have more children and are followed up for longer which explains why the people with child loss are more likely to be women, younger at first birth, with more children and to have longer follow-up.

I remain concerned as to whether comparing people who survived at least one child with those who did not gives the effect of child loss on survival. Essentially the people who have child loss may have to live longer, which may obscure the effect of losing a child. Of course, it is possible that these issues do not matter, but it would be helpful if the authors could explain why not, or amend accordingly.

The reviewer raises several important issues.

The reviewer remains concerned for a potential bias based on that “people who have child loss may have to live longer” than the comparison cohort of parents who do not experience a child loss. However, such type of bias, if any, caused by a selection of “survivors” to the exposed population would rather yield conservative estimates of the increase in parent mortality after a child loss. Also, as explained in more details below, we carefully designed our study to compare the two groups of parents, with and without a child loss, from the same age (i.e. the age when the loss of the child in one group occurred) whenever possible. Thus, both groups were at the same age at the start of follow-up (i.e., at the time of matching).

Thus, we remain convinced that our study design entails few, if any, possible conditions giving rise to survivor bias as described in the paper cited by the reviewer in the last round of review (2). We followed parents who lose a child starts immediately after child loss and, due to the fact that we use parental age as the underlying time scale, the follow-up of the matched parents not yet exposed to child loss starts at the same age/time point as for the exposed parents yielding identical conditions for ascertaining mortalities among these two compared groups. The sibling controls are followed from the same age as the parent that lost a child, or from the age when becoming a parent if they are not yet parents when their sibling lost a child; the latter, according to the definition above, may give rise to survivor bias. For this reason, we reran our sibling-based comparison, limiting it to unexposed siblings who were already parents at the age as the index sibling who lost a child. The results remained largely similar (if anything, the point estimates indicated even stronger associations) and we have added this sensitivity analysis to the revised manuscript (Figure 2—figure supplement 8).

Thus, we believe that the reviewer’s concern of survivor bias should be eliminated through our study design only considering person time of unexposed siblings who have reached the exact age at the index parent’s loss. In other words, both exposed and unexposed parent’s need to survive through the same age in order to be compared. Therefore, the amount of “immortal time” is similar across these groups. We have now clarified this issue by adding the following statements to the manuscript:

Modified manuscript text:

“We further restricted unexposed siblings to those who were already parents at the age when the index sibling lost a child, which also yielded largely similar results (Figure 2—figure supplement 8).”

“Indeed, parents who lost a child may have to live longer to experience such an event, however, such survivor bias (if any) due to a selection of survivors to the exposed group would rather yield conservative estimates of parent mortality after a child loss. […] However, we observed similar results when we in a sensitivity analysis restricted to unexposed siblings who were already parents at the same age as when the index parents lost a child.”

“To reduce the potential differences between siblings we, in sensitivity analyses restricted to parent-sibling pairs born within five years from the parents who lost a child (N=96,041, 55.2%) as well as to unexposed siblings who were already parents at the same age as when the index parents lost a child (N=135,719, 78.0%).”

In contrast to the abovementioned concern for selection (survivor) bias induced by the study design where certain individuals (or person time) are selected to or excluded from a study (3), we believe that the remaining points raised by the reviewer represent conditions of confounding but not selection bias (4). Our study is indeed a nationwide complete follow-up of all parents during a 199-year period; no individual is excluded from the study on the basis of sex, age at first childbirth or total number of children. As the reviewer correctly points out, these factors may be associated with parental survival and some of them are unequally distributed across the compared parental groups and therefore need to be accounted for in the analysis. We have indeed made extensive efforts to account for these factors in our analysis.

Sex is accounted for in the population analysis by matching on sex, i.e. the rates of mortality are only compared between pairs of exposed and unexposed mothers and then fathers, separately. To preserve statistical power, we did not match on sex in the sibling analysis while controlling for sex in all models. We have now performed additional analysis confined to sister- and then brother-pairs (discordant on child loss) yielding similar estimates in the primary analysis albeit with lower precision. These results have been added to the manuscript (Figure 2—figure supplement 6).

Modified manuscript text:

“In a sensitivity analysis, restricting to parent-sibling pairs born within five years from the parents who lost a child, comparing same-sex siblings, or adjusting for the number of children alive at time of matching, we observed similar estimates of parental mortality rates after child loss (Figure 2—figure supplements 5-7).”

“To further address the potential influences of parent sex and number of children, in another sensitivity analyses, we only compared same-sex siblings and additionally adjusted for number of alive children at the time of loss/matching, respectively.”

Other factors mentioned by the reviewer have appropriately, as potential confounders, been addressed in Table 2 in subgroup analyses. There we e.g. demonstrate that maternal mortality rates below the age of 50 are considerably elevated across parent groups with varying number of children, with the greatest excess observed after loss of an only child. In the last round of revisions, we further performed additional analyses controlling for the number of children at the time of matching revealing virtually identical, if anything stronger, results (Figure 2—figure supplement 7).

In accordance with the reviewer’s comment, we have now added a subgroup analysis of age at first child birth in Table 2, revealing that the association between child loss and maternal mortality rates below the age of 50 is present in every subgroup of parental age at first child birth.

Modified manuscript text:

“We further performed subgroup analyses by the age of parents at loss (or matching), parent age at first childbirth, follow-up time and number of children alive at loss, and demonstrated that the association between child loss and parent mortality was observed across the range of these covariates.”

“Parent mortality before and after age of 50 years was further analyzed in subgroups of age of the deceased child (0, 1-5, 6-17, or ≥18 years old), number of children alive at loss (0, 1-3, or ≥4), sex of the deceased child, parental age at loss (13-30, 31-40, or 41-50 years old), and parental age at first childbirth (13-21, 22-24, 25-27, or ≥28 years old).”

Finally, as we explained in the last round of revisions, the discrepancy between the parental groups in follow-up time is driven by the high incidence of child loss, particularly during 20th century. Our study design mimics a prospective study, where the investigators have at every time point of the observation no information on which parent will live or die, or lose a child. Therefore, parents were randomly selected as matched “control parents” while they had not lost a child and then censored from that control group at the time of child loss (if it occurred). Due to the high incidence of child loss, many parents were later censored from the control group at the time of child loss and assigned to a new set (now as index-individuals) yielding the relatively shorter follow-up time of parents not exposed to child loss relative to the exposed. Thus, the Cox model does not generate estimates of parental mortality rates based on comparison between two distinct groups of parents, exposed and unexposed to child loss, but on the basis of hazards during exposed and unexposed person time.

1) Levesque LE, Hanley JA, Kezouh A, et al. Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ 2010;340:b5087. doi: 10.1136/bmj.b5087 published Online First: 2010/03/17]

2) Lévesque LE et al. Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ, 2010;340:b5087.

3) Rothman KJ. Epidemiology: an introduction. Oxford university press; 2012: p126.

4) Newton J. Confounding is mistakenly called selection bias. BMJ, 1998. Accessed on https://www.bmj.com/rapid-response/2011/10/27/confounding-mistakenly-called-selection-bias

Author details

1. Unnur A Valdimarsdóttir

   1. Center of Public Health Sciences, Faculty of Medicine, University of Iceland, Reykjavik, Iceland
   2. Department of Epidemiology, Harvard TH Chan School of Public Health, Boston, United States
   3. Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Solna, Sweden

   Contribution

   Conceptualization, Supervision, Funding acquisition, Methodology, Writing—original draft, Writing—review and editing

   Contributed equally with

   Donghao Lu

   For correspondence

   unnurav@hi.is

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-5382-946X

2. Donghao Lu

   1. Center of Public Health Sciences, Faculty of Medicine, University of Iceland, Reykjavik, Iceland
   2. Department of Epidemiology, Harvard TH Chan School of Public Health, Boston, United States
   3. Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Solna, Sweden
   4. Channing Division of Network Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, United States

   Contribution

   Formal analysis, Methodology, Writing—original draft, Writing—review and editing

   Contributed equally with

   Unnur A Valdimarsdóttir

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-4186-8661

3. Sigrún H Lund

   deCODE Genetics, Reykjavik, Iceland

   Contribution

   Formal analysis, Writing—review and editing

   Competing interests

   No competing interests declared

4. Katja Fall

   1. Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Solna, Sweden
   2. Clinical Epidemiology and Biostatistics, School of Medical Sciences, Örebro University, Örebro, Sweden

   Contribution

   Writing—review and editing

   Competing interests

   No competing interests declared

5. Fang Fang

   Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Solna, Sweden

   Contribution

   Writing—review and editing

   Competing interests

   No competing interests declared

6. Þórður Kristjánsson

   deCODE Genetics, Reykjavik, Iceland

   Contribution

   Project administration, Writing—review and editing

   Competing interests

   No competing interests declared

7. Daníel Guðbjartsson

   1. deCODE Genetics, Reykjavik, Iceland
   2. School of Engineering and Natural, Sciences, University of Iceland, Reykjavik, Iceland

   Contribution

   Supervision, Methodology, Writing—review and editing

   Competing interests

   No competing interests declared

8. Agnar Helgason

   1. deCODE Genetics, Reykjavik, Iceland
   2. Department of Anthropology, University of Iceland, Reykjavik, Iceland

   Contribution

   Conceptualization, Data curation, Supervision, Writing—original draft, Writing—review and editing

   Contributed equally with

   Kári Stefánsson

   Competing interests

   No competing interests declared

9. Kári Stefánsson

   1. deCODE Genetics, Reykjavik, Iceland
   2. Faculty of Medicine, School of Health Sciences, University of Iceland, Reykjavik, Iceland

   Contribution

   Conceptualization, Data curation, Writing—review and editing

   Contributed equally with

   Agnar Helgason

   For correspondence

   kstefans@decode.is

   Competing interests

   No competing interests declared

• 8,339

  Page views

• 142

  Downloads

• 6

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.