Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Exploring the risk of hypospadias in children born from mothers living close to a vineyard

Abstract

Hypospadias (H) is a common birth defect affecting the male urinary tract. It has been suggested that exposure to endocrine disrupting chemicals might increase the risk of H by altering urethral development. However, whether H risk is increased in places heavily exposed to agricultural pesticides, such as vineyards, remains debated and difficult to ascertain. The objective of the work is to test the possible association of H with residential proximity to vineyards. Residential address at birth of 8,766 H cases born 1980–2011 was taken from 17 specialized surgery centers. The geographical distribution of vineyards was obtained from the European Land Parcel Identification System (LPIS) and the distance of address to the nearest vineyard was computed. A first estimate of the variation of H relative risk with distance to vineyards was obtained using as controls 13,105 cryptorchidism (C) cases operated during the same period in the same centers. A separate estimate was obtained from a case-control study using “virtual controls” (VC) defined as points of the map sampled to match the demographic distribution of births within the recruitment territories of the study centers. Non-exposed patients were defined as those with a residence between 5,000 and 10,000 m from the closest vineyard. The residential distance to vineyard was smaller for H than for C cases (p<10−4). We found 42/8766 H cases (0.48%) and 50/13,105 C cases (0.38%) born to mothers living within 20 m of a vineyard. The odds ratios for H were 2.48 (CI: 1.0 to 5.1) and 2.4 (CI: 1.3 to 4.4), vs C or vs VC, respectively, when pregnant mothers lived 10–20 m from a vineyard. In conclusion, our study supports that children born to mothers living close to a vineyard have a two-fold increased risk of H. For environmental research, the use of VC provides an alternative to classical case control technique.

Introduction

Hypospadias (H) is a defect in the urethral development occurring between the 8th and the 16th week of pregnancy causing the meatus to be ill-located on the ventral side of the penis. While H is a frequent congenital malformation that needs early correction by urethroplasty [1], its epidemiology has remained imprecise. Since the early 2000s, H incidence has shown large and unexplained differences across populations and time periods [25]. In 2003, H incidence was estimated to be 8/10,000 live male births in the US [2] and 18.6/ 10,000 in a study of 23 European national registries [3]. H incidence had increased until 1999 according to 36 years of Swedish data, [6], but no variation was observed between 2004 and 2010 across European registries [3].

Since less than 1% of H cases are due to rare monogenic syndromes, the vast majority are still called “idiopathic” and are presumably due to unknown developmental interaction between gene variants and environmental factors. Genetic predisposition is important, as shown by the 12 to 20-fold augmentation of H) among first-degree relatives [7]. Genome-wide studies have identified a few low-risk genomic variants associated with H [8], while epigenetic studies have started to explore the association of H with DNA methylation [9, 10]. The contribution of environmental factors is likely to be prominent, but has been difficult to assess. A link of H with decreasing male fertility suggests the existence of a “testicular dysgenetic syndrome” that has emerged in developed countries [11].

Agricultural pesticides—prominently suspected endocrine disruptors—remain the main environmental suspects, but studies have yielded variable results across countries, regions and time periods [12], Case-control studies of parental occupation did not show an increased risk of H in sons of women working in farming and gardening [1315], but none had specifically focused on people working in vineyards.

A direct way to assess the exposures of cases and controls to agricultural pesticides is offered by Geographical Information Systems (GIS) and the increasing availability of environmental data bases informing on land use [16]. For example, the comparison of exposures of 354 cases with 727 controls with regard to 38 pesticides within 500 m of pregnant mother’s address in Arkansas found a slight association of H with diclofop-methyl (OR = 1.08), and surprisingly a stronger protective effect of other pesticides (e.g. permethrin: OR = 0.37) [17]. Among 304,906 singletons born in North Carolina, the estimated amounts of pesticides used within 500 meters of maternal residence yielded only a marginal association with H in 856 cases [18]. In France, an increase of H incidence at the proximity of Montpellier vineyards has been debated [19] [20]. The fact that many vineyards are located near houses in many countries calls for robust epidemiological evidence based on reliable controls. Notably, a reliable control is”someone who–if he had developed the disease- would have been recruited among the cases of the study” [21].

Surgical centers do not usually record epidemiological information on the exposure history of their patients, but all of them record basic demographic information for billing purposes, including mother’s address at childbirth. This low-cost hospital information crossed with available geographic environmental databases was used for the current study. We have used two independent approaches to define appropriate controls to which the cases should be compared. First we used as controls the cases of C recruited in the same network of surgical centers than those used for H. These C patients could be considered as H controls since all included cases of C had a normal penis and meatus, while none of the studied H cases had mal-descended testes. We did a second independent case-control analysis, where controls were “virtual” controls (VC) defined through an algorithmic approach described below.

Material and methods

Patients

To be included in this study, H cases were the clinically significant cases that had reconstructive surgery in the line of the current practice [22], and C cases had to have orchidopexy [23]. Minimal glandular hypospadias cases that are considered clinically not significant, thus not operated, were not part of our study. Under coordination by L.Esterle and Pr P.Mouriquand (Lyon), the following 17 French centers of urological pediatric surgery (Table 1) collected retrospectively the date of birth and address at birth of the parents of all boys born after January 1st 1980 who had H (coded as ICD Q54) and C (coded as ICD Q53): Angers (Pr G. Podevin), Bordeaux (Pr B. Frémond), Caen (Pr P. Ravasse), Colmar (Pr S. Geiss), Grenoble (Pr B. Boillot), Lyon (Pr P. Mouriquand), Marseille and Toulon (Pr P. Alessandrini), Vandoeuvre-lès-Nancy (Pr JL. Lemelle), Nantes (Pr MD. Leclair), Reims (Pr PJ. Lefebvre), Rennes (Pr E. Dobremez), Saint Étienne (Pr F. Varlet), Strasbourg (Pr R. Moog), Saint Vincent de Paul at Paris (Pr F. Bargy), Toulouse (Pr J. Moscovici), Valence (Pr B.Defauw).

thumbnail
Table 1. Numbers (and %) of cases and controls analyzed in each collaborative center.

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

All children operated for hypospadias in the University hospitals participating to this study had all been carefully examined at a pediatric endocrinology unit belonging to the Centre National for Disorders of Sex Development (DSD). Cases resulting from a genetic disorder of sex differentiation or hypogonadotropic hypogonadism were excluded. At total, 8,766 H and 13,105 C cases were recorded by May 2013. (see Table 1). We computed the distance of each case to the nearest vineyard (see below), and constituted the database which is analyzed in this paper.

The distribution of the year of birth of the H and C cases is shown in Fig 1.

thumbnail
Fig 1. Distribution of the year of birth distribution cases (H: Left, C: Right).

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

Ethics, data confidentiality

The research protocol was approved by the Ethics committee of Ile de France (DC-2009-1041). The parents of children provided written informed consent. To protect the confidentiality and data of participants, all were assigned a unique identification number without identifying information. The security conditions of the database received the agreement of the “national information science and liberties commission” (ref CNIL 912529). Under this agreement, CNIL can allow external researchers to get access to the original databases. Requests must be sent to the board of the HCFCG (HypoCrypto French Collaborative Group, Inserm, Batiment Pincus, Hopital de Bicêtre, 94276- Le Kremlin Bicêtre).

Assessment of the residential proximity of people to vineyards

We used the Land Parcel Identification System (LPIS) to assess the residential proximity of the homes of the children to vineyards. The origin of LPIS is a decision of the European Union Council Regulations in 2003 to create an information system on the land use of European agricultural parcels.

In LPIS, an agricultural parcel is defined as a continuous area of land on which an individual farmer cultivates a single crop. The data is encoded using the shapefile format, which is a GIS file format describing vector features such as points, lines, and polygons. Polygons are used in order to specify the borders of the parcels. The geographic projection system used to specify the edges of these polygons is Lambert-93. Each item element in shapefile has attributes. The attribute of a polygon is a number that identifies the crop being grown on that parcel. LPIS information was made available for years 2010 to 2014. We used in this work the 2012 database that was freely downloadable from the French data.gouv.fr web portal (https://www.data.gouv.fr/fr/datasets/registre-parcellaire-graphique-2012-contours-des-ilots-culturaux-et-leur-groupe-de-cultures-majorita/). At total the current study is based on the 340,336 LPIS parcels with vineyards, for a total surface of 381,395 Ha.

The exposure metric that we use in this work is the individual’s distance to the nearest parcel of vineyard. In practice, we had to calculate these distances for more than 50,000 individuals (cases plus controls) to 340,336 parcels of vineyard. Each parcel was modeled as a polygon. To avoid calculating exact distances to polygons that are sufficiently far away for having any chance to induce a risk, we chose a cutoff at 10,000 meters and pooled all distances greater than 10,000 m in a single class.

Statistical analysis of the risk associated with distance to the nearest vineyard

Two analyses were performed. The first one was based on patients whose home was located in three large classes of distances to a vineyard: 0–100 m, 100–200 m and 200–500 m. The second analysis studied the risk at distances of less than 50 m from a vineyard, despite the fact that only 82 cases were from homes located there. The three classes of “small” distances studied were 0–10 m, 10–20 m and 20–50 m.

In both analyses we studied H risk in two independent case-control studies, taking the 5,000–10,000 m class of distance as reference. In the first analysis, we compared the distances to the nearest vineyard in the two groups of H and C cases. For this comparison, we computed the variation of the odds ratio of the H risk, by reference to the C risk, in the different zones of distance to a vineyard, taking the 5,000 to 10,000 m distance as the reference.

In the second analysis, we computed the variation of the odds ratio (OR) for the H risk of people living within zones close to vineyards, taking as reference the risk of children born to mothers living between 5,000 and 10,000 m, i.e. roughly 21% of the patients (Table 2). The specificity of our approach is that the controls of this case-control study were “virtual controls” (VC). A “virtual control” is a geographical point of the map which represents a real physical person that would have been chosen as control, if it was practically possible to sample randomly a subject within the same source population as the “case”. We modeled the source population of each hospital as composed of the population living in the circle around this center, in which there was a proportion p of the cases. In the analyses presented here, p was taken at 80%. The algorithm used to sample the VC is the PPS algorithm implemented by Jack G. Gambino as an R package (PPS stands for probability proportional to size). https://rdrr.io/cran/pps/. This algorithm allows to sample people in a territory proportionally to the local density of population. The French National Institute of Statistics (INSEE) 2009 database (https://www.insee.fr/fr/statistiques/2520034) provided us with an estimate the population density. This database provides the 2009 French density of population in each age-class within a 200 m resolution grid. From this 2009 information, we can estimate the spatial distribution of children born, say, in 1995: assuming that mobility was low (which is the case in France), this distribution is simply the distribution of the children aged 14 (= 2009–1995) in the 2009 grid. This spatial distribution is then used to randomly chose–thanks to PPS- a cell within all the cells that constitute the source population. The VC is then set at the center of the 200m x 200m cell. We took 5 virtual controls per case by repeating the algorithm 5 times.

thumbnail
Table 2. Proportions of hypospadias and cryptorchidism in different classes of distance to the nearest vineyard.

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

Socioeconomic environment

The 2009 version of the geographic information system of the French National Institute of Statistics (INSEE) (https://www.insee.fr/fr/statistiques/2520034) provided us with an estimate of the density of population, by age class, in each square of a 200m x 200m grid, repeat and an estimate of the local socioeconomic environment by using the Townsend deprivation index (TDI) [24] at the place of birth. A higher TDI score implies more severe social deprivation.

Statistical packages

The odds ratios quantifying the association of disease and distance to the nearest vineyard were estimated using logistic regression adjusted for the coordinating center, to account for difference in the distance distribution between centers. Additional analyses were adjusted for TDI. Subgroup analyses by birth period were also carried out. Analyses were carried out using the R statistical software version 3.6.3, with the glm function of the stats package.

Results

The distribution of the distance of the children affected with H or C to the nearest vineyard is shown in Table 2. The two distributions were significantly different (chi square = 38, 9 df, p<10−4). H cases were more frequent than C in all classes of distances below 5,000 m to vineyards, with the exception of those located between 0 and 10 m.

We first analyzed wide classes of distance (0–100 m, 100–200 m, 200–500 m, 500–1000 m), where the number of H patients was large (between 145 and 376). We found that odds ratios were > 1 in all these classes of distance, and were significantly > 1 in the 0–100m class (Table 3) when controls were patients with cryptorchidism (Column 1 of Table 3), and when controls were “virtual controls” after adjustment on the Townsend index (column 2b of Table 3).

thumbnail
Table 3. Odds ratios, with 95% confidence intervals, of the risk of hypospadias by classes of distance to the closest vineyard.

(large distances).

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

In the second part of the analysis, we studied the patients whose houses were located at close distances of vineyards (less than 50 m), although their numbers were small (see Table 2). Again, we estimated the risk of H possibly associated with a short distance from the nearest vineyard by computing the odds ratios of H patients versus C patients, and the odds ratios of patients versus VC. The two analyses gave comparable results. The odds ratio of H patients in the 0–50 m was significantly > 1 (column 1, where controls are patients with C and column 2b, where controls were “virtual controls”). When we focused on classes of addresses at very short distances of a vineyard (second part of Table 4), we found that the odds ratio of H patients belonging to a home located between 10 and 20 m from the closest vineyard was significantly > 1 (p = 0.029 when controls are C patients, p = 0.007 when controls are “virtual controls”). The results were comparable after adjustment on the Townsend Index.

thumbnail
Table 4. Odds ratios, with 95% confidence intervals, of the risk of hypospadias by classes of distance to the closest vineyard.

(short distances).

https://doi.org/10.1371/journal.pone.0249800.t004

Discussion-conclusion

Common “idiopathic” H results from unknown mechanisms, among which endocrine disruptors are suspected to interfere with gonadal and genital development of the male fetus.

We found in this work that the risk of H was higher in children born to mothers living close to a vineyard.

The strength of our observation is that it was achieved through two independent methodological approaches. Both methods guarantee that controls were from the same hospital recruitment territory than cases. This avoids the bias of controls made up of children with common pediatric diseases that usually live close to the center, whereas cases of H come from far away for a very specialized intervention. Such bias could create the false impression that patients from rural areas would be more at risk than those from urban areas. The controls who were used in the two approaches match the ideal definition of what must be a control group: "individuals who, if they had had the disease (H), would have been recruited in the same centers as the (H) cases were" [21]. In the first approach, the controls were children with C and no H having been operated in the same clinical centers as H. A few studies found a possible positive association between C and exposure to pesticides. However, even if this link is true, it would not invalidate our results. Indeed, the consequence of such a link would be that the “true” odd-ratio is even greater than the one we found. In the second approach, the “virtual controls” were not real persons. However the algorithm that we used guarantees that they are identical to the population of children that would be sampled in the “territory” of the cases. To realize this, we modeled the recruitment “territory” of an hospital as the population living in a virtual circle around this hospital containing a proportion p = 80% of the cases. We also performed the analyses using other models to define hospital territory, and obtained similar results.

The increased H prevalence at short distances from vineyards is an indirect indication supporting, but not proving, a harmful effect of pesticides since our study did not qualify nor quantify the loads and nature of pesticides spread in the vineyards. Indeed, this was an impossible task given the 31-year duration of our H case collection. Not only pesticides, such as fungicides, insecticides, herbicides, bactericides, rodenticides, fumigants, but also fertilizers and other toxic chemicals were and still are commonly used for viticulture. During the long studied period providing our data, the complexity, variety and dynamic mix of chemicals used in vineyards were considerable, and no database was available to quantify precisely their nature and load. The found association of vineyard proximity with H risk could also be due to other environmental factors present close to vineyards. Known endocrine disruptors commonly used in viticulture at the time of the study are vinclozolin, procymidone, and linuron, all forbidden since 2007–2008, but still present in some soils close to vines. The fungicide vinclozolin and its metabolites act as androgen receptor antagonists [25]. Procymidone is another anti-androgenic fungicide that can induce hypospadias in rodents [26]. Linuron, a widely used herbicide, is a weaker androgen receptor antagonist than procymidone and vinclozolin [27]. Nitrates are another category of inorganic pollutants that can disrupt gonadal steroidogenesis [28, 29]. It was interesting in this respect to search whether the association of H with vineyard proximity was stronger in the years before 1995, but the small number of cases did not allow a reliable analysis.

Agricultural substances may diffuse via surface runoff, leaching to field drains, and atmospheric spray drift. Pesticide drift may occur during ground application and after it, since applied substances have different volatility and local persistence. Our choice of focusing on the narrow zone bordering vineyards was influenced by methods used to apply pesticides, which never include aerial sources such as in the US but only land-based devices.

Vineyards are a privileged place for pesticide spraying, undoubtedly one of the most important of all French agriculture. A weakness of our study is that it did not get occupational data for the mothers of cases, although one can expect that a much higher percentage of those living in the vicinity of vineyards are occupationally involved in wine producing activity, thus possibly exposed to pesticide stocks, bisphenol A, phthalates, or other potential endocrine disruptors [30].

It must be noted that the possible risk associated to being a fœtus developed within a short distance of a vineyard concerns only a small proportion of the H cases: H cases born in a house at less than 20 m from a vineyard represent only 0.5% of the total number of H patients (those born at less than 50 m are 1%, and at less than 100 m are 1.8% of the overall H prevalence in France).

The VC approach makes the best possible use, at almost no cost, of the huge data resource that is sleeping in hospitals instead of serving epidemiology. It does not require filling questionnaires or track families to get controls, but instead relies on the basic medical information system. In addition, the method avoids some selection bias, in that we did not have to obtain subjects’ consent to participate, as we used only information in the public domain. For all of these reasons, we believe that the VC method could be used increasingly for well-defined diseases using more detailed, and frequently refreshed environmental databases. A limitation of the method is that VC cannot be studied to compare biological exposome with cases, including pesticide and metabolite concentrations in biological fluids or cells [31].

In conclusion, we have found a significant statistical link between the occurrence of H and the residence of pregnant mothers close to vineyards. However we are conscious that generalizations that focus on ‘positive’ findings rather than more comprehensive views that incorporate null findings and study limitations should be avoided. Future environmental research will have to account for the changing nature and loads of pesticides in modern viticulture. This could be done with the method presented, if high-resolution public reliable databases on the geographical use of different pollutants become available to public health research. In addition, further identification of precise hormone-disruptors will require a deep, systematic and specific chemical investigation of vineyards and their close neighborhood.

More generally, the technique of “virtual controls” that was used in this work can be applied to search for candidate environmental factors in fetal diseases where the only available information is the location of the patient at birth.

Acknowledgments

We gratefully acknowledge the administrative and medical staff of the surgical centers for providing the address of patients used in this study.

References

  1. 1. Springer A, Tekgul S, Subramaniam R. An Update of Current Practice in Hypospadias Surgery. European Urology Supplements. 2017;16(1):8–15.
  2. 2. Manson JM, Carr MC. Molecular epidemiology of hypospadias: review of genetic and environmental risk factors. Birth Defects Res A Clin Mol Teratol. 2003;67(10):825–36. pmid:14745936.
  3. 3. Bergman JE, Loane M, Vrijheid M, Pierini A, Nijman RJ, Addor MC, et al. Epidemiology of hypospadias in Europe: a registry-based study. World journal of urology. 2015. pmid:25712311.
  4. 4. Fernandez N, Perez J, Monterrey P, Poletta FA, Bagli DJ, Lorenzo AJ, et al. ECLAMC Study: Prevalence patterns of hypospadias in South America: Multi-national analysis over a 24-year period. Int Braz J Urol. 2017;43(2):325–34. pmid:27802003.
  5. 5. Yu X, Nassar N, Mastroiacovo P, Canfield M, Groisman B, Bermejo-Sánchez E, et al. Hypospadias Prevalence and Trends in International Birth Defect Surveillance Systems, 1980–2010. Eur Urol. 2019;76(4):482–90. Epub 2019/07/14. pmid:31300237
  6. 6. Nordenvall AS, Frisen L, Nordenstrom A, Lichtenstein P, Nordenskjold A. Population based nationwide study of hypospadias in Sweden, 1973 to 2009: incidence and risk factors. J Urol. 2014;191(3):783–9. pmid:24096117.
  7. 7. Schnack TH, Zdravkovic S, Myrup C, Westergaard T, Christensen K, Wohlfahrt J, et al. Familial aggregation of hypospadias: a cohort study. Am J Epidemiol. 2008;167(3):251–6. Epub 2007/11/29. pmid:18042671.
  8. 8. Geller F, Feenstra B, Carstensen L, Pers TH, van Rooij IA, Körberg IB, et al. Genome-wide association analyses identify variants in developmental genes associated with hypospadias. Nat Genet. 2014;46(9):957–63. Epub 2014/08/12. pmid:25108383.
  9. 9. Vottero A, Minari R, Viani I, Tassi F, Bonatti F, Neri TM, et al. Evidence for epigenetic abnormalities of the androgen receptor gene in foreskin from children with hypospadias. J Clin Endocrinol Metab. 2011;96(12):E1953–62. Epub 2011/09/23. pmid:21937623.
  10. 10. Richard MA, Sok P, Canon S, Nembhard WN, Brown AL, Peckham-Gregory EC, et al. Altered mechanisms of genital development identified through integration of DNA methylation and genomic measures in hypospadias. Sci Rep. 2020;10(1):12715. Epub 2020/07/31. pmid:32728162
  11. 11. Skakkebaek NE, Rajpert-De Meyts E, Main KM. Testicular dysgenesis syndrome: an increasingly common developmental disorder with environmental aspects. Human reproduction. 2001;16(5):972–8. Epub 2001/05/02. pmid:11331648.
  12. 12. Carbone P, Giordano F, Nori F, Mantovani A, Taruscio D, Lauria L, et al. Cryptorchidism and hypospadias in the Sicilian district of Ragusa and the use of pesticides. Reproductive toxicology. 2006;22(1):8–12. pmid:16530380.
  13. 13. Weidner IS, Møller H, Jensen TK, Skakkebaek NE. Cryptorchidism and hypospadias in sons of gardeners and farmers. Environ Health Perspect. 1998;106(12):793–6. Epub 1998/11/30. pmid:9831539
  14. 14. Jorgensen KT, Jensen MS, Toft GV, Larsen AD, Bonde JP, Hougaard KS. Risk of cryptorchidism among sons of horticultural workers and farmers in Denmark. Scandinavian journal of work, environment & health. 2014;40(3):323–30. pmid:24220013.
  15. 15. Rocheleau CM, Romitti PA, Sanderson WT, Sun L, Lawson CC, Waters MA, et al. Maternal occupational pesticide exposure and risk of hypospadias in the National Birth Defects Prevention Study. Birth Defects Res A Clin Mol Teratol. 2011;91(11):927–36. Epub 2011/09/29. pmid:21954192
  16. 16. Nuckols JR, Ward MH, Jarup L. Using geographic information systems for exposure assessment in environmental epidemiology studies. Environ Health Perspect. 2004;112(9):1007–15. Epub 2004/06/17. pmid:15198921
  17. 17. Meyer KJ, Reif JS, Veeramachaneni DN, Luben TJ, Mosley BS, Nuckols JR. Agricultural pesticide use and hypospadias in eastern Arkansas. Environ Health Perspect. 2006;114(10):1589–95. pmid:17035148
  18. 18. Rappazzo KM, Warren JL, Meyer RE, Herring AH, Sanders AP, Brownstein NC, et al. Maternal residential exposure to agricultural pesticides and birth defects in a 2003 to 2005 North Carolina birth cohort. Birth Defects Res A Clin Mol Teratol. 2016;106(4):240–9. pmid:26970546
  19. 19. Blandin M. Les risques et dangers pour la santé humaine de substances chimiques d’usage courant: éthers de glycol et polluants de l’air intérieur. Évaluation de l’expertise publique et des choix opérés. Office parlementaire d’évaluation des choix scientifiques et technologiques, 2008.
  20. 20. Rambourg M, Pascal L, Lasalle J. Les malformations congénitales du petit garçon en Languedoc-Roussillon—Étude des cryptorchidies et des hypospadias opérés à partir des données du PMSI de 1998 à 2001. Institut de Veille Sanitaire, 2004.
  21. 21. Miettinen OS. The "case-control" study: valid selection of subjects. J Chronic Dis. 1985;38(7):543–48. pmid:4008595.
  22. 22. Steven L, Cherian A, Yankovic F, Mathur A, Kulkarni M, Cuckow P. Current practice in paediatric hypospadias surgery; a specialist survey. J Pediatr Urol. 2013;9(6 Pt B):1126–30. pmid:23683539.
  23. 23. Schneuer FJ, Holland AJ, Pereira G, Jamieson S, Bower C, Nassar N. Age at Surgery and Outcomes of an Undescended Testis. Pediatrics. 2016;137(2):e20152768. pmid:26801912.
  24. 24. Townsend P, Phillimore P, Beattie A. Health and Deprivation: Inequality and the North. London: Routledge; 1988.
  25. 25. Gray LE Jr., Wilson VS, Stoker T, Lambright C, Furr J, Noriega N, et al. Adverse effects of environmental antiandrogens and androgens on reproductive development in mammals. International journal of andrology. 2006;29(1):96–104; discussion 5–8. Epub 2006/02/10. pmid:16466529.
  26. 26. Ostby J, Kelce WR, Lambright C, Wolf CJ, Mann P, Gray LE Jr.. The fungicide procymidone alters sexual differentiation in the male rat by acting as an androgen-receptor antagonist in vivo and in vitro. Toxicology and industrial health. 1999;15(1–2):80–93. Epub 1999/04/03. pmid:10188193.
  27. 27. Gray LE Jr. Xenoendocrine disrupters: laboratory studies on male reproductive effects. Toxicol Lett. 1998;102–103:331–5. Epub 1999/02/18. pmid:10022274.
  28. 28. Edwards TM, McCoy KA, Barbeau T, McCoy MW, Thro JM, Guillette LJ Jr.. Environmental context determines nitrate toxicity in Southern toad (Bufo terrestris) tadpoles. Aquatic toxicology (Amsterdam, Netherlands). 2006;78(1):50–8. Epub 2006/03/30. pmid:16569448.
  29. 29. Guillette LJ Jr., Edwards TM. Is nitrate an ecologically relevant endocrine disruptor in vertebrates? Integrative and comparative biology. 2005;45(1):19–27. Epub 2005/01/01. pmid:21676740.
  30. 30. Plank CM, Trela BC. A Review of Plastics Use in Winemaking: HACCP Considerations. American Journal of Enology and Viticulture. 2018;69:307–20.
  31. 31. Vrijheid M, Slama R, Robinson O, Chatzi L, Coen M, van den Hazel P, et al. The human early-life exposome (HELIX): project rationale and design. Environ Health Perspect. 2014;122(6):535–44. pmid:24610234