Visual Overview
Abstract
Background Epidemic levels of CKD of undetermined cause, termed Mesoamerican nephropathy in Central America, have been found in low- and middle-income countries. We investigated the natural history of, and factors associated with, loss of kidney function in a population at high risk for this disease.
Methods We conducted a 2-year prospective, longitudinal study with follow-up every 6 months in nine rural communities in northwestern Nicaragua and included all men (n=263) and a random sample of women (n=87) ages 18–30 years old without self-reported CKD, diabetes, or hypertension. We used growth mixture modeling to identify subgroups of eGFR trajectory and weighted multinomial logistic regression to examine associations with proposed risk factors.
Results Among men, we identified three subpopulations of eGFR trajectory (mean baseline eGFR; mean eGFR change over follow-up): 81% remained stable (116 ml/min per 1.73 m2; −0.6 ml/min per 1.73 m2 per year), 9.5% experienced rapid decline despite normal baseline function (112 ml/min per 1.73 m2; −18.2 ml/min per 1.73 m2 per year), and 9.5% had baseline dysfunction (58 ml/min per 1.73 m2; −3.8 ml/min per 1.73 m2 per year). Among women: 96.6% remained stable (121 ml/min per 1.73 m2; −0.6 ml/min per 1.73 m2 per year), and 3.4% experienced rapid decline (132 ml/min per 1.73 m2; −14.6 ml/min per 1.73 m2 per year; n=3 women). Among men, outdoor and agricultural work and lack of shade availability during work breaks, reported at baseline, were associated with rapid decline.
Conclusions Although Mesoamerican nephropathy is associated with agricultural work, other factors may also contribute to this disease.
CKD of undetermined cause (CKDu), also termed Mesoamerican nephropathy (MeN) in Central America, has led to the deaths of tens of thousands of young adults in rural Nicaragua and El Salvador.1,2 Cross-sectional studies have shown low (<60 ml/min per 1.73 m2) eGFR at a prevalence between 2% and 50% among the population of lowland agricultural communities in the region.3–6 Forms of CKDu occur in other tropical climates, with reports of high prevalence in Sri Lanka7,8 (where similar but not identical histopathologic findings have been reported9), India,10 and Egypt,11 although whether this represents the same disease entity remains unclear.
Men from communities affected by MeN predominantly work in agriculture, primarily sugar production from cane. Agricultural activity in this industry is concentrated in the dry season, which in Nicaragua, occurs between November and May. Although a leading hypothesis in Mesoamerica is that the disease relates to heat stress, a number of other causes, including agrichemicals, infection, and heavy metals, have been proposed.1,12–14
Empirical evidence for causes of CKDu has to date been limited to identification of factors associated with either reduced eGFR in cross-sectional studies3,15,16 or loss of eGFR across the harvest season in two workplace-based follow-up studies.17,18 Given the potential for reverse causation (i.e., reduced eGFR resulting in changes in exposure) and misclassification of exposures and/or outcome in the cross-sectional designs along with the nongeneralizability and the substantial loss to follow-up that occurred in the longitudinal workplace studies, evidence on risk factors for and evolution of CKDu is extremely limited.19
Our aim was to investigate the natural history of disease, specifically early loss of kidney function, along with risk factors and urinary markers (albumin-to-creatinine ratio [ACR] and neutrophil gelatinase–associated lipocalin [NGAL]) associated with decline in eGFR. Therefore, we conducted a community-based longitudinal study of an initially apparently healthy young rural population in northwest Nicaragua.
Methods
Cohort
Both local and United Kingdom–based institutional review boards approved the study, and participants provided written informed consent. The rationale and description of the study design have been published elsewhere.20 Briefly, this was a 2-year longitudinal, community-based study following 350 participants ages 18–30 years old in the Leon and Chinandega regions of Nicaragua (Figure 1). After engagement work, we performed a census of all adults ages 18–30 years old in nine rural communities. Because we were specifically interested in associations with early kidney injury in MeN, all potential participants with a self-reported diagnosis of CKD, diabetes, or hypertension were excluded. All remaining men (because men have been reported to suffer more CKDu) and women selected at random (in numbers leading to a men-to-women ratio of 3:1) were invited to take part. Participants were predominantly recruited in November 2014, with an additional 7% recruited in May 2015, because recruitment targets had not been met in November.
Study participants were recruited from nine communities in Northwest Nicaragua and study retention rates were high. (A) Location of the nine study communities in Nicaragua. (B) Cartoon showing the study timeline along with population, recruitment, and follow-up numbers. Two participants died from ESRD.
Procedures
Questionnaire data, clinical measurements, and biosamples were collected at baseline and then every 6 months until November 2016. Participants were asked to respond to questions on demography, occupational history and current job, lifestyle factors, and symptoms. Urinary tract infection was recorded where participants reported a clinical diagnosis (which is common in this part of Nicaragua), typically without urinalysis or microbiologic confirmation. Body weight was measured with minimal clothes using electronic scales (Seca, Birmingham, United Kingdom), and height was measured using a portable stadiometer (Seca). BP and heart rate were measured in a sitting position using a calibrated digital sphygmomanometer (Omron, Kyoto, Japan) after 5 minutes of quiet seated rest. A mean of three measurements was recorded. Participants were asked to attend fasted first thing in the morning (before work) in an attempt to reduce within- and between-person variation in serum creatinine.
Biochemical Methods
Serum creatinine and cystatin C were both measured in a single batch using quality control referenced to international standards (for creatinine, isotope dilution mass spectrometry–quantified National Institute of Standards and Technology Standard Reference Material 967). eGFR was calculated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) formula combining creatinine and cystatin C.21 ACR along with semiquantitative protein and specific gravity by test stick were performed in baseline urine samples thawed for the first time. In addition, 55 samples (thawed for a second time) selected using a nested patient-control approach were analyzed for NGAL.
Statistical Methods
The collection and categorization of exposure variables are described in Supplemental Material. Because eGFR trajectories clustered in discrete subgroups (Supplemental Figures 1 and 2) and differently between sexes, we used growth mixture modeling (GMM) separately in men and women to empirically derive latent classes of eGFR trajectory.22 The GMM is a longitudinal finite mixture model that allows identification of unobserved latent classes of individuals following similar progression of the outcome over time without imposing a priori constraints on the levels of eGFR or rates of eGFR change (or the proportion of participants experiencing any class of change). Each individual’s probability of belonging to a particular latent class is derived entirely from the observed eGFR measurements, with individual departures from the mean trajectory within each class represented by random effects. We primarily used the Bayesian Information Criterion to determine the optimal number of classes as suggested in this setting.23 The GMM was estimated by maximum likelihood using an expectation maximization algorithm, with 95% confidence intervals (95% CIs) for the mean rate of eGFR decline derived using conventional SEM.
Each individual was assigned a probability of each class (eGFR trajectory) and then for the purposes of the descriptive figure, tables, and urinary findings, allocated to the highest probability group.
To test whether proposed causal exposures (alcohol or nonsteroidal anti-inflammatory drug use, occupational factors, heat stress, agrochemical exposure, fever, dysuria, or water quantity/quality/source in men only) were associated with rapidly declining eGFR trajectory, we conducted age- and educational level–adjusted analyses using probability-weighted logistic regression (with weighting according to the participant’s probability of each eGFR trajectory as per the GMM), examining exposures individually using stable with preserved eGFR trajectory as a reference. Associations where the 95% CI of the odds ratio (OR) did not include unity were interpreted as significant. We also performed a sensitivity analysis using exposures assessed at visit 2 (only in those men recruited at visit 1) and rapid decline given the seasonal variation in occupational exposures. Those with baseline dysfunction were not the primary focus of this study, but another analysis additionally exploring associations between risk factors and this eGFR trajectory was also performed using probability-weighted logistic regression (Supplemental Material).
Differences in urinary markers in each eGFR trajectory group (defined on the basis of the highest probability as above) were investigated either in the whole population for ACR or using a nested patient-control approach in the case of NGAL. Differences between groups were explored using ANOVA with the Dunnett post hoc test. Positive and negative predictive values were calculated for urinary NGAL for the rapid decline versus stable group.
Results
Cohort and Follow-Up
Five hundred twenty adults ages 18–30 years old were identified in the study communities. After exclusion of 4% of the potential participants because of self-reported CKD, diabetes, or hypertension, 350 participants (of the 360 invited after random selection of eligible women; 97%) were included in the study.20 Overall, participants attended a total of 1581 study visits over the 2-year follow-up (92% of planned visits). Two participants died from ESRD during the study period. The cohort is described in Figure 1 and Table 1.
Selected demographic, lifestyle, and occupational characteristics of the study cohort
The median eGFR in men was 116.2 ml/min per 1.73 m2 (interquartile range [IQR], 102.4–124.6) at baseline, and 110.4 ml/min per 1.73 m2 (IQR, 92.5–120.5) at the end of follow-up. The corresponding figures for women were 122.0 ml/min per 1.73 m2 (IQR, 116.3–127.2) at baseline and 120.2 ml/min per 1.73 m2 (IQR, 110.6–126.6) at the end of follow-up. The eGFR varied by season (Figure 2), with a median of 116.0 ml/min per 1.73 m2 (IQR, 102.7–123.8) at the end of the rainy season (November; i.e., before sugarcane harvest, all years combined) compared with 113.4 ml/min per 1.73 m2 (IQR, 100.8–122.4) at the end of the dry season (May; i.e., after sugarcane harvest, all years combined). This effect was greatest in those participants with lower eGFRs, but it was also present in those with stable kidney function (Supplemental Table 2).
Median eGFR was lower postharvest than preharvest. Box and whisker plot of eGFR across the study population. The dashed line represents the median of all eGFR values in the population across the study.
eGFR Trajectory Groups
Using GMM, we identified three different subgroups in men and two subgroups in women on the basis of the model intercept (baseline eGFR) and slope (change in eGFR over time). Among men (Figure 3A), the majority (81%) of men had preserved and stable eGFR; however, 9.5% (n=25) had baseline kidney dysfunction (eGFR of approximately 60 ml/min at recruitment), and another 9.5% experienced rapid decline in eGFR (with a mean loss of 18 ml/min per 1.72 m2 per year) despite preserved eGFR at baseline. Almost all of the women (Figure 3B) had preserved and stable eGFR, but 3.4% (n=3) also experienced rapid decline (with a mean loss of 14 ml/min per 1.72 m2 per year). No differences were seen between communities in the proportions of participants in these subgroups.
A substantial proportion of men and a small number of women experience rapid decline in eGFR. Individual eGFR values over time stratified by trajectory subgroup. (A) Three subgroups were identified in 263 men, and (B) two subgroups were identified in 87 women. Each line represents the individual eGFR of a single participant. Each participant was allocated to the group of highest probability derived from the growth mixture model. Coefficients for the three groups of men: preserved and stable eGFR (n=213; intercept [mean eGFR at baseline], 113.3 ml/min per 1.73 m2; 95% confidence interval [95% CI], 111.3 to 115.3; slope [mean eGFR decline over time], −0.6 ml/min per 1.73 m2 per year; 95% CI, 0.0 to −0.9); rapid decline in eGFR (n=25; intercept, 109.5 ml/min per 1.73 m2; 95% CI, 99.1 to 119.9; slope, −18.2 ml/min per 1.73 m2 per year; 95% CI, −13.5 to −22.9); and baseline dysfunction (n=25; intercept, 55.6 ml/min per 1.73 m2; 95% CI, 48.5 to 62.7; slope, −3.8 ml/min per 1.73 m2 per year; 95% CI, −0.7 to −6.9). Coefficients for the two groups women: preserved and stable eGFR (n=84; intercept, 120.5 ml/min per 1.73 m2; 95% CI, 118.1 to 122.9; slope, −0.6 ml/min per 1.73 m2 per year; 95% CI, 0.2 to −1.4). We also identified a small number with rapid decline in kidney function (n=3; intercept, 127.5 ml/min per 1.73 m2; 95% CI, 119.3 to 135.7; slope, −14.6 ml/min per /1.73 m2 per year; 95% CI,−7.5 to −21.7).
Baseline sociodemographics, occupational history, occupational exposures, lifestyle factors, and symptoms stratified by the assigned kidney trajectory groups are presented in Supplemental Tables 2 and 3. The frequencies of indoor work and availability of shade were both lower in the rapidly declining subgroup. Of the three women who fell into the rapid decline group, one had worked in (nonsugarcane) agriculture, and two worked exclusively at home.
Adjusted Associations with Rapid Decline Trajectory
Baseline age– and educational level–adjusted probability-weighted associations with the rapid decline in eGFR trajectory in men using the preserved and stable trajectory as the reference are presented in Table 2. Outdoor work (OR, 10.35; 95% CI, 1.35 to 79.24), (nonsugarcane) agricultural work (OR, 3.57; 95% CI, 1.14 to 11.13), and lack of shade available during work breaks (OR, 3.74; 95% CI, 1.59 to 8.76) were associated with this outcome. However, we found no evidence for associations between rapid decline and years of work in sugarcane or agriculture; self-reported physical effort in the last week at work; self-reported occupational heat or agrochemical exposure over last 6 months; alcohol consumption, self-reported fluid consumption, or water quality or source; heat/dehydration-related symptoms; or use of nonsteroidal anti-inflammatory drugs.
Age- and education level–adjusted associations of rapid decline in eGFR by baseline exposure in study participants who were men
We were concerned that the questionnaire administered at baseline might fail to capture important occupational exposures, because for most participants, it was conducted 6 months after the harvest season. Therefore, we conducted a sensitivity analysis (men recruited at the November visit only; n=213) examining the association with the same rapid decline eGFR trajectory as above and occupational exposures, hydration variables, and heat-related symptoms captured at the second study visit (May 2015; immediately after harvest) (Supplemental Table 4). At this time point, no associations were detected between working outdoors or lack of shade and rapid decline in eGFR trajectory (although very few participants were not exposed). There was an association between both those working in a sugarcane cutting role (OR, 3.84; 95% CI, 1.17 to 12.58) and those reporting fever over the last 6 months (OR, 5.77; 95% CI, 2.03 to 16.33) and rapid decline trajectory, but in line with the baseline exposure analysis, no associations were observed between self-reported measures of heat exposure, combined heat-related symptoms, or fluid intake and outcome (Supplemental Tables 5 and 6).
Urinary Findings
No associations were found between dipstick proteinuria, specific gravity, or ACR and eGFR trajectory subgroups (Tables 3 and 4). Urinary NGAL levels among men differed between the three groups tested (Figure 4). The positive and negative predictive values of NGAL≥5.5 pg/mmol for rapid decline were 28.5% and 62.5%, respectively.
Description of urinary findings at baseline by assigned eGFR trajectory groups in men
Description of urinary findings at baseline by assigned eGFR trajectory groups in women
Urinary NGAL concentrations were higher in men with baseline kidney dysfunction and those who experienced rapid decline in eGFR. Box and whisker plot of urinary neutrophil gelatinase–associated lipocalin (NGAL)/creatinine concentrations by assigned eGFR trajectory group among male study participants. Lines indicate medians. Boxes are interquartile ranges. Whiskers indicate 1.5× interquartile ranges. Dots are outlying values. Stable group, n=55; rapid decline group, n=25; baseline dysfunction, n=24. *P=0.03 using ANOVA followed by the Dunnett multiple comparisons test (using the stable and preserved eGFR group as the reference); ****P≤0.001 using ANOVA followed by the Dunnett multiple comparisons test (using the stable and preserved eGFR group as the reference).
Discussion
This is the first community-based cohort study from an area with high reported prevalence of MeN and the first longitudinal study of at least moderate size with follow-up of >6 months in an area at high risk of disease. Even after excluding those with self-reported CKD, 9.5% of the apparently healthy men (but no women) in the study had evidence of baseline renal dysfunction. Rapid loss of eGFR from normal baseline levels was found in another 9.5% of men and 3.4% of women. Among men, risk factors at baseline for rapid decline included working outdoors, agricultural work, and lack of shade availability, but none of the other questionnaire responses aimed at capturing heat stress, time-accumulated occupation, or other proposed causes of MeN were associated with the outcome at baseline. Because of small numbers, we were unable to examine associations in women.
Other important findings from our study include the cyclical annual fluctuation in renal function across the entire population, with the average eGFR approximately 2.5-ml/min per 1.73 m2 lower after the dry (harvest) season compared with 6 months earlier. Furthermore, although there were no differences in albuminuria between those with different kidney function trajectories, urinary NGAL was substantially higher among those with baseline dysfunction and marginally elevated in the rapid decline group.
Although CKDu has been anecdotally reported as an aggressive disease,1 the rate of loss of kidney function in those in the rapid decline group who make up almost 10% of the unselected population of young men in our study is, to our knowledge, without precedent. Even compared with eGFR decline in other forms of CKD seen in clinic populations, the observed loss of kidney function is alarming. Although a recent biopsy study that enrolled patients with established CKDu reported a rate of decline in eGFR of 7.0 ml/min per 1.73 m2 per year among men with a history of work in the sugarcane,24 there have been no longitudinal studies that have examined medium- or long-term (>1-year) changes in kidney function in the at-risk population. The rate of eGFR decline has been explored in more detail in other forms of CKD; for example, a longitudinal study in 55 clinic patients with diabetic nephropathy from Belgium reported that approximately 15% of patients suffered severe decline in kidney function (defined as eGFR loss >4 ml/min per 1.73 m2 per year).25 Most recently, Boucquemont et al.26 examined eGFR decline in a patient population with CKD in France using a similar latent class-based modeling approach to that used in this analysis. This study reported severe eGFR decline in only 0.6% of patients (approximately 50 ml/min per 1.73 m2 over almost 6 years). Therefore, our study findings underline the unique and severe nature of kidney disease in this region.
The associations with rapid decline trajectory in men suggest that occupation (outdoor agricultural work) is an important risk factor for loss of kidney function, consistent with previous reports.18 The temporary nature of work in this population makes distinguishing relationships between specific occupations and eGFR loss challenging; however, it is interesting to note that neither time-accumulated sugarcane work nor agricultural work were associated with outcome. Furthermore, the association between lack of available shade at baseline and rapid decline trajectory suggests that working environment may play an important role in disease evolution either by (not reducing) solar exposure or as a surrogate for generally poor occupational conditions. Consistent with this and in line with previous crossharvest studies,17 we identified an association between rapid decline and a cane/seed cutting role (a job role that has been associated with particularly hot working conditions) in a sensitivity analysis examining associations with exposures assessed postharvest.
The abscence of an association between variables aimed at capturing heat stress (self-reported physical effort the previous week at work and both work in very hot environment and combined dehydration/heat stress symptoms in the last 6 months) and the outcome measure, both at baseline and in the sensitivity analysis with exposures assessed at visit 2 raises further questions. Although self-reported measures of thermal sensation and physical exertion have been shown to robustly capture acute physiologic heat stress,27 our (similar) instruments (and/or our combined measure of heat symptoms) may not be valid in the rural Nicaraguan population, or they may not reflect time-accumulated heat stress. Alternatively, we may have had inadequate power to detect heat stress as a partial contributor to eGFR decline, or otherwise, it may be that nonheat-related occupational exposures promote the development of CKDu. Finally, the association between self-reported fever over the previous 6 months at the second study visit and the rapid decline trajectory might support a proposed infective/inflammatory contributor to MeN,28 although this finding was from a sensitivity analysis and should be treated with caution.
In summary, our data do not provide clear evidence for a cause of the disease. Along with occupation, the importance of nonoccupational factors is supported by (1) the range of jobs undertaken by the men experiencing rapid decline and (2) the 3.4% of women in our study who also showed a rapid loss of eGFR. As others have suggested,2 separate initiating and exacerbating factors should be considered, such as in other forms of CKD. For example, the progression of kidney disease due to known causes (e.g., diabetes or GN) can be exacerbated by episodes of volume depletion. Therefore, the possibility of an initial (currently unknown) subclinical insult, which is then exacerbated by the harsh working conditions, might explain the increased rates of eGFR loss and excess of advanced disease in men.
Although other studies have identified changes in urinary biomarkers in sugarcane workers over the harvest season in Mesoamerica,29 none have examined associations with subsequent eGFR loss over the medium term. There were no associations between dipstick proteinuria or ACR and eGFR trajectory group. Although albuminuria is a strong risk factor for renal decline in most populations, this is consistent with previous reports from Mesoamerica, where patients with established CKDu show only low-grade proteinuria.6,24,30 Urinary NGAL levels were substantially raised in those with baseline dysfunction, but levels in the rapid decline group overlapped with the stable group, making this test poorly predictive at an individual level.
Finally, it is worth noting the seasonal variation of eGFR in the population. Other studies (unrelated to CKDu) have described similar seasonal differences in renal function31,32; whether this variation is in any way related to the factors that cause MeN is unclear, but this finding does need to be considered when interpreting the change in eGFR reported in crossharvest studies.5,18 Ideally, any future longitudinal biomarker study should be of >1 year in duration to ensure that small falls in eGFR do not reflect cyclical seasonal changes.
Our study has several strengths. Overall response rates were high, and the eGFR was estimated using robust methods. We excluded those with self-report of diabetes and hypertension in an attempt to focus our study on eGFR decline due to MeN, and the prospective nature of our study enabled us to identify those with aggressive disease without necessarily meeting definitions for CKD. Furthermore, we excluded those with established renal disease (either by self-report from the study as a whole or by examining only those with preserved eGFR at baseline for the risk factor analysis), and hence, we could overcome issues associated with reverse causation.
Our study also has limitations. We did not formally exclude diabetes in our participants. Although often undiagnosed,33 the prevalence of diabetes is low in Nicaraguans of this age group,34 and none of those in the rapid decline group showed albuminuria (or glycosuria; data not presented), making an underlying diabetic lesion highly unlikely. We also relied on self-report to quantify the majority of occupational and environmental exposures. Although questionnaire-based assessments are useful instruments, none of them have been validated in the Nicaraguan population; therefore, some exposures may be prone to misclassification. The study took place in a confined geographical area, which limits generalizability. Resources restricted our study to a moderate sample size, and we had to alter our statistical approach. We were nonetheless able to detect a number of strong associations with eGFR trajectory, but the analytical change did lead to a reduction in power. Therefore, we would have expected to identify associations with a primary cause of disease that had been reliably captured by questionnaire but may have missed weaker associations, particularly with contributing exposures. The baseline dysfunction group is unrepresentative due to selection criteria (those with established CKD were intentionally excluded at recruitment) and possibly survivor bias (due to the small number of deaths in this group), and the nature of the study design means that the relationship between rapid decline in eGFR and hard outcomes could not be described. However, we hope to perform extended follow-up to investigate the longer-term outcomes in the cohort. Finally, the CKD-EPI formula has not been validated for this population, although because we were interested in within-person change in eGFR, this is unlikely to be of major importance.
In conclusion, this is the first community-based cohort study that describes the natural history of eGFR in those at risk of MeN. Almost 10% of apparently healthy young men and 3.4% of young women showed a marked decline in kidney function. Additional studies with at least 1-year of follow-up are needed to understand the causes of this decline, including the risks associated with outdoor (agricultural) work. Efforts to identify biomarkers of this early loss of eGFR rather than established disease are essential to gain a better understanding of etiology as well as to identify the population(s) that would benefit from interventions. A combined multidisciplinary approach is called for in partnership with the affected communities and local employers to address this devastating disease.
Disclosures
None.
Acknowledgments
The authors would like to thank the participants and each of the community leaders for their support during the data collections across study visits. We would also like to thank the interview team, drivers, phlebotomists, and staff of the Research Centre for Health, Work and Environment, National Autonomous University of Nicaragua at Leon for their assistance during the study.
The study was supported by a grant from the Colt Foundation, UK. In addition, the Dutch National Postcode Lottery provided funding to Solidaridad to support a proportion of the fieldwork costs.
No funding source was involved in any part of the study design or the decision to submit the manuscript for publication.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
This article contains supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2018020151/-/DCSupplemental.
- Copyright © 2018 by the American Society of Nephrology