Glyphosate Use Predicts ADHD Hospital Discharges in the Healthcare Cost and Utilization Project Net (HCUPnet): A Two-Way Fixed-Effects Analysis

There has been considerable international study on the etiology of rising mental disorders, such as attention-deficit hyperactivity disorder (ADHD), in human populations. As glyphosate is the most commonly used herbicide in the world, we sought to test the hypothesis that glyphosate use in agriculture may be a contributing environmental factor to the rise of ADHD in human populations. State estimates for glyphosate use and nitrogen fertilizer use were obtained from the U.S. Geological Survey (USGS). We queried the Healthcare Cost and Utilization Project net (HCUPNET) for state-level hospitalization discharge data in all patients for all-listed ADHD from 2007 to 2010. We used rural-urban continuum codes from the USDA-Economic Research Service when exploring the effect of urbanization on the relationship between herbicide use and ADHD. Least squares dummy variable (LSDV) method and within method using two-way fixed effects was used to elucidate the relationship between glyphosate use and all-listed ADHD hospital discharges. We show that a one kilogram increase in glyphosate use, in particular, in one year significantly positively predicts state-level all-listed ADHD discharges, expressed as a percent of total mental disorders, the following year (coefficient = 5.54E-08, p<.01). A study on the effect of urbanization on the relationship between glyphosate and ADHD indicates that the relationship is marginally significantly positive after multiple comparison correction only in urban U.S. counties (p<.025). Furthermore, total glyphosate use is strongly positively associated with total farm use of nitrogen fertilizers from 1992 to 2006 (p<.001). We present evidence from the biomedical research literature of a plausible link among glyphosate, nitrogen dysbiosis and ADHD. Glyphosate use is a significant predictor of state hospitalizations for all-listed ADHD hospital discharges, with the effect concentrated in urban U.S. counties. This effect is seen even after controlling for individual state characteristics, strong correlations over time, and other significant associations with ADHD in the literature. We draw upon the econometric results to propose unique mechanisms, borrowing principles from soil and atmospheric PLOS ONE | DOI:10.1371/journal.pone.0133525 August 19, 2015 1 / 31 OPEN ACCESS Citation: Fluegge KR, Fluegge KR (2015) Glyphosate Use Predicts ADHD Hospital Discharges in the Healthcare Cost and Utilization Project Net (HCUPnet): A Two-Way Fixed-Effects Analysis. PLoS ONE 10(8): e0133525. doi:10.1371/journal. pone.0133525 Editor: Valsamma Eapen, University of New South Wales, AUSTRALIA Received: January 20, 2015 Accepted: June 29, 2015 Published: August 19, 2015 Copyright: © 2015 Fluegge, Fluegge. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability Statement: Glyphosate data is available from the United States Geological Survey (USGS). All mental health discharge data are available at (HCUPnet). The links to access these public data are contained within the manuscript. Funding: This research was supported by the NIH National Heart, Lung, and Blood Institute (Grants T32HL007567 and HL007567-29 (T32)). Competing Interests: The authors have declared that no competing interests exist. sciences, for how glyphosate-based herbicides may be contributing to the rise of ADHD in all populations.

Introduction cascade that could eventually adversely impact dopaminergic activity [18], and dopamine has been implicated in attention-related neural networking [19]. In addition to Bt protein, glyphosate could also be a contributing factor to the increasing prevalence of ADHD.
Glyphosate has been shown to disrupt cytochrome P450 (CYP) enzymes [20], which, among many other downstream effects, can inhibit activation of vitamin D3, which depends on CYP enzymes in both the liver and the kidney [21]. Vitamin D regulates serotonin synthesis in the brain [22], and reduced central nervous system serotonergic activity has been implicated as a factor in ADHD [23].
Our preliminary investigations led us to consider a complicit role of nitrogen dysbiosis in the disease process in addition to the more direct mechanisms presented previously. We became aware that glyphosate usage on crops, such as maize, necessitates an increased application rate of nitrogen fertilizers, because glyphosate disrupts the uptake of nitrogen by plants [24]. Even in glyphosate resistant soybean crops, glyphosate application has been shown to disrupt nitrogen fixation [25]. We also propose that glyphosate disrupts the soil bacteria, leading to a major change in the way nitrogen is handled in the soil. This can cause an accumulation of various toxic nitrogen oxides in rain water exposed to soil, that then causes ammonia and nitrogen oxides in the soil to seep into the major waterways. The action of microbes in the soil and water can cause release of nitrous oxide (N 2 O) into the air, or the accumulation of nitrates in the water. Both of these are known to have toxic effects on human physiology: the dopaminergic system was shown to mediate the antinociceptive effects of N 2 O [26], and levels of central neurotransmitters, dopamine (DA) and norepinephrine (NE), were significantly elevated after repeated N 2 O exposure in CD-1 mice [27]. Dysregulation of both catecholamine systems has been implicated in ADHD [28]. Furthermore, BTBR mice-an inbred mouse strain that exhibits many behaviors typical of autism [29]-have severe heparan sulfate deficiency in the cerebrospinal fluid [30], a shared feature with autism [31], and the reduction of sulfate by Desulfovibrio may be a contributing factor to sulfate insufficiency [32].
Our analysis has focused on state-provided data on ADHD hospital discharge rates and corresponding information about glyphosate usage, nitrogen fertilizer usage, industrial and food usage of maize and associated N 2 O emissions. We use herbicide resistant weed event statistics for further clarification on glyphosate use. Our objectives are shown as an empirical schematic in Fig 1. We have been somewhat limited by the lack of available state-level hospital discharge data on ADHD prior to 2007. However, we have been able to make use of data on other mood and behavioral disorders to help inform our analysis.
We have uncovered a pattern that is challenging to explain: glyphosate use significantly positively predicts glyphosate resistance events in weeds prior to 2007, but not subsequently. We suspect that a modification in the formulation and the use of maize as a residual cover crop in agriculture caused a widespread suppression of weed resistance development, while at the same time causing an increase in the rate of ADHD in the population. We expect that the two are related, and we propose hypotheses to explain this relationship. We have also investigated the relationship between ADHD and population density, by dividing data into metropolitan, urban, and rural subsets.

Herbicide Use
Pesticide use was determined using data from the United States Geological Survey [33]. The USGS uses the estimated pesticide use rate (known as EPest rate) to determine pesticide usage rates for 39 pesticides among 76 crops and in 304 crop reporting districts (CRD) in the United States from 1992-2009. Pesticide usage was defined as pounds applied per harvested-crop acre in the CRD. The proprietary pesticide use survey data for each CRD was obtained from GfK Kynetec, Inc. These data were normalized to the harvested acreage in nearby counties to arrive at county estimates of pesticide use rate. Harvested acreage is non-specific and includes the crops themselves, the soil on which the crops were grown, as well as the air above the acreage. In our analysis, we have summed EPest high county estimates (we also use low county estimates for purposes of replication) to generate state-level pesticide use for each year 2006-2009, a practice recommended by the USGS [34]. The difference between high and low county estimates is that high county estimates include more counties. We have also included use of other herbicides-arbitrarily selecting atrazine, halosulfuron, and clethodim-as referent herbicides from the most prevalent groups of herbicide resistance, including Photosystem II inhibitors, acetolactate synthase (Als) inhibitors, and acetyl-Coenzyme A carboxylase (ACCase) inhibitors, respectively.
In a subset analysis, we have assigned an urbanization score to each of the counties included in the USGS glyphosate usage rate. We have used the 1993 and 2003 urban-rural continuum codes from the USDA Economic Research Services [35]. All counties receiving a code of 0, 1, 2, or 3 were deemed metropolitan counties, and we assigned an urbanization score of 1 (metropolitan) or 2 (fringe metropolitan) to these counties. Counties receiving a USDA continuum code of 4, 5, 6, or 7 were given an urbanization score of 3 (urban). Counties receiving a USDA continuum code of 8 or 9 were given an urbanization score of 4 (rural). We have used the 1993 year codes for all analysis in the years between 1992 and 2002 and the 2003 codes for all analysis in 2003 and later.

Nitrogen Inputs
We used county-level estimates of nitrogen fertilizer use to compile state statistics for farm use of nitrogen fertilizer between the years 1992 and 2006. The selection protocol can be found in [36]. Briefly, the authors gathered annual fertilizer sales data from the Association of American Plant Food Control Officials. This state data was used to determine total farm and non-farm proportions of nitrogen fertilizer use. The farm nitrogen tonnage was then estimated by multiplying the farm nitrogen rate by the nitrogen percentage from N-P-K values and the total nitrogen product tonnage.  Since the county continuum codes changed from 1993 to 2003, we have standardized the codes so that a comparison among counties could be made. To that effort, all counties with a continuum code 0 between 1992 and 2006 were given a code 1. All other codes remained the same.
State-level estimates of total food and industrial use of maize was gathered from data published by the University of Nebraska-Lincoln [37]. The data was derived from the Feed Grain Database: Yearbook Tables, Corn: Food, Seed, and Industrial Use, Table 31 from the Economic Research Service at the USDA [38]. State data were expressed as a percent of the total food and industrial use of maize grain in the U.S.

Nitrous Oxide (N 2 O) Emissions
Agricultural N 2 O emissions were obtained from the U.S. Energy Information Administration [39]. Emissions are expressed as million metric tons carbon dioxide equivalent. The EIA classifies emissions according to the following sources: agricultural, energy, industry, and waste management. To account for all sources, we expressed agricultural emissions as a percentage of total N 2 O emissions.

Herbicide Resistance Events
Statewide data for herbicide-resistant weeds were obtained from the International Survey of Herbicide Resistant Weeds (ISHRW) in August 2014 and replicated in April 2015 [40]. We used four herbicide resistant categories, according to the sites of action resisted. To measure the effect of weeds resistant to EPSP synthase inhibitors (which account for 9.6%, or 14 out of 146, of total herbicide tolerant weeds in the United States at the time of this publication), we chose to include in our model the herbicides that accounted for at least ten percent of the total cumulative resistant weeds in the United States at the time of this publication as reported by ISHRW. Therefore, we included the following herbicide groups, showing the percent of the total resistant weeds in the United States in parentheses: ALS inhibitors (31.5%, or 46 out of 146), Photosystem II inhibitors (17.8%, or 26 out of 146), and ACCase inhibitors (10.3%, or 15 out of 146). However, given that one weed can possess multiple sites of resisted action, we used the number of resisted sites of action for each of these groups in each state. To account for groups with less than ten percent of the total cumulative resistant weeds, we expressed data for herbicide weed resistance in each group of interest as a percentage of the total cumulative weed resistance events in each state for the period 2000-2012. Data before 2000 was incorporated into the cumulative analysis.

Hospitalization Data
The Healthcare Cost and Utilization Project (HCUPNET) State Inpatient Database was queried to identify state-level hospital discharge data trends for all-listed diagnoses of mental disorders [41]. The number of states reporting such data to HCUPNET was 37, although not all states reporting provided data for every year. The time period of interest was 2007 to 2010. Records for attention-deficit, conduct, and disruptive behavior disorders in HCUPNET were begun in 2007. For the all-listed diagnostic category, we searched for the number of discharges that received a diagnosis of attention-deficit, conduct, and disruptive behavior disorder (ADHD). HCUPNET defines all-listed diagnoses as being ". . .the principal diagnosis plus additional conditions that coexist at the time of admission, or that develop during the stay, and which have an effect on the treatment or length of stay in the hospital." All diagnoses that patients receive while admitted to the hospital are assigned to an ICD-9 (International Classification of Disease, 9 th Revision-Clinical Modification) code. The ICD-9 codes that are included in the clinical classification software (CCS) diagnosis of attention-deficit, conduct, and disruptive behavior disorder are shown in S1 Table, along with a translation to respective codes in the  Diagnostic and Statistical Manual for Mental Disorders IV, Text Revision. Because we have previously hypothesized that glyphosate may exert a deleterious effect on many mental conditions, we have normalized the number of all-listed ADHD hospital discharges to a percentage of total discharges for all-listed mental disorders recorded in HCUPNET.
Hospitalization discharge data for all-listed ADHD diagnoses were also categorized according to the location of the patients' residence, as a percentage of the total all-listed ADHD hospital discharges. HCUPNET provided four categories in this respect: large central metro, large fringe metro (suburbs), medium and small metro, micropolitan and noncore (rural). We have created four categories, called metropolitan, fringe metropolitan, urban, and rural, to mirror the glyphosate and nitrogen fertilizer usage urbanization coding system mentioned previously. Our metropolitan category included the large central metro heading in HCUPNET. Fringe metro was the second category. Medium/small metro defined our urban category, and our rural category included the micropolitan and noncore (rural) data in HCUPNET. Hospitalization data for each of these categories were expressed in HCUPNET as a percentage of the total hospital discharges for all-listed ADHD diagnoses in each state, and we retained this unit of measure. We also gathered data on age and gender of patients with an all-listed ADHD discharge.
To verify the discharge data among the reporting states in HCUPNET, we have also gathered CDC data on the number of hospital discharges for all-listed ADHD (ICD-9-CM code: 314.01) between 1993 to 2005 [42]. We compared these data to prior year's estimated glyphosate using ordinary least squares regression.

Least Squares Dummy Variable (LSDV) Regression
Panel regressions were performed using the LSDV in plm package in R, version 3.1.3 [43]. Two-ways within estimation was also used to confirm the LSDV estimations and all model diagnostics and was used solely when time (T) for each subject was greater than 4 years. The LSDV method creates dummy variables for both state and time, while the two-ways within effect estimations use deviations from group means. When noted, simple ordinary least squares regression modeling was performed in R using the fit argument. Annual county input data from the USGS was summed and organized into study categories by Microsoft Access 2010. Fixed effects data were checked for heteroskedasticity in R using the Breusch-Pagan (BP) test (package lmtest) [44]. Heteroskedasticity is the phenomenon that a dependent variable (such as hospitalization discharges) shows variability during the time period in which the independent variable (such as herbicide use) predicts its occurrence. Data were also tested for serial correlation using the Durbin-Watson (DW) statistic [43,44]. We chose the Fisher-type test in Stata 11.1 [45] to check for unit roots in unbalanced panels for each model since the test uses Newey-West standard errors to account for any serial correlation present among the residuals in each model. To our knowledge, an equivalent test was not available in R. The Augmented Dickey-Fuller (ADF) test (package tseries) [46] was used to test for unit roots in Model 3 panel sets. We used the Arellano variance covariance (vcov) HC function covariance estimator [47] to control for heteroskedasticity and serial correlation in all models. We used student's paired t-test when making comparisons between age and urban subsets of ADHD hospital discharge among available HCUPNET states during our study period (2007)(2008)(2009)(2010)). Data were first tested for normality using Shapiro-Wilk test. Unless otherwise specified, we define significance at the .05 level for all statistical testing performed, while a .10 level was used for model diagnostic testing (i.e., BP, DW, ADF, Fisher test). Graphical output was performed using R and GraphPad Prism 6 Demo for Windows, GraphPad Software, La Jolla California USA, www.graphpad. com.
Model 1. Our first model tested our hypothesis that glyphosate use in one year predicts all-listed ADHD hospital discharges the following year. In addition to glyphosate and maize use, we have identified four time-variant covariates, or regressors, including other societal and health care related variables that have been associated with the disorder and could explain the increase in ADHD diagnoses. Furthermore, we have included glyphosate resistance to clarify the herbicide's effect on human health. We subjected each covariate, excluding our variables of interest, glyphosate use and food/industrial maize use, to a single linear ordinary least squares model for each year from 2007 to 2010, inclusive. If a significant relationship was noted in any year at the .05 alpha level, we included the covariate in the final fixed effects model. Model 1 is expressed below where Y it is the proportion of all mental disorders for which ADHD was a discharging condition in state i and year t, c 1 is a constant, α i represents a state fixed effect to control for permanent differences between states and v t denotes a time fixed effect, X m i;tÀ1 is a (N 1 x 1) column vector representing an individual lagged covariate m, with m = 1,. . ., M, and N 1 is the total number of state-time observations, X p i;t is a (N 1 x 1) column vector representing an individual contemporaneous covariate p, with p = 1,. . ., P, and ε it captures unobserved heterogeneity in state i and year t. The parameters to be estimated include b m i;tÀ1 and b p i;t . The states of Nevada and Hawaii were excluded in this model since there were no weed resistance or glyphosate use data, respectively, for the years of interest. All covariates in model 1 are explained in Table 1.
Model 2. Our second model-an ad hoc analysis-sought to test the hypothesis that, if glyphosate is significantly contributing to ADHD, there would be regional influences underpinning this relationship. That is to say that glyphosate is a main input in the agricultural industry, which largely has existed in rural (i.e., less populated) parts of the country, although inroads have been made in sustaining and expanding adaptive metropolitan agriculture systems [56]. To study this, we have categorized all counties from the USGS glyphosate use estimates into four codes (1 = metropolitan, 2 = fringe metropolitan, 3 = urban, 4 = rural) to generate four data points for each HCUPNET state in each year. The state of Hawaii was excluded since no estimates of glyphosate use were available. These data were then matched to the ADHD data in HCUPNET, according to the model specifications in Table 2. The other three predictors were unchanged from Model 1. Hospitalization data by location of patient residence and patient age was masked ( Ã ) in HCUPNET if less than or equal to ten discharges or fewer than two hospitals were reported, and these data were removed from the analysis. Beginning in 2007, HCUPNET began revising data on location of patients' residence to accommodate the specifications provided by the Economic Research Service at the USDA. Our four models were therefore where each model corresponds to one of four unique k values, Y k it is the proportion of total alllisted ADHD disorders in state i and year t that occurred among patients residing in a particular county classification code (k = 1, 2, 3, 4), c 2 , α i and v t are as defined in the previous sub-section for each model k, X it is a contemporaneous continuous variable (i.e., population percent, hence the p superscript equals 1 and is therefore omitted for simplicity) in state i and year t as in model 1, X m i;tÀ1 denotes one of two lagged variables predicting ADHD discharges and ε k it captures unobserved heterogeneity in state i and year t for model k. Our main predictor of interest, glyphosate use (G k i;tÀ1 ), is expressed in kilograms as the previous year's total use in counties categorized only by urbanization code k in state i. The sixteen parameters to be estimated include g k it , g 1;k i;tÀ1 , g 2;k i;tÀ1 and b k i;tÀ1 for k = 1, 2, 3, and 4. Model 3. Third, we tested our hypothesis that glyphosate use could be indirectly altering the land biota to impact the rise of ADHD. We assessed the influence of county level estimates of glyphosate application on county level estimates of farm nitrogen fertilizer use between 1992 and 2006. The hypothesis is that increasing amounts of glyphosate use would perturb soil such that greater farm use of nitrogen fertilizer would be needed. We summed glyphosate and farm nitrogen use for all USDA county continuum codes between 1992 and 2006. We then organized our panel sets by study categories set forth in Table 2. We focused our attention to glyphosate use in urban counties given our results with the Model 2 set. For comparison, we also analyzed results for fringe metro counties. Our two models were where each model corresponds to one of two unique k values, Z k j;t is a (N 3 x 1) column vector (n is the number of RUCC, and N 3 is represented by (n Ã t)) and denotes the farm use of nitrogenbased fertilizers (in kilograms) in the RUCC county code j with classification category k in year t, c k 3 is a constant that varies by model, a k j is a RUCC fixed effect, and v k t denotes a time fixed effect, X k j;t indicates the glyphosate use estimate (in kilograms) in the RUCC county code j with category k in year t, and k j;t captures unobserved heterogeneity in county code j in year t. The two parameters to be estimated and reported include b k j;t , where k = 2 and 3. These estimates generate a comparison of the effect of glyphosate's use in fringe metro and urban counties on application of nitrogen-based fertilizers.  Data on estimated glyphosate use by USDA rural-urban continuum code between 1992 and 2006 indicate that the greatest relative increase in use occurred in county codes 4 and 5 (1523 and 1240%, respectively). The largest relative increase in estimated glyphosate use during the period of study (2006)(2007)(2008)(2009) occurred in county code 9 (55%).
USGS county sums of farm use of nitrogen fertilizer show that the greatest relative percentage increase occurred in those counties with a USDA rural-urban continuum code of either 3 or 4 (64 and 49%, respectively).
Total U.S. food and industrial use of maize grain increased by 2 billion bushels, or roughly

Agricultural N 2 O Emissions
According to the United States Environmental Protection Agency, nitrous oxide (N 2 O) accounted for 6% of all U.S. greenhouse gas emissions. Man-made sources of N 2 O include agriculture, energy, industry, and waste management practices, with agricultural soil management accounting for 75% of U.S. N 2 O emissions. Because of its long half-life, N 2 O is thought to be a source with a global warming potential over 300 times that of other greenhouse gases, like carbon dioxide [57]. Total N 2 O agricultural emissions have been increasing since 1990. Total agricultural emissions have increased, as a percentage of total N 2 O emissions, 9.20% or about 12.3 million metric tons from the period 1990-2009 [39].
Direct agricultural soil N 2 O emissions, including those from the use of synthetic nitrogen fertilizers, animal manure, sewage sludge and other non-manure organics, above and below ground crop residues, and soil mineralization, increased 9.1% as a percentage of total N 2 O emissions, or about 9 million metric tons carbon dioxide equivalent, from 1990-2009. Almost half of that normalized increase (48.77%) came during the period of the current study (2006)(2007)(2008)(2009). On an absolute basis, direct agricultural nitrous oxide emissions have increased 8.23% between 1990 and 2009, with almost a third (30.50%) of that increase occurring between 2006 and 2009. Most of that increase between 2006 and 2009 was attributable to above and below ground crop residues. Use of high residue crops in soil management is a growing trend in agriculture [58].
Indirect agricultural soil N 2 O emissions, including those from soil leaching, runoff and atmospheric deposition, increased 11.28% as a percentage of total N 2 O emissions, or about 2 million metric tons carbon dioxide equivalent, between 1990 and 2009. Most of this increase is attributable to soil leaching and runoff. It should also be noted that significant N 2 O emissions may be coming from indirect sources such as bodies of freshwater [59]. Although formal estimations of N 2 O emissions from bodies of water are limited, there is growing evidence to suggest that these comprise a significant source of emissions and must be included in future estimates, especially considering the increasing agricultural runoff into outlying drainage channels [60].

Herbicide-Resistance Events
Since it is the policy of the ISHRW that data derived from their database shall not be used in a significant way in a scientific paper without the expressed permission of the director, we have chosen to include the results of this analysis as an unpublished supplement to this paper that can be viewed at the non-profit's website [61]. Our data leads us to believe that the development of herbicide weed resistance is not a significant component in our hypothesized relationship between glyphosate use and all-listed ADHD hospital discharges.
We can validate, using two-way fixed effects, that use of herbicides with contact modalities of action does significantly predict class resistance between 2000 and 2010. Glyphosate use in states (N = 470), as a lag indicator, does significantly positively predict EPSP synthase inhibitor weed resistance (adjusted estimate: 3.05e-08±1.14e-08, p-value < .01). That is, for a one kilogram increase in glyphosate used, there is a 3.05e-08% increase in glyphosate weed resistance, as a percentage of total weed resistance events in each state, the following year. Similarly, halosulfuron use in states (N = 470) in one year significantly positively predicts acetolactate synthase (ALS) inhibitor weed resistance, as a percentage of total resistance events, the following year (adjusted estimate: 2.09e-05±8.42e-06, p-value: 01). Atrazine use in states (N = 470) in one year also does significantly positively predict photosystem inhibitor weed resistance the following year as a percentage of all resistance events (adjusted estimate: 3.62e-08 ±1.75e-08, p-value: 0.04).
We, therefore, feel confident that, if use predicts resistance, the phenomenon of weed resistance can be a possible mechanism (through consumption of treated vegetation) to explain the hypothesized health effects from the use of similarly purposed herbicides (i.e., same site of targeted action). However, in the current study period, glyphosate use does not significantly positively predict EPSP synthase inhibitor weed resistance in HCUPNET states (adjusted estimate: 7.55e-09±1.24e-08, p-value: 0.54) (S1 Fig). Similarly, halosulfuron use does not significantly positively predict acetolactate synthase (Als) inhibitor weed resistance (adjusted estimate: 6.83e-06±7.27e-06, p-value: 0.35). Atrazine does not significantly positively predict photosystem II weed resistance events either (adjusted estimate: 1.20e-08±3.44e-08, p-value: 0.73). Therefore, the development of herbicide resistance (i.e., herbicide application on vegetation) among contact herbicides, like glyphosate, is insignificant and unrelated to any link between use and disorders of interest between 2007 and 2010, and this data is confirmed by other reports [62]. We find distinct trends emerging from a breakdown of all-listed ADHD discharges by patient age in HCUPNET reporting states from 2007 to 2010. Persons aged 1-17 years accounted for an average of 48.8% of the all-listed ADHD discharges in 2007 but that percentage decreased to 40.9% by 2010, a statistically significant drop (p-value = 0.01, 95% CI: 1.67-14.01), while adults aged 18-44 gained as a percentage of total all-listed ADHD discharges, climbing from 32.7% to 38.0% across all HCUPNET reporting states (p-value < 0.05, 95% CI: -9.90--0.74) (Fig 4).

Hospitalization Data
The state-level ADHD data by patient age comports with national hospital discharge statistics in HCUPNET. There were 174,746±11,765 discharges with a diagnosis of ADHD in Hospital ADHD discharge diagnoses by gender revealed that majority were males, although there was a significant drop in the percent (58.6% to 56.4%, p-value < 0.01, 95% CI: 0.81-3.50) of discharges among males and a significant rise in the percent (41.4% to 43.6%, p-value < 0.01, 95% CI: -3.51--0.82) of discharges among females during the study period. This data matches the national HCUPNET statistics on gender and ADHD, with 58.65% of all-listed ADHD cases being male and 41.22% being female in 2007. In 2010, 56.29% of all-listed ADHD discharges diagnoses were male, while 43.55% were female.
Urbanization data among both HCUPNET reporting states and national HCUPNET statistics indicate that the greatest percentage of all-listed ADHD diagnoses occur among persons living in medium and small metro (urban populations) areas. The average percentage of alllisted ADHD hospital discharges occurring in persons living in medium and/or small metro areas (urban) across all HCUPNET reporting states used in this analysis, excluding Hawaii, was significantly higher than other urbanization categories (Fig 4). This finding matches the national HCUPNET discharge statistics, with 32.47% (56,736±7,774 discharges) of total alllisted ADHD hospital discharges occurring in persons residing in urban counties in 2007. In 2010, there were 86,081±9,608 discharges with an all-listed ADHD diagnosis among persons living in urban counties, comprising 33.35% of the national total.
The national all-listed ADHD hospital discharge data from the CDC reaffirms the normalized increase seen in HCUPNET reporting states. The all-listed ADHD hospital discharges (ICD-9-CM code, 314.01) increased 84,000 from 21,000 in 1993 to over 100,000 in 2005, representing a 400% increase.

Fixed Effects
Model 1. The results from the simple ordinary least squares model screening process are presented in Table 3. For year 2007, the only covariate to pass the screening at the .05 level was population (estimate: -0.053±0.02, p-value: 0.01). Population retained a significant association with all-listed ADHD hospital discharges in all subsequent years. Precipitation as a lagged indicator was correlated with all-listed ADHD in 2008. Therefore, population and precipitation were added to the two-way fixed effects final model. Specifically, P = 1 and M = 3 given these screening results. None of the arbitrarily selected herbicides, including halosulfuron, clethodim, and atrazine, passed the initial screening at either USGS estimate.   For every kilogram increase in glyphosate use (high estimate) the prior year, there is a 5.54E-08% increase (p<0.01) in all-listed ADHD hospital discharges as a percent of all mental health disorder discharges in each HCUPNET reporting state the next year (2007)(2008)(2009)(2010), and this effect is seen in the absence of glyphosate use significantly predicting EPSP synthase inhibitor resistance events (Table 3). Fig 5 shows the heteroskedastic trends in this model. These findings were replicated when low USGS herbicide estimates were used for this model (data not shown): the estimate for glyphosate use jumps 2%.
To put this result in context, the state of Nebraska increased glyphosate use by 3,252,428 kilograms between 2006 and 2009. That translates to a 0.18% increase in all-listed ADHD hospital discharges, representing 42.9% of the total increase in all-listed ADHD hospital discharges between 2007 and 2010 in that state.
For every percentage increase in maize for food and /or industrial use in the prior year, as a percentage of total U.S. maize for food and/or industrial use, there is a 0.027% increase (p<0.05) in all-listed ADHD hospital discharges as a percent of all mental health disorder discharges in each HCUPNET reporting state the following year (2007-2010).
To further help inform our analysis on glyphosate as a single predictor for other mental disorders, we used the same HCUPNET states and years and found a marginally significantly negative relationship after data adjustment between glyphosate use and mood disorders as a percent of total all-listed mental health disorder discharges (adjusted estimate:-1.32e-07± 7.37e-08, p-value: 0.076). During the period from 2001 to 2006, however, we find the estimate to be positive, although the estimate is not statistically significant after data adjustment: glyphosate use predicts affective disorders (adjusted estimate: 1.69e-07±1.07e-07, p-value: 0.12), and this occurs during a period when glyphosate use significantly predicts EPSP synthase inhibitor resistance events (adjusted estimate: 4.58e-08±2.30e-08, p-value: .049).
Although patient age was not a focus in this model, we did find unique associations between glyphosate and the ADHD patient age. When we replace the dependent variable with the percentage of total all-listed ADHD diagnoses that occurred within a certain age cohort (see Fig 4), we found the association to be only positive, although not significant, in young adults. Interestingly, there was a strong negative association between glyphosate use and middle-aged adults, 45-64, (unadjusted estimate: -9.34e-07±2.50e-07, p-value: 0.0003) and seniors (not significant). With respect to patient gender, when we replace the dependent variable with the percent of total all-listed ADHD diagnoses that occurred within a gender (male or female), we found the association to be positive, although not significant (p < .05), with males and negative, although not significant, with females (data not shown).
As an additional supplemental note, we have been able to replicate this model in a simple linear regression with USGS estimated national glyphosate use as the predictor variable and non-normalized national all-listed ADHD hospital discharges from the CDC Vital and Health Statistics Reports (from 1993 to 2005) (unadjusted estimate: 1.31E-06±1.04E-07, p-value: 7.23E-08).
Model 2`. The study significance was adjusted for multiple comparisons across model categories; the revised p-value of 0.0125 was used to determine significance. The urbanization data suggests that glyphosate's significantly positive relationship with all-listed ADHD hospital discharges is marginally significant in urban counties. As shown in Table 4, for every kilogram increase in glyphosate used in urban counties, the model predicts a 2.40e-06 percent increase in the percentage of all-listed ADHD hospital discharges that are attributable to patients living in these areas. The Durbin-Watson test indicated no serial correlation among the residuals in this sample (DW = 2.30, p-value = 0.96), and the data was determined to be heteroskedastic using the Breusch-Pagan test (BP = 9.90, df = 4, p-value = 0.042). The Fisher type test indicated the presence of stationarity in at least one panel for all covariates. S2, S3, S4, and S5 Figs provide model diagnostics for each panel set included in Model 2.
In urban counties in Missouri, for example, there was an increase of 372,246.10 kilograms of glyphosate used (representing an 18.7% increase) between 2006 and 2009. The percentage of the total all-listed ADHD diagnoses among patients living in urban areas, as a percent of total all-listed ADHD discharges in the state, rose 3.18% from 2007 to 2010. Using this econometric output, about 28% of the increase in ADHD discharges can be attributed to the increase in glyphosate used in these areas in the prior year.
We have performed additional iterations of this model to accommodate the flexible interpretations of the rural and urban continuum codes among government agencies. For example, the USDA ERS categorizes all counties in the state of Rhode Island as major metropolitan (code 1), while HCUPNET provides data for both metropolitan and fringe metropolitan (codes 2 and 3). We, therefore, combined these categories in an additional analysis and found the association between estimated glyphosate use and all-listed ADHD discharges to be negative in Table 4. The effect of urbanization on the relationships between glyphosate use and ADHD hospital discharges, using two ways fixed effects within model in R and adjusting for heteroskedasticity and serial correlation and including all significant covariates that passed screening in Model 1.
Glyphosate use (2006)(2007)(2008)(2009)  the combined metropolitan code, although it was not statistically significant. In a second iteration, we revised our urban category to codes 4 and 5 only and relocated codes 6 and 7 to the rural category. The same results were noted. These results confirm the first model iteration and suggest that the model is robust to flexible continuum interpretations. Model 3. The study significance was adjusted for two comparisons across model categories; the revised p-value of 0.025 was used to determine significance for the glyphosate estimates. The two-way within model fixed effects indicate that glyphosate use is associated with farm use of nitrogen fertilizers. In fringe metropolitan counties, the association is significantly positive (unadjusted estimate: 188.38±16.52, p-value < .001). Serial correlation was not present among the residuals (DW = 2.78, p-value = 0.99), and the data was homoskedastic (BP = 0.25, df = 1, p-value = 0.62). We used second-order differencing of both variables and confirmed no unit roots with the ADF test in this longer time series.
In the urban counties as we have defined them, the association is significantly positive. (adjusted estimate: 150.53±24.91, p-value < .001). Serial correlation was not present (DW = 3.41, p-value = 1) and data were homoskedastic (BP = 0.21, df = 1, p-value = 0.64). We differenced both variables twice to ensure stationarity in the urban time series. S6 and S7 Figs present model diagnostics for the panel sets used in this model.

Discussion
The current study suggests that glyphosate use may be contributing to all-listed ADHD hospital discharges in the United States. The main findings of the study include the following: 1. Normalized glyphosate use, as estimated by the USGS and used as a lag indicator, significantly positively predicts normalized ADHD in the 33 HCUPNET reporting states from 2007-2010. Normalized total food and industrial use of maize as a lagged indicator also significantly positively predicts normalized ADHD in the same data sample. These findings are highly statistically significant even after controlling for individual state differences, strong correlations over time, and other significant associations reported in the ADHD literature.
2. Glyphosate use strongly positively predicts glyphosate weed resistance events in states from 2001 to 2010 but not in the 33 HCUPNET reporting states or all available states (N = 47) from 2007 to 2010. Therefore, the development of weed resistance is not a significant component in the relationship between glyphosate and ADHD during our study period. Our working hypothesis is that the formulation was altered significantly and marketed heavily around the 2006 time frame and maize was increasingly used as a residual cover crop in agriculture. The data on glyphosate and other mental disorders further supports this hypothesis. The changes in these agricultural inputs and perhaps others contributed to the impairment of weeds to develop resistance to herbicide use in the subsequent years. We further hypothesize that these input changes may be contributing to ADHD in humansespecially young adults who tend to be more adventurous outdoors and therefore exposedthrough the release of N 2 O from altered microbial constitutions in soil and water systems. Since the ingredients in the herbicide are under patent control, they are not required to be listed, and thus it is difficult to determine the contributing factors.
3. Our urbanization coding system revealed that glyphosate's positive effect on ADHD remained marginally significant only in urban counties between 2007 and 2010, suggesting that urbanization may be a factor in the relationship. Recent historic population migration from nonmetro counties might explain the negative association with glyphosate in rural areas [63].
4. Glyphosate use is significantly positively associated with farm use of nitrogen fertilizers between 1992 and 2006 in urban and fringe metro counties, suggesting that farm use of nitrogen fertilizers may be contributing significantly to the association between glyphosate use and all-listed ADHD hospital discharges prior to 2007, but there may be other land interactions, such as glyphosate's direct impact on soil microbiota and nitrification inhibition especially after patent reformulation, that could be driving the marked increase in ADHD seen after 2006.

The link between glyphosate, nitrogen inputs, and soil microbiota
Glyphosate's role as a mineral chelator in soil [64] could cumulatively dampen manganese bioavailability in soil. Manganese is important for nitrogen assimilation in plants. If less manganese is available, less nitrogen assimilation occurs. This would explain the strong positive association we have found between glyphosate and nitrogen-based fertilizers. This strong association could contribute to indirect agricultural N 2 O emissions from leached anthropogenic nitrogen into urban drainage networks [59,60]. The findings in this paper suggest some kind of modification in the use of glyphosate around 2006. The ability for estimated glyphosate use alone to predict mental health disorders changed significantly from 2001 to 2010. For example, estimated glyphosate use positively predicts affective (i.e., mood) disorders from 2001 to 2006, but the estimate for mood disorders turns negative during 2007 to 2010. Estimated glyphosate use does not significantly predict psychotic disorders early in the decade but then does in the period defined in the current study, 2007 to 2010. ADHD was not a recorded discharge in HCUPNET until 2007, although data is available for national all-listed ADHD hospital discharges. These changes suggest some sort of change in glyphosate use in agriculture, and the inability of estimated glyphosate use to predict EPSP synthase inhibitor resistance events in the latter part of the decade, whereas it had in the earlier part of the decade, is evidence of a formulaic inconsistency [65].
Egamberdiyeva et al. [66] have shown that adding nitrification inhibitors such as oxalates inhibits net nitrification in calcareous soils. Furthermore, Wan et al. [67] demonstrated that ammonia-oxidizing archaea (AOA) were more responsive to a simazine herbicide application to agricultural soils. Denitrifying bacterial populations play a critical role in nitrous oxide emissions from amended soil [68]. Therefore, increasing use of herbicides and adding nitrification inhibitors like oxalates to fertilization programs in the United States would theoretically result in an attenuation of the nitrification process that is the conventionally regarded rate-limiting reaction contributing to N 2 O release from agricultural soils [69].
However, this inhibition of nitrification could potentially lead to nitrogen mismanagement in the soil, including an accumulation and potential release of soil ammonia (NH 3 ) after ureabased fertilizer application [70]. Compared to ammonia-oxidizing (AO) bacteria (AOB), AOA strains have a much higher substrate affinity for ammonia and oxygen and can therefore thrive in the crust of low nutrient, microaerophilic soils with varying ammonia concentrations [71][72][73]. The growth of archaea in agricultural soils could, therefore, reasonably serve to remediate soil from intensive agricultural herbicide and fertilization programs worldwide. Microbes, including the documented emergence of archaeal strains in 2007, have been shown to be integral decomposers of residual inputs and biosolids in corn agroecosystems [74][75][76]. In this process, these archaeal strains are also known to be very potent contributors to soil N 2 O emissions and may be able to contribute even amid the presence of any nitrification inhibition [77], indicating a widespread genomic adaptability distinct from AOB populations [78][79][80]. It has also been suggested that N 2 O is emitted as a spontaneous metabolic AO intermediate [77]. Recent research reveals that ammonia-oxidizing archaea are underestimated contributors to N 2 O emissions in soil [81], and this discovery is consistent with changing glyphosate dynamics.

N 2 O and ADHD
N 2 O has been used as an anesthetic agent in health care settings, primarily dentistry. However, there is a growing reappraisal of the compound's potential adverse health effects. We propose that N 2 O may be a potential nexus between estimated glyphosate use and increasing all-listed ADHD hospital discharges, especially during the period of study.
Yagiela [82] explains that many of the adverse effects seen with N 2 O inhalation are attributable to its reaction with the reduced form of vitamin B12 in the body [83]. Vitamin B12, or cobalamin, is a critical co-factor for methionine synthase, an enzyme necessary for nucleic acid synthesis and methylation reactions. Exposure of infants to N 2 O during surgery resulted in a significant increase in plasma levels of homocysteine, due to an impaired ability to convert it to methionine [84].
Methionine synthase inactivity interferes with methylation of the homocysteine subunit of the dopamine receptor 4 [85], leading to an inactivation of dopamine-induced phospholipid methylation (PLM) through dopamine receptor 4 (DR4). NE has been shown to activate DR4 in rat lateral habenula [86], so it is plausible for there to be excess NE amid DR4 inactivation. Another mode of action of N 2 O is thought to be inhibition of NMDA receptors [87,88], which has been implicated in ADHD [89]. Since the United States Food and Drug Administration has recently associated ADHD medications with adverse symptoms of psychosis [90], we propose that environmental exposure to N 2 O -likely from a potent reshaping of soil and water dynamics, especially in urban areas, attributable to intensive agricultural operations, principally including reliance on glyphosate-is the mechanism behind the positive association between glyphosate and all-listed ADHD and psychotic disorders during the period 2007 to 2010.

N 2 O and Autism Spectrum Disorders (ASD)
Musser et al. [91] have shown a strong familial transmission of autism spectrum disorder (ASD) and ADHD. Mothers who have ADHD have a significantly higher risk of having a child with ADHD or ASD. Furthermore, Hawi et al. [92] demonstrated a preferential paternal transmission of several ADHD risk genes, and N 2 O has been implicated in causing genetic damage in occupational workers [93]. This data suggests an etiological link between ADHD and ASD, and we propose that link could be inhalational N 2 O exposure.
This increase has been mitigated, at least in part, by large-scale fossil fuel replacement efforts. However, Nobel Prize laureate Crutzen and others [94] have noted that increasing N 2 O from agro-biofuel activities may counter the savings seen from fossil fuel replacement strategies, and this revealing finding did not thoroughly consider mode and effect of pesticide delivery in agricultural systems. The fixed effects analyses and the theoretical implications surrounding it suggest that these data ought to be considered when exploring total agricultural N 2 O emission burden.
In rats, gestational N 2 O exposure has been shown to affect growth and reflex responses in pups between postnatal days 14 and 21 [95] as well as body laterality in cultured rat embryos [96], and these are consistent phenotypes found in subjects with high functioning autism and Asperger's syndrome [97,98]. It has been shown that cardiac parasympathetic activation is dependent upon nitric oxide (NO) synthesis [99], and higher NO production is associated with ASD [100], suggesting parasympathetic dominance in ASD. Inhalational nitrous oxide exposure in human subjects contributed to a suppression in sympathetic cardiac activity [101].
Aerobic exercise, a marker of sympathetic activity, for at least 30 minutes three times a week among pregnant women induced a significantly lower autonomic fetal heart rate (HR) and greater heart rate variability (HRV) and this continued in infancy [102,103], suggesting that parasympathetic domination in utero, perhaps attributable to chronic inhalational N 2 O exposure during fetal organogenesis, may predispose countering effects that constitute the baseline physiology seen in ASD. Toichi and Kamio [104] found parasympathetic activity remained activated in autistic subjects during a mental task. Similarly, Casanova et al. [105] saw a significant increase in cardiac vagal control in autistic subjects treated with low-frequency (LF) repetitive transcranial magnetic stimulation (rTMS). Daluwatte et al. [106] demonstrated constricted pupillary light reflexes (PLR) in ASD subjects, suggesting a lowered parasympathetic modulation during sensory stimulation. It should also be mentioned that heavy metal exposures may exacerbate parasympathetic dominant states [107,108]; interestingly, it may be the case that levels of these heavy metals may be elevated in ASD subjects out of a consequence of this heightened parasympathetic state induced in utero. In other words, the ASD phenotype, through reduced baseline sympathetic activation, may be better able to tolerate higher endogenous levels of such toxins [109], whereas sympathetic activation in non-ASD subjects may aggravate the toxicity of heavy metal exposures [110,111]. The current N 2 O hypothesis therefore supports the notion of autonomic dysregulation in ASD and suggests that the origin of such dysregulation may occur during gestation, although more comprehensive intergenerational research is needed.
It is likely that the sympathetic hyperarousal seen in ASD is the physiological manifestation to contain this baseline parasympathetic domination arising from N 2 O exposure in utero and any other xenobiotic exposure after birth. That is to say that many of the behavioral phenotypes that are noted in ASD, including anxiety, hyperactivity, stereotypical compulsivity and compromised social interactions, may arise out of a compensating, homeostatic preference to curtail the dominating parasympathetic influences and may not be explicitly indicative of any genuine etiological link to the disorder. A few animal studies on the effect of prenatal N 2 O exposure indicate behavioral differences, although care should be taken to avoid direct human ASD comparisons [112,113].
It is advised that policy makers pass greater funding measures to institute sustained and comprehensive N 2 O monitoring. Aside from its implications as a potent greenhouse gas, we have presented limited existing evidence, using animal and human models, on the complementarity between the fetal/infant health effects from inhalational N 2 O exposure in utero and ASD. We note that N 2 O monitoring from agricultural sources discontinued in 2009 at the EIA. Moreover, there is a worrisome lack of data on N 2 O emissions from bodies of water, even though government agencies like the United States Geological Survey and others have pointed out the potential significance of such sources in overall environmental N 2 O burden. Emerging ASD advocacy points to a wandering desire to seek bodies of water in autistic subjects, and this may arise from a low-grade addictive phenotype to N 2 O developed in utero [114]. Future epidemiological studies must consider the influence of exposure to large bodies of water and risk of ADHD and/or autism, and case studies may be needed to ascertain whether a preference for bodies of water in autistic subjects is a significant underlying feature in ASD or not. As mentioned, studies on soil science need be conducted to clarify the emerging role of archaeal strains in agriculture soils and how modern agricultural practices might be contributing to their emergence. Furthermore, data from HCUPNET suggest that adult ADHD is gaining share as a percentage of total ADHD hospital discharges and given our conjecturing here and the findings of others, this may forewarn a greater increase in ASD is still yet to come.

Study Limitations
There exist several limitations to our study. These limitations include, principally, the data itself. We were limited in the number of states we could include due to HCUPNET availability and the availability of regressor data, although there is strong agreement between data sources used in this study. Normalized estimated lagged glyphosate use does positively predict both national (non-normalized CDC all-listed ADHD hospital discharges: 1993-2005, OLS) and state-level (HCUPNET normalized state-level estimates: 2007-2010, two-way fixed effects) alllisted ADHD hospital discharges. Between 2007 and 2010, national all-listed ADHD hospital discharges in HCUPNET, as a percent of total all-listed mental health discharges increased 18.5%, while the average cross-sectional increase across available HCUPNET states between 2007 and 2010 was more modest at 9.6%, suggesting our findings may be conservative. Furthermore, the USGS data on estimated glyphosate use is normalized according to harvested agricultural acreage within a county and does not overtly account for non-farm use. It is possible that non-farm use of glyphosate and nitrogen fertilizers could account for some of the increase in ADHD hospital discharges, as our data cannot exclude this as a possibility. The USGS report indicates that the ratio of non-farm to total nitrogen input estimates from 1987 to 2006 was higher in several states reporting large increases in all-listed ADHD hospital discharges from 2007 to 2010 (i.e., New Hampshire, Massachusetts, and Rhode Island). Our data, however, suggest that agricultural use of glyphosate is a significant predictor of all-listed ADHD hospital discharges in HCUPNET reporting states from 2007 to 2010, even after controlling for state fixed effects, strong correlations over time, and other significant associations with ADHD. While we have shown that glyphosate herbicide resistance events and farm use of nitrogen fertilizers are strongly associated with estimated glyphosate use, neither phenomenon can fully explain the increase in all-listed ADHD hospital discharges in our time period of study.
A second limitation of our study is the lack of genetic markers in our model. Genetic risk factors are being studied for ADHD disorders, in addition to other mental disorders. More specifically, LaHoste et al. [115] have implicated genetic polymorphisms within the DR4 gene as a potential contributor to ADHD, and there are many other studies that point to a genetic predisposition for substance and alcohol addiction [116,117]. In our analysis of national HCUPNET hospitalization data, we find that, as a percentage of total discharges for all mental disorders, discharges for substance/alcohol abuse have actually plateaued (8.03 percent in 2007 to 7.38 percent in 2010 for substance abuse and 8.36 percent in 2007 to 7.66 percent in 2010 for alcohol-related disorders), suggesting that genetic polymorphisms contributing to these conditions have remained static in the health-care seeking population. A significant increase in genetic polymorphisms associated with the mental conditions under study is therefore probably not a plausible explanation for the results obtained, although we cannot say unequivocally that this is the case.
A third limitation is that we have not fully considered the interaction of other variables. For example, we have included maize grain in the current study and found that it is a significant predictor of all-listed ADHD hospital discharges even after controlling for glyphosate, suggesting the mechanism could be related to the genetic composition of the corn as discussed earlier, its use as an agricultural residual cover crop, or another interaction yet to be identified. There are other environmental toxicants that have been associated with ADHD, as well, and we have not accounted for these variables in our models [118,119], principally because we were unable to obtain reliable annual state-level estimates for the years of interest. Furthermore, our study period was subject to influence by the worldwide economic recession. As socioeconomic status has been implicated in ADHD [50], it is possible that our variables of interest (i.e., home price index) were not adequate enough in capturing the potential effect of the economic recession, and we find this to be a limitation. It should also be mentioned that, while the ADHD population in this study was based on ICD-9 diagnoses, there have been significant diagnostic improvements and greater awareness of mental conditions like ADHD in recent years. It is therefore difficult to know for certain whether such conditions are actually increasing or are the result of major diagnostic advancements and greater public awareness [120]. This is a significant quandary, especially in light of the economic difficulty of the period [121]. Therefore, attempts should be made to replicate these results using a diverse set of patient populations, diagnostic environments, and economic conditions.

Conclusions
We have shown that agricultural estimated glyphosate use and total food and industrial use of maize grain are significant lagged predictors of all-listed ADHD hospital discharges in HCUP-NET reporting states from 2007 to 2010 even after controlling for state fixed effects, strong correlations over time, and other documented associations with ADHD in the literature. The association appears to be most significant in urban counties. Fig 6 is a global schematic that demonstrates the indirect mechanism that may underlie the increase in ADHD seen during our study period. These significant associations occur amid the absence of estimated glyphosate use predicting resistance events and marked swings in the ability of estimated glyphosate use to predict other mental health discharges. We propose that a patented reformulation in glyphosate could be enhancing both its herbicidal and bactericidal properties, enabling the dominance of the more genetically adaptable archaea in agricultural soils. Maize residuals may also be used as an herbicide control, and this may be contributing to the increase in direct agricultural N 2 O emissions, most especially given the changing soil dynamics possibly induced by glyphosate use. Impairment of several physiological mechanisms, including NMDA receptors and methionine synthase, is a likely explanation for how glyphosate may significantly contribute to an increase in all-listed ADHD hospital discharges. The precise mechanisms mentioned here require further empirical support, and this ought to include a systematic reevaluation of the crucial intermediates involved. Histograms were produced in R using the hist(residual(x)) argument. Residual plots were created in R by subtracting residuals from the observed ADHD percent and plotting these fitted values (yhat) against model residuals. (TIF) S1 Table. The ICD-9 codes that are included in the HCUPNET CCS diagnosis category (652) of attention-deficit, conduct, and disruptive behavior disorder, along with a DSM-IV-TR translation provided by the American Psychological Association. (DOCX)