Volume 24, Number 10—October 2018
Mapping Histoplasma capsulatum Exposure, United States
Maps of Histoplasma capsulatum infection prevalence were created 50 years ago; since then, the environment, climate, and anthropogenic land use have changed drastically. Recent outbreaks of acute disease in Montana and Nebraska, USA, suggest shifts in geographic distribution, necessitating updated prevalence maps. To create a weighted overlay geographic suitability model for Histoplasma, we used a geographic information system to combine satellite imagery integrating land cover use (70%), distance to water (20%), and soil pH (10%). We used logistic regression modeling to compare our map with state-level histoplasmosis incidence data from a 5% sample from the Centers for Medicare and Medicaid Services. When compared with the state-based Centers data, the predictive accuracy of the suitability score–predicted states with high and mid-to-high histoplasmosis incidence was moderate. Preferred soil environments for Histoplasma have migrated into the upper Missouri River basin. Suitability score mapping may be applicable to other geographically specific infectious vectors.
Histoplasmosis is a regionally endemic mycosis caused by inhalation of the spores of the fungus Histoplasma capsulatum from the soil, leading occasionally to symptomatic disease and asymptomatic pulmonary nodules (1). Soils act as reservoirs for Histoplasma, especially where temperatures are 22°C–29°C and annual rainfall is 35–50 inches (2). According to a 1969 map by Edwards et al. of histoplasmin skin test reactivity among single-county-resident Navy recruits in the United States, H. capsulatum has been most prevalent in the Mississippi and Ohio River basins (3). Since then, much has changed with the environment, climate, and anthropogenic land use. Indeed, the advent of HIV/AIDS and the use of immunosuppressive medications for rheumatologic, dermatologic, and gastrointestinal conditions has unmasked areas of previously hidden histoplasmosis endemicity, nationally and globally (4). Recent outbreaks in Montana and Nebraska have led the Centers for Disease Control and Prevention to suggest that histoplasmosis is now endemic to these regions (5,6). State-level incidence rates of histoplasmosis among older patients also demonstrate a shift toward Nebraska and the northern Great Lakes area (7). In addition, a recent publication of surveillance data from 2011 through 2014 in 12 states suggests that histoplasmosis also occurs in Minnesota, Wisconsin, and Michigan, areas where histoplasmosis was not thought to be endemic (8).
Accurate maps of H. capsulatum prevalence assist in the early recognition, diagnosis, and treatment of acute infections. Pulmonary histoplasmosis may generate cancer-mimicking lung granulomas, causing false-positive radiographic images on high-resolution computed tomography and fluorodeoxyglycose positron emission tomography (FDG-PET) (9–11). Thus, an accurate understanding of the conditions for H. capsulatum is of practical epidemiologic and clinical value.
Edwards’ 1969 survey of histoplasmosis endemicity is impractical to repeat now or in the foreseeable future (3). It would be prohibitively expensive to survey a sufficiently large sample of the population and, in our mobile society, nearly impossible to identify an appropriate number of lifetime-single-county residents who are representative of the broader population. Furthermore, the histoplasmin skin test is no longer available, and another approach is needed. Our study objective was to model known environmental factors preferred by H. capsulatum and produce a large-scale map to enable risk stratification of patients on the basis of their geographic history.
We created a site suitability model for H. capsulatum weighted by 3 environmental soil characteristics measured across the continental United States within a geographic information system (ArcGIS 10.3.1; ESRI, Redlands, CA, USA). Suitability scores ranged from 1 through 9 within each soil characteristic; higher scores signified more favorable conditions for H. capsulatum. Selection of raster datasets was based on known criteria for histoplasmosis behavior—land cover, water, and soil pH—because cultivation and anthropogenic land use increase the spread of H. capsulatum from the soil to the air (1).
We obtained land cover data from the 2006 National Land Cover Database (12). This publicly available dataset applies 16 classes of land cover or land use types over the conterminous United States at a spatial resolution of 30 × 30 meters. The risk for histoplasmosis growth is highest in areas of high soil humidity (e.g., near rivers), and soil moisture affects the temperature at which H. capsulatum can survive (13–15).
Distance from Water
Because no continuous flood risk map exists for the United States, we defined water presence by attributes extracted from the 2006 National Land Cover Database (e.g., open water, woody wetlands, or emergent herbaceous woodlands). Suitability for H. capsulatum was modeled to decrease with Euclidean distance from these water types across 9 equally distributed categories. Suitability scores ranged from a best-suited value of 9, which occurred from 0 through 222 meters, to a low value of 1, at >1,777 meters (Table). Open water was assumed to preclude presence of histoplasmosis.
We obtained soil pH raster data from the University of Wisconsin Atlas of the Biosphere, compiled by the International Geosphere-Biosphere Programme Global Soils Data Task in 1998 (16). Original data projection of the Environmental Systems Research Institute ArcView gridded format with a spatial resolution of 10,000 meters, or 10 km, was reclassified into 9 categories. The highest score of 9 was assigned to pH values near 7.4, the pH of normal lungs, the presumed ideal environment for the H. capsulatum exposure of interest. Decreasing pH values were modeled in a descending fashion to a score of 4 for pH >8 and a score of 1 for pH <4.5 (Table).
We calculated the overall score for each ZIP code area as a weighted sum of assigned values from each of the 3 environmental factors. Each layer was preprocessed for analysis by mosaicking native coordinate systems across a common geographic coordinate system (World Geodetic System 1984, https://www.nga.mil/ProductsServices/GeodesyandGeophysics/Pages/WorldGeodeticSystem.aspx) and then clipping the data for each raster layer to the conterminous US border. Pixels for each layer were reclassified to a common scale of 1–9, ranging from least to most suitable for H. capsulatum. The unique map layers were then weighted and combined. Weighting was an iterative process to approximate the historical distribution of Histoplasma and recent histoplasmosis infection patterns, especially in the lower Mississippi River valley and, given the limitations of computational efficiency, to interpolate the final risk surface (17). The final suitability model included 3 layers of data: the 2006 National Land Cover Database with 16 classes of land cover or land use types (weighted 70%), Euclidean distance from the nearest open water source (20%), and soil pH (10%) (Table). A granular risk score was computed for each pixel. Next, the pixel values from the suitability map were extracted as points, and those points were joined to state and census tract polygons to provide a mean, median, and SD for exposure suitability scores aggregated to each spatial unit. To create a continuous exposure risk surface, we applied the inverse distance weights to interpolate a risk surface, which was then clipped to the US border.
We determined the 70%, 20%, and 10% weights through an iterative process based on Edwards’ 1969 survey county-level data for a few states (Tennessee with high within-state variability, Mississippi and Kansas with low variability, and Florida as a coastal state). During the developmental phase of the suitability score, we compared different proportions with state-level Medicare data, but there were no better models. Because of the lack of granularity of the state data, our sensitivity analyses may not have been very informative for detecting minor differences in state prevalence.
We summarized our data to the state level to permit comparison with previously published research (7). The reference standard was state-level histoplasmosis incidence data from a random 5% sample of the Centers for Medicare and Medicaid Services (CMS) claims data for 1999–2008 (7). Incident cases of endemic histoplasmosis required 1 inpatient claim or 2 outpatient claims over 7–90 days. We used 2 logistic regression models to assess the ability of our suitability histoplasmosis score to predict the binary outcome of correctly identifying states with high histoplasmosis incidence rates, defined as >5 cases/100,000 person-years, and states with mid-to-high incidence rates, defined as >3 cases/100,000 person-years. We generated the appropriate area under the receiver operator curves (AUCs) and calculated the Wilcoxon-Mann-Whitney test by using Stata version 12.0 (StataCorp LLC, College Station, TX, USA).
Figure 1 depicts the granular suitability score map with mean scores by US ZIP code. The data and model are applicable to estimates of Histoplasma in soils east of the Rocky Mountains because of lack of suitable soil habitats and limits in data availability. We calculated suitability score estimates west of the Rocky Mountains, but many geographic blocks had no or limited water data and thus may overestimate Histoplasma suitability (e.g., distances >3,000 meters from substantial amounts of water may preclude the presence of Histoplasma irrespective of land use). Figure 2 compares the state-level suitability scores with 1999–2008 state-level incidence data from Baddley et al. (7). Our suitability score model identified the states with high and mid-to-high histoplasmosis incidence rates with AUCs of 0.72 and 0.74, respectively, when compared with state-based diagnosis rates among elderly persons during 1999–2008.
Our histoplasmosis exposure suitability score predicted state-level histoplasmosis incidence rates from recent CMS data with an AUC >0.72. Our model confirms previous reports that preferred soil environments for Histoplasma have expanded into the upper Missouri River basin. This progression probably reflects changes in climate and human land use with increased cultivation in that region and urbanization of other regions of the country.
Our suitability score has limitations. We compared our score to state-level histoplasmosis incident cases as derived from administrative data and were thus unable to assess the validity of our score to predict disease incidence for a more granular geographic region. However, CMS is not likely to release county-level histoplasmosis infection rates because of the risk of identifying individual patients. The administrative data we used for our comparative analysis have their own limitations, namely reliance on codes from the International Classification of Disease, Ninth Revision for histoplasmosis rate estimates, with unknown sensitivity and specificity.
We also did not include known specific vectors of histoplasmosis exposure (e.g., droppings in areas of bird or bat roosts). This exclusion is particularly relevant, given the microfocus concept that Histoplasma exist in small environmental foci, resulting in localized outbreaks or other exposures in non–histoplasmosis-endemic areas (17,18). It is also worth noting that the coarse resolution of the soil pH pixels limited the sensitivity of the weighted overlay result. Local studies that use soil pH from US Department of Agriculture data within each state at a smaller resolution may be preferred to upsampling the coarse International Geosphere-Biosphere Programme data. The accuracy of our approach should be confirmed with soil pH measurements. In addition, as with any aggregate measure, our score is subject to the ecologic fallacy, whereby inferences about individual risk are estimated by an overall estimate from the geographic region in which the person resides. This suitability score could be combined with occupational and residential histories to risk-stratify patients for fungal lung disease risk and to improve interpretation of chest images.
The accurate diagnosis of histoplasmosis is challenging and often requires a battery of antigen detection, serology, and histopathology studies (19,20). Serologic tests are imperfect and, although often specific, can have low sensitivity, particularly for immunosuppressed patients. Knowledge of the distribution of H. capsulatum exposure and histoplasmosis infection is also critical for the appropriate regional use of imaging to diagnose fungal lung disease in immunocompromised persons and persons with lung malignancies. Investigators have questioned the utility of FDG-PET for diagnosing lung nodules in regions where granulomatous disease is prevalent (9,10,21). A recent meta-analysis concluded that FDG-PET was less specific for diagnosing malignancy in populations where lung mycosis is endemic than in those where it is not (11).
In conclusion, we generated a detailed suitability map for Histoplasma exposure east of the Rocky Mountains, an approach that can be updated as land use changes. Our suitability map may provide a detailed prediction of the risk for acute histoplasmosis infection, including areas of new endemicity. Accurate mapping of Histoplamosis exposure may also aid in the interpretation of chest images, particularly in the context of lung cancer screening programs (22). Future work might add elevation, soil moisture, water-holding capacity, and geologic data into the suitability model. Suitability score models have the potential to be applied to other infectious agents strongly associated with geographic-specific vectors, with the promise to inform healthcare and improve public health assessments and interventions.
Dr. Maiga is a resident in General Surgery at Vanderbilt University Medical Center and a former Veterans Affairs Quality Scholar at the Tennessee Valley Healthcare System. Her research interests include clinician decision-making and predictive modeling, primarily related to lung cancer.
A.W.M. is supported by the Office of Academic Affiliations, Department of Veterans Affairs National Quality Scholars program. E.L.G. is a recipient of the Department of Veterans Affairs, Veterans Health Administration, Health Services Research and Development Service Career Development award (10-024). M.C.A. is supported by National Institutes of Health National Cancer Institute grant K07 CA172294.
- Christie A. Histoplasmosis and pulmonary calcification; geographic distribution. Am J Trop Med Hyg. 1951;31:742–52.
- Sykes J, Taboada J. Histoplasmosis. In: Sykes J, editor. Canine and feline infectious disease. St. Louis: Elsevier Inc; 2014. p. 587–98.
- Edwards LB, Acquaviva FA, Livesay VT, Cross FW, Palmer CE. An atlas of sensitivity to tuberculin, PPD-B, and histoplasmin in the United States. Am Rev Respir Dis. 1969;99(Suppl):1–132.
- Bahr NC, Antinori S, Wheat LJ, Sarosi GA. Histoplasmosis infections worldwide: thinking outside of the Ohio River valley. Curr Trop Med Rep. 2015;2:70–80.
- Centers for Disease Control and Prevention. Histoplasmosis in a state where it is not known to be endemic—Montana, 2012-2013. MMWR Morb Mortal Wkly Rep. 2013;62:834–7.
- Centers for Disease Control and Prevention (CDC). Outbreak of histoplasmosis among industrial plant workers—Nebraska, 2004. MMWR Morb Mortal Wkly Rep. 2004;53:1020–2.
- Baddley JW, Winthrop KL, Patkar NM, Delzell E, Beukelman T, Xie F, et al. Geographic distribution of endemic fungal infections among older persons, United States. Emerg Infect Dis. 2011;17:1664–9.
- Armstrong PA, Jackson BR, Haselow D, Fields V, Ireland M, Austin C, et al. Multistate epidemiology of histoplasmosis, United States, 2011–2014. Emerg Infect Dis. 2018;24:425–31.
- Grogan E, Deppen SA, Ballman K, Andrade GM, Verdail FC, Aldrich MC, et al. Accuracy of FDG-PET to diagnose lung cancer in the ACOSOG Z4031 trial. Alliance for Clinical Trials in Oncology; 2012 June 29–July 1; Chicago, IL.
- Deppen S, Putnam JB Jr, Andrade G, Speroff T, Nesbitt JC, Lambright ES, et al. Accuracy of FDG-PET to diagnose lung cancer in a region of endemic granulomatous disease. Ann Thorac Surg. 2011;92:428–32, discussion 433.
- Deppen SA, Blume JD, Kensinger CD, Morgan AM, Aldrich MC, Massion PP, et al. Accuracy of FDG-PET to diagnose lung cancer in areas with infectious lung disease: a meta-analysis. JAMA. 2014;312:1227–36.
- Fry J, Xian G, Jin S, Dewitz J, Homer C, Yang L, et al. Completion of the 2006 National Land Cover Database for the conterminous United States. Photogramm Eng Remote Sensing. 2011;77:858–64.
- Mahvi TA. Factors governing the epidemiology of Histoplasma capsulatum in soil. Mycopathol Mycol Appl. 1970;41:167–76.
- Goodman NL, Larsh HW. Environmental factors and growth of Histoplasma capsulatum in soil. Mycopathol Mycol Appl. 1967;33:145–56.
- Menges RW, Furcolow ML, Larsh HW, Hinton A. Laboratory studies on histoplasmosis. I. The effect of humidity and temperature on the growth of Histoplasma capsulatum. J Infect Dis. 1952;90:67–70.
- University of Wisconsin-Madison Center for Sustainability and the Global Environment. International Geosphere-Biosphere Programme-DIS. SoilData, version 0. A program for creating global soil-property databases, IGBP Global Soils Data Task, France. 1998 [cited 2018 Jul 22]. https://nelson.wisc.edu/sage/data-and-models/datasets.php
- Benedict K, Mody RK. Epidemiology of histoplasmosis outbreaks, United States, 1938–2013. Emerg Infect Dis. 2016;22:370–8.
- Nett RJ, Skillman D, Riek L, Davis B, Blue SR, Sundberg EE, et al. Histoplasmosis in Idaho and Montana, USA, 2012-2013. Emerg Infect Dis. 2015;21:1071–2.
- McKinsey DS, McKinsey JP. Pulmonary histoplasmosis. Semin Respir Crit Care Med. 2011;32:735–44.
- Wheat LJ, Azar MM, Bahr NC, Spec A, Relich RF, Hage C. Histoplasmosis. Infect Dis Clin North Am. 2016;30:207–27.
- Croft DR, Trapp J, Kernstine K, Kirchner P, Mullan B, Galvin J, et al. FDG-PET imaging and the diagnosis of non-small cell lung cancer in a region of high histoplasmosis prevalence. Lung Cancer. 2002;36:297–301.
- Aberle DR, Adams AM, Berg CD, Black WC, Clapp JD, Fagerstrom RM, et al.; National Lung Screening Trial Research Team. Reduced lung-cancer mortality with low-dose computed tomographic screening. N Engl J Med. 2011;365:395–409.
TableCite This Article
Original Publication Date: 9/12/2018