Key PointsQuestion
Is early-life mental health associated with biomarkers in midlife and premature mortality?
Findings
In a birth cohort study of 17 415 participants, early-life mental health was associated with less favorable levels of biomarkers 28 years later and elevated risk of premature mortality. The associations with biomarkers in midlife were more commonly seen in female individuals.
Meaning
These findings, if causal and generalizable to younger cohorts, may have implications for overall population health and age-related disease as effective interventions on early-life mental health have the potential to shift the distribution of risk and improve population health.
Importance
Early-life mental health is known to be associated with socioeconomic adversity and psychological distress in adulthood, but less is known about potential associations with biomarkers and mortality.
Objective
To investigate the association between early-life mental health trajectories with biomarkers in midlife and premature mortality.
Design, Setting, and Participants
This study used data from the British National Child Development Study, a population-based birth cohort. The initial sample of 17 415 individuals consisted of all infants born in Great Britain in a single week in 1958. Analysis began Feburary 2017 and ended May 2020.
Main Outcomes and Measures
Biomarkers collected at age 44 to 45 years were fibrinogen, C-reactive protein, glycated hemoglobin, high- and low-density lipoprotein cholesterol, forced expiratory volume, blood pressure, and waist-to-hip ratio. Information on all-cause mortality was available up to age 58 years and cause-specific mortality up to age 50 years.
Results
Biomarkers were analyzed from 9377 participants (age, 44-45 years, 4712 female [50.3%]) and mortality data from 15 067 participants (age, 58 years; 7379 female [49.0%]). A 4-group longitudinal typology of early-life conduct problems and affective symptoms was identified: (1) stable low, (2) teacher-identified adolescent onset, (3) moderate, and (4) stable high. Compared with the stable-low group, the stable-high (B, 2.308; 95% CI, 0.213-4.404) and teacher-identified adolescent-onset (B, 2.114; 95% CI, 0.725-3.503) groups had less favorable levels of fibrinogen in middle age. Effect modification was observed by sex (P < .005) such that compared with the stable-low group, differences in high-density lipoprotein and abdominal obesity were only observed in female individuals. The stable-high and teacher-identified adolescent-onset groups had elevated risk for all-cause mortality (hazard ratio, 1.878; 95% CI, 1.501-2.350 and hazard ratio, 1.412; 95% CI, 1.174-1.698). Psychopathology-associated mortality was higher in both groups but unintentional injuries–associated mortality only in the stable-high group.
Conclusions and Relevance
Experiencing affective symptoms and conduct problems in childhood and adolescence is associated with less favorable levels of biomarkers 28 years later and elevated risk of premature mortality. These findings, if causal and generalizable to younger cohorts, may imply that effective interventions on early-life mental health have the potential to shift the distribution of risk and improve population health.
Mental health problems in childhood and adolescence are known to be associated with an array of unfavorable outcomes in later life, including psychological distress,1 low educational attainment, unemployment, unstable family formation, and criminal offending.2,3 There are also emerging links with premature all-cause mortality4,5 but a paucity of evidence concerning the prospective association with physical health in adulthood. In the few studies that have investigated this link,6,7 the developmental perspective in the emergence of mental health symptomatology in childhood and adolescence8 has been neglected.
There are various plausible mechanisms of action through which early-life mental health may be linked with midlife health and premature mortality. For example, early-life mental health has been shown to be associated with lower levels of physical activity and a greater prevalence of smoking, alcohol use,9-11 and socioeconomic adversity in adulthood,3,11,12 while specific links to mortality include suicide, drug overdose, unintentional injuries, and homicide.13 We capitalize on the availability of 3 assessments of mental health spanning childhood and adolescence in a population-based prospective birth cohort to characterize trajectories of early-life affective symptoms and conduct problems and to investigate their association with biomarkers in midlife and premature mortality.
Participants in the British National Child Development Study14 have been resurveyed on 10 occasions since birth. The initial sample of 17 415 individuals, consisting of all infants born in Great Britain in a single week in 1958 provide prospective data on social, biological, physical, and psychological phenotypes. In 2002 and 2003, when cohort members were aged 44 to 45 years, a biomedical survey was conducted in 9377 respondents. Ethical approval for this study was obtained by the University College London Institute of Education Research Ethics Committee. Written consent for all measures and samples was collected during the interview.
Conduct problems and affective symptoms in childhood and adolescence were assessed using the modified version of the Rutter Child Scale A.15 This version of the scale was completed by mothers of participants at ages 7 and 11 years and by mothers and teachers at age 16 years. Mother and teacher reports were used to capture symptoms at home and school, as is well known that maternal and teacher reports are weakly correlated and triangulating information from multiple informants may bring unique insights from different environments into children’s behavior.16,17 More information on the Rutter Child Scale is available in eTable 1 and eAppendix 4 in the Supplement.
We used biomarkers as collected at age 44 to 45 years: fibrinogen, a marker of inflammation and cardiovascular disease; C-reactive protein (CRP), an indicator of inflammation and cardiovascular disease; glycated hemoglobin, an index of glucose metabolism throughout the previous 30 to 90 days, which is used as a marker of diabetes; high-density lipoprotein (HDL) and low-density lipoprotein cholesterol as markers of cardiovascular risk; and high blood pressure (3 measures of systolic and diastolic blood pressure were taken). The mean of valid readings was used, and an individual was recorded as having high blood pressure if the mean value was above 140/90 mm Hg. Participants who reported use of blood pressure–related medication were classified as having high blood pressure. Waist and hip circumferences were measured and the ratio of waist over hip was used as a measure of abdominal obesity (≥0.94 in male individuals and ≥0.85 in female individuals18,19). Results with continuous versions of both dichotomized, for comparison purposes with other studies, outcomes were similar to those presented here (eAppendix 11 and eFigures 5-6 in the Supplement). To assess respiratory function, we used forced expiratory volume, with the highest measurement used. Forced expiratory volume is a measure of how much air a person can exhale during forced breath during the first second. Details on the measurement procedures and equipment are available on the biomedical sweep technical report.20 With the exception of glycated hemoglobin, which represents the percentage of total hemoglobin, all other continuous biomarkers were transformed with the natural logarithm and multiplied by 100, so all ordinary least squares regression coefficients can be interpreted as percentage differences in means.21 Appropriate corrections for medication use were also applied22,23 (eAppendix 3 and eTable 3 in the Supplement).
Mortality was recorded by National Health Service digital notifications combined with information during fieldwork and the address database held at the Centre for Longitudinal Studies (eAppendix 2 in the Supplement). All-cause mortality data were available up to age 58 years and cause-specific mortality up to age 50 years, both with information on year of death. We analyzed 3 causes of death that have been hypothesized to be associated with early-life mental health5,13,24: (1) unintentional injuries (transport-related injuries and other external causes of unintentional injury); (2) deaths due to drug and alcohol poisonings, suicide, chronic liver diseases, and cirrhosis, collectively known as death of despair25; and (3) deaths related to the circulatory system. Further information on the classification of deaths and results with suicide risk as a stand-alone outcome are available in eAppendix 3 and eAppendix 7 in the Supplement.
We included a rich set of early-life potential confounders such as birth characteristics, parental characteristics, indicators of early-life socioeconomic position, and characteristics of the cohort member from birth to age 7 years. Detailed information about potential confounders is presented in eAppendix 1 and eTable 2 in the Supplement.
We derived separate latent summaries of conduct problems and affective symptoms at ages 7, 11, and 16 years by modeling the probability of response to the Rutter Child Scale items with a 2-parameter probit latent variable measurement model.26,27 The conduct problems and affective symptoms scores were entered in a latent profile analysis28 to derive a longitudinal typology of mental health from ages 7 to 16 years that was used as a predictor of biomarkers in midlife and mortality in models that included confounders up to age 7 years (eAppendix 4 in the Supplement). For continuous outcomes, we estimated ordinary least squares regression and for binary outcomes (abdominal obesity and high blood pressure), a log binomial model with robust standard errors that returns risk ratios as both outcomes were not rare (>20%) to avoid bias due to noncollapsibility of the odds ratio.29 The association between early-life mental health trajectories and mortality was analyzed with the standard Cox regression model. Considering that female individuals live longer than male individuals and that female individuals are on average less healthy, we investigated modification of all hypothesized associations by sex. We used Stata software versions 15 and 16 (StataCorp) and Mplus software, version 8 (Muthén and Muthén). Two-sided P values with varying levels of significance due to corrections for multiple testing were used. Analysis began in February 2017 and ended in May 2020.
We used multiple imputation with chained equations30 and generated 50 imputed data sets to reduce bias and maintain power. More information on our missing data analysis strategy is available in eAppendix 6 in the Supplement.
Early-Life Mental Health Trajectories
A 4-group solution was selected based on various criteria and quality of allocation/misclassification indices (eAppendix 4 and eTable 4 in the Supplement). Descriptive statistics of all summary variables are presented in Table 1. The 4 longitudinal groups included (1) the stable-low group (5779 [46.3%]), which was characterized by absence of conduct problems and affective symptoms at all ages; (2) the teacher-identified adolescent-onset (TIAO) group (2236 [17.9%]) (absence of affective symptoms and conduct problems at ages 7 and 11 years, followed by sharp increase in levels of both dimensions, but predominantly conduct problems, at age 16 years as reported by teachers); (3) the moderate group (3472 [27.8%]), which experienced moderate levels of conduct problems and affective symptoms in childhood (ages 7 and 11 years), followed by a decline, which signifies improvement, in both dimensions at age 16 years as reported by teachers but a mild increase in both dimensions reported by mothers; and (4) the stable-high group (994 [8%]), in which persistent experience of both conduct problems and affective symptoms from age 7 to 16 years were reported consistently by both mothers and teachers. More information on selection to the 4 groups is available in eAppendix 5 and eTables 6 to 9 in the Supplement.
The analytic sample ranged from 7466 (CRP) to 9298 (abdominal obesity) individuals depending on the outcome. A systematic pattern emerged with respect to all outcomes (Table 2), with the 4 groups being consistently ordered. The stable-high group had the least favorable levels on all biomarkers and the highest mortality rates. The TIAO group had the second worst biomarker levels and mortality rates, followed by the moderate group. The stable-low group had the most favorable biomarker levels and the lowest mortality rates. In Table 3, we present ordinary least squares regression coefficients, risk ratios, and hazard ratios for the association between early-life mental health trajectories and all outcomes after multivariable adjustment. With a conventional 5% α level, early-life mental health was associated with less favorable levels of fibrinogen, CRP, HDL cholesterol, and abdominal obesity. The associations with fibrinogen, HDL, and abdominal obesity were owing to the stable-high group having different levels in these biomarkers compared with the stable-low group. Furthermore, compared with the stable-low group, the TIAO group had higher levels of fibrinogen and CRP and had increased risk for abdominal obesity.
We observed evidence of modification by sex in the association between early-life mental health, HDL, and abdominal obesity. With respect to HDL, there was strong evidence of interaction with sex in the TIAO group (B for interaction, −6.203; 95% CI, −9.333 to −3.074; P < .001) and considering that interaction terms require more power, there was also suggestive evidence for the stable-high group (B for interaction, −4.214; 95% CI, −8.795 to 0.367; P = .07). For abdominal obesity, there was strong interaction with sex in the stable-high group (risk ratio for interaction, 1.378; 95% CI, 1.109-1.711; P = .004). In Figure 1 and Figure 2, we present results (point estimates and 95% CIs) from models stratified by sex, where early-life mental health is associated with both outcomes only in female individuals. With respect to HDL, the TIAO and stable-high groups had less favorable levels compared with the stable-low group (B, −3.404; 95% CI, −5.699 to −1.110; P = .004 and B, −3.986; 95% CI, −7.351 to −0.621; P = .02, respectively). Similarly, the TIAO and stable-high groups had increased risk for abdominal obesity, compared with the stable-low group (risk ratio, 1.190; 95% CI, 1.019-1.389; P = .03 and risk ratio, 1.522; 95% CI, 1.278-1.813; P < .001).
From a null hypothesis significance testing perspective, carrying out multiple tests increases the probability of type I error. We further evaluated our findings with the Bonferroni correction31 that controls for familywise error rate. Considering each biomarker as a family of tests, with 3 tests carried out, the α level was .05 / 3 = .016. Using this level, we were able to detect associations with fibrinogen, CRP, HDL (female individuals only), and abdominal obesity. In a more conservative scenario, where all biomarkers were considered a family of tests, the Bonferroni-adjusted α level would be .05 / (3 × 8) = .002. With this very conservative correction, we were able to detect an association with abdominal obesity. Going beyond the conventional null hypothesis significance testing dichotomy and considering the magnitude and patterns of point estimates as well as coverage of CIs as presented in Table 3, we observed a systematic pattern of associations with fibrinogen, CRP, HDL, glycated hemoglobin, and abdominal obesity.
Having excluded participants who had died prior to age 16 years, our analytic sample for the mortality analyses was n = 15 067. With no evidence of effect modification by sex in the mental health–mortality association, we present results from the pooled sample. For those who participated in the data collection at age 16 years, there were 1068 deaths and we analyzed a total of 765 209 person-years, with age 58 years being the last year of observation. More deaths were observed in the stable-high and TIAO groups compared with the moderate and stable-low groups (Table 2). The early-life mental health trajectories were associated with premature all-cause mortality, with stable-high and TIAO groups having elevated risk of mortality compared with the stable-low group (stable-high group: hazard ratio, 1.878; 95% CI, 1.501-2.350; P < .001; TIAO group: hazard ratio, 1.412; 95% CI, 1.174-1.698; P < .001). We were able to detect an association with all-cause mortality under any familywise error rate correction.
With respect to cause-specific mortality, we analyzed a total of 718 254 person-years, with age 50 years being the last year of observation. For those who participated in the data collection at age 16 years, we observed 82 unintentional injuries, 82 deaths due to drug and alcohol poisonings, suicide, chronic liver diseases, and cirrhosis (with 55 suicides among those), and 154 circulatory system–related deaths. Considering the low statistical power owing to the small number of events, the stable-high group had elevated risk for unintentional injuries and deaths due to drug and alcohol poisonings, suicide, chronic liver diseases, and cirrhosis, as well as suicide when analyzed separately (results for suicide analyzed separately available in eAppendix 7 in the Supplement). The TIAO group had elevated risk for deaths due to drug and alcohol poisonings, suicide, chronic liver diseases, and cirrhosis. The elevated risk for these deaths in the stable-high compared with the stable-low group was detectable under any family wise error rate correction (hazard ratio, 3.167; 95% CI, 1.529-6.558; P = .002).
Our main finding was that we observed associations between early-life mental health with biomarkers in midlife as well as premature mortality. The associations with biomarkers were more commonly seen in female individuals, whereas as expected from earlier follow-ups of the present birth cohort and 1 other study,4,5,24,32 associations with mortality were observed in both male and female individuals. Persistent experience of high levels of conduct problems and affective symptoms from age 7 to 16 years was associated with higher fibrinogen and CRP levels, lower HDL levels, and increased risk for abdominal obesity 28 years later. It was also associated with elevated risk for premature all-cause mortality and when separately examined, unintentional injuries and deaths due to drug and alcohol poisonings, suicide, chronic liver diseases, and cirrhosis. Experiencing the onset of high levels of conduct problems and affective symptoms in adolescence was associated with CRP and abdominal obesity in midlife, as well as increased risk for premature all-cause mortality and suicide (results for suicide risk available in eAppendix 7 in the Supplement). However, we did not observe an association with unintentional injuries nor cardiovascular mortality, which indicates that trajectory-specific mechanisms underlie the association with premature mortality. In this group, it appears to be primarily driven by direct results of psychopathology, whereas in the stable-high group was also driven by unintentional injuries and perhaps by circulatory system deaths, but our estimates were imprecise.
Our findings, if causal and generalizable to younger cohorts, may have implications for overall population health but also for age-related disease as prevention of morbidity and disability in old age could be achieved by improving early-life mental health.33 The rising levels of early-life and midlife psychological distress in more recently born cohorts34-36 point to early-life mental health being one of the drivers of the increases in midlife mortality caused by deaths due to drug and alcohol poisonings, suicide, chronic liver diseases, and cirrhosis37 and perhaps to the recent stagnation of life expectancy in the UK.38 The observed link with markers of health and mortality implies that generational differences in early-life mental health may partly underlie the observed expansion of morbidity in the UK.39 As a target of interventions, early-life mental health has the potential to mitigate the effect that adverse socioeconomic circumstances in childhood, adverse childhood experiences, and early-life low cognitive ability have in adulthood as it can be plausibly assumed to lie on the pathway that links these with adult health, mortality, and social and economic outcomes.40-42 Furthermore, interventions to improve early-life mental health could prevent future excess morbidity and mortality due to the possible (at the time of writing) psychological effect on children and adolescents that coronavirus disease 2019–related measures might have.43,44
Experiencing the onset of mental health symptoms in the transition from childhood to adolescence was found to be detrimental with respect to both biomarkers and premature mortality. As the assessment of individuals prior to 16 years was carried out at age 11 years, it appears that the transition from childhood to adolescence may be a sensitive period that warrants further research to elucidate the mechanisms that links the onset of mental health problems during this period to adult health and premature mortality.
The mechanisms of action that underlie the observed modification by sex in the association with HDL and abdominal obesity that were only observed in female individuals may include societal pressures and/or sex-related inequalities that may exacerbate the association of conduct problems and affective symptoms with midlife health in female individuals, pointing to sex-specific pathways, in a similar manner as to those suggested for adult depression.45 It appears that at least in this cohort, early-life mental health may be one of the factors that underlie the so-called sex survival paradox, a term that describes the observation that female individuals are on average less healthy but live longer than male individuals.46,47
Strengths and Limitations
Strengths of this study include the availability of a large population-based and representative prospective study, the 3 assessments of early-life mental health, which is of major importance given the time-varying nature of mental health in childhood and adolescence,48 the inclusion of both conduct problems and affective symptoms in the assessment of early-life mental health, and the wealth of information on potential confounders. However, our findings can only be generalized to those born in Britain in 1958 or close to this year. Furthermore, our data are derived from an observational longitudinal study and bias due to unmeasured confounding cannot be ruled out. However, sensitivity analyses with the E value49,50 suggest that for the observed associations with abdominal obesity and all-cause mortality, strong confounding, stronger than that observed in our data, would be needed to completely explain away our findings (eAppendix 8 and eFigures 1-4 in the Supplement). An alternative specification of our models where early-life mental health trajectories where restricted to data from ages 7 and 11 years, returned similar results (eAppendix 10 and eTables 13-15 in the Supplement). This analysis considerably shortened the temporal gap from the measurement of the confounders to when early-life mental health was assessed and provided further reassurance against unmeasured confounding. As in any longitudinal survey, missing data due to attrition are unavoidable. We used multiple imputation, augmenting our models with auxiliary variables in the imputation phase to maximize the plausibility of the missing-at-random assumption and restore sample representativeness, but bias due to a nonignorable missing data-generating mechanism cannot be ruled out.
We did not include hyperactivity symptoms in our analysis, since building on previous work, we focused on affective symptoms and conduct problems as potential mechanisms of action with our outcomes have been shown to exist in the National Child Development Study and other cohorts. Another limitation is that reports from both mothers and teachers were only available at age 16 years. While they agreed in their assessment of those who experienced persistent mental health symptoms (stable-high group) and those who did not (stable-low group), there was stark disagreement in the TIAO group. The similar assessments for 2 of the 4 groups makes the possibility of teachers misreporting unlikely. A more plausible explanation is that symptoms were context specific and lack of insight into the other context might lead to misreporting. However, results from sensitivity analysis using only maternal reports were similar to our main analysis, indicating that our findings are not reporter driven (eAppendix 9 and eTables 10-12 in the Supplement). Further sensitivity analysis taking into account measurement error in the mental health trajectories returned similar results to the ones presented here (eAppendix 4 and eTable 5 in the Supplement). Despite these efforts, the extent to which undetected systematic measurement error may have biased our results is unknown.
Our findings, if causal and generalizable to younger cohorts, may have implications for public health policy, especially if mental health is worse in more recently born cohorts.34,35 Effective interventions on early-life mental health have the potential to shift the distribution of risk51 and may have additional benefits, reducing future risk of premature mortality and improving multiple physical health outcomes.
Corresponding Author: George B. Ploubidis, PhD, Centre for Longitudinal Studies Department of Social Science, University College London, 55-59 Gordon Square, Room 212, London WC1H 0NU, England (g.ploubidis@ucl.ac.uk).
Accepted for Publication: July 19, 2020.
Published Online: September 30, 2020. doi:10.1001/jamapsychiatry.2020.2893
Author Contributions: Dr Ploubidis had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: Ploubidis, Batty, Bann, Goodman.
Acquisition, analysis, or interpretation of data: Ploubidis, Batty, Patalay, Goodman.
Drafting of the manuscript: Ploubidis.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Ploubidis.
Obtained funding: Goodman.
Conflict of Interest Disclosures: None reported.
Funding/Support: This study was supported by the Economic and Social Research Council (grants ES/M008584/1 and ES/M001660/1).
Role of the Funder/Sponsor: The funder 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.
Additional Contributions: We thank all participating cohort members for their continued commitment to this study.
1.Clark
C, Rodgers
B, Caldwell
T, Power
C, Stansfeld
S. Childhood and adulthood psychological ill health as predictors of midlife affective and anxiety disorders: the 1958 British Birth Cohort.
Arch Gen Psychiatry. 2007;64(6):668-678. doi:
10.1001/archpsyc.64.6.668
PubMedGoogle ScholarCrossref 7.Winning
A, Glymour
MM, McCormick
MC, Gilsanz
P, Kubzansky
LD. Psychological distress across the life course and cardiometabolic risk: findings from the 1958 British Birth Cohort Study.
J Am Coll Cardiol. 2015;66(14):1577-1586. doi:
10.1016/j.jacc.2015.08.021
PubMedGoogle ScholarCrossref 8.Patalay
P, Moulton
V, Goodman
A, Ploubidis
GB. Cross-domain symptom development typologies and their antecedents: results from the UK Millennium Cohort Study.
J Am Acad Child Adolesc Psychiatry. 2017;56(9):765-776.e2. doi:
10.1016/j.jaac.2017.06.009
PubMedGoogle ScholarCrossref 15.Rutter
M, Tizard
J, Whitmore
K. Education, Health and Behaviour. Longman Publishing Group; 1970.
16.De Los Reyes
A. Introduction to the special section: More than measurement error: Discovering meaning behind informant discrepancies in clinical assessments of children and adolescents.
J Clin Child Adolesc Psychol. 2011;40(1):1-9. doi:
10.1080/15374416.2011.533405PubMedGoogle ScholarCrossref 17.Van Roy
B, Groholt
B, Heyerdahl
S, Clench-Aas
J. Understanding discrepancies in parent-child reporting of emotional and behavioural problems: effects of relational and socio-demographic factors.
BMC Psychiatry. 2010;10:56-56. doi:
10.1186/1471-244X-10-56
PubMedGoogle ScholarCrossref 22.Pinto Pereira
SM, Ki
M, Power
C. Sedentary behaviour and biomarkers for cardiovascular disease and diabetes in mid-life: the role of television-viewing and sitting at work.
PLoS One. 2012;7(2):e31132. doi:
10.1371/journal.pone.0031132
PubMedGoogle Scholar 26.Muthén
B. A general structural equation model with dichotomous, ordered categorical, and continuous latent variable indicators.
Psychometrika. 1984;49(1):115-132. doi:
10.1007/BF02294210
Google ScholarCrossref 28.Oberski
D. Mixture models: Latent profile and latent class analysis. In: Robertson J, Kaptein M, eds.
Modern Statistical Methods for HCI. Springer; 2016:275-287. doi:
10.1007/978-3-319-26633-6_12
36.Patalay
P, Gage
SH. Changes in millennial adolescent mental health and health-related behaviours over 10 years: a population cohort comparison study.
Int J Epidemiol. 2019;48(5):1650-1664. doi:
10.1093/ije/dyz006
PubMedGoogle ScholarCrossref 38.Leon
DA, Jdanov
DA, Shkolnikov
VM. Trends in life expectancy and age-specific mortality in England and Wales, 1970-2016, in comparison with a set of 22 high-income countries: an analysis of vital statistics data.
Lancet Public Health. 2019;4(11):e575-e582. doi:
10.1016/S2468-2667(19)30177-X
PubMedGoogle ScholarCrossref 40.Birnie
K, Cooper
R, Martin
RM,
et al; HALCyon study team. Childhood socioeconomic position and objectively measured physical capability levels in adulthood: a systematic review and meta-analysis.
PLoS One. 2011;6(1):e15564. doi:
10.1371/journal.pone.0015564
PubMedGoogle Scholar