Abstract
Right ventricular systolic dysfunction (RVSD) at baseline (pre-treatment) predicts early death in patients with pulmonary hypertension (PH). However, RVSD can only be detected reliably by prohibitively invasive or expensive techniques. N-terminal B-type natriuretic peptide concentration ([NT-proBNP]) correlates with RV function in PH; however, an [NT-proBNP] threshold that indicates RVSD in individual patients has not previously been determined.
Twenty-five patients with PH (pulmonary arterial hypertension (n = 19) or chronic thromboembolic PH (n = 6)) underwent cardiovascular magnetic resonance (CMR) imaging and NT-proBNP measurement at baseline. [NT-proBNP] was correlated against RV dimensions and ejection fraction (RVEF) measured directly by CMR imaging. The ability of NT-proBNP to detect RVSD (defined as a CMR-derived RVEF >2 sds below control values) was tested and predictors of [NT-proBNP] identified.
[NT-proBNP] correlated negatively with RVEF. RVSD was present in nine out of 25 patients. An [NT-proBNP] threshold of 1,685 pg·mL−1 was sensitive (100%) and specific (94%) in detecting RVSD. RVEF and RV mass index independently predicted [NT-proBNP].
In pulmonary hypertension, a baseline N-terminal B-type natriuretic peptide concentration of >1,685 ng·L−1 suggests right ventricular systolic dysfunction, and thus an increased risk of early death. N-terminal B-type natriuretic peptide could prove useful as an objective, noninvasive means of identifying patients with pulmonary hypertension who have right ventricular systolic dysfunction at presentation.
Pulmonary hypertension (PH) is a rare condition characterised by increased resistance to pulmonary blood flow, progressive right ventricular (RV) systolic dysfunction (RVSD) and early death. Without specific treatment, up to 50% of patients diagnosed with PH will die within 2 yrs, usually as a result of RVSD 1. In earlier studies that identified RVSD either by invasive haemodynamic measurements 1–5 or by the detection of overt clinical signs (i.e. patients within class IV of the World Health Organization classification of disease severity) 1, evidence of RVSD at baseline (i.e. pre-treatment) was associated with an increased mortality risk.
Accurate assessment of baseline RV function is therefore desirable in all patients with PH. This information can help identify individuals in whom prompt invasive assessment and the initiation of disease-targeted therapy, is appropriate. Early RV dysfunction might be difficult to detect on clinical examination, therefore transthoracic echocardiography is often used for this purpose 6. Although inexpensive and widely available, RV ejection fraction (RVEF) can only be estimated indirectly from measurements acquired in two dimensions by echocardiography. Variables derived by this method that reflect RV function and provide useful prognostic information in PH are either complicated to define and therefore reliant upon an experienced operator (e.g. the Tei index) 7, or most frequently detected in patients in whom RV failure is already clinically overt (e.g. a pericardial effusion) 2.
Cardiovascular magnetic resonance (CMR) imaging has recently emerged as the new gold standard for assessment of the RV in PH and other diseases 8, 9; it allows image acquisition in any plane and at almost any angle, allowing direct measurement of RV volumes and RVEF. At present, CMR imaging remains too expensive for routine clinical use and is therefore not available in many PH centres. An alternative noninvasive method of identifying baseline RVSD in PH that is objective, operator independent and widely available would, therefore, be clinically useful.
B-type natriuretic peptide (BNP) and N-terminal-proBNP (NT-proBNP) are released from both cardiac ventricles in response to raised chamber pressure or volume overload 10. Previous studies have shown that circulating BNP concentration correlates with invasive pulmonary haemodynamics 11, 12, RVEF measured by computed tomography (CT) 12, survival 13, 14, and the response to therapy 15, 16 in PH patients. NT-proBNP, which is biologically inactive, is more stable in vivo and in vitro than BNP 17 and can therefore be used instead of BNP as a cardiovascular biomarker in various disease states. Circulating NT-proBNP concentration ([NT-proBNP]) has recently been shown to correlate with invasive pulmonary haemodynamics 18, survival- and echocardiography-derived measures of RV function in PH 19. However, a threshold [NT-proBNP] indicating that RVSD is present in an individual patient has yet to be determined, perhaps because of insufficient precision in the previous reference method of RV systolic function measurement (echocardiography).
The present authors aimed to relate [NT-proBNP] to RV volumes and systolic function measured directly by CMR imaging and to identify an [NT-proBNP] threshold for the detection of RVSD in patients with PH at baseline.
METHODS
Patients
A prospective cross-sectional study was performed. Thirty consecutive patients were recruited who attended for diagnosis and baseline assessment at the Scottish Pulmonary Vascular Unit (Glasgow, UK). All subjects underwent routine clinical assessment, including right heart catheterisation to establish a diagnosis of PH. The following were used as inclusion criteria for the study: 1) a diagnosis of either pulmonary arterial hypertension (PAH) or chronic thromboembolic PH (CTEPH) reached by conventional diagnostic methods 6; and 2) informed written consent to a study protocol approved by the local institutional review committee.
Patients were excluded from the study (after review of their past medical history, admission blood tests and transthoracic echocardiogram and before CMR imaging and NT-proBNP measurement) if they met any of the following exclusion criteria. 1) Evidence of congenital heart disease (CHD), either at echocardiography or as a previous diagnosis. Patients with CHD were excluded because differing patterns of natriuretic peptide release have been demonstrated in conditions that cause RV volume overload, such as CHD, in comparison with those that cause RV pressure overload (e.g. PAH and CTEPH) 20. 2) A history of any other condition that might cause elevation of [NT-proBNP]. These include ischaemic heart disease 21; myocardial infarction 22; left ventricular (LV) systolic dysfunction 23; systemic hypertension 24; left-sided valvular heart disease 25; chronic renal impairment 26 or an abnormal serum creatinine on admission to the Scottish Pulmonary Vascular Unit (normal was defined as <120 μmol·L−1); and diabetes mellitus 27 or a fasting blood glucose on admission to the Scottish Pulmonary Vascular Unit of >6.1 mmol·L−1. 3) CMR-related exclusion criteria included an indwelling cardiac pacemaker and claustrophobia.
On the basis of these criteria, five out of 30 patients were excluded from the final analysis (atrial septal defect (two out of five), chronic renal impairment (one out of five), significant mitral regurgitation (one out of five) and systemic hypertension (one out of five)). The remaining 25 subjects underwent CMR imaging and venous blood sampling for [NT-proBNP] within 24 h of invasive assessment. The final diagnoses reached in these 25 patients were as follows: idiopathic PAH in 11 out of 25, CTEPH in six out of 25 and PAH associated with connective tissue disease (PAH-CTD) in eight out of 25. At the time of the study, no patient had received any disease-targeted therapy (i.e. prostacyclin analogues, bosentan, sildenafil or calcium-channel blockers) for PH.
Twelve control subjects, with no history of any cardiorespiratory disease, underwent CMR imaging to provide local control values for ventricular dimensions and function. All control subjects were staff of the Western Infirmary, Glasgow, UK, and were subject to the same exclusion criteria as PH patients. Controls were approached individually and asked to participate. All controls gave informed written consent but did not undergo right heart catheterisation or NT-proBNP measurement.
There were no significant differences between PH patients and controls in mean±sd age (PH: 60±12, controls: 57±12, p = 0.455) or mean systemic blood pressure (PH: 95±13, controls: 92±17, p = 0.536). PH (17 males and eight females) and control populations (eight males and four females) were well sex matched.
Venous blood sampling and measurement of NT-proBNP
After ≥20 min of recumbent rest, 5 mL of venous blood was drawn from the antecubital fossa of all 25 PH patients immediately prior to CMR imaging. Blood samples were drawn into vacuum-sealed containers containing ethylenediamine tetraacetic acid and immediately centrifuged using a PK110 centrifuge (ALC, Winchester, VA, USA). Supernatant was then removed and individual samples were coded, labelled and stored at -80°C until analysis.
NT-proBNP analysis was performed on an ELECSYS 2010 bench top analyser utilising a chemiluminescent assay with a coefficient of variation <5% (Roche Diagnostics, Lewes, UK). All biochemical analyses were performed by a single operator (J.J. Morton) who was blinded to haemodynamic, CMR and clinical results. NT-proBNP measurements are subsequently quoted as plasma concentration (ng·L−1).
CMR image acquisition
CMR imaging was undertaken on a 1.5 T MRI scanner (Sonata Magnetom; Siemens, Erlangen, Germany) using a protocol recently described in detail 28. An identical CMR protocol was followed in the control subjects. Fast imaging with steady-state free precession cines (TrueFISP; Siemens) were utilised throughout. Methodological details that are of particular importance include standardised imaging parameters of repetition time, echo time, flip angle, voxel size, field of view = 3.14 ms, 1.6 ms, 60°, 2.2×1.3×8.0 mm, 340 mm; 8-mm imaging slices were used with a 2-mm interslice gap. Short-axis imaging was initiated at the atrioventricular valve plane identified on a horizontal long-axis view of the heart, and then propagated sequentially to the cardiac apex, providing complete coverage of left and right ventricles.
CMR image analysis
CMR images in PH patients and controls were analysed by a single operator (K.G. Blyth) using the Argus analysis software (Siemens). Individual scans were coded by number and analysed in batches by K.G. Blyth, who was blinded to the identity and haemodynamic results of any given subject at the time of analysis. CMR analysis was performed before blood samples were analysed for NT-proBNP levels. RV and LV volumes (RV end-diastolic, RV end-systolic, LV end-diastolic and LV end-systolic volumes) were determined by manual planimetry of selected short-axis images, as described previously 29, 30. RV and LV ventricular stroke volume (RVSV and LVSV), RVEF and LV ejection fraction, and RV and LV mass (RVM and LVM) were determined as previously described 29, 30. RVM was determined as RV free wall mass. The interventricular septum was considered part of the LV. RVM index was determined as RVM/LVM 31. LV and RV volumes were corrected for body surface area and reported as indexed measurements (RV end-diastolic, RV end-systolic, LV end-diastolic and LV end-systolic volume indices).
Clinical assessment and right heart catheterisation
A standard diagnostic algorithm 6 was employed in the clinical assessment of all PH patients studied. This included right heart catheterisation using established techniques 6, isotope perfusion lung scanning and CT pulmonary angiography in all patients. If CTEPH was suspected following these initial investigations, selective pulmonary angiography was performed in the catheterisation laboratory (CTEPH was confirmed using this algorithm in six out of 25 PH patents).
Statistical analysis
For all variables, a normal distribution was verified using histograms and Kolmogorov–Smirnov tests. Non-normally distributed values ([NT-proBNP] and right atrial pressure) were logarithmically transformed and all correlation analyses were performed using Pearson's method. Comparisons between PH patients and controls (and between subgroups of PH patients) were made using an independent samples (unpaired) t-test (using Welch's correction if appropriate).
Contingency tables (2×2) were used to calculate the sensitivity and specificity of various plasma concentrations of NT-proBNP as a means of detecting RVSD. RVSD was defined as an RVEF <42% (this constituted an RVEF >2 sds below the mean RVEF of the control population (66±7%)). The plasma [NT-proBNP] with the greatest diagnostic accuracy (that which produced the lowest number of false-positive and false-negative results) was chosen as a “threshold” value for the detection of RVSD. Negative and positive predictive values for results either side of this NT-proBNP threshold were then calculated. The implications of an above-threshold result were quantified by likelihood ratios and Fisher's exact test.
To identify independent predictor(s) of circulating [NT-proBNP], variables significantly associated with [NT-proBNP] (RV end-diastolic volume and RV end-systolic volume indices, RVEF, RVM index, LVSV index, mean pulmonary arterial pressure (Ppa), systolic Ppa, diastolic Ppa, cardiac index, pulmonary vascular resistance and mixed venous oxygen saturation (Sv,O2)) were entered into a multivariable stepwise linear regression model, in which NT-proBNP was the dependent variable. Covariables with a correlation coefficient ≥0.8 were tested within separate regression models to avoid the effects of co-linearity. All statistical analyses were two-sided and assumed a significance level of 5%. All values are presented as mean±1 sd, unless otherwise stated.
RESULTS
Right heart catheterisation
Invasive measurements acquired in the PH patients are summarised in table 1⇓.
Results of right heart catheterisation in 25 patients with pulmonary hypertension
CMR imaging
Ventricular dimensions and systolic function of the 25 PH patients and 12 controls are presented in table 2⇓. In summary, RV end-diastolic volume, RV end-systolic volume, RVM and RVM index were significantly increased, and RVEF was significantly depressed in PH patients relative to control subjects. LV end-diastolic volume and LVSV were significantly lower in PH patients compared with controls. All LV measurements, in both PH patients and controls, were within previously published normal limits 29, 32. Where available, CMR normal ranges determined by steady-state free precession imaging sequences similar to those employed in an earlier study have been used 32.
Ventricular dimensions and function by cardiovascular magnetic resonance imaging in patients with pulmonary hypertension(PH) and controls
Correlation analyses
[NT-proBNP] was non-normally distributed in the population studied. Median (interquartile range) [NT-proBNP] was 669 (94–2615) ng·L−1. After logarithmic transformation, mean±1 sd log10 [NT-proBNP] was 2.71±0.77. A complete list of correlation coefficients and accompanying p-values for variables compared with [NT-proBNP] is presented in table 3⇓. There was a powerful negative logarithmic correlation between [NT-proBNP] and RVEF (fig. 1⇓). With the exception of LVSV, there was no significant correlation between [NT-proBNP] and LV measurements (see Discussion).
Cardiovascular magnetic resonance imaging and N-terminal B-type natriuretic peptide (NT-proBNP) measurements were acquired in 25 patients with pulmonary hypertension. After logarithmic transformation, NT-proBNP concentration ([NT-proBNP]) correlated powerfully with right ventricular ejection fraction (RVEF) (r = -0.66, p<0.001). – – –: the [NT-proBNP] threshold chosen (1,685 ng·L−1) for the detection of RV systolic dysfunction (RVSD); ·········: the RVEF (42%) used to define RVSD. Only one false-negative result (▴) was returned.
Correlation coefficients between log10 [NT-proBNP] (ng·L−1) and measurements acquired in 25 pulmonary hypertension patients at cardiovascular magnetic resonance imaging and right heart catheterisation
Differences in log10 [NT-proBNP] in subgroups of PH patients
There were no statistically significant differences in [NT-proBNP] based on sex (mean log10 [NT-proBNP] was 2.7±0.6 in males and 2.7±0.8 in females) or aetiology of PH (mean log10 [NT-proBNP] was 2.6±0.8 in PAH, 2.4±0.9 in idiopathic PAH, 3.0±0.7 in PAH-CTD and 3.0±0.5 in patients with CTEPH).
Right ventricular systolic dysfunction
RVSD was present in nine out of 25 PH patients (36%). The distribution of [NT-proBNP] values for the nine patients with RVSD and the 16 patients without RVSD were normal. Mean [NT-proBNP] was significantly higher in patients with RVSD (4,127±2,951 ng·L−1) in comparison with those without RVSD (354±435 ng·L−1; t = 3.813, degrees of freedom = 8, p = 0.005).
Diagnostic performance of NT-proBNP
NT-proBNP performed best at a plasma concentration of 1,685 ng·L−1, detecting RVSD with 100% sensitivity (95% confidence interval (CI) 63.0–100%) and 94% specificity (95% CI 71.3–99.9%). A value of 1,685 ng·L−1 was therefore chosen as the threshold value for RVSD detection.
Seventeen out of 25 patients had an [NT-proBNP] below the 1,685 ng·L−1 threshold. Sixteen out of 17 of these patients were true negatives (i.e. RVEF was normal (≥42%)) and one was a false negative. In the false-negative subject, RVEF measured 34% and [NT-proBNP] was 669 ng·L−1. Thus, the negative predictive value of a low [NT-proBNP] was 94%.
Eight out of 25 patients had an [NT-proBNP] >1,685 ng·L−1. RVSD was demonstrated by CMR imaging (RVEF <42%) in all of these individuals. Thus, the positive predictive value of a high [NT-proBNP] was 100%. These results are presented graphically in figure 1⇑. Based on likelihood ratio testing, individuals with an [NT-proBNP] >1,685 ng·L−1 were 17 times more likely to have RVSD than those with below threshold results.
Predictors of [NT-proBNP]
Since RVEF correlated powerfully with RVM index (r = -0.854, p<0.001), these variables were entered into two separate multivariable linear regression models in which [NT-proBNP] was the dependent variable. In these models, RVEF (odds ratio (OR) = -0.035, 95% CI -0.052– -0.018, p<0.001) and RVM index (OR = 2.254, 95% CI 1.403–3.106, p<0.001) proved to be the only independent predictors of circulating [NT-proBNP].
DISCUSSION
In the current study, a close correlation between RVEF, measured directly by CMR imaging, and [NT-proBNP], is reported in patients with PH, assessed before the initiation of any disease-targeted therapy. A previous study by Nagaya et al. 12 demonstrated a similar relationship between RVEF, measured by electron beam CT imaging and BNP concentration, in a similar cohort of patients. The present study is of additional interest to this previous report for three important and clinically relevant reasons.
1) In routine clinical practice, most centres use NT-proBNP, not BNP, to aid clinical decision making in heart failure and a variety of other conditions. As previously discussed, this is because NT-proBNP is more stable both in vivo and in vitro than the active BNP molecule itself 17.
2) The study by Nagaya et al. 12 used electron beam CT (EBCT) imaging to measure RVEF, while the present authors have utilised CMR imaging as their reference method. Temporal resolution is better with CMR imaging 33, without the need for ionising radiation or intravenous contrast 34. As a result, CMR has emerged over recent years as the gold standard for detailed study of the right ventricle 9, 35–37, and has become an established modality for the physiological assessment of PH patients in cross-sectional studies 31, 38–40, longitudinal follow-up studies 41 and clinical trials of therapy 16. EBCT, in contrast, is not routinely used to assess PH patients. The current study is therefore useful, as it defines the relationship between [NT-proBNP] and RV volumes and ejection fraction measured directly by CMR imaging.
3) In addition to a correlation between [NT-proBNP] and RVEF, the present authors found that an [NT-proBNP] >1,685 ng·L−1 (detected in eight out of 25 patients) proved a sensitive marker of RVSD (100% sensitivity). Below this threshold, only one false-negative result was returned, suggesting that RVSD is very unlikely when circulating [NT-proBNP] is <1,685 ng·L−1, at least in patients assessed at baseline, before treatment initiation (negative predictive value 94%). Stratifying patients on either side of this NT-proBNP threshold may prove useful as an operator-independent, objective means of classifying individuals at presentation based upon their probability of having RVSD and, thus, an increased risk of early death.
Physiological determinants of [NT-proBNP]
Several variables correlated closely with circulating [NT-proBNP] in the present study (table 3⇑). Of these variables, only RVEF and RVM index, analysed within separate regression models to avoid the effects of colinearity, proved independent predictors of [NT-proBNP]. Although the RVM index relates linearly to mean Ppa in PH 31, it has not been proven to predict early death in patients with this condition. For this reason, the current authors sought to detect RVSD, not RV hypertrophy, using [NT-proBNP], however the influence of RV hypertrophy (as described by RVM index) should be acknowledged as a potential confounder of the present study’s results. The vast majority (seven out of eight (88%)) of patients in the present study with above-threshold [NT-proBNP] levels, indicative of RVSD, also had RV hypertrophy (defined as an RVM >84 g (>2 sds above the mean RV mass of the control population (48±18 g)), reflecting the covariation between RVM index, RVEF and [NT-proBNP] that was found. Despite this, RVEF appeared to be an important determinant of circulating [NT-proBNP] in the current study.
Clinical implications
A recent study reported by Quaife et al. 42 demonstrated a close correlation between RVEF and RV wall stress, which is the principal stimulus to NT-proBNP release in PH. This relationship may explain the usefulness of NT-proBNP as a means of detecting RVSD reported herein. In interpreting the present results, it is worth considering the findings of Fijalkowska et al. 19, who recently examined the prognostic significance of NT-proBNP in PH. Nonsurvivors in this study were identified by an [NT-proBNP] >1,400 ng·L−1. Since nonsurvivors with PH are likely to die from RV failure, and [NT-proBNP] correlated powerfully with echocardiographic measures of RV systolic function in the study by Fijalkowska et al. 19, it seems plausible that the 1,400 ng·L−1 threshold reported by Fijalkowska et al. 19 and the present study’s threshold of 1,685 ng·L−1 reflect genuine RVSD in the patients studied.
Clearly, the true clinical significance of the [NT-proBNP] threshold for RVSD (1,685 ng·L−1) identified in the present study can only be determined in a longitudinal survival analysis. For this reason, it is difficult to justify detailed proposals for the future use of NT-proBNP in clinical practice at this time. However, our findings indicate that NT-proBNP, measured during a patient's first consultation at a specialist PH centre, could, in conjunction with existing means of assessing RV function (i.e. clinical examination and echocardiography), help to identify patients with prognostically significant RVSD. Assuming other reasons for a high NT-proBNP were excluded, patients with a high [NT-proBNP] could then be invasively assessed at an early stage and treatment initiated if appropriate. Although there is as yet no evidence that early treatment of PH alters prognosis per se, it is known that a history of overt RV failure before treatment initiation is associated with an increased mortality risk 43. Prompt initiation of therapy in patients with high NT-proBNP levels and subclinical RVSD may, therefore, avoid the development of overt RV failure and reduce early mortality.
RV ejection fraction
In the current study, RVEF was utilised to define RVSD. Although historical data have suggested that RVEF does not reflect RV systolic function alone (since it is also dependant upon right ventricular afterload) 44–46, RVEF determined by echocardiography and radionuclide ventriculography has recently been shown to predict adverse outcomes (cardiac death or deterioration) in patients with PH 47. However, until longitudinal survival studies are reported in which CMR-derived variables are assessed, the true prognostic significance of CMR-derived RVEF, and the cut-off value of 42% used here to define RVSD cannot be known.
NT-proBNP and LV stroke volume
A statistically significant relationship between [NT-proBNP] and LVSV was found, but there was no correlation between [NT-proBNP] and RVSV. This apparent contradiction reflects the inaccuracy of CMR planimetry in determining true (i.e. forward) stroke volume in patients with significant valvular regurgitation. CMR planimetry of a dilated RV with significant tricuspid valve regurgitation will always overestimate true forward RVSV because a proportion of the ejected blood volume will move backwards through the incompetent tricuspid valve. Since forward LVSV and forward RVSV will, over several heart beats, be equivalent in the absence of any communicating shunt, and since no PH patient in the present study had any evidence of mitral regurgitation, the statistically significant relationship that was found between LVSV and [NT-proBNP] is likely to reflect an underlying correlation between [NT-proBNP] and true forward RVSV that is obscured by this predictable measurement error. Although CMR planimetry may not always accurately determine forward RVSV in PH patents, RVEF still reflects the systolic contractile function of the ventricle, assessment of which was the primary aim of the present study. An alternative confounder to the CMR planimetry results might be the trabeculation of the RV. Despite this concern, acceptable intra-observer variation has been found for ventricular mass and ejection fraction (3.8–7.4 and 3.4–4.4%, respectively) measured by this method in the Glasgow Cardiac Magnetic Resonance Unit (unpublished data).
Study limitations
The relatively small number of patients involved in the present study is of concern in the proposal of a specific cut-off level for use in the detection of RVSD in clinical practice. However, CMR imaging represents an extremely accurate and reproducible gold standard of RVEF, limiting the number of patients required for such studies. [NT-proBNP] can be affected by a number of factors, including renal function, left-sided heart disease, age, sex 48 and recent exertion. In the current study, patients with other reasons for a high [NT-proBNP] were excluded, and the timing and the conditions surrounding NT-proBNP sampling were standardised. A recent study has reported a close relationship between Sv,O2 and [NT-proBNP] in patients with cyanotic congenital heart disease and PH 49; however, Sv,O2 did not predict [NT-proBNP] in the present study, and patients with CHD were specifically excluded. The present authors wish to emphasise that their conclusions can only be applied to selected PH patients assessed in the described manner. In addition, it is not known whether therapies such as sildenafil, which acts directly upon the cardiac myocyte 50, will affect the relationship reported herein between RVEF and [NT-pro BNP], or whether any observed fall in [NT-pro BNP] during treatment will prove proportional to an improvement in RV systolic function. Further insight into these issues is necessary to understand the clinical significance of changes in [NT-proBNP] during treatment.
Conclusion
The present study is the first to utilise a specific N-terminal B-type natriuretic peptide threshold concentration (1,685 ng·L−1) as a means of detecting right ventricular systolic dysfunction in patients with pulmonary hypertension. The method that is described for classifying patients either side of such a threshold is simple, inexpensive and, most importantly, non-invasive and it may help to inform the largely subjective process of stratifying new patients for invasive assessment and the initiation of disease-targeted therapy.
Footnotes
-
For editorial comments see page 622.
- Received July 20, 2006.
- Accepted November 20, 2006.
- © ERS Journals Ltd