Environmental predictors of pulmonary nontuberculous mycobacteria (NTM) sputum positivity among persons with cystic fibrosis in the state of Florida

Abstract

Nontuberculous mycobacteria (NTM) are opportunistic human pathogens that are commonly found in soil and water, and exposure to these organisms may cause pulmonary nontuberculous mycobacterial disease. Persons with cystic fibrosis (CF) are at high risk for developing pulmonary NTM infections, and studies have shown that prolonged exposure to certain environments can increase the risk of pulmonary NTM. It is therefore important to determine the risk associated with different geographic areas. Using annualized registry data obtained from the Cystic Fibrosis Foundation Patient Registry for 2010 through 2017, we conducted a geospatial analysis of NTM infections among persons with CF in Florida. A Bernoulli model in SaTScan was used to identify clustering of ZIP codes with higher than expected numbers of NTM culture positive individuals. Generalized linear mixed models with a binomial distribution were used to test the association of environmental variables and NTM culture positivity. We identified a significant cluster of M. abscessus and predictors of NTM sputum positivity, including annual precipitation and soil mineral levels.

About the Authors:

Sydney L. Foote

Roles Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Office of Data Science and Emerging Technologies, Office of Science Management and Operations, National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH), Rockville, MD, United States of America

ORCID logo https://orcid.org/0000-0002-2351-6910

Ettie M. Lipner

Roles Methodology, Writing – original draft, Writing – review & editing

Affiliations Center for Genes, Environment, and Health, National Jewish Health, Denver, CO, United States of America, Department of Epidemiology, Colorado School of Public Health, Aurora, CO, United States of America

D. Rebecca Prevots

Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

Affiliation: Division of Intramural Research, Epidemiology and Population Studies Unit, NIAID, NIH, Rockville, MD, United States of America

Emily E. Ricotta

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Division of Intramural Research, Epidemiology and Population Studies Unit, NIAID, NIH, Rockville, MD, United States of America

ORCID logo https://orcid.org/0000-0001-9928-8275

Introduction

Nontuberculous mycobacteria (NTM) are opportunistic human pathogens that reside in the environment and are commonly found in soil and water [1]. Exposure to these organisms may cause pulmonary nontuberculous mycobacterial (PNTM) disease, which poses a threat to high-risk groups including older adults and individuals with chronic lung conditions; persons with CF are particularly vulnerable [2]. Infection likely results from a combination of behavioral and environmental exposure which jointly increase the risk of NTM infection [3]. Studies have tried to evaluate common exposure sources including household plumbing such as showerheads [4,5], water heating units, and dust, as well as external sources such as soil, watersheds, and climatic factors [6–8]. These factors play a role in determining environmental suitability for the pathogen, with higher abundance of NTM associated with increased rainfall, humidity, certain watersheds, and soil composed of particular elements. Studies have also shown that prolonged exposure to certain environments can increase the risk of PNTM [9]. It is therefore important to determine the risk associated with different geographic areas, so high risk individuals can either avoid an area or focus prevention efforts to reduce their risk of exposure. Previous studies have shown significant clustering of NTM in multiple regions of the United States, including Florida, with geographical heterogeneity in overall and species-specific risk of pulmonary disease. In Florida, the 5-year NTM sputum positivity prevalence was 31% from 2010 through 2014 [2], and this state has the highest prevalence of NTM in the contiguous US [10,11]. Identifying risk factors for exposure and subsequent infection with NTM are central to prevention efforts in the CF community [12]. In this study we describe spatial clusters of NTM infections in Florida and identify environmental predictors of NTM sputum positivity.

Materials and methods

Data sources

We conducted a nested case-control study, using annualized registry data obtained from the Cystic Fibrosis Foundation Patient Registry (CFFPR) for 2010 through 2017 [13]. Our study population comprised patients aged ≥ 12 years residing ≥ 2 consecutive years in Florida. Patients with a history of lung transplant or Mycobacterium tuberculosis infection were excluded. Incident NTM cases were persons with a positive pulmonary culture after ≥ 1 negative culture(s) and with residence in Florida the year before and the year of their first positive NTM culture. Controls were defined as persons with ≥1 negative NTM culture(s) during the study period, residence in Florida the year of and the year before a negative culture, and no positive cultures during Florida residency. Because the majority of individuals had multiple cultures, we used the first culture and residential ZIP code associated with the year of first culture for each person meeting these criteria for analysis.

We selected environmental variables for analysis based on prior findings. Variables that have been previously found to be predictive of sputum positivity include evapotranspiration [10], saturated vapor pressure [14], vapor pressure [15], temperature [6], and rainfall [6], as well as soil or water mineral concentration including copper [10], sodium [10], manganese, [8,10], calcium [7], and molybdenum [7]. Environmental data sources used in this study are described in Table 1. Soil geochemistry collected from 2007 through 2010 included data on calcium, copper, molybdenum, manganese, and sodium content from samples measured in the top 5 cm of soil. Annual temperature and rainfall, and eight-day evapotranspiration data were extracted for the years 2010 through 2017. Evapotranspiration was averaged to create annual estimates. Ordinary kriging was performed to estimate the broader spatial distribution over Florida from the original sampling sites or weather stations for all soil geochemistry variables, temperature, and rainfall. To adapt county-level census data to the ZIP-code level, we determined the county each ZIP code primarily fell within by overlaying ZIP code polygons on a map of Florida counties using the R packages “rgdal” and “sp” [16–18]. The ZIP code-level mean of each environmental variable was calculated, scaled, and mean-centered by year, when data were available, for analysis. Driving distance from ZIP code centroids to CF clinics in Florida were calculated using a Distance Matrix API via Google Cloud Services and the R package “gmapsdistance” [19] to control for potential spatial clustering near CF clinics. The closest clinic for individuals was assumed to be the pediatric or adult clinic with the shortest driving distance for all patients in each ZIP code under or at least 18 years of age, respectively.

[Figure omitted. See PDF.]

Table 1. Data sources for variables used in analysis.

https://doi.org/10.1371/journal.pone.0259964.t001

Spatial and statistical analysis

We used a Bernoulli model comparing individuals with CF who were NTM culture positive (cases) to individuals with CF who were NTM culture negative (controls) to identify clustering of ZIP codes with higher than expected numbers of NTM culture positive individuals from 2011 through 2017, using SaTScan version 9.6 [20] with default setting, limiting cluster radius to 50km. This was done as a way to assess the presence of geographic locations in Florida where people with CF might be at increased risk of NTM. We repeated the analysis for overall NTM as well as each species of NTM separately. To test the association of environmental variables and NTM culture positivity, we used generalized linear mixed models with a binomial distribution. The variance inflation factor was used to assess collinearity and model fit was assessed via Akaike Information Criterion. The final environmental variables included were copper, manganese, molybdenum, sodium, annual precipitation, and evapotranspiration from the year in which an individual had their first NTM culture. Additionally, we included patient gender, age, number of years receiving chronic macrolides through the year of culture, the closest CF clinic in miles, whether the ZIP code was considered metropolitan, ZIP code median household income, and a random intercept for ZIP code. Statistical analyses were conducted using R versions 3.6.1–4.0.2 [21]. The study was determined to be not Human Subjects Research by the NIH Office of Human Subjects Research Protection.

Results

Of the 1293 patients in the CFFPR residing ≥ 2 consecutive years in Florida from 2010 through 2017, 979 patients met inclusion criteria; 261 (26.7%) were classified as cases and 718 (73.3%) as controls (Fig 1). Species identified were distributed as follows: 109 (41.8%) M. avium complex (MAC), 118 (45.2%) M. abscessus and its subspecies, and 62 (23.8%) other species. The proportions of sex, age, and years of chronic macrolide use were similar between cases and controls, while years of Florida residency from 2010 to the year of the first culture were generally greater for cases (Table 2). We found one statistically significant high-risk cluster in southeast Florida (p-value: 0.035, radius: 48.8 km), which included Broward, Miami-Dade, and Palm Beach counties (Fig 2). This high risk-cluster was associated with M. abscessus; no significant clustering was observed for other species.

[Figure omitted. See PDF.]

Fig 1. Study population for incident NTM cases and negative controls in the CFFPR from 2010 through 2017.

https://doi.org/10.1371/journal.pone.0259964.g001

[Figure omitted. See PDF.]

Fig 2. Distributions of persons residing in Florida for ≥ 2 years from 2011 through 2017.

A. Persons with cystic fibrosis (CF), n = 979. B. Nontuberculous mycobacteria (NTM) cases among persons with CF, n = 261. C. ZIP codes within a high-risk cluster of NTM cases among persons with CF.

https://doi.org/10.1371/journal.pone.0259964.g002

[Figure omitted. See PDF.]

Table 2. Baseline demographics of the study population by NTM culture result, Florida, 2011–2017.

https://doi.org/10.1371/journal.pone.0259964.t002

Sputum-positive patients were more likely to live in a ZIP code with higher average yearly precipitation, with a 34% increase in the odds of a positive culture for each standard deviation (SD) increase in average annual precipitation (adjusted odds ratio [aOR]: 1.34, 95% confidence interval [95% CI]: 1.13–1.58). Soil geochemistry was also associated with NTM positivity; a one SD increase in levels of sodium in the soil was associated with a 92% increased risk of culture positivity (aOR: 1.92, 95% CI: 1.46–2.52), and a one SD increase in soil manganese was associated with a 40.7% decreased risk (aOR: 0.59, 95% CI, 0.46–0.77). Species-specific analysis showed the same associations for M. abscessus (Table 3).

[Figure omitted. See PDF.]

Table 3. Generalized linear mixed models with a binomial distribution and a random intercept for zip code of demographic and environmental risk factors for NTM culture positivity in persons with cystic fibrosis, Florida, 2011–2017.

https://doi.org/10.1371/journal.pone.0259964.t003

Discussion

Because the prevalence of PNTM in Florida is so high [10,11], understanding whether there are environmental predictors or geospatial clustering of NTM is of interest to the CF community and public health. We found a significant cluster of NTM culture positivity, specifically M. abscessus, among persons with CF living in Florida in the southeast part of the state. M. abscessus is one of the most commonly isolated species of NTM from persons with CF, found in up to 16–68% of NTM-positive sputum cultures [22,23] and is considered difficult to treat due to high levels of inducible resistance to macrolides, the typical first-line therapy for infections with other NTM species [24]. It is therefore important to understand geographical areas of higher risk to acquiring this pathogen.

In addition to spatial clustering, we also found an association between sputum positivity and annual precipitation, soil sodium levels, and levels of soil manganese. The risk associated with higher sodium and lower manganese levels in soil is consistent with two other studies in the US, however these studies were not able to examine species-specific associations [8,10]. An association with rainfall has been recently identified for the province of Queensland, Australia; this relationship varied by species and geographic region [6]. One of the challenges of studying the environmental risks for PNTM disease is that the incubation period for PNTM is unknown, which creates uncertainty about the appropriate timescale for measuring exposure. While we studied incident infections to limit this potential bias, it will be important to conduct further analyses varying the timescale of possible exposure to environmental variables. Additionally, longitudinal, granular data on soil geochemistry is currently unavailable; the soil mineral concentrations throughout Florida likely vary more than represented in our study, limiting our ability to adjust these variables based on time and location. These are particularly important considerations as studies have shown that prolonged exposure to certain environments increases the risk of PNTM [9], so cumulative exposures to a variety of high risk sources may increase the risk of infection.

Climatic and environmental factors that contribute to increased mycobacterial abundance likely vary by region, making identification of a uniform set of determinants contributing to PNTM disease across national and global geographic areas challenging. Factors related to the built environment also likely interact with soil and water sources to affect the presence of mycobacteria in the environment; a recent study quantifying mycobacteria in showerheads found that both type of chlorination and showerhead material influenced their abundance [5]. Future studies which estimate risk related to mycobacterial abundance as well as components of the natural and built environmental will allow more complete and precise elucidation of these factors. Because persons with CF are at increased risk for NTM infection, continued studies to determine high-risk geographical areas and specific predictors of disease are critical so that precautions can be taken to reduce risk of exposure.

Acknowledgments

The authors thank the patients, care providers, and clinic coordinators at CF centers throughout the United States for their contributions to the Cystic Fibrosis Foundation Patient Registry.

Citation: Foote SL, Lipner EM, Prevots DR, Ricotta EE (2021) Environmental predictors of pulmonary nontuberculous mycobacteria (NTM) sputum positivity among persons with cystic fibrosis in the state of Florida. PLoS ONE 16(12): e0259964. https://doi.org/10.1371/journal.pone.0259964

References

1. Falkinham JO. Ecology of nontuberculous mycobacteria-where do human infections come from? Seminars in Respiratory and Critical Care Medicine. 2013;34: 95–102. pmid:23460009

2. Adjemian J, Olivier KN, Prevots DR. Epidemiology of pulmonary nontuberculous mycobacterial sputum positivity in patients with cystic fibrosis in the United States, 2010–2014. Annals of the American Thoracic Society. 2018;15: 817–825. pmid:29897781

3. Halstrom S, Price P, Thomson R. Review: Environmental mycobacteria as a cause of human infection. International Journal of Mycobacteriology. 2015;4: 81–91. pmid:26972876

4. Thomson R, Tolson C, Carter R, Coulter C, Huygens F, Hargreaves M. Isolation of nontuberculous mycobacteria (NTM) from household water and shower aerosols in patients with pulmonary disease caused by NTM. Journal of Clinical Microbiology. 2013;51: 3006–3011. pmid:23843489

5. Gebert MJ, Delgado-Baquerizo M, Oliverio AM, Webster TM, Nichols LM, Honda JR, et al. Ecological Analyses of Mycobacteria in Showerhead Biofilms and Their Relevance to Human Health. mBio. 2018;9: 2020. pmid:30377276

6. Thomson RM, Furuya-Kanamori L, Coffey C, Bell SC, Knibbs LD, Lau CL. Influence of climate variables on the rising incidence of nontuberculous mycobacterial (NTM) infections in Queensland, Australia 2001–2016. Science of the Total Environment. 2020;740: 139796. pmid:32563864

7. Lipner EM, French J, Bern CR, Walton-Day K, Knox D, Strong M, et al. Nontuberculous Mycobacterial Disease and Molybdenum in Colorado Watersheds. International Journal of Environmental Research and Public Health. 2020;17: 3854. pmid:32485845

8. Lipner EM, Knox D, French J, Rudman J, Strong M, Crooks JL. A Geospatial Epidemiologic Analysis of Nontuberculous Mycobacterial Infection: An Ecological Study in Colorado. Annals of the American Thoracic Society. 2017;14: 1–36. pmid:27870577

9. Adjemian J, Frankland TB, Daida YG, Honda JR, Olivier KN, Zelazny A, et al. Epidemiology of Nontuberculous Mycobacterial Lung Disease and Tuberculosis, Hawaii, USA. Emerging Infectious Diseases. 2017;23: 439–447. pmid:28221128

10. Adjemian J, Olivier KN, Seitz AE, Falkinham JO, Holland SM, Prevots DR. Spatial clusters of nontuberculous mycobacterial lung disease in the United States. American Journal of Respiratory and Critical Care Medicine. 2012;186: 553–558. pmid:22773732

11. Adjemian J, Olivier KN, Prevots DR. Epidemiology of Pulmonary Nontuberculous Mycobacterial Sputum Positivity in Patients with Cystic Fibrosis in the United Adjemian. Ann Am Thorac Soc. 2018;15: 817–825. pmid:29897781

12. Gross JE, Martiniano SL, Nick JA. Prevention of transmission of Mycobacterium abscessus among patients with cystic fibrosis. Current Opinion in Pulmonary Medicine. 2019;25: 646–653. pmid:31436542

13. Cystic Fibrosis Foundation. Patient Registry. [cited 6 Nov 2020]. Available: https://www.cff.org/Research/Researcher-Resources/Patient-Registry/.

14. Adjemian J, Olivier KN, Prevots DR. Nontuberculous Mycobacteria among Cystic Fibrosis Patients in the United States: Screening Practices and Environmental Risk. American Journal of Respiratory and Critical Care Medicine. 2014;190. pmid:25068291

15. Prevots DR, Adjemian J, Fernandez AG, Knowles MR, Olivier KN. Environmental risks for nontuberculous mycobacteria: Individual exposures and climatic factors in the cystic fibrosis population. Annals of the American Thoracic Society. 2014;11: 1032–1038. pmid:25068620

16. Bivand R, Keitt T, Rowlingson B. rgdal: Bindings for the “Geospatial” Data Abstraction Library. 2021. Available: https://proj.org.

17. Pebesma EJ, Bivand R. Classes and Methods for Spatial Data. 2005.

18. Bivand RS, Pebesma E, Gomez-Rubio V. Applied spatial data analysis with R, Second edition. Second. Springer, NY; 2013. Available: https://asdar-book.org/.

19. Melo RA, Rodriguez D, Zarruk D. Package “gmapsdistance”: Distance and Travel Time Between Two Points from Google Maps. 2018.

20. Kuldorff M, Information Management Services Inc. SaTScan—Software for the spatial, temporal, and space-time scan statistics. 2018.

21. R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2020. Available: https://www.r-project.org/.

22. Martiniano SL, Nick JA, Daley CL. Nontuberculous Mycobacterial Infections in Cystic Fibrosis. Thoracic surgery clinics. 2019;29: 95–108. pmid:30454926

23. Adjemian J, Prevots DR, Gallagher J, Heap K, Gupta R, Griffith D. Lack of adherence to evidence-based treatment guidelines for nontuberculous mycobacterial lung disease. Annals of the American Thoracic Society. 2014;11: 9–16. pmid:24236749

24. Daley CL, Iaccarino JM, Lange C, Cambau E, Wallace RJ, Andrejak C, et al. Treatment of Nontuberculous Mycobacterial Pulmonary Disease: An Official ATS/ERS/ESCMID/IDSA Clinical Practice Guideline. Clinical Infectious Diseases. 2020;71: 905–913. pmid:32797222