http://orcid.org/0000-0002-3785-7417
Figures
Abstract
Current knowledge on the pathogenesis of diastolic heart failure predominantly rests on case-control studies involving symptomatic patients with preserved ejection fraction and relying on invasive diagnostic procedures including endomyocardial biopsy. Our objective was to gain insight in serum and urinary biomarkers reflecting collagen turnover and associated with asymptomatic diastolic LV dysfunction. We randomly recruited 782 Flemish (51.3% women; 50.5 years). We assessed diastolic LV function from the early and late diastolic peak velocities of the transmitral blood flow and of the mitral annulus. By sequencing urinary peptides, we identified 70 urinary collagen fragments. In serum, we measured carboxyterminal propeptide of procollagen type 1 (PICP) as marker of collagen I synthesis and tissue inhibitor of matrix metalloproteinase type 1 (TIMP-1), an inhibitor of collagen-degrading enzymes. In multivariable-adjusted analyses with Bonferroni correction, we expressed effect sizes per 1-SD in urinary collagen I (uCI) or collagen III (uCIII) fragments. In relation to uCI fragments, e’ decreased by 0.183 cm/s (95% confidence interval, 0.017 to 0.350; p = 0.025), whereas E/e’ increased by 0.210 (0.067 to 0.353; p = 0.0012). E/e’ decreased with uCIII by 0.168 (0.021 to 0.316; p = 0.018). Based on age-specific echocardiographic criteria, 182 participants (23.3%) had subclinical diastolic LV dysfunction. Partial least squares discriminant analysis contrasting normal vs. diastolic LV dysfunction confirmed the aforementioned associations with the uCI and uCIII fragments. PICP and TIMP-1 increased in relation to uCI (p<0.0001), whereas these serum markers decreased with uCIII (p≤0.0006). Diastolic LV dysfunction was associated with higher levels of TIMP-1 (653 vs. 696 ng/mL; p = 0.013). In a general population, the non-invasively assessed diastolic LV function correlated inversely with uCI and serum markers of collagen I deposition, but positively with uCIII. These observations generalise previous studies in patients to randomly recruited people, in whom diastolic LV function ranged from normal to subclinical impairment, but did not encompass overt diastolic heart failure.
Citation: Zhang Z-Y, Ravassa S, Yang W-Y, Petit T, Pejchinovski M, Zürbig P, et al. (2016) Diastolic Left Ventricular Function in Relation to Urinary and Serum Collagen Biomarkers in a General Population. PLoS ONE 11(12): e0167582. https://doi.org/10.1371/journal.pone.0167582
Editor: Elena Cavarretta, Universita degli Studi di Roma La Sapienza, ITALY
Received: July 20, 2016; Accepted: November 16, 2016; Published: December 13, 2016
Copyright: © 2016 Zhang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Consent given by study participants did not include data sharing with third parties. Anonymized data can be made available to investigators for research based on a motivated request to be addressed to the corresponding author via jan.staessen@med.kuleuven.be.
Funding: The European Union (HEALTH-FP7-278249-EUMASCARA, HEALTH-F7-305507 HOMAGE) and the European Research Council (Advanced Researcher Grant 2011-294713-EPLORE and Proof-of-Concept Grant 713601-uPROPHET) and the Fonds voor Wetenschappelijk Onderzoek Vlaanderen, Ministry of the Flemish Community, Brussels, Belgium (G.0881.13, G.088013, and 11Z0916N) currently support the Studies Coordinating Centre in Leuven. The funders provide support in the form of salaries for authors (ZYZ, WYY, TP, FFW, LT, LJ), but did not have any additional role in the study design, data collection or analysis, decision to publish, or preparation of the manuscript. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Competing interests: M Pejchinovski, P Zürbig and T Koeck are employees of Mosaiques-Diagnostics GmbH. H Mischak is the co-founder and co-owner of Mosaiques Diagnostics GmbH (http://mosaiques-diagnostics.de/diapatpcms/mosaiquescms/front_content.php?idcat=161) and Diapat GmbH (http://www.diapat.com/DiaPat-Diagnostik/diapat-en), Hannover, Germany. This commercial affiliation did not alter our adherence to PLOS ONE policies on sharing data. The other authors did not declare a conflict of interest.
Introduction
Diastolic heart failure, also known as heart failure with preserved ejection fraction (HFpEF) represents half of all heart failure cases. It has a prognosis as dire as heart failure with reduced ejection fraction with a rate of cardiovascular mortality of over 30% within one year of the first hospital admission. In patients with diastolic heart failure, the left ventricular (LV) wall undergoes fibrosis characterised by increased interstitial deposition [1] and cross-linking of collagen I at the detriment of collagen III [2,3]. Small increases in the collagen I/III ratio augment myocardial stiffness, thereby reducing early diastolic LV filling and increasing LV filling pressure [4,5]. However, information on the asymptomatic phases of the disease remains scarce. This knowledge gap is particularly relevant, because the prevalence of asymptomatic diastolic LV dysfunction in the general population is as high as 25%[6] with a 10% risk of deterioration over 5 years [7].
Capillary electrophoresis coupled with high-resolution mass spectrometry (CE-MS) enables detection of over 5000 distinct peptide fragments in urine samples [8]. Our previous studies revealed a unique urinary proteomic signature, which in case-control studies [9] and in the general population [10] was reproducibly associated with subclinical diastolic LV dysfunction and which predicted the incidence of adverse cardiovascular outcomes over and beyond traditional cardiovascular risk factors [11]. We hypothesised that jointly linking diastolic LV dysfunction to urinary and serum markers of collagen turnover might increase our understanding of LV dysfunction. With the aim to generalise observations in patients with diastolic heart failure [2–5] to the early still asymptomatic stage of the disease in the population at large, we analysed the Flemish Study on Environment, Genes and Health Outcomes (FLEMENGHO) [10]. First, we searched for association of diastolic LV function with individual urinary peptides with known amino-acid sequence, thereby identifying collagen types as the parent proteins. Next, we correlated single urinary peptides significantly associated with diastolic LV function with circulating markers of cardiac collagen turnover, which in previous studies [5,12] predicted mortality and cardiovascular events.
Materials and Methods
Study population
FLEMENGHO complies with the Helsinki declaration [13] for research in human subjects. The Ethics Committee of the University of Leuven approved the study [10]. At each contact, participants gave informed written consent. Recruitment for FLEMENGHO started in 1985 [10,14–17]. From August 1985 until November 1990, a random sample of the households living in a geographically defined area of Northern Belgium was investigated with the goal to recruit an equal number of participants in each of six subgroups by sex and age (20–39, 40–59, and ≥60 years). All household members with a minimum age of 20 years were invited to take part, provided that the quota of their sex-age group had not yet been satisfied. From June 1996 until January 2004 recruitment of families continued using the former participants (1985−1990) as index persons and also including teenagers. The initial participation rate was 78·0%. The participants were repeatedly followed up at the field centre in the catchment area (North Limburg, Belgium). From May 2005 until May 2010, an invitation letter was mailed to 1208 former participants for a follow-up examination. However, 153 were unavailable, because they had died (n = 26), had been institutionalised or were too ill (n = 27), or had moved out of the area (n = 100). Of the remaining 1055 former participants, 828 renewed informed consent [10]. The participation rate at re-examination was 78.5%. We excluded 46 participants from analysis, because serum samples were unavailable (n = 22), or because arrhythmia (n = 8), paced heart rhythm (n = 3) or poor echocardiographic image quality (n = 13) rendered assessment of diastolic LV function difficult. Thus, the number of participants statistically analysed totalled 782.
Echocardiography
The acquisition, offline analysis, staging of diastolic LV function and quality control are extensively described in previous publications [6,7,18] and summarised below.
Data acquisition.
One observer (TKu) did the ultrasound examination [6], using a Vivid7 Pro (GE Vingmed, Horten, Norway) device interfaced with a 2.5- to 3.5-MHz phased-array probe. For off-line analysis, she recorded at least five heart cycles according to the recommendations of the American Society of Echocardiography [19]. M-mode echocardiograms of the left ventricle (LV) were recorded from the parasternal long-axis under control of the two-dimensional image. The ultrasound beam was positioned just below the mitral valve at the level of the posterior tendinous chords. To record the mitral and pulmonary vein (PV) flow velocities from the apical window, the observer positioned the Doppler sample volume at the mitral valve tips, in the right superior PV, and between the LV outflow and mitral inflow, respectively. From the apical window, the observer placed a 5-mm Doppler sample at the septal, lateral, inferior and posterior sites of the mitral annulus to record low-velocity, high-intensity myocardial signals at a high frame rate (>190 frames per second), while ensuring parallel alignment of the ultrasound beam with the myocardial segment of interest.
Off-line analysis.
One reader (TKu) analysed the digitally stored images, averaging three heart cycles, using a workstation running EchoPac software, version 4.0.4 (GE Vingmed, Horten, Norway). LV internal diameter and interventricular septal and posterior wall thickness were measured at end-diastole from the 2-dimensionally guided M-mode tracing. When optimal orientation of M-mode ultrasound beam could not be obtained, the reader performed linear measurements on correctly oriented 2-dimensional images. End-diastolic LV dimensions were used to calculate LV mass [19]. Left atrial (LA) volume was calculated using the prolate-elipsoid method from the LA dimensions in three orthogonal planes and indexed to body surface area [19]. From the transmitral flow signal, the reader determined peak early diastolic velocity (E), peak late diastolic velocity (A), the E/A ratio, and the transmitral A flow duration. From the PV flow signal, she measured the duration of PV reversal flow during atrial systole. From the tissue Doppler recordings, the observer measured peak early (e') and peak late (a') diastolic mitral annular velocities, and the e'/a' ratio at the four acquisition sites (septal, lateral, inferior, and posterior).
Staging of diastolic LV function.
S1 Table summarises the physiological interpretation of the echocardiographic measurements reflecting diastolic LV function. As shown in Fig 1, guideline-driven echocardiographic criteria to stage patients with advanced diastolic LV dysfunction [20] leave a large proportion of people unclassified in population studies [6,21]. We therefore developed age-specific criteria in a healthy reference sample drawn from FLEMENGHO [6] and replicated these criteria in a European population study [21]. To stage diastolic LV function in our current study, as previously described [6,21], we combined the velocities of the transmitral blood flow and the mitral annular movement. Group 1 included patients with an abnormally low age-specific transmitral E/A ratio indicative of impaired relaxation, but without evidence of increased LV filling pressures (E/e' ≤ 8.5). Group 2 had mildly-to-moderately elevated LV filling pressure (E/e' > 8.5) and an E/A ratio within the normal age-specific range. Differences in durations between the transmitral A flow (Ad) and the reverse flow in the pulmonary veins (ARd) during atrial systole (Ad < ARd + 10) and the left atrial volume index (≥ 29 mL/m2; the 97.5th percentile of the distribution in a reference sample of 239 healthy FLEMENGHO participants [6]) were checked to confirm elevation of LV filling pressure (S1 Table). Group 3 had an elevated E/e' ratio and an abnormally low age-specific E/A ratio (combined dysfunction). We combined these three groups for analysis.
Distributions of: (A) the peak velocity of the mitral annular movement during early diastole (e’); (B) the ratio of the peak velocities of transmitral blood flow to mitral annular movement in early diastole; and (C) the left atrial volume index (LAVI). In line with the US guideline (J Am Soc Echocardiogr 2016;29:277–314), e’ was measured at septal wall in panel A, but was the average of the septal and lateral sample sites in panel B. Arrows indicate the cut-off values for diastolic LV dysfunction recommended for clinical use in patients with advanced diastolic dysfunction.
Quality control.
Intra-observer reproducibility was the 2-SD interval about the mean of the relative differences across pairwise readings. The intra-observer reproducibility for the tissue Doppler peak velocities across the four sampling sites ranged from 4.5% to 5.3% for e’ and from 4.0% to 4.5% for a’ [6]. Reproducibility was 2.2% for internal end-diastolic LV diameter, 4.6% for LV wall thickness, and 4.3% for LV mass [18].
Urinary peptides
Participants collected 24-h urine samples within one week of the echocardiographic examinations. Aliquots (0.7 mL) were stored at –80°C and thawed immediately before analysis. Previous publications provide a detailed account of the methods applied for preparation of the urine samples and running the CE-MS analysis [22,23]. For targeted sequencing, urine samples were analysed on a Dionex Ultimate 3000 RSLS nano flow system (Dionex, Camberly, UK) or on a Beckman CE, coupled to an Orbitrap Velos MS instrument (Thermo Scientific, Waltham, Massachusetts, US) [24]. The data files were analysed using Proteome Discoverer 1.2 (precursor mass tolerance, 10 ppm; fragment mass tolerance, 0.05 Da) and were searched against the UniProt human non-redundant database without enzyme specificity. No fixed modifications were selected. Oxidation of methionine and proline were considered as variable modifications. The criteria for accepting sequences were high confidence (Xcorr ≥ 1.9) and absence of unmodified cysteine. A strong correlation between peptide charge at the working pH of 2 and capillary electrophoresis migration time was used to avoid falsely characterised peptides [25].
We identified 88 sequenced peptides with mean signal amplitude different from undetectable in over 95% of participants. Proteomics and peptidomics data usually display missing values, which result from either the absence of a particular peptide in certain datasets or from the fact that this particular peptide may be below of the limit of detection of the technological platform. For the current analysis, the values of peptides undetectable in less than 5% of randomly varying study participants were set to the minimum of the distribution of each sequenced peptide. This strategy has been validated by other research groups [26,27]. Furthermore, the cardiac extracellular matrix (ECM) is predominantly composed of fibrillar collagen I (85%) and III (11%) and small amounts of collagen IV and V co-distributed with collagen I. We therefore confined our analyses to 70 urinary peptide fragments, which could be traced back to the aforementioned collagen types (S2 Table).
Circulating biomarkers
On the day of echocardiography, venous blood samples were drawn after at least 8 hours of fasting. Serum was analysed for carboxyterminal propeptide of procollagen I (PICP), a marker of collagen I synthesis; carboxyterminal telopeptide of collagen I (CITP) and tissue inhibitor of the matrix metalloproteinase type 1 (TIMP-1) as markers of degradation of collagen I; and amino terminal propeptide of procollagen type III (PIIINP), a marker of synthesis and degradation of collagen III (S2 Fig). Serum markers were measured in 740 of 782 participants (94.6%). As described previously [3], PICP (Quidel Corporation, San Diego, CA), TIMP-1 (GE Healthcare Life Sciences, Buckinghamshire, UK) and PIIINP (MyBioSource, San Diego, CA) were quantified by sandwich enzyme linked immunosorbent assay and CITP by a quantitative enzyme immunoassay (Orion Diagnostica, Espoo, Finland). The detection limits were 0.2 μg/L for PICP (inter- and intra-assay coefficients of variability, 6.4% and 4.5%), 0.3 μg/L for CITP (13.1% and 10.0%), 1.25 ng/mL for TIMP-1 (12.8% and 2.6%), and 31.3 pg/mL for PIIINP (<15.0%) [3].
Other measurements
Blood pressure was the average of five consecutive auscultatory readings obtained according to European guidelines [28] with a standard mercury sphygmomanometer with the participant resting in the seated position for at least 10 minutes. Hypertension was a blood pressure of at least 140 mm Hg systolic or 90 mm Hg diastolic or use of antihypertensive drugs. Body mass index was weight in kilograms divided by height in meters squared. Plasma glucose and serum total and high-density lipoprotein (HDL) cholesterol, creatinine, γ-glutamyltransferase (biomarker of alcohol intake) and insulin were measured by automated methods in a single certified laboratory. Glomerular filtration rate was derived by the Chronic Kidney Disease Epidemiology Collaboration equation [29]. We graded chronic kidney disease according to the National Kidney Foundation K/DOQI guideline [30] into stages 1, 2, 3A, 3B, 4 and 5 based on eGFR values of ≥90, 60–89, 45–59, 30–44, 15–29 and <15 mL/min/1.73 m2. Diabetes mellitus was a self-reported diagnosis, a fasting glucose level of 7 mmol/L or higher, or use of antidiabetic agents [31].
Statistical Analysis
For database management and statistical analysis, we used the SAS system, version 9.4 (SAS Institute Inc., Cary, NC). Means were compared using the large-sample z-test and proportions by Fisher’s exact test. We normalised the distributions of γ-glutamyltransferase, insulin and PIIINP by a logarithmic transformation.
We identified covariables to be retained in the analyses by a stepwise regression procedure with p-values for covariables to enter and stay in the models set at 0.15. In continuous analyses, we standardised diastolic LV function for the average in the whole study population (mean or ratio) of the covariables identified by stepwise regression. While accounting for covariables, we regressed the indexes of diastolic LV function on the urinary peptide markers and constructed–log10 probability plots. Based on the number of parent collagen proteins (types I, III, IV and V), we adjusted significance levels by the Bonferroni method. In the next step of our analyses, we applied partial least squares discriminant analysis (PLS-DA). PLS-DA is a statistical technique that constructs predictive models for categorical outcomes in relation to correlated high dimensional predictors. PLS-DA allowed us to identify a set of independent latent factors that were linear combinations of the urinary peptides and that maximised the covariance between the urinary peptides and the variables describing diastolic LV function. We studied normal vs. diastolic dysfunction in relation to the latent factors. We retained the smallest number of latent factors for which the predicted residual sums of squares (PRESS, calculated using leave-one-out cross-validation) did not differ significantly (p>0.10) from the model with the minimum PRESS value as assessed by the van der Voet T2 statistic. The importance of each urinary peptide in the construction of the PLS factors and in the association of diastolic function was assessed from the Variable Importance in Projection (VIP) scores of Wold. Using principal component analysis, we summarized the urinary peptides with a VIP score higher than 1.0 and a correlation coefficient of less than –0.04 or higher than 0.04 into a single factor. Finally, we constructed Receiver Operating Characteristic (ROC) plots and evaluated the area under the ROC curve (AUC).
Results
Characteristics of participants
Of 782 participants, 401 (51.3%) were women. All were White Europeans. Mean values (SD) in the 782 participants were 50.5 (15.5) years for age, 26.4 (4.2) kg/m2 for body mass index, 129.1 (17.5) mm Hg and 79.6 (9.5) mm Hg for systolic and diastolic blood pressure, and 5.25 (0.97) mmol/L for total cholesterol. Among all participants, 327 (41.8%) had hypertension, of whom 199 (60.9%) were on antihypertensive drug treatment, and 9 (1.2%) had diabetes.
The prevalence of diastolic LV dysfunction amounted to 182 (23.3%). Two participants had tricuspid regurgitation and 10 had an E/e’ ratio exceeding 14. Among participants with normal or impaired diastolic LV function only two and seven had an ejection fraction of less than 50%. There was no difference in the ejection fraction between these two groups (68.2 vs. 69.0%, p = 0.29). Table 1 shows that cardiovascular risk factors were more prevalent or elevated in participants with diastolic dysfunction compared with those with normal LV function, with exception of the proportion of women, smokers and patients with diabetes. Fig 1 shows the distributions of echocardiographic indexes, e’ and E/e’. Diastolic LV dysfunction was characterised by impaired relaxation in 68 patients (37.4%) or an elevated filling pressure in the presence of a normal (90 [49.5%]) or low (24 [13.2%]) age-specific E/A ratio [6,21]. The echocardiographic characteristics of participants by these categories of diastolic LV function appear in Table 2. Compared with participants with normal LV function, patients with diastolic LV dysfunction more frequently (p≤0.003) used medications: any diuretic (5.3% vs. 25.3%), loop diuretics (0.2% vs. 2.2%), spironolactone (1.0% vs. 7.7%), β-blockers (9.0% vs. 36.3%), angiotensin-converting enzyme inhibitors (3.0% vs. 9.3%), angiotensin I type-1 receptor blockers (2.2% vs. 10.4%) and vasodilators including calcium-channel blockers and α-blockers (2.8% vs. 11.0%).
S3 Table lists the covariables considered for entry and retained for multivariable adjustment of diastolic LV function. Based on the results of the stepwise regression analysis, we standardised all echocardiographic indexes of diastolic function for sex, age, body mass index, mean arterial pressure, heart rate, serum total cholesterol, γ−glutamyltransferase and creatinine, plasma glucose, LV mass index and treatment with diuretics, β−blockers and inhibitors of the renin-angiotensin system.
Urinary biomarkers
In general, urinary peptides derived from collagen I correlated with one another. For instance, the correlation coefficients relating p70635 to the other collagen I fragments ranged from 0.14 to 0.30 (p<0.0001 for all). The correlation between the two collagen III fragments, p107460 and p112106, was 0.19 (p<0.0001). There was a weak inverse correlation (r = –0.10; p = 0.0054) between p70635 (collagen I fragment) and p107460 (collagen III fragment).
Among all 782 participants, eGFR averaged 80.7 (16.4) mL/min/1.73 m2. The prevalence of renal function stages 1, 2, 3A, 3B, 4 and 5 amounted to 175 (22.4%), 552 (70.8%), 49 (6.3%), 5 (0.6%), 1 (0.1%) and 0, respectively. With adjustments applied for mean arterial pressure, waist-to-hip ratio, smoking, γ-glutamyltransferase, the total-to-HDL cholesterol ratio, plasma glucose, and use of antihypertensive medications by drug class, none of the associations of eGFR with the urinary collagen fragments reached significance (p≥0.082).
Continuous analysis
S1 Fig. shows the–log10(p) probability plot of the multivariable-adjusted associations of various indexes of diastolic LV function with the urinary peptides. With Bonferroni correction applied, the urinary peptides that remained significantly associated with the Doppler indexes of diastolic LV function included six fragments of collagen I (p70635, p72896, p73697, p77018, p77952 and p115491) and two fragments of collagen III (p107460 and p112106). Focusing on collagen I (Table 3), e’ peak velocity and the e’/a’ ratio decreased respectively with p70635 (effect size per 1-SD increment, –0.183; p = 0.025) and p77952 (–0.041; p = 0.006), whereas the E/e’ ratio increased with p72896 (0.164; p = 0.020), p77018 (0.210; p = 0.0012) and p115491 (0.162; p = 0.019). The a’ peak velocity declined with p72896 (–0.192; p = 0.0024) and p73697 (–0.160; p = 0.016). In relation to collagen III fragments, A peak velocity (–1.450; p = 0.0024) and the E/e’ ratio (–0.168; p = 0.018) declined with p107460 and the A peak also with p112106 (–1.334; p = 0.006). Sensitivity analyses with additional adjustment for HDL cholesterol and insulin produced confirmatory results (Table 3). None of the indexes of diastolic LV function was significantly associated with the sequenced collagen IV or V fragments (S2 Table).
Categorical analysis.
We dichotomised the study population in 600 participants with normal LV function and 182 with diastolic LV dysfunction. The PLS-DA procedure yielded two latent factors accounting for 10.2% and 6.3% of the variance in the urinary peptides and 16.5% in total.
Using a VIP score of 1.5 and a correlation coefficient of –0.04 as cut-offs, normal diastolic LV function (Fig 2, left top side of the V-plot) was associated with the collagen I fragments p35339, p57531, and p91542. Using a VIP score of 1.5 and a correlation coefficient of 0.04 as cut-offs, diastolic dysfunction (Fig 2, right top side of the V-plot) was associated with collagen I fragment p77763, collagen III fragments p50840 and p105352, and collagen V fragment p104786. p77763 is a collagen I fragment with an amino-acid sequence very similar to p77018 (one proline residue being replaced by hydroxyproline; S2 Table), which was related to the indexes of diastolic LV function in the continuous analyses (Table 3 and S1 Fig). In general, the sequences of collagen I fragments associated with normal diastolic function were shorter than those associated with dysfunction (S2 Table). The level of the collagen I fragment p77018 was higher in patients with diastolic LV dysfunction than in participants with normal diastolic function, whereas the opposite was true for the collagen III fragment p107460 (Table 4).
We dichotomised the study population in 600 participants with normal LV function and 182 with subclinical diastolic LV dysfunction. VIP is the importance of each urinary fragments in the construction of the PLS factors. The correlation coefficients reflect the association of diastolic LV dysfunction with the urinary collagen fragments. Fragments associated with normal and diastolic LV dysfunction were sitting on the left and right arms, respectively (see Table 3). Fragments derived from collagen I, III, IV and V are labelled blue, red, yellow and green, respectively.
To evaluate diagnostic accuracy, we combined the urinary peptides with a VIP score higher than 1.5 and a correlation coefficient lower than –0.04 or higher than 0.04 into a single factor (p35339, p50840, p57531, p77763, p91542, p104786, p105352), using principal component analysis. The AUC for diagnosing diastolic LV dysfunction was 0.92 (95% CI, 0.89–0.94; p<0.0001) yielding a sensitivity, specificity and positive and negative predictive values of 67.6%, 93.3%, 75.5% and 90.5%, respectively.
Circulating biomarkers
With adjustments applied for metabolic confounders (body mass index, serum total cholesterol, γ-glutamyltransferase and creatinine, and plasma glucose) and concurrent treatment by drug class (Table 4), serum levels of CITP (6.26 vs. 5.34 μg/L), TIMP-1 (696 vs. 653 ng/mL) and PIIINP (562 vs. 447 pg/mL) were significantly higher (p≤0.020) in patients with diastolic LV dysfunction (n = 175) than in people with normal function (n = 565), with no between-group differences in serum PICP.
Table 5 lists the associations of circulating collagen biomarkers with the urinary collagen I and III fragments that reached a significance level of 1% or less. With association sizes expressed per 1-SD increment in collagen I fragments, PICP, CITP and PIIINP increased (p≤0.0016) by 14.9 μg/L, 0.31 μg/L and 4.63% in relation to p73697. TIMP-1 increased by 21.4 ng/mL (p = 0.0013) and by 39.4 ng/mL (p<0.0001) in relation to p77018 and p77763, respectively. The association sizes of PICP, CITP and TIMP-1 with collagen III fragments amounted to –5.39 μg/L, –0.62 μg/L and –30.9 ng/mL for p107460 (p≤0.0006) and to –0.44 μg/L for CITP.
Discussion
The novel findings in our article, all obtained in a general population, can be summarised as follows: (i) the correlations between urinary collagen I and III fragments were inverse; (ii) in continuous analyses, e’ peak and e’/a’ decreased and E/e’ increased with urinary collagen I fragments, whereas the A peak and E/e’ decreased with urinary collagen III fragments (S1 Fig and Table 3); (iii) the PLS-DA analysis contrasting normal vs. diastolic LV dysfunction confirmed the associations with urinary collagen I and III fragments (Fig 2); (iv) PICP, CITP and TIMP-1 increased in relation to urinary collagen I fragments, whereas these serum markers decreased in relation to urinary collagen III (Table 5); and (v) in categorical analyses, diastolic LV dysfunction was associated with higher levels of urinary collagen I fragments, lower levels of urinary collagen III degradation products, and higher levels of CITP and TIMP-1, but not PICP (Table 4).
Collagen I is a stiff fibrillar protein providing tensile strength, whereas collagen III forms an elastic network storing kinetic energy that is released during elastic recoil [32]. Histopathological [33–35] and expression [34] studies of endomyocardial biopsies suggested that in humans chronic heart failure with preserved ejection fraction is characterised by myocardial fibrosis with a predominant increase in collagen I. Zile and colleagues measured myocardial stiffness directly in myocardial biopsies of 70 patients undergoing coronary artery bypass grafting [36]. In comparison with controls without comorbidity and controls with hypertension, patients with hypertension and diastolic heart failure had an increased end-diastolic LV pressure, left atrial volume, and collagen-dependent passive stiffness [36]. In keeping with experimental studies in rats [37], among patients with dilated cardiomyopathy [32], tissue samples of patients with heart failure compared with those from controls with mild global LV dysfunction had a 2- to 6-fold increase in collagen I mRNA, a 2-fold increase in collagen III mRNA, resulting in a higher collagen I/III expression ratio (8.6 vs. 6.4). Our current findings are novel, because they show in a general population that echocardiographic indexes of diastolic LV function are associated with sequenced urinary collagen I and III fragments (S1 Fig and Table 3). In addition, the I/III ratio of urinary collagen fragments was higher in patients with diastolic LV dysfunction compared to participants with normal LV function (Table 4). Thus, sequencing of the urinary peptide fragments allowed us to translate previous observations in endomyocardial biopsies [32–35] to people randomly recruited from the general population, in whom diastolic LV function ranged from normal to subclinical dysfunction, but did not encompass overt diastolic heart failure.
During the cardiac cycle, the left atrium acts as a reservoir, receiving pulmonary venous return during LV systole; as a conduit, passively transferring blood to the LV during early diastole; and as a pump, actively priming the LV in late diastole [38]. Stiffening of the LV requires a greater contribution of the atrial contraction to late diastolic LV filling and is associated with a higher a’ peak velocity. Next, as diastolic LV function deteriorates, the a’ peak velocity decreases [39]. This so-called pseudo-normalisation (moderate diastolic LV dysfunction in S1 Table) might therefore underlie the inverse association between a’ peak velocity and urinary collagen I fragments, as observed in our current study (Table 3). However, excluding 13 patients with an e’/a’ ratio higher than unity or all patients with an E/e’ ratio exceeding 8.5 [6] did not confirm this interpretation. An alternative explanation is that worsening of diastolic LV function leads to collagen deposition in the atria [40] with higher collagen I/III ratio [40], impairment of the atrial reservoir function [41], increase in the left atrial volume [36], deterioration of electromechanical coupling [40], and therefore to lower a’ peak velocity.
PICP is released in a 1:1 stoichiometric ratio during conversion of procollagen I to collagen I (S2 Fig) and therefore its serum concentration is a direct indicator of concurrent collagen I synthesis [42]. In patients with hypertensive heart disease [43], circulating PICP, at least in part, originates from the heart, because there is a positive concentration gradient from the coronary sinus towards the antecubital vein with a high correlation between coronary and peripheral levels. This association was not present in normotensive controls [42]. Moreover, serum PICP concentration correlates well with histologically proven myocardial collagen type I deposition [44] and in response to pharmacological intervention changes in serum PICP associate with changes in myocardial collagen type I deposition [45]. Our current study moves the field forward by showing that in the general population PICP increased in relation to p73697, a urinary collagen I fragment, but decreased in relation to p107460, a marker of the more elastic [32] collagen III (Table 5).
Metalloproteinases catalyse the degradation of collagen I resulting in the release of CITP in a 1:1 stoichiometric ratio (S2 Fig). The role of CITP as a reliable biomarker of collagen breakdown is not firmly established, because its association with myocardial fibrosis was inconsistently reported as negative [46] or positive [47]. Notwithstanding the uncertainty in the clinical interpretation of circulating CITP levels, our study (Table 5) revealed that serum CITP correlated positively with p73697, a urinary collagen I fragment associated with worse diastolic LV function, and inversely with p107460, a urinary marker of collagen III breakdown associated with more performant diastolic LV function.
Circulating TIMP-1 inhibits the metalloproteinases and is a pro-fibrotic stimulus. Similar to PICP [43], a positive gradient and a direct correlation exist between the TIMP-1 concentrations in coronary sinus and antecubital vein blood in patients with hypertensive heart disease, but not in normotensive controls [33]. In hypertensive patients with heart failure but normal ejection fraction, elevated estimated capillary wedge pressure compared with normal LV filling pressure was associated with higher TIMP-1 levels and a lower metalloproteinase-1 to TIMP-1 ratio, indicative of lower breakdown of collagen [48]. Zile and coworkers confirmed that in patients with hypertension with or without diastolic heart failure, circulating TIMP-1 levels, but not metalloproteinases, were elevated compared to normotensive controls [36]. Our study moves current knowledge forward by demonstrating that in a general population diastolic LV dysfunction was associated with higher levels of TIMP-1 (Table 4) and that TIMP-1 increased in relation to urinary collagen I fragments (Table 5). By linking circulating TIMP-1 to urinary collagen I fragments, our observation support the hypothesis that an excess of TIMP-1 inhibits collagen degradation, thereby promoting collagen deposition in the myocardium and diastolic LV dysfunction characterised by higher LV filling pressure [42]. On the other hand, the urinary collagen III fragment p107460 was associated with better diastolic LV function and lower LV filling pressure (Table 3 and Fig 2) and lower levels of TIMP-1 (Table 5). In patients with heart failure due to ischaemic heart disease or dilated cardiomyopathy [37], serum PIIINP levels are highly correlated with the myocardial collagen III volume fraction. The positive association between the urinary collagen I fragment p73697 and circulating PIIINP (Table 5), formed in a 1:2 stoichiometric ratio during the conversion of procollagen III to mature collagen III (S2 Fig), probably reflects the joint increase in both collagen subtypes [49] during myocardial fibrosis.
Strong points of our study are the availability of Doppler indexes of early subclinical diastolic LV dysfunction measured on a continuous scale, the application of two approaches in the statistical analysis, and the demonstration of a pathophysiologically plausible correlation between sequenced urinary collagen fragments and the serum biomarkers of collagen turnover. The epidemiological angle enhances the relevance of our findings over and beyond that of case-control studies involving selected heart failure patients, who represent the end stage of a long pathogenetic process confounded by multiple comorbidities and poly-medication. However, our present study must also be interpreted within the context of its limitations. First, our findings originate from a cross-sectional analysis and therefore reflect a snapshot in each individual participant. From this point of view our results should be considered as hypothesis generating. Whether or not, the urinary proteomic and serum biomarkers can predict the course over time of diastolic LV dysfunction remains to be confirmed in longitudinal studies. Second, the pathogenetic drivers leading to diastolic LV dysfunction are multifaceted each with different contributions among people at risk. Whether or not, the urinary collagen markers can predict the course over time of diastolic LV dysfunction remains to be proven in longitudinal studies. Third, we could not apply the simplified US criteria for the diagnosis of diastolic LV dysfunction in clinical practice (Fig 1) for the simple reason that they do not align with the gradation from normal to impaired diastolic function in the general population. However, our classification system passed expert peer review [7,9–11,18,21]. Finally, epidemiological studies demonstrate association and the causal interpretation of associations between traits of interest and biomarkers rests on a careful interpretation of the literature.
In conclusion, by sequencing urinary collagen I and III fragments and by linking diastolic LV function with urinary and serum collagen biomarkers, our current findings generalize previous observations in patients to the population at large. Our current observations support the concept of porting the use of multidimensional biomarkers measured on diverse platforms or in different media, e.g. urine and serum, to clinical practice to enable a personalised approach to the diagnosis, prevention and treatment of diastolic LV dysfunction, a high-risk condition [7] that affects 25% of the general population [6]. Furthermore, our current study highlights research tracks to be pursued in the future, such as showing parallelism between concomitant changes in diastolic LV function and biomarkers in longitudinal studies and proving concordance in the proteomic profiles of urine, serum or plasma and the myocardium.
Supporting Information
S1 Table. Glossary of echocardiographic measurement reflecting diastolic left ventricular function
https://doi.org/10.1371/journal.pone.0167582.s001
(DOC)
S3 Table. Covariables selected by stepwise regression
https://doi.org/10.1371/journal.pone.0167582.s003
(DOC)
S4 Table. Multivariable-adjusted associations of tissue Doppler indexes with urinary peptides
https://doi.org/10.1371/journal.pone.0167582.s004
(DOC)
S5 Table. Urinary biomarkers by category of diastolic LV function
https://doi.org/10.1371/journal.pone.0167582.s005
(DOC)
S1 Fig. –Log10(p) probability plot of the multivariable-adjusted associations of various indexes of diastolic left ventricular function with the urinary peptides
https://doi.org/10.1371/journal.pone.0167582.s006
(DOC)
S2 Fig. Circulating biomarkers of collagen turnover
https://doi.org/10.1371/journal.pone.0167582.s007
(DOC)
Acknowledgments
The authors gratefully acknowledge the clerical assistance of V De Leebeeck and R Wolfs and the technical support of L Custers, MJ Jehoul, D Thijs and H Truyens in data collection.
Author Contributions
- Conceptualization: HM JAS.
- Data curation: ZYZ MP LT LJ T Kuznetsova JAS.
- Formal analysis: ZYZ MP LT LJ JAS.
- Funding acquisition: JAS.
- Investigation: ZYZ SR WYY TP MP PZ BL FFW AG T Kuznetsova JD.
- Methodology: SR BL CP AG T Kuznetsova JD HM JAS.
- Project administration: JAS.
- Resources: JAS.
- Software: LT LJ HM.
- Supervision: JUV PV JD HM JAS.
- Visualization: ZYZ.
- Writing – original draft: ZYZ JAS.
- Writing – review & editing: ZYZ SR WYY TP MP PZ BL FFW CP LT LJ AG T Koeck CD JUV PV T Kuznetsova JD HM JAS.
References
- 1. Hess OM, Schneider J, Koch R, Bamert C, Grimm J, Krayenbuehl HP. Diastolic function and myocardial structure in patients with myocardial hypertrophy. Circulation. 1981; 63: 360–371. pmid:6450003
- 2. Weber KT. Cardiac interstitium in health and disease: the fibrillar collagen network. J Am Coll Cardiol. 1989; 13: 1637–1652. pmid:2656824
- 3. López B, Ravassa S, González A, Zubillaga E, Bonavila C, Bergés M, et al. Myocardial collagen cross-linking is associated with heart failure hospitalization in patients with hypertensive heart failure. J Am Coll Cardiol. 2016; 67: 251–260. pmid:26796388
- 4. Pauschinger M, Knopf D, Petschauer S, Doerner A, Poller W, Schwimmbeck PL, et al. Dilated cardiomyopathy is associated with significant changes in collagen type I/III ratio. Circulation. 1999; 99: 2750–2756. pmid:10351968
- 5. Westermann D, Lindner D, Kasner M, Zietsch C, Sawatis K, Escher F, et al. Cardiac inflammation contributes to changes in the extracellular matrix in patients with heart failure and normal ejection fraction. Circ Heart Fail. 2011; 4: 44–52. pmid:21075869
- 6. Kuznetsova T, Herbots L, López B, Jin Y, Richart T, Thijs L, et al. Prevalence of left ventricular diastolic dysfunction in a general population. Circ Heart Fail. 2009; 2: 105–112. pmid:19808325
- 7. Kuznetsova T, Thijs L, Knez J, Cauwenberghs N, Petit T, Gu YM, et al. Longitudinal changes in left ventricular diastolic dysfunction in a general population. Circ Cardiovasc Imaging. 2015; 8: e002882. pmid:25873723
- 8. Mischak H, Vlahou A, Ioannidis JP. Technical aspects and inter-laboratory variability in native peptide profiling: the CE-MS experience. Clin Biochem. 2013; 46: 432–443. pmid:23041249
- 9. Kuznetsova T, Mischak H, Mullen W, Staessen JA. Urinary proteome analysis in hypertensive patients with left ventricular diastolic dysfunction. Eur Heart J. 2012; 33: 2342–2350. pmid:22789915
- 10. Zhang Z, Staessen JA, Thijs L, Gu Y, Liu Y, Jacobs L, et al. Left ventricular diastolic dysfunction in relation to the urinary proteome: a proof-of-concept study in a general population. Int J Cardiol. 2014; 176: 158–165. pmid:25065337
- 11. Zhang ZY, Thijs L, Petit T, Gu YM, Jacobs L, Yang WY, et al. The urinary proteome and systolic blood pressure as predictors of 5-year cardiovascular and cardiac outcomes in a general population. Hypertension. 2015; 66: 52–60. pmid:26063667
- 12. Velagaleti RS, Gona P, Sundstrom J, Larson MG, Siwik D, Colucci WS, et al. Relations of biomarkers of extracellular matrix remodeling to incident cardiovascular events and mortality. Arterioscler Thromb Vasc Biol. 2012; 30: 2283–2288.
- 13. World Medical Association declaration of Helsinki. Recommendations guiding physicians in biomedical research involving human subjects. JAMA. 1997; 277: 925–926. pmid:9062334
- 14. Staessen JA, Wang JG, Brand E, Barlassina C, Birkenhäger WH, Herrmann SM, et al. Effects of three candidate genes on prevalence and incidence of hypertension in a Caucasian population. J Hypertens. 2011; 19: 1349–1358.
- 15. Li Y, Zagato L, Kuznetsova T, Tripodi G, Zerbini G, Richart T, et al. Angiotensin-converting enzyme I/D and α-adducin Gly460Trp polymorphisms. From angiotensin-converting enzyme activity to cardiovascular outcome. Hypertension. 2007; 49: 1291–1297. pmid:17452507
- 16. Wei FF, Drummen NEA, Schutte AE, Thijs L, Jacobs L, Petit T, et al. Vitamin K dependent protection of renal function in multi-ethnic population studies. EBioMed. 2016; 4: 162–169.
- 17. Yang WY, Petit T, Thijs L, Zhang ZY, Jacobs L, Hara A, et al. Coronary risk in relation to genetic variation in MEOX2 and TCF15 in a Flemish population. BMC Genet. 2015; 16: 116. pmid:26428460
- 18. Kuznetsova T, Codd V, Brouilette S, Thijs L, González A, Jin Y, et al. Association between left ventricular mass and telomere length in a population study. Am J Epidemiol. 2010; 172: 440–450. pmid:20660518
- 19. Gottdiener JS, Bednarz J, Devereux R, Gardin J, Klein A, Manning WJ, et al. American Society of Echocardiography recommendations for use of echocardiography in clinical trials. A report from the American Society of Echocardiography's Guidelines and Standard Committee and the Task Force on Echocardiography in Clinical Trials. J Am Soc Echocardiogr. 2004; 17: 1086–1119. pmid:15452478
- 20. Jha V, Garcia-Garcia G, Iseki K, Li Z, Naicker S, Plattner B, et al. Chronic kidney disease: global dimension and perspectives. Lancet. 2013; 382: 260–272. pmid:23727169
- 21. Kloch-Badelek M, Kuznetsova T, Sakiewicz W, Tikhonoff V, Ryabikov A, González A, et al. Prevalence of diastolic left ventricular dysfunction in European populations based on cross-validated diagnostic thresholds. Cardiovascular Ultrasound. 2012; 10: 10. pmid:22429658
- 22. Jantos-Siwy J, Schiffer E, Brand K, Schumann G, Rossing K, Delles C, et al. Quantitative urinary proteome analysis for biomarker evaluation in chronic kidney disease. J Proteome Res. 2009; 8: 268–281. pmid:19012428
- 23. Mischak H, Kolch W, Aivalotis M, Bouyssie D, Court M, Dihazi H, et al. Comprehensive human urine standards for comparability and standardization in clinical proteome analysis. Proteomics Clin Appl. 2010; 4: 464–478. pmid:21137064
- 24. Klein J, Papadopoulos T, Mischak H, Mullen W. Comparison of CE-MS/MS and LC-MS/MS sequencing demonstrates significant complementarity in natural peptide identification in human urine. Electrophoresis. 2014; 35: 1060–1064. pmid:24254231
- 25. Zürbig P, Renfrow MB, Schiffer E, Novak J, Walden M, Wittke S, et al. Biomarker discovery by CE-MS enables sequence analysis via MS/MS with platform-independent separation. Electrophoresis. 2006; 27: 2111–2125. pmid:16645980
- 26. Meleth S, Deshane J, Kim H. The case for well-conducted experiments to validate statistical protocols for 2D gels: different pre-processing = different lists of significant proteins. BMC Biotechnol. 2005; 5: 7. pmid:15707480
- 27. Lazar C, Gatto L, Ferro M, Bruley C, Burger T. Accounting for the multiple natures of missing values in label-free quantitative proteomics data sets to compare imputation strategies. J Proteome Res. 2016; 15: 1116–1125. pmid:26906401
- 28. O'Brien E, Asmar R, Beilin L, Imai Y, Mallion JM, Mancia G, et al. European Society of Hypertension recommendations for conventional, ambulatory and home blood pressure measurement. J Hypertens. 2003; 21: 821–848. pmid:12714851
- 29. Levey AS, Stevens LA, Schmid CH, Zhang Y, Castro AF III, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009; 150: 604–612. pmid:19414839
- 30. National Kidney Foundation. K/DOQI Clinical Practice Guidelines for Chronic Kidney Disease: evaluation, classification, and stratification. Am J Kidney Dis. 2002; 39: S1–S266. pmid:11904577
- 31. Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Report of the expert committee on the diagnosis and classification of diabetes mellitus. Diabet Care. 2003; 26 (suppl 1): S5–S20.
- 32. Lapiere MC, Nusgens B, Pierard GE. Interaction between collagen type I and type III in conditioning bundles organization. Connect Tissue Res. 1977; 5: 21–29. pmid:141359
- 33. López B, González A, Querejeta R, Larman M, Díez J. Alterations in the pattern of collagen deposition may contribute to the deterioration of systolic function in hypertensive patients with heart failure. J Am Coll Cardiol. 2006; 48: 89–96. pmid:16814653
- 34. Kasner M, Westermann D, Lopez B, Gaub R, Escher F, Kühl U, et al. Diastolic tissue Doppler indexes correlate with the degree of collagen expression and cross-linking in heart failure and normal ejection fraction. J Am Coll Cardiol. 2011; 57: 977–985. pmid:21329845
- 35. Aoki T, Fukumoto Y, Sugimura K, Oikawa M, Satoh K, Nakano M, et al. Prognostic impact of myocardial interstitial fibrosis in non-ischemic heart failure. Comparison between preserved and reduced ejection fraction heart failure. Circ J. 2011; 75: 2605–2613. pmid:21821961
- 36. Zile MR, Baicu CF, Ikonomidis JS, Stroud RE, Nietert PJ, Bradshaw AD, et al. Myocardial stiffness in patients with heart failure and a preserved ejection fraction: contributions of collagen and titin. Circulation. 2015; 131: 1247–1259. pmid:25637629
- 37. Wei S, Chow LT, Shum IO, Qin L, Sanderson JE. Left and right ventricular collagen type I/III ratios and remodeling post-myocardial infarction. J Card Fail. 1999; 5: 117–126. pmid:10404351
- 38. Cianciulli TF, Saccheri MC, Lax JA, Bermann AM, Ferreiro DE. Two-dimensional speckle tracking echocardiography for the assessment of atrial function. World J Cardiol. 2010; 2: 163–170. pmid:21160748
- 39. Redfield MM, Jacobsen SJ, Burnett JC Jr, Mahoney DW, Bailey KR, Rodeheffer RJ. Burden of systolic and diastolic ventricular dysfunction in the community. Appreciating the scope of the heart failure epidemic. JAMA. 2003; 289: 194–202. pmid:12517230
- 40. Gramme B, Böhm J, Dufour A, Benz M, Lange R, Bauernschmitt R. Atrial fibrosis in heart surgery patients. Decreased collagen III/I ratio in postoperative atrial fibrillation. Basic Res Cardiol. 2005; 100: 288–294. pmid:15690103
- 41. Longobardo L, Todaro MC, Zito C, Piccione MC, Di Bella G, Oreto L, et al. Role of imaging in assessment of atrial fibrosis in patients with atrial fibrillation: state-of-the-art review. Eur Heart J Cardiovasc Imaging. 2014; 15: 1–5. pmid:23798579
- 42. López B, González A, Ravassa S, Beaumont J, Moreno MU, San JG, et al. Circulating biomarkers of myocardial fibrosis: the need for a reappraisal. J Am Coll Cardiol. 2015; 65: 2449–2456. pmid:26046739
- 43. Querejeta R, López B, González A, Sánchez E, Larman M, Martínez Ubago JL, et al. Increased collagen type I synthesis in patients with heart failure of hypertensive origin. Relation to myocardial fibrosis. Circulation. 2004; 110: 1263–1268. pmid:15313958
- 44. López B, Querejeta R, González A, Larman M, Díez J. Collagen cross-linking but not collagen amount associates with elevated filling pressures in hypertensive patients with stage C heart failure. Potential role of lysyl oxidase. Hypertension. 2012; 60: 677–683. pmid:22824984
- 45. López B, Querejeta R, Varo N, González A, Larman M, Martínez Ubago JL, et al. Usefulness of serum carboxy-terminal propeptide of procollagen type I in assessment of the cardioreparative ability of antihypertensive treatment in hypertensive patients. Circulation. 2001; 104: 286–291. pmid:11457746
- 46. Izawa H, Murohara T, Nagata K, Isobe S, Asano H, Amano T, et al. Mineralocorticoid receptor antagonism ameliorates left ventricular diastolic dysfunction and myocardial fibrosis in mildly symptomatic patients with idiopathic dilated cardiomyopathy: a pilot study. Circulation. 2005; 112: 2940–2945. pmid:16275882
- 47. Klappacher G, Franzen P, Haab D, Mehrabi M, Binder M, Plesch K, et al. Measuring extracellular matrix turnover in the serum of patients with idiopathic or ischemic dilated cardiomyopathy and impact on diagnosis and prognosis. Am J Cardiol. 1995; 75: 913–918. pmid:7733000
- 48. González A, López B, Querejeta R, Zubillaga E, Echeverria T, Díez J. Filling pressures and collagen metabolism in hypertensive patients with heart failure and normal ejection fraction. Hypertension. 2010; 55: 1418–1424. pmid:20404218
- 49. Bishop JE, Greenbaum R, Gibson DG, Yacoub M, Laurent GJ. Enhanced deposition of predominatly type I collagen in myocardial disease. J Mol Cell Cardiol. 1990; 22: 1157–1165. pmid:2095438