Elevated rheumatoid factor and long term risk of rheumatoid arthritis: a prospective cohort studyBMJ 2012; 345 doi: https://doi.org/10.1136/bmj.e5244 (Published 06 September 2012) Cite this as: BMJ 2012;345:e5244
- Sune F Nielsen, senior scientist12,
- Stig E Bojesen, consultant123,
- Peter Schnohr, consultant3,
- Børge G Nordestgaard, professor123
- 1Department of Clinical Biochemistry, 54M1, Herlev Hospital, Copenhagen University Hospital, Herlev Ringvej 75, DK-2730 Herlev, Denmark
- 2Faculty of Health Sciences, University of Copenhagen, Denmark
- 3The Copenhagen City Heart Study, Bispebjerg Hospital, Copenhagen University Hospital
- Correspondence to: B G Nordestgaard
- Accepted 12 July 2012
Objective To test whether elevated concentration of rheumatoid factor is associated with long term development of rheumatoid arthritis.
Design A prospective cohort study, the Copenhagen City Heart Study. Blood was drawn in 1981-83, and participants were followed until 10 August 2010.
Setting Copenhagen general population.
Participants 9712 white Danish individuals from the general population aged 20-100 years without rheumatoid arthritis at study entry.
Main outcome measures Rheumatoid arthritis according to baseline plasma IgM rheumatoid factor level categories of 25-50, 50.1-100, and >100, versus <25 IU/mL.
Results Rheumatoid factor levels were similar from age 20 to 100 years. During 187 659 person years, 183 individuals developed rheumatoid arthritis. In healthy individuals, a doubling in levels of rheumatoid factor was associated with a 3.3-fold (95% confidence interval 2.7 to 4.0) increased risk of developing rheumatoid arthritis, with a similar trend for most other autoimmune rheumatic diseases. The cumulative incidence of rheumatoid arthritis increased with increasing rheumatoid factor category (Ptrend<0.0001). Multivariable adjusted hazard ratios for rheumatoid arthritis were 3.6 (95% confidence interval 1.7 to 7.3) for rheumatoid factor levels of 25-50 IU/mL, 6.0 (3.4 to 10) for 50.1-100 IU/mL, and 26 (15 to 46) for >100 IU/mL, compared with <25 IU/mL (Ptrend<0.0001). The highest absolute 10 year risk of rheumatoid arthritis of 32% was observed in 50-69 years old women who smoked with rheumatoid factor levels >100 IU/mL.
Conclusion Individuals in the general population with elevated rheumatoid factor have up to 26-fold greater long term risk of rheumatoid arthritis, and up to 32% 10 year absolute risk of rheumatoid arthritis. These novel findings may lead to revision of guidelines for early referral to a rheumatologist and early arthritis clinics based on rheumatoid factor testing.
Rheumatoid arthritis is an autoimmune disease affecting 0.5-2% of the population.1 2 Although modern treatments for rheumatoid arthritis can induce remission in many patients, diagnosis of rheumatoid arthritis in early disease stages is important for preventing irreversible damage to the synovial lining and cartilage of diseased joints and for preventing progression into later disease stages.3 4 5 6 7 At present, there is no good clinical available indicator for long term development of rheumatoid arthritis.8 9 10
Rheumatoid factor is an autoantibody targeting the Fc region of IgG antibodies.1 Testing for rheumatoid factor is the most widely used blood test in the classification of rheumatoid arthritis.3 11 In the current classification criteria for rheumatoid arthritis3 “definite rheumatoid arthritis” is based on the confirmed presence of synovitis in at least one joint, absence of an alternate diagnosis better explaining the synovitis, and achievement of a total score of ≥6 (of a possible 10) on a scoring system. The score is derived from four criteria: the number and site of affected joints (range 0-5), serological abnormality (elevated levels of rheumatoid factor or anti-citrullinated protein antibody; range 0-3), elevated acute phase response (range 0-1), and symptom duration (<6 v ≥6 weeks; range 0-1). It is often stated that levels of rheumatoid factor increase with age,1 but convincing data for this statement is difficult to find. About 80% of all patients with rheumatoid arthritis will eventually be seropositive for rheumatoid factor, while only 40% are positive at clinical onset of the disease.10 12 13 However, it is unknown whether elevated levels of rheumatoid factor in individuals in the general population without rheumatoid arthritis is associated with later development of rheumatoid arthritis.
We tested the hypothesis that elevated levels of rheumatoid factor is associated with long term development of rheumatoid arthritis. For this purpose, we measured baseline plasma levels of IgM rheumatoid factor in 9712 white Danish individuals without rheumatoid arthritis from the general population, the Copenhagen City Heart Study, and followed them for up to 28 years, during which time 183 developed rheumatoid arthritis.
The studies were approved by Herlev Hospital, Copenhagen University Hospital, and the Danish ethics committees for Copenhagen and Frederiksberg. Participants provided written informed consent.
The national Danish Civil Registration System records all births, immigrations, emigrations, and deaths in Denmark through the civil registration number, which uniquely identifies all inhabitants of Denmark and provides information on age and sex.14 The national Danish Civil Registration System is 100% complete—that is, for the present study all people were accounted for during the entire follow-up period.
The Copenhagen City Heart Study is a prospective study of a random sample of the Danish general population drawn using the Danish Civil Registration System and initiated in 1976-78.15 16 We studied white participants of Danish descent attending the 1981-83 examination: 19 698 individuals aged 20-100 years were invited, 12 698 (64%) attended, and 9712 (49%) had plasma available for measurement of rheumatoid factor in 2009-10. We excluded 52 individuals with rheumatoid arthritis diagnosed before plasma collection.
Participants filled in a questionnaire that was reviewed together with an investigator on the day of attendance. Subsequently, a physical examination was conducted and blood samples were drawn.
Participants diagnosed clinically with incident rheumatoid arthritis11 from 1 January 1977 to 10 August 2010 were identified by means of the national Danish Patient Registry,17 requiring only a single inpatient or outpatient hospital diagnosis. However, in an attempt to exclude patients with a misdiagnosis of rheumatoid arthritis, we also examined patients with at least two hospitalisations for rheumatoid arthritis at least six weeks apart in accordance with current classification criteria.3 Rheumatoid arthritis was identified using ICD-8 ((international classification of diseases, eighth revision) diagnostic codes 712.1, 712.2, 712.3 (for 1977-94) and ICD-10 codes M05-M06 (for 1995-2010). Rheumatoid arthritis end points were actively reported by hospitals nationwide immediately after discharge of patients from hospitals in the entire period of follow-up.
Other autoimmune rheumatic diseases
Participants diagnosed clinically with other incident autoimmune rheumatic diseases were identified similarly using the national Danish Patient Registry and ICD-8 and ICD-10 codes for 1977-94 and 1995-2010 respectively: systemic lupus erythematosus was coded 734.1 and L93, M32; Sjögren’s syndrome was 734.90 and M32.0; systemic sclerosis was 731.0 and L94; and polymyositis or dermatomyositis was 716.0 and M33.0-M33.1.
Turbidity were used to measure concentrations of rheumatoid factor of IgM type in plasma (measuring range 15-440 IU/mL) (Konelab, Thermo Fischer Scientific, Helsinki, Finland). Plasma samples were drawn in 1981-83 and frozen at −20°C until measurement in 2009-10. Investigators were blinded to rheumatoid arthritis diagnoses when measuring rheumatoid factor, and blinded to rheumatoid factor level when diagnosing rheumatoid arthritis.
Tobacco smoking (cigarettes, pipe, cigarillos, cigars) was the daily amount of cigarettes or equivalent smoked (1 cigarette or equivalent = 20 g tobacco) at the examination, together with the cumulative amount of cigarettes or equivalent smoked up until examination in pack years (1 pack year is 20 cigarettes or equivalent smoked daily for 1 year). Body mass index (weight (kg)/(height (m)2) was measured. Alcohol use was based on self reported number of drinks per week of beer, wine, and spirits (1 drink ≈ 12 g alcohol). We categorised parity for women (0, 1-2, or >2 children), marital status (single, married, separated or divorced, or widowed), and level of education (elementary (1-9 years’ schooling), high school (10-12 years), or academic (>12 years)).
All values are as reported by individuals at study entry in 1981-83, but these covariates were also re-examined in the 1991-94 and 2001-03 examinations and used in the multivariable adjustment as time varying covariates. Information on age and sex was 100% complete, while information on other covariates was 99% complete (1% of participants each missed one or more covariates; see table⇓). Missing information was multivariable imputed before categorisation, thus all statistical analyses were complete for all 9712 participants.
Statistical analyses were performed with Stata 12.1 SE software. All tests were two sided. We used log transformation of rheumatoid factor levels to calculate geometric means for each age group. Plasma levels of rheumatoid factor were illustrated using box plots with 1%, 25%, 50%, 75%, and 99% on the log scale.
Receiver operator characteristics curves of plasma levels of rheumatoid factor at baseline and all future events of rheumatoid arthritis were used to determine an optimal threshold concentration of rheumatoid factor of 25 IU/mL—that is, the level above and below which individuals are best separated into those with and without future rheumatoid arthritis (see supplementary fig 1 in online data supplement). From this threshold we chose categories of rheumatoid factor of <25, 25-50, 50.1-100, and >100 IU/mL—cut-off values were doublings starting from 25 IU/mL. For trend tests across these categories, we used the logarithm of individuals’ rheumatoid factor values on a continuous scale.
First, the association between rheumatoid factor and risk of rheumatoid arthritis and other autoimmune rheumatic diseases were investigated on a continuous scale of doublings of rheumatoid factor as the predictor variable. Second, the association between rheumatoid factor and risk of rheumatoid arthritis was studied by comparing participants by categories of rheumatoid factor.
Cumulative incidence curves were estimated by the method of Kaplan-Meier and Fine-Gray, and log rank trend tests examined differences between rheumatoid factor categories. We calculated hazard and subhazard ratios with 95% confidence intervals by means of Cox and Fine-Gray regression models with age as the time scale, left truncation, and delayed entry at age of entering the Copenhagen City Heart Study. The Cox and Fine-Gray models were adjusted for known risk factors and markers of lifestyle and social status—that is, for age (automatic adjustment as age is time scale), sex, alcohol intake, body mass index, current daily tobacco use, cumulative tobacco use, marital status, parity, and years of education. the covariates were at study entry, but were also used as time-varying covariates at follow-up examinations in 1991-94 and 2001-03. Missing values in covariates were imputed using multivariable normal regression imputation (mi impute mvn), where age at examination, sex and birth year were independent variables, and alcohol use, body mass index, current daily smoking, cumulative tobacco use, marital status, parity and years of education were dependent variables in the model. In sensitivity analyses we examined those individuals with complete information on all covariates for the 1981-1983 examination. Supplementary table 1 in the online data supplement compares individuals with complete information against the entire cohort with respect to variables used for multivariable adjusted analysis, and lists the fraction with missing information for each covariate.
For Cox proportional hazards regression analyses, we assessed the assumption of proportional hazards graphically by plotting log(cumulative hazard) for different rheumatoid factor categories as a function of log(age). We detected no major violations of the proportional hazards assumption. Interaction was tested for by introducing a two factor interaction term in the Cox models.
All individuals were followed from study entry and censored at the date of emigration (n=52), 10 August 2010, or diagnosis of rheumatoid arthritis (n=183), whichever came first. In the Cox regression models death (n=6380) also led to censoring, whereas in the Fine-Gray regression models death acted as a competing event. Given that 66% of participants died during follow-up, competing risks from death could potentially bias the findings. Cox regression appropriately stops allowing individuals from contributing person-time when they died and then no longer are at risk, but this means that failure from other causes is unobservable. Fine-Gray regression is a proportional hazards model for the cumulative incidence function, where failure from other causes is observable.18 If findings from both regressions are similar, then the influence of competing events on the association between rheumatoid factor and rheumatoid arthritis is limited, which is the case in the present study. Thus, competing risk of death (information on which is 100% complete in the Danish registries) was accounted for in the analysis. The five main causes of death for the cohort (as assessed from the national Danish Causes of Death Registry) are shown in supplementary table 2 in the online data supplement, distributed according to categories of rheumatoid factor. For evaluation of risk of rheumatoid arthritis developed within 10 years of blood sampling, all surviving participants were censored after 10 years for follow-up; this was done as plasma biomarkers may have larger risk estimates for shorter follow-up than longer follow-up time.
Poisson regression models were used for estimation of absolute 10 year risks of rheumatoid arthritis. These risks were estimated for each rheumatoid factor category as the percentage within 10 years from baseline who developed rheumatoid arthritis, stratified for age, sex, and smoking. Such 10 year absolute risks are used in prediction tools for other major diseases, such as cardiovascular disease,19 breast cancer,20 and prostate cancer.21
We included 9712 individuals from the general population aged 20-100 years with a measured rheumatoid factor and without a prior diagnosis of rheumatoid arthritis at study entry in a 28 year follow-up from 1981-83 to 10 August 2010. During these 187 654 person years of follow-up, 183 participants developed rheumatoid arthritis. Baseline characteristics are shown in table 1⇓. Median age at diagnosis of rheumatoid arthritis was 70 years (interquartile range 63-76). The median time from providing the blood sample to developing rheumatoid arthritis was 15 years for those with rheumatoid factor levels of <25 IU/mL, 12 years for those with 25-50 UI/mL, seven years for 50.1-100 IU/mL, and seven years for >100 IU/mL.
Plasma levels of rheumatoid factor by age
The median plasma level of rheumatoid factor was 18 IU/mL (supplementary fig 2 in online data supplement). Distributions of rheumatoid factor in individuals in the general population stratified in 10 year age groups showed similar levels of rheumatoid factor across all age categories (fig 1⇓). Plasma levels were measured up to 28 years after blood sampling and were not reported to participants or their doctors, and consequently could not have influenced the risk estimates described below.
Risk of autoimmune rheumatic diseases
Risk of rheumatoid arthritis during the 28 years’ follow-up increased by a factor of 3.3 (95% confidence interval 2.7% to 4.0%) for a doubling of rheumatoid factor level (fig 2⇓). The risk of the other autoimmune rheumatic diseases Sjögren’s syndrome, systemic lupus erythematosus, and systemic sclerosis showed similar trends for a doubling of rheumatoid factor level, but the number of events for individual diseases were insufficient to reach statistical significance. None of the participants developed polymyositis or dermatomyositis. Sensitivity analyses restricted to participants with complete information on all covariates at baseline showed similar results (supplementary fig 3 in online data supplement).
Risk of rheumatoid arthritis
The cumulative incidence of rheumatoid arthritis as a function of age for the four categories of baseline level of rheumatoid factor (<25, 25-50, 50.1-100, and >100 IU/mL) is shown in fig 3⇓ as Kaplan-Meier estimates (and in supplementary fig 4 in the data supplement as Fine-Gray estimates). The cumulative incidence of rheumatoid arthritis increased with increasing rheumatoid factor (log rank trend P<0.0001).
During the full 28 years of follow-up, the multivariable adjusted hazard ratios for rheumatoid arthritis were 3.6 (95% confidence interval 1.7 to 7.3) for rheumatoid factor levels of 25-50 IU/mL, 6.0 (3.4 to 10) for 50.1-100 IU/mL, and 26 (15 to 46) for >100 IU/mL compared with <25 IU/mL (Ptrend<0.0001) (fig 4⇓, top left panel). During the first 10 years of follow-up, the corresponding hazard ratios were 6.0 (2.1 to 17) for 25-50 IU/mL, 14 (6.7 to 28) for 50.1-100 IU/mL, and 39 (18 to 85) for >100 IU/mL (Ptrend<0.0001) (fig 4⇓, bottom left panel). Corresponding risk estimates that took account of the competing risk of death were similar (supplementary fig 5, left hand panels, in the data supplement). There was no interaction between sex and rheumatoid factor categories on risk of rheumatoid arthritis.
In an attempt to exclude patients misclassified with rheumatoid arthritis, we also examined patients with at least two hospitalisations for rheumatoid arthritis at least six weeks apart. This reduced the number of patients with rheumatoid arthritis but strengthened the association between elevated rheumatoid factor and risk of rheumatoid arthritis (fig 4⇑, right hand panels). These results were also not influenced by competing risk of death (supplementary fig 5, right hand panels). Sensitivity analyses restricted to participants with complete information on covariates at baseline also showed similar results (supplementary fig 6).
Absolute 10 year risk of rheumatoid arthritis
The highest absolute 10 year risk of rheumatoid arthritis of 32% was observed for 50-69 year old women who smoked and had rheumatoid factor levels >100 IU/mL (fig 5⇓, far right middle panel). This means that one out of three of such women will develop rheumatoid arthritis within 10 years from blood sampling. The lowest absolute 10 year risk of rheumatoid arthritis of 0.1% was observed for men ≥70 years old with rheumatoid factor levels <25 IU/mL irrespective of smoking status (fig 5⇓, bottom panels far left and inner left).
The principal findings in this study of 9712 individuals without rheumatoid arthritis recruited from the general population of Copenhagen are those with elevated levels of rheumatoid factor had up to 26-fold higher long term risk of developing rheumatoid arthritis and up to 32% 10 year absolute risk of developing rheumatoid arthritis. These findings are novel. Importantly, these data do not serve as evidence that rheumatoid factor plays a causal role in the pathogenesis of rheumatoid arthritis.
Mechanism and comparison with other studies
Development of rheumatoid arthritis is thought to be an inflammatory process from early arthritis through rheumatoid arthritis and possibly to severe extra-articular rheumatoid arthritis.10 13 22 23 Among patients with early arthritis, only 40% are seropositive for rheumatoid factor, whereas in the final stage of rheumatoid arthritis 80% of patients are seropositive.10 23 A current debate is whether elevated levels of rheumatoid factor, elevated levels of anti-citrullinated protein antibody, variations in the PTPN22 gene, or some combination of these offer the best means of predicting short term (2-3 years) risk of rheumatoid arthritis.8 9 24 25 26 27 28 29 30 31 However, the debate does not include long term risk prediction, simply because good evidence is lacking. Our finding that elevated rheumatoid factor levels are associated with an increased long term (up to 28 years) risk of rheumatoid arthritis provides such data. Unfortunately we do not have information available for elevated anti-citrullinated protein antibody or genetic variations in PTPN22. Our data show that elevated levels of rheumatoid factor can be present many years before the clinical manifestation of arthritis, supported by a few earlier studies.32 33 34
Rheumatoid factor levels are believed to increase with age in the general population,1 but we could not confirm this. Although the age of individuals with rheumatoid factor >100 IU/mL was 62 years compared with 58 years for the reference group with rheumatoid factor <25 IU/mL, the corresponding median times to diagnosis of rheumatoid arthritis were seven years and 15 years, indicating that the higher rheumatoid factor level >100 IU/mL group can be explained by a short time to diagnosis rather than by greater age. On direct investigation, we found that rheumatoid factor levels were constant across age groups from 20 years to 100 years old.
Tobacco smoking has long been known to play a role in the pathogenesis of rheumatoid arthritis.35 Smoking related alterations to the cytokine balance, stress to the immune system, and modifications of autoantibodies are strongly associated with rheumatoid arthritis.36 37 However, we adjusted extensively for smoking—including daily tobacco use at examination and cumulative tobacco use in pack years at baseline and as time varying covariates at follow-up examinations—so smoking is unlikely as a confounder for our observation of increased long term risk of rheumatoid arthritis in those with elevated rheumatoid factor levels.
Strengths and limitations of study
Our study has several strengths. The large study population was homogeneous and representative of the general population, the study population was well characterized, follow-up was more than 28 years, and we had 100% follow-up. Therefore, with a median age of entry to the study above 50 years, the 28 years of follow-up is well beyond the age of 70 years, when the incidence of rheumatoid arthritis peaks.1 Also, because rheumatoid factor concentrations were measured up to 28 years after blood sampling and thus not reported to participants or their doctors, these measurements did not influence ascertainment of rheumatoid arthritis during follow-up. In other words, we here study the natural course from an elevated rheumatoid factor to a clinical diagnosis of rheumatoid arthritis.
Potential limitations include selection bias; however, as we enrolled a sample from the general population randomly selected without regard to disease status, selection bias is unlikely as an explanation for our findings. Another potential limitation is the validity of the diagnostic information; however, most of the rheumatoid arthritis diagnoses in the Danish registries were confirmed,38 and random misclassification would tend to bias the results towards the null hypothesis and therefore cannot explain the present results (as illustrated by our sensitivity analysis). Also, as rheumatoid factor enters into the classification of rheumatoid arthritis, individuals with lifelong elevated rheumatoid factor may have a heightened probability of being diagnosed with rheumatoid arthritis. Hospitalisation with a discharge diagnosis of rheumatoid arthritis can happen many years after the first diagnosis of rheumatoid arthritis while a patient is under the care of a rheumatologist or being cared for in some other outpatient setting; thus, those cases identified from a hospital discharge diagnosis may have been classified as not having rheumatoid arthritis when they were simply not hospitalised with it. Similarly, those with rheumatoid arthritis who eventually were hospitalised might have contributed substantially more person years than they should. Likewise important, the patients in this study were presumably at a later or more severe disease stage than patients consulting private rheumatologists; however, this implies that time to diagnosis in the outpatient setting would be shorter and risk estimates correspondingly higher than observed in this study. Given that the categorisation of rheumatoid factor was devised to maximise the relation between rheumatoid factor and rheumatoid arthritis, the risk estimates based on these categories could be exaggerated and need to be replicated in another study. Finally, as we studied only white people, our results may not apply to other races.
Implications for clinicians and future work
Our finding of a particularly high absolute risk of developing rheumatoid arthritis in women with elevated rheumatoid factor who smoked is yet another reason to promote cessation of smoking among such women. In addition, our finding of high risks of developing rheumatoid arthritis based on elevated levels of rheumatoid factor alone suggests the need for early referral to a rheumatologist or to early arthritis clinics for examination on the basis of a positive rheumatoid factor test—even in the absence of the typical arthritic joint symptoms—because of the better response to therapy the earlier it is initiated in rheumatoid arthritis.
The superior specificity of anti-citrullinated protein antibody tests has somewhat diminished interest in rheumatoid factor at time of diagnosis of rheumatoid arthritis, but not necessarily when it comes to prediction of long term risk. Future work should include a comparison of rheumatoid factor and anti-citrullinated protein antibody tests for the long term prediction of development of rheumatoid arthritis. Also, larger studies might be able to confirm or refute that elevated rheumatoid factor is associated with future risk for other autoimmune rheumatic diseases.
Finally, our absolute risk estimates could be used for designing a randomised controlled trial studying early intervention in high risk people, by designing the inclusion criteria around patient groups with, say, a ≥20% elevated 10 year absolute risk of developing rheumatoid arthritis.
Individuals in the general population without rheumatoid arthritis but with an elevated plasma level of rheumatoid factor have up to 26-fold greater long term risk of developing rheumatoid arthritis, and up to 32% 10 year absolute risk of rheumatoid arthritis.
What is already known on this topic
Rheumatoid arthritis is an autoimmune disease affecting 0.5-2% of the population
At present, there is no good clinical available indicator for long term risk of developing rheumatoid arthritis
What this study adds
Individuals in the general population without rheumatoid arthritis, but with an elevated plasma level of rheumatoid factor have up to 26-fold greater long term risk of developing rheumatoid arthritis, and up to 32% 10 year absolute risk of rheumatoid arthritis
These findings may lead to revision of guidelines for early referral to a rheumatologist and early arthritis clinics based on rheumatoid factor testing
Cite this as: BMJ 2012;345:e5244
We thank staff and participants of the Copenhagen City Heart Study for their contribution.
Contributors: SFN and BGN had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. SFN, SEB, and BGN designed the study. PS, SEB, and BGN collected data. SFN performed all statistical analyses. SEB and BGN oversaw statistical analyses and contributed to interpretation of data. SFN wrote the first draft of the paper. SEB, PS, and BGN edited the paper critically, and all authors approved this paper in its final form. SFN and BGN are the guarantors.
Funding: This study was funded by Herlev Hospital, Copenhagen University Hospital and the Danish Heart Foundation. The sponsors had no role in the design of the study; collection, analysis, and interpretation of data; writing of the report; or decision to submit the article for publication. The views expressed in this paper are those of the authors and not those of any funding body or others whose support is acknowledged.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: The studies were approved by Herlev Hospital, Copenhagen University Hospital, and a Danish ethical committee (the Copenhagen and Frederiksberg committee Nos. KF-100.2039/9 and KF-01-144/01), and were conducted according to the Declaration of Helsinki. Participants provided written informed consent.
Data sharing: No additional data available.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.