Nearly 80% of antibiotics in the United States are sold for use in livestock feed.

Back

Read Article on line

High-Density Livestock Operations, Crop Field Application of Manure, and Risk of Community-Associated Methicillin-Resistant Staphylococcus aureus Infection in Pennsylvania FREE

Joan A. Casey, MA1,2; Frank C. Curriero, PhD, MA1,3; Sara E. Cosgrove, MD, MS2,4; Keeve E. Nachman, PhD, MHS1,5,6; Brian S. Schwartz, MD, MS1,2,7

ABSTRACT

ABSTRACT | METHODS | RESULTS | DISCUSSION | ARTICLE INFORMATION | REFERENCES

Importance  Nearly 80% of antibiotics in the United States are sold for use in livestock feeds. The manure produced by these animals contains antibiotic-resistant bacteria, resistance genes, and antibiotics and is subsequently applied to crop fields, where it may put community members at risk for antibiotic-resistant infections.

Objective  To assess the association between individual exposure to swine and dairy/veal industrial agriculture and risk of methicillin-resistant Staphylococcus aureus (MRSA) infection.

Design, Setting, and Participants  A population-based, nested case-control study of primary care patients from a single health care system in Pennsylvania from 2005 to 2010. Incident MRSA cases were identified using electronic health records, classified as community-associated MRSA or health care–associated MRSA, and frequency matched to randomly selected controls and patients with skin and soft-tissue infection. Nutrient management plans were used to create 2 exposure variables: seasonal crop field manure application and number of livestock animals at the operation. In a substudy, we collected 200 isolates from patients stratified by location of diagnosis and proximity to livestock operations.

Main Outcomes and Measures  Community-associated MRSA, health care–associated MRSA, and skin and soft-tissue infection status (with no history of MRSA) compared with controls.

Results  From a total population of 446 480 patients, 1539 community-associated MRSA, 1335 health care-associated MRSA, 2895 skin and soft-tissue infection cases, and 2914 controls were included. After adjustment for MRSA risk factors, the highest quartile of swine crop field exposure was significantly associated with community-associated MRSA, health care-associated MRSA, and skin and soft-tissue infection case status (adjusted odds ratios, 1.38 [95% CI, 1.13-1.69], 1.30 [95% CI, 1.05-1.61], and 1.37 [95% CI, 1.18-1.60], respectively); and there was a trend of increasing odds across quartiles for each outcome (P ≤ .01 for trend in all comparisons). There were similar but weaker associations of swine operations with community-associated MRSA and skin and soft-tissue infection. Molecular testing of 200 isolates identified 31 unique spa types, none of which corresponded to CC398 (clonal complex 398), but some have been previously found in swine.

Conclusions and Relevance  Proximity to swine manure application to crop fields and livestock operations each was associated with MRSA and skin and soft-tissue infection. These findings contribute to the growing concern about the potential public health impacts of high-density livestock production.

The US Food and Drug Administration reports that 80% of antibiotics in the United States are used in food animals, mainly to promote growth of livestock in high-density production.1,2 In this industrial model, thousands of livestock animals are housed together and fed subtherapeutic doses of antimicrobial agents.3 This can select for antibiotic-resistant bacteria like methicillin-resistant Staphylococcus aureus (MRSA).4 About 75% of administered antibiotics are not absorbed by the animal and end up in manure.5 In addition to antibiotics, antibiotic-resistant bacteria and resistance genes have been isolated from manure and are known to persist in the environment.6 Application of manure to crop fields close to human dwellings may lead to increased risk of antibiotic-resistant infections.7


High prevalence of MRSA colonization among pig farmers was first noted in the Netherlands in 2005.8 Transmission of MRSA among cows, calves, and humans has also been reported.9,10 A US study reported that 45% of swine workers were colonized, which is 30 times the national average.11,12 Few studies have assessed risk of MRSA infection. A Danish case-control study (21 cases) found that living or working on a farm was a risk factor for MRSA infection.13 Other studies of MRSA infection have described the spread of infection within pig farming families.14,15 To our knowledge, the impact of proximity to high-density livestock production on risk of MRSA infection has not been evaluated. Two studies assessed community livestock density and found no association with MRSA colonization.16,17 Most studies have implicated clonal complex 398 (CC398) as the most common variant harbored by swine.9,12- 15,17,18


In the United States, community-associated MRSA (CA-MRSA) now accounts for more than half of all noninvasive MRSA infections.19- 22 CA-MRSA often occurs as a skin and soft-tissue infection (SSTI) in young, otherwise healthy individuals who lack common health care risk factors.23- 25 In contrast, health care–associated MRSA (HA-MRSA) generally occurs in older individuals with comorbidities. Now, the epidemiology and microbiology of 2 epidemics previously considered distinct are merging.26- 28


The possible contribution of high-density livestock production to these epidemics has been inadequately studied.29 US data are needed to evaluate which livestock species pose risk, what aspects of the operations and manure handling are involved, and the magnitude of the risk. We conducted a population-based, nested case-control study using electronic health record (EHR) data from the Geisinger Health System, an integrated health services organization with 4 hospitals and 41 outpatient clinics. We evaluated residential proximity to swine and dairy/veal high-density livestock operations and manure-applied crop fields in relation to CA-MRSA, HA-MRSA, and SSTI (without a history of MRSA infection or colonization) case status.


METHODS

ABSTRACT | METHODS | RESULTS | DISCUSSION | ARTICLE INFORMATION | REFERENCES

Data Sources

To select cases and controls, we used EHR data from January 1, 2005, to February 9, 2010, from 446 480 patients with a Geisinger primary care provider. Geisinger’s primary care patients represent the general population in the region.30 The study area comprised the health system’s primary care market and bordering counties totaling 38 counties in Pennsylvania (Figure 1); approximately 3.8 million people resided in these counties according to the 2000 US Census. The study was approved by institutional review boards at the Geisinger Health System and the Johns Hopkins Bloomberg School of Public Health.


Figure 1.

Crop Field Locations and Methicillin-Resistant Staphylococcus aureus (MRSA) Rates


Crop field locations and rates of MRSA per 1000 Geisinger Health System (GHS) primary care patients in townships, boroughs, and cities. Rates in communities with fewer than 50 GHS patients were not estimated. The map demonstrates that crop fields were often located in areas with a range of human population densities.


 Image not available.

View Large  |  Save Figure  |  Download Slide (.ppt)

Pennsylvania Act 38 of 2005 requires high-density livestock operations to develop and implement nutrient management plans (NMPs) for manure handling. We defined high-density livestock operations as agricultural operations where animal density exceeded 2 animal equivalent units (AEUs, 1000 pounds of live weight) per acre and where total AEUs exceeded 8; or operations that exceeded 300 AEUs. When livestock operations have insufficient land for manure application, manure is exported to other areas and applied to crop fields. We focused on swine and dairy/veal operations because of previously identified links with MRSA.8- 10,13- 15 The NMPs provided data on livestock operation location, animal type, livestock quantity (AEUs), amount of manure applied to crop fields by season, and crop field area and location. We obtained NMPs from the Pennsylvania Department of Environmental Protection and County Conservation Districts.


Case Ascertainment and Control Selection

Incident MRSA cases were identified primarily using laboratory cultures and secondarily by diagnosis codes (eg, International Classification of Diseases, Ninth Revision [ICD-9]) that indicated MRSA infection, as previously described.22 Cases were then classified as either CA-MRSA or HA-MRSA based on presence of health care risk factors (eg, hospitalization, surgery, dialysis, nursing home residence, indwelling device)22,31 or diagnosis more than 2 days after hospital admission using ICD-9 codes21,23,32 and Current Procedural Terminology codes. We then randomly selected patients with SSTI but no history of MRSA using 29 ICD-9 codes (eg, carbuncle, furuncle, abscess)22 and controls with no history of MRSA, and we frequency matched both groups with case patients by age (0-6, 7-18, 19-45, 46-62, 62-74, ≥75 years), sex, and diagnosis or an outpatient encounter in the same year as MRSA diagnosis. The SSTI cases were evaluated as a separate case group because some SSTIs occurring during the study period were likely to have been caused by MRSA but not diagnosed as such, and high-density livestock production could cause SSTIs from other bacteria. Therefore, we selected patients with SSTIs without reference to any specific pathogen. If a control had multiple outpatient encounters during the year, a single encounter was randomly selected as the date for exposure assignment.


Geographic Location of Patients, Livestock Operation, and Crop Fields

We identified the latitude and longitude of patient addresses using ArcGIS, version 10 (Esri),33 and calculated MRSA rates by community (Figure 1). The latitude and longitude of livestock operation addresses were located using Google Earth with visual confirmation of the presence of a barn at the address. Three different methods were used to locate crop fields because 131 fields (17.3%) were missing address data (Figure 2). In the absence of an address, we used the county and community of the field from the NMP to locate the crop field.


Figure 2.

Three Methods Used to Identify and Locate Crop Fields


A, Aerial photograph (top) or map (bottom) were located using Google Earth (n = 135). B, Operation addresses known and located using ArcGIS, version 10 (Esri) (n = 420). C, County and township known and addresses were located by identifying cropland, hay land, and pastureland on a land use map and randomly selecting a point within the eligible land use types (n = 131).


 Image not available.

View Large  |  Save Figure  |  Download Slide (.ppt)

Exposure Assessment

We estimated each individual’s exposure to livestock operations and to manure-applied crop fields, for swine and dairy/veal, using an inverse distance-squared approach34 in R, version 2.14.2 (R Foundation for Statistical Computing). For livestock operations, we used the equation

Image not available.

where n is the number of operations, ai is AEUs of livestock at operation I, and dij2 is the squared distance (in meters) between operation centroid i and patient j. For crop field manure application during the season of diagnosis or visit, we used the equation

Image not available.

where n is the number of crop fields, ci is the concentration of manure (gallons per square meter) applied during the season of diagnosis or visit at field i, and dij2 is the squared distance (in meters) between the crop field centroid i and patient j. This resulted in exposure units of gallons per quartic meter (m4). For missing seasonal application data (40% of fields), we used the regional seasonal average. Exposure variables were modeled as quartiles because they were not normally distributed.


Genetic Typing of MRSA Isolates

In a prospective substudy that occurred from January to December 2012, we identified all primary care patients with a culture-confirmed MRSA infection. We then randomly selected isolates from patients stratified by location of diagnosis (inpatient or outpatient) and presence or absence of a livestock operation or crop field in their zip code of residence. We oversampled outpatient isolates from zip codes with high-density livestock production. All isolates were characterized by DNA sequence analysis of the protein A gene variable repeat region (spa typing).35 The spa types were also assigned a Ridom StaphType, using the Ridom SpaServer database (available at http://www.spaserver.ridom.de).36 Panton-Valentine leukocidin (PVL) genes were detected by polymerase chain reaction.37


Statistical Analysis

To compare the 4 study groups, we used Kruskal-Wallis tests for comparisons of means and χ2 tests for comparisons of proportions. We used multilevel multinomial (2 MRSA groups and controls) and multilevel logistic regression (SSTI cases and controls) to assess the association between exposure and case status. We adjusted models for predictors obtained from the EHR: age, sex, race/ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, and other); ever-smoking status using ICD-9 codes22; antibiotic order in the 2 years preceding diagnosis; and Medical Assistance for health insurance, residential community (city, borough, or township), and community socioeconomic deprivation.33Medical Assistance is a means-tested program, which we used as a surrogate indicator for low socioeconomic status.38,39 It was identified for each encounter by 1 of 24 codes and modeled as ever received (if received for >2 encounters) vs never received. The livestock operation model was additionally adjusted for season of infection or visit. Final models were selected based on associations reported in our group’s earlier article,22 and adequate model fit was confirmed using Pearson goodness-of-fit tests and likelihood ratio tests. Semivariograms were used to determine that spatial autocorrelation did not account for the results.40 We used a 2-sided type 1 error rate of 0.05 as the threshold for statistical significance and made no adjustments for multiple comparisons. Logistic and multinomial logistic regression analyses were completed using Stata statistical software version 11.2 (StataCorp Inc) and the GLLAMM program (www.gllamm.org).


We calculated population attributable fraction (PAF) for CA-MRSA and SSTI cases by using the adjusted odds ratios (AORs) for the third and fourth quartile of crop field swine manure exposure obtained from the final multivariable logistic regression model using the formula PAF = (AOR – 1)/AOR × (proportion of cases exposed to risk factor).41 The PAF for the third and fourth quartiles were summed to calculate the final PAF. Confidence intervals were computed using the delta method with the postestimation command nlcom in Stata.


We performed several sensitivity analyses. First, we evaluated 2 other time windows for crop field manure exposures: a 1-year and a 30-day window before either the date of MRSA diagnosis or the date of a preceding SSTI diagnosis. Second, models were repeated using different methods to identify crop fields (Figure 2). Third, we evaluated effect modification by community type and season. Fourth, analysis was performed to evaluate whether distance to nearest crop field, independent of manure application, was associated with disease outcomes. Finally, we repeated analyses using only culture-confirmed MRSA cases.


RESULTS

ABSTRACT | METHODS | RESULTS | DISCUSSION | ARTICLE INFORMATION | REFERENCES

Patients, Farms, and Crop Fields Included in the Study

A total of 1734 CA-MRSA and 1519 HA-MRSA cases were identified between January 2005 and February 2010. Most cases (72.3%) were identified by a positive MRSA culture. We frequency matched 3336 SSTIs and 3336 controls to the cases. Of the 9925 patients, we were able to identify latitude and longitude coordinates for home addresses for 87.5%. Thirty-one patients (0.35%) were excluded because they did not reside in counties for which we collected NMPs. This resulted in a total of 1539 CA-MRSA, 1335 HA-MRSA, and 2895 SSTI cases and 2914 controls with residences in 32 counties and 574 communities defined by minor civil divisions, ie, census-designated divisions of a county.


We collected NMPs and identified the location of 326 high-density livestock operations, 123 swine and 203 dairy/veal operations, in 27 counties and 168 townships. Ninety-eight of the swine operations (79.7%) and 71 of the dairy/veal farms (35.0%) exported at least a portion of their manure to a total of 424 crop fields. Crop fields were located in 29 counties, 8 boroughs, and 218 townships, for an average of 3.4 crop fields per township. A reported 637 266 595 gallons of manure was applied to crop fields annually. Among those in the highest quartiles of swine livestock operation exposure and swine crop field manure exposure, the median (IQR) distance to the nearest farm or field was 4.6 km (3.1-6.7 km) and 4.0 km (2.6-7.5 km), respectively. The Pearson r correlation between the individual-level exposure assignments from the 2 models was 0.59 (both log transformed).


Demographic Characteristics of Cases and Controls

There were no differences in demographic characteristics between patients included in the study and those excluded because we could not map their home address (n = 1211), except that patients with SSTIs who were included were significantly more likely to have ever smoked than were those who were excluded (eTable 1 in the Supplement). Compared with controls, patients in the CA-MRSA group were significantly younger; they were more likely to be black, to smoke, and to receive Medical Assistance; and they lived in communities with greater socioeconomic deprivation (Table 1). White race/ethnicity (95.3% of the study population) was associated with higher levels of swine and dairy/veal crop field manure exposure (Table 2). Exposures to manure were also the lowest among those who had not received an antibiotic prescription in the 2 years prior to diagnosis or visit, those who had received Medical Assistance, and those residing in cities.


Table 1.  Patient Demographic and Clinical Characteristics

Image not available.

View Large  |  Save Table  |  Download Slide (.ppt)

Table 2.  Characteristics of CA-MRSA Group by Crop Field Manure Exposure Quartiles (n = 1539)

Image not available.

View Large  |  Save Table  |  Download Slide (.ppt)

Association of Crop Field Manure Exposure With Case Status

After adjusting for potential confounding variables, we found a significantly increased odds of CA-MRSA, HA-MRSA, and SSTI with higher swine manure exposure (fourth vs first quartile AORs, 1.38 [95% CI, 1.13-1.68], 1.30 [95% CI, 1.05-1.61], and 1.37 [95% CI, 1.18-1.60], respectively; P = .01, P < .001, and P < .001 for trend, respectively) (Table 3). The fourth quartile (vs first quartile) of dairy/veal exposure was also associated with increased odds of CA-MRSA (AOR, 1.24 [95% CI, 1.01-1.52]). The PAFs (95% CI) for the highest 2 quartiles of swine crop field manure exposure for CA-MRSA infection and SSTI were 10.7% (5.0-16.4%) and 11.5% (7.0-16.0%), respectively.


Table 3.  Association of Seasonal Crop Field Manure Exposure and Livestock Operation Exposure With CA-MRSA, HA-MRSA, and SSTI From Full Multilevel Multivariate Modela

Image not available.

View Large  |  Save Table  |  Download Slide (.ppt)

Association of High-Density Livestock Operations With Case Status

The fourth quartile vs first quartile of swine livestock operation exposure was associated with increased odds of CA-MRSA and SSTI (AORs, 1.25 [95% CI, 0.99-1.58] and 1.27 [95% CI, 1.08-1.50], respectively; P = .04 and P = .002, respectively) but not HA-MRSA (Table 3). No associations were seen with dairy/veal operations. Receipt of Medical Assistance remained independently associated with MRSA and SSTI outcomes in all adjusted analyses.


MRSA Isolate Substudy

A total of 200 isolates were randomly selected from 1128 patients with isolates (Table 4). Of these, 133 (66.5%) were common community strains (spa types t008, t024, t064, t206, t211); 33 (16.5%) were common health care strains (t002, t010, t062, t105); and the remainder were divided among 22 additional types. Notably, none was a spa type associated with CC398. A total of 27 community-onset isolates were PVL-negative.


Table 4.  Summary of MRSA Genetics

Image not available.

View Large  |  Save Table  |  Download Slide (.ppt)

Sensitivity Analyses

Evaluation of the 2 other time windows for manure application revealed slightly attenuated but otherwise similar results, except in the case of annual dairy/veal exposure, where the association strengthened for both CA-MRSA and SSTI cases (fourth vs first quartile AORs, 1.41 [95% CI, 1.13-1.77] and 1.26 [95% CI 1.08-1.48], respectively; P = .01 and P = .02, respectively) (eTable 2 in the Supplement). Exclusion of crop fields without aerial photographs and/or addresses, including interaction terms for community type or season, and restricting analysis to culture-confirmed MRSA cases did not alter results. Distance to nearest crop field, independent of manure application, was not associated with MRSA or SSTI.


DISCUSSION

ABSTRACT | METHODS | RESULTS | DISCUSSION | ARTICLE INFORMATION | REFERENCES

High-density swine production was an independent risk factor for CA-MRSA and HA-MRSA infection and SSTI in patients without a history of MRSA. There was evidence that both a crop field manure application model that incorporated distance, swine manure volume, season of application, and field area, as well as a livestock operation model that incorporated distance and swine count, were each associated with increased risk of infection. Exposure assignments from these 2 models were moderately correlated, suggesting 2 independent sources of risk. Furthermore, while livestock operations are easily identifiable and fixed, manure-applied crop fields are not necessarily easily identifiable, are scattered, and may change from year to year. Associations with dairy/veal operations were less consistent and weaker than those for swine operations, which was not surprising, given the limited number studies on the topic. Our data also suggest that approximately 11% of CA-MRSA and SSTI cases in the study population could be attributed to crop field application of swine manure.


The association with SSTIs is notable because there were more than 50 000 such cases, compared with 4000 MRSA cases, in the region since 2001.22 Crop field manure application may lead to SSTI either by causing MRSA infection that goes undiagnosed as such or by exposure to other bacteria that cause SSTI. We cannot attribute these infections to any single pathogen. Subanalysis of 258 SSTI cases (9%) identified from culture data as methicillin-susceptible S aureus (MSSA) revealed an association with swine livestock exposure (AOR, 1.70 [95% CI, 1.15-2.52]) that was stronger than for all SSTIs together. The public health burden of SSTIs25,42 warrants future study of links to high-density livestock production. Concerning the association with HA-MRSA, this finding does not imply that livestock-associated strains are in hospitals. Rather, these older patients may have been colonized in the community with livestock strains and then had health care contact (eg, hospitalization, surgery, dialysis, indwelling device) before infection.


To our knowledge, no studies have examined the association between high-density livestock operations and MRSA infection in the community, and only 2 studies have assessed risk of colonization.16,17 Study of colonization may not be useful when evaluating the risk of MRSA infection associated with industrial agriculture because, first, environmental exposure or direct contact may more often lead to CA-MRSA infection than does colonization43 and second, MRSA colonization may not persist for more than 24 hours after livestock exposure.44 Prior studies limited their ability to identify risk from high-density livestock exposure by including only patients with CC398,13,45 by studying only people who lived or worked on farms,9,12,46- 48 and by not including crop field manure exposure.


A unique aspect of this study is the consideration of risk from manure application to crop fields. Sensitivity analyses across time (eg, annual, seasonal, and 30-day exposure) and across space (eg, restricting analysis to fields with aerial photographs or maps) did not substantively change results. Our findings are likely not an artifact of rural residents in the study population because, while manure-applied crop fields were associated with risk, residence in a township (where 98.5% of crop fields are located) and distance to nearest crop field independent of manure application were not associated with the 3 outcomes.


Previous studies have used community-level livestock density as a marker of livestock exposure.16,17,45 We estimated exposure at the individual level, incorporating several features of these operations and fields. Operation and field sites were confirmed using Google Earth rather than physical inspection of the area, which may have led to some misclassification. Additional exposure misclassification may have occurred since the exposure assessment did not capture all smaller farms. In Pennsylvania, only large operations, specifically concentrated animal feeding operations, or operations with greater than 2 AEUs per acre are required to produce NMPs. The 2007 US Census of Agriculture revealed that our exposure models included 89% of swine AEUs present in the study area. Because characteristics of high-density livestock practices differ by and within states, our results may not apply in full to other areas with these operations.


A limitation of this study is that we did not measure how livestock-associated pathogens may reach individuals in a community, but several ways have been considered in prior studies. Aerosolized MRSA has been isolated from the air up to 150 m downwind and from the soil up to 300 m downwind of swine operations.49,50 Studies have also reported that tetracycline, the most commonly used antimicrobial agent used in food-producing animals in the United States,1 was found at concentrations sufficient to select for resistant bacteria in both swine manure and manure-fertilized soil.51 These studies provide biologic plausibility for our findings, but more research is needed to establish sources, media, routes, and behaviors that may lead to infection. Other animals,52 slaughterhouses,53 and meat consumption54 may be involved.


We did not find any spa types associated with clonal complex CC398. Most MRSA found on farms or in farmers has been characterized as CC398, though other types have been identified.9,12,13,45,55,56 However, few studies have assessed MRSA infection, instead focusing on MRSA carriage, and those that do come from Europe14,15,57,58 and/or deal exclusively with CC398 MRSA.58- 60 The most common strains on US farms are not well known. In Ohio, the closest state to Pennsylvania for which there are data, common health care strains were most often isolated from swine; CC398 was secondary.53 We were limited in that genetic data did not come from the same time period as the main study, and isolate strains could have changed over time.


We used ICD-9 codes to classify MRSA cases as CA or HA and to derive the predictor variables used in the analysis. The primary HA epidemiologic factors31 are well captured by the EHR. Ever-smoking status can also be effectively obtained using ICD-9 codes.61 Unfortunately, the EHR did not directly capture data on individual-level socioeconomic status or occupation, so our results could have been due to these variables. To evaluate this potential confounding, we adjusted the final models for 2 surrogates for socioeconomic status: Medical Assistance at the individual level and community socioeconomic deprivation at the community level. As expected, patients with CA-MRSA were more likely to receive Medical Assistance and more likely to live in more deprived communities than controls, associations that persisted in the adjusted models.62,63 Previous literature has suggested increased prevalence of MRSA colonization in livestock workers.8 However, given the low prevalence of this employment (1% of population in any agricultural work, according to US Census data) even in townships and boroughs in the study area, it seems unlikely that farming occupation could account for the reported associations. Moreover, inclusion of community-level percentage of agricultural workers in the multilevel model did not change associations.


In conclusion, proximity to and size of high-density livestock production were associated MRSA infection and SSTI, and the population-attributable fraction of crop field manure application exceeded 10% for CA-MRSA and SSTI. The findings contribute to the growing concern about the potential public health impacts of high-density livestock production.


ARTICLE INFORMATION

ABSTRACT | METHODS | RESULTS | DISCUSSION | ARTICLE INFORMATION | REFERENCES

Corresponding Author: Brian S. Schwartz, MD, MS, Department of Environmental Health Sciences, Johns Hopkins Bloomberg School of Public Health, 615 N Wolfe St, Room W7041, Baltimore, MD 21205 (bschwart@jhsph.edu).


Accepted for Publication: June 11, 2013.


Published Online: September 16, 2013. doi:10.1001/jamainternmed.2013.10408.


Author Contributions: Ms Casey and Dr Schwartz had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.


Study concept and design: Casey, Curriero, Nachman, Schwartz.


Acquisition of data: Casey, Schwartz.


Analysis and interpretation of data: Casey, Curriero, Cosgrove, Nachman, Schwartz.


Drafting of the manuscript: Casey.


Critical revision of the manuscript for important intellectual content: Casey, Curriero, Cosgrove, Nachman, Schwartz.


Statistical analysis: Casey, Curriero.


Obtained funding: Casey, Nachman, Schwartz.


Administrative, technical, or material support: Schwartz.


Study supervision: Cosgrove.


Conflict of Interest Disclosures: None reported.


Funding/Support: This study was jointly funded by the New York University–Geisinger Seed Grant Program and by the Johns Hopkins Center for a Livable Future. Additional support was provided by the Johns Hopkins Sommer Scholarship and the National Institute of Environmental Health Sciences Training Grant ES07141 (J Casey).


Role of the Sponsor: The sponsors had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.


Additional Contributions: We thank Dione G. Mercer, BS, BA, for study management; Joseph J. DeWalle, BS, for assistance with patient geocoding; Sy H. Brandau, AAS, for isolate acquisition, handling, and shipping; and Bo Shopsin, MD, PhD, and Hannah R. Rose, BA, for microbiologic testing of isolates. These individuals received compensation from the study sponsors for their contributions.


REFERENCES

ABSTRACT | METHODS | RESULTS | DISCUSSION | ARTICLE INFORMATION | REFERENCES

1

Love  DC, Davis  MF, Bassett  A, Gunther  A, Nachman  KE.  Dose imprecision and resistance: free-choice medicated feeds in industrial food animal production in the United States. Environ Health Perspect. 2011;119(3):279-283.

PubMed   |  Link to Article

2

Bottemiller H. Slaughter reintroduces bill to curb Ag antibiotics. http://www.foodsafetynews.com/2011/03/rep-slaughter-reintroduces-bill-to-limit-antibiotic-use-in-ag/#.UoJpQOLNld0. Accessed November 12, 2013.

3

Marshall  BM, Levy  SB.  Food animals and antimicrobials: impacts on human health. Clin Microbiol Rev. 2011;24(4):718-733.

PubMed   |  Link to Article

4

Silbergeld  EK, Graham  J, Price  LB.  Industrial food animal production, antimicrobial resistance, and human health. Annu Rev Public Health. 2008;29(1):151-169.

PubMed   |  Link to Article

5

Elmund  GK, Morrison  SM, Grant  DW, Nevins  SM.  Role of excreted chlortetracycline in modifying the decomposition process in feedlot waste. Bull Environ Contam Toxicol. 1971;6(2):129-132.

PubMed   |  Link to Article

6

Chee-Sanford  JC, Mackie  RI, Koike  S,  et al.  Fate and transport of antibiotic residues and antibiotic resistance genes following land application of manure waste. J Environ Qual. 2009;38(3):1086-1108.

PubMed   |  Link to Article

7

Kumar  K, Gupta  SC, Chander  Y, Singh  AK.  Antibiotic use in agriculture and its impact on the terrestrial environment. Adv Agron. 2005;87:1-54.

8

Voss  A, Loeffen  F, Bakker  J, Klaassen  C, Wulf  M.  Methicillin-resistant Staphylococcus aureus in pig farming. Emerg Infect Dis. 2005;11(12):1965-1966.

PubMed   |  Link to Article

9

Graveland  H, Wagenaar  JA, Heesterbeek  H, Mevius  D, van Duijkeren  E, Heederik  D.  Methicillin-resistant Staphylococcus aureus ST398 in veal calf farming: human MRSA carriage related with animal antimicrobial usage and farm hygiene. PLoS One. 2010;5(6):e10990.

PubMed   |  Link to Article

10

Juhász-Kaszanyitzky  E, Jánosi  S, Somogyi  P,  et al.  MRSA transmission between cows and humans. Emerg Infect Dis. 2007;13(4):630-632.

PubMed   |  Link to Article

11

Gorwitz  RJ, Kruszon-Moran  D, McAllister  SK,  et al.  Changes in the prevalence of nasal colonization with Staphylococcus aureus in the United States, 2001-2004. J Infect Dis. 2008;197(9):1226-1234.

PubMed   |  Link to Article

12

Smith  TC, Male  MJ, Harper  AL,  et al.  Methicillin-resistant Staphylococcus aureus (MRSA) strain ST398 is present in midwestern U.S. swine and swine workers. PLoS One. 2009;4(1):e4258.

PubMed   |  Link to Article

13

Lewis  HC, Mølbak  K, Reese  C,  et al.  Pigs as source of methicillin-resistant Staphylococcus aureus CC398 infections in humans, Denmark. Emerg Infect Dis. 2008;14(9):1383-1389.

PubMed   |  Link to Article

14

Hartmeyer  GN, Gahrn-Hansen  B, Skov  RL, Kolmos  HJ.  Pig-associated methicillin-resistant Staphylococcus aureus: family transmission and severe pneumonia in a newborn. Scand J Infect Dis. 2010;42(4):318-320.

PubMed   |  Link to Article

15

Huijsdens  XW, van Dijke  BJ, Spalburg  E,  et al.  Community-acquired MRSA and pig farming. Ann Clin Microbiol Antimicrob. 2006;5:26.

PubMed   |  Link to Article

16

Bisdorff  B, Scholhölter  JL, Claußen  K, Pulz  M, Nowak  D, Radon  K.  MRSA-ST398 in livestock farmers and neighbouring residents in a rural area in Germany. Epidemiol Infect. 2012;140(10):1800-1808.

PubMed   |  Link to Article

17

van Cleef  BA, Verkade  EJ, Wulf  MW,  et al.  Prevalence of livestock-associated MRSA in communities with high pig-densities in the Netherlands. PLoS One. 2010;5(2):e9385.

PubMed   |  Link to Article

18

Price  LB, Stegger  M, Hasman  H,  et al Staphylococcus aureus CC398: host adaptation and emergence of methicillin resistance in livestock. MBio.2012;3(1):e00305-e00311.

Link to Article

19

Crum  NF, Lee  RU, Thornton  SA,  et al.  Fifteen-year study of the changing epidemiology of methicillin-resistant Staphylococcus aureus. Am J Med. 2006;119(11):943-951.

PubMed   |  Link to Article

20

Liu  C, Graber  CJ, Karr  M,  et al.  A population-based study of the incidence and molecular epidemiology of methicillin-resistant Staphylococcus aureus disease in San Francisco, 2004-2005. Clin Infect Dis. 2008;46(11):1637-1646.

PubMed   |  Link to Article

21

Tracy  LA, Furuno  JP, Harris  AD, Singer  M, Langenberg  P, Roghmann  MC.  Staphylococcus aureus infections in US veterans, Maryland, USA, 1999-2008. Emerg Infect Dis. 2011;17(3):441-448.

PubMed   |  Link to Article

22

Casey  JA, Cosgrove  SE, Stewart  WF, Pollak  J, Schwartz  BS.  A population-based study of the epidemiology and clinical features of methicillin-resistant Staphylococcus aureus infection in Pennsylvania, 2001-2010. Epidemiol Infect. 2013;141(6):1166-1179.

PubMed   |  Link to Article

23

Fridkin  SK, Hageman  JC, Morrison  M,  et al; Active Bacterial Core Surveillance Program of the Emerging Infections Program Network.  Methicillin-resistant Staphylococcus aureus disease in three communities. N Engl J Med. 2005;352(14):1436-1444.

PubMed   |  Link to Article

24

David  MZ, Daum  RS.  Community-associated methicillin-resistant Staphylococcus aureus: epidemiology and clinical consequences of an emerging epidemic. Clin Microbiol Rev. 2010;23(3):616-687.

PubMed   |  Link to Article

25

Landrum  ML, Neumann  C, Cook  C,  et al.  Epidemiology of Staphylococcus aureus blood and skin and soft tissue infections in the US military health system, 2005-2010. JAMA. 2012;308(1):50-59.

PubMed   |  Link to Article

26

Klevens  RM, Morrison  MA, Fridkin  SK,  et al; Active Bacterial Core Surveillance of the Emerging Infections Program Network.  Community-associated methicillin-resistant Staphylococcus aureus and healthcare risk factors. Emerg Infect Dis. 2006;12(12):1991-1993.

PubMed   |  Link to Article

27

Maree  CL, Daum  RS, Boyle-Vavra  S, Matayoshi  K, Miller  LG.  Community-associated methicillin-resistant Staphylococcus aureus isolates causing healthcare-associated infections. Emerg Infect Dis. 2007;13(2):236-242.

PubMed   |  Link to Article

28

Tattevin  P, Schwartz  BS, Graber  CJ,  et al.  Concurrent epidemics of skin and soft tissue infection and bloodstream infection due to community-associated methicillin-resistant Staphylococcus aureus. Clin Infect Dis. 2012;55(6):781-788.

PubMed   |  Link to Article

29

Fitzgerald  JR.  Livestock-associated Staphylococcus aureus: origin, evolution and public health threat. Trends Microbiol. 2012;20(4):192-198.

PubMed   |  Link to Article

30

Liu  AY, Curriero  FC, Glass  TA, Stewart  WF, Schwartz  BS.  The contextual influence of coal abandoned mine lands in communities and type 2 diabetes in Pennsylvania. Health Place. 2013;22:115-122.

PubMed   |  Link to Article

31

Morrison  MA, Hageman  JC, Klevens  RM.  Case definition for community-associated methicillin-resistant Staphylococcus aureus. J Hosp Infect. 2006;62(2):241.

PubMed   |  Link to Article

32

Klevens  RM, Morrison  MA, Nadle  J,  et al; Active Bacterial Core surveillance (ABCs) MRSA Investigators.  Invasive methicillin-resistant Staphylococcus aureus infections in the United States. JAMA. 2007;298(15):1763-1771.

PubMed   |  Link to Article

33

Schwartz  BS, Stewart  WF, Godby  S,  et al.  Body mass index and the built and social environments in children and adolescents using electronic health records. Am J Prev Med. 2011;41(4):e17-e28.

PubMed   |  Link to Article

34

Talen  E, Anselin  L.  Assessing spatial equity: an evaluation of measures of accessibility to public playgrounds. Environ Plan A. 1998;30(4):595-613.

Link to Article

35

Mathema  B, Mediavilla  J, Kreiswirth  BN. Sequence analysis of the variable number tandem repeat in Staphylococcus aureus protein A gene: spa typing. In: DeLeo  F, Otto  M, eds. Methods in Molecular Biology: Bacterial Pathogenesis.Vol 431. Totowa, NJ: Humana Press; 2008:285-305.

36

Harmsen  D, Claus  H, Witte  W,  et al.  Typing of methicillin-resistant Staphylococcus aureus in a university hospital setting by using novel software for spa repeat determination and database management. J Clin Microbiol. 2003;41(12):5442-5448.

PubMed   |  Link to Article

37

Lina  G, Piémont  Y, Godail-Gamot  F,  et al.  Involvement of Panton-Valentine leukocidin-producing Staphylococcus aureus in primary skin infections and pneumonia. Clin Infect Dis. 1999;29(5):1128-1132.

PubMed   |  Link to Article

38

Grineski  SE, Staniswalis  JG, Bulathsinhala  P, Peng  Y, Gill  TE.  Hospital admissions for asthma and acute bronchitis in El Paso, Texas: do age, sex, and insurance status modify the effects of dust and low wind events? Environ Res. 2011;111(8):1148-1155.

PubMed   |  Link to Article

39

Ayanian  JZ, Kohler  BA, Abe  T, Epstein  AM.  The relation between health insurance coverage and clinical outcomes among women with breast cancer. N Engl J Med. 1993;329(5):326-331.

PubMed   |  Link to Article

40

Waller  L, Gotway  C. Applied Spatial Statistics for Public Health Data. Hoboken, NJ: John Wiley & Sons Inc; 2004.

41

Rockhill  B, Newman  B, Weinberg  C.  Use and misuse of population attributable fractions. Am J Public Health. 1998;88(1):15-19.

PubMed   |  Link to Article

42

Pallin  DJ, Espinola  JA, Leung  DY, Hooper  DC, Camargo  CA  Jr.  Epidemiology of dermatitis and skin infections in United States physicians’ offices, 1993-2005. Clin Infect Dis. 2009;49(6):901-907.

PubMed   |  Link to Article

43

Yang  ES, Tan  J, Eells  S, Rieg  G, Tagudar  G, Miller  LG.  Body site colonization in patients with community-associated methicillin-resistant Staphylococcus aureus and other types of S. aureus skin infections. Clin Microbiol Infect. 2010;16(5):425-431.

PubMed   |  Link to Article

44

van Cleef  BA, Graveland  H, Haenen  AP,  et al.  Persistence of livestock-associated methicillin-resistant Staphylococcus aureus in field workers after short-term occupational exposure to pigs and veal calves. J Clin Microbiol. 2011;49(3):1030-1033.

PubMed   |  Link to Article

45

Feingold  BJ, Silbergeld  EK, Curriero  FC, van Cleef  BA, Heck  ME, Kluytmans  JA.  Livestock density as risk factor for livestock-associated methicillin-resistant Staphylococcus aureus, the Netherlands. Emerg Infect Dis. 2012;18(11):1841-1849.

PubMed   |  Link to Article

46

Cuny  C, Nathaus  R, Layer  F, Strommenger  B, Altmann  D, Witte  W.  Nasal colonization of humans with methicillin-resistant Staphylococcus aureus (MRSA) CC398 with and without exposure to pigs. PLoS One. 2009;4(8):e6800.

PubMed   |  Link to Article

47

Köck  R, Harlizius  J, Bressan  N,  et al.  Prevalence and molecular characteristics of methicillin-resistant Staphylococcus aureus (MRSA) among pigs on German farms and import of livestock-related MRSA into hospitals. Eur J Clin Microbiol Infect Dis. 2009;28(11):1375-1382.

PubMed   |  Link to Article

48

Neela  V, Mohd Zafrul  A, Mariana  NS, van Belkum  A, Liew  YK, Rad  EG.  Prevalence of ST9 methicillin-resistant Staphylococcus aureus among pigs and pig handlers in Malaysia. J Clin Microbiol. 2009;47(12):4138-4140.

PubMed   |  Link to Article

49

Gibbs  SG, Green  CF, Tarwater  PM, Mota  LC, Mena  KD, Scarpino  PV.  Isolation of antibiotic-resistant bacteria from the air plume downwind of a swine confined or concentrated animal feeding operation. Environ Health Perspect. 2006;114(7):1032-1037.

PubMed   |  Link to Article

50

Schulz  J, Friese  A, Klees  S,  et al.  Longitudinal study of the contamination of air and of soil surfaces in the vicinity of pig barns by livestock-associated methicillin-resistant Staphylococcus aureus. Appl Environ Microbiol. 2012;78(16):5666-5671.

PubMed   |  Link to Article

51

Hamscher  G, Sczesny  S, Höper  H, Nau  H.  Determination of persistent tetracycline residues in soil fertilized with liquid manure by high-performance liquid chromatography with electrospray ionization tandem mass spectrometry. Anal Chem. 2002;74(7):1509-1518.

PubMed   |  Link to Article

52

Geenen  PL, Graat  EA, Haenen  A,  et al.  Prevalence of livestock-associated MRSA on Dutch broiler farms and in people living and/or working on these farms. Epidemiol Infect. 2013;141(5):1099-1108.

PubMed   |  Link to Article

53

Molla  B, Byrne  M, Abley  M,  et al.  Epidemiology and genotypic characteristics of methicillin-resistant Staphylococcus aureus strains of porcine origin. J Clin Microbiol. 2012;50(11):3687-3693.

PubMed   |  Link to Article

54

Waters  AE, Contente-Cuomo  T, Buchhagen  J,  et al.  Multidrug-resistant Staphylococcus aureus in US meat and poultry. Clin Infect Dis. 2011;52(10):1227-1230.

PubMed   |  Link to Article

55

Khanna  T, Friendship  R, Dewey  C, Weese  JS.  Methicillin resistant Staphylococcus aureus colonization in pigs and pig farmers. Vet Microbiol. 2008;128(3-4):298-303.

PubMed   |  Link to Article

56

Wulf  MW, Sørum  M, van Nes  A,  et al.  Prevalence of methicillin-resistant Staphylococcus aureus among veterinarians: an international study. Clin Microbiol Infect. 2008;14(1):29-34.

PubMed   |  Link to Article

57

Wulf  MW, Verduin  CM, van Nes  A, Huijsdens  X, Voss  A.  Infection and colonization with methicillin resistant Staphylococcus aureus ST398 versus other MRSA in an area with a high density of pig farms. Eur J Clin Microbiol Infect Dis. 2012;31(1):61-65.

PubMed   |  Link to Article

58

Salmenlinna  S, Lyytikäinen  O, Vainio  A,  et al.  Human cases of methicillin-resistant Staphylococcus aureus CC398, Finland. Emerg Infect Dis. 2010;16(10):1626-1629.

PubMed   |  Link to Article

59

Golding  GR, Bryden  L, Levett  PN,  et al.  Livestock-associated methicillin-resistant Staphylococcus aureus sequence type 398 in humans, Canada. Emerg Infect Dis. 2010;16(4):587-594.

PubMed   |  Link to Article

60

Witte  W, Strommenger  B, Stanek  C, Cuny  C.  Methicillin-resistant Staphylococcus aureus ST398 in humans and animals, Central Europe. Emerg Infect Dis. 2007;13(2):255-258.

PubMed   |  Link to Article

61

Wiley  LK, Shah  A, Xu  H, Bush  WS.  ICD-9 tobacco use codes are effective identifiers of smoking status. J Am Med Inform Assoc. 2013;20(4):652-658.

PubMed   |  Link to Article

62

Bagger  JP, Zindrou  D, Taylor  KM.  Postoperative infection with methicillin-resistant Staphylococcus aureus and socioeconomic background. Lancet. 2004;363(9410):706-708.

PubMed   |  Link to Article

63

Bratu  S, Landman  D, Gupta  J, Trehan  M, Panwar  M, Quale  J.  A population-based study examining the emergence of community-associated methicillin-resistant Staphylococcus aureus USA300 in New York City. Ann Clin Microbiol Antimicrob. 2006;5:29.

PubMed   |  Link to Article