BACKGROUND
identifying mechanisms of pathogen transmission, including potential environmental sources, is critical to control disease  <cit> .
molecular epidemiology integrates conventional epidemiological approaches with molecular techniques to track specific strains of pathogens in order to understand the distribution of pathogens in populations and environments  <cit> .
this can be used to elucidate inter- and intra-specific transmission pathways and environmental risk factors, from individual to population, and from local to broader spatial scales.
the genus mycobacterium comprises over  <dig> species and several subspecies.
over  <dig> of these can cause disease in livestock, wildlife and humans, occurring worldwide.
mycobacterial diseases such as bovine tuberculosis  have become a major sanitary and conservation problem even in relatively unmanaged natural areas across the world.
similar to other shared diseases, the existence of wildlife reservoirs is limiting the effectiveness of eradication schemes in livestock  <cit> .
in btb, known risk factors for wild ungulates include age, gender, density, spatial aggregation, intra and inter-specific contact, fencing and other habitat features as well as genetic factors  <cit> .
however, most data derive from large scale studies , while detailed information at small spatial scales is still very scarce , and usually, fine associations with spatial and environmental factors are not addressed.
even less information is available regarding the effects on wildlife and livestock of mycobacteria other than tuberculosis .
environmental mycobacteria or mott include a large number of species that can cause serious illnesses in humans, particularly in immunocompromised patients  <cit> .
for example, mycobacterium interjectum has been identified as a causative agent of cervical lymphadenitis in children  <cit> , and of cutaneous infections in immunosuppressed patients  <cit> .
m. xenopi may cause pulmonary disease in humans  <cit> , and m. scrofulaceum may cause cutaneous infections and lymphadenitis  <cit> .
in humans, risk factors for mott infections include immunosuppression, contaminated water and aerosol exposure, and short or old age  <cit> .
mott are widely distributed in the environment, particularly in wet soil, marshland, streams, rivers and estuaries, but each species shows different preferences  <cit> .
because of its habitat characteristics, extension and their sizeable wild and domestic animal populations, doñana national park  in southern spain has been proposed as a good natural laboratory for studying wildlife mycobacteriosis  <cit> .
molecular typing of m. bovis isolates for the period 1998- <dig> showed that wildlife species in dnp were infected only with those m. bovis typing patterns  that were more prevalent in local cattle.
furthermore, the results were suggestive of micro-evolutionary events in the local m. bovis population  <cit> .
in the same period, m. bovis infection prevalence in dnp was 33% in european wild boar , 21% in red deer , and 26% in fallow deer   <cit> .
in a more recent study, we confirmed infection with m. bovis in 52% wild boar, 27% red deer and 18% fallow deer from dnp in 2006- <dig>  and evidenced that m. bovis prevalence decreased from north to south in wild boar and red deer, whereas no clear spatial pattern was observed for fallow deer  <cit> .
three wild ungulates coexist in dnp, wild boar, fallow deer and red deer, along with domestic cattle subjected to btb eradication programs.
we included the wild species as our study models as all are highly susceptible to btb and are known to show high prevalence in the area  <cit> .
in addition, their different ecology and behavior peculiarities  <cit>  can play a role in the epidemiology of mycobacteria, for example, variations in sociability or gregariousness, and scavenging habits.
in addition, different habitats could provide variable environmental suitability for m. bovis persistence  <cit> .
in this sense, scrublands and woodlands are preferably used by red deer and wild boar compared with fallow deer  <cit> .
in this study we used molecular epidemiological techniques to establish the extent of m. bovis strain richness and other environmental mycobacterial species in isolates collected in wildlife and cattle from the dnp, so as the association with social, spatial and environmental factors in this multi-host and multi-pathogen situation.
we hypothesized that infection by mycobacteria would differ among hosts and sites depending on host, pathogen and local environment ecology.
fine-tuning mycobacterial epidemiology in dnp allowed rising a number of relevant questions:  do hosts get infected twice by m. bovis and mott, and can this interfere in m. bovis infection or vice versa?
have new m. bovis types appeared or have any changes in type composition taken place in recent years?
is there an effect of the social group on infection risk?
is there a spatial structure in mycobacteria distribution?
are there species-specific variants of mycobacteria that could be attributed to species-specific behavior patterns  and/or to advanced host species-pathogen interactions?
methods
study area
the study was carried out in dnp, located in south-western spain  and covering  <dig>  ha.
this is a flat region of sandy soils bordering the atlantic ocean, with a maximum elevation of  <dig> m. the climate is mediterranean sub-humid with marked seasons.
in the wet season , most of the marshlands are flooded and wildlife and cattle tend to graze in the more elevated scrublands  <cit> .
in summer, the wetter and more productive ecotone between the scrublands and the marshes supports aggregations of wild and domestic ungulates.
human access is restricted and management is carried out by park authorities.
limited traditional exploitation of some natural resources, such as logging, and cattle and horse rising are allowed.
after  <dig>  when btb in wildlife was first diagnosed in dnp a government-sponsored program was initiated to eradicate btb-positive cattle.
ungulate populations have been culled by shooting .
animal sampling
from april  <dig> to april  <dig>   <dig> european wild boar,  <dig> red deer, and  <dig> fallow deer were sampled within the park by shooting.
the culling of wild ungulates was approved by the research commission of doñana national park in accordance with management rules established by the autonomous government of andalucía.
for each animal we recorded the exact position with gps.
sex and age, based on tooth eruption patterns , were recorded in the field.
a necropsy was performed on site and the presence of tuberculosis-like lesions recorded by macroscopic inspection of lymph nodes and abdominal and thoracic organs  <cit> .
this protocol included the examination of the lungs for the presence of tb-compatible macroscopic lesions during field inspection and a sample was collected.
a tonsil and a head lymph node sample from each individual were collected for culture .
in wild boar, one piece of the tonsils and one from both mandibular lymph nodes were submitted for microbiological studies.
in red deer and fallow deer, one piece of the tonsils and head lymphnode samples, always containing at least half left and half right medial retropharyngeal lymph node, were submitted for culture.
due to logistic and budget constraints, no thoracic or abdominal lymphoid tissues were cultured except when tb-compatible macroscopic lesions were evidenced.
m. scr.
= mycobacterium scrofulaceum; m. int.
= mycobacterium interjectum, m. xen.
= mycobacterium xenopi, m. int.
= mycobacterium intracellulare
social groups were defined as animals sampled the same day at the same site, and with characteristics that were compatible with forming a stable  or seasonal  group.
only part of the individuals belonging to a given social group was sampled.
sampling was performed according to european  and spanish laws , and current guidelines for ethical use of animals in research  and uclm animal experimentation committee.
microbiological procedures
culture: the specimens were liquefied and decontaminated with an equal volume of n-acetyl-l-cysteine and 2% naoh, mixed by centrifugal swirling, and then incubated for  <dig> min .
the reaction was neutralized by adding  <dig> m phosphate-buffered saline , to a final volume of  <dig> ml.
the specimens were concentrated by centrifugation at  <dig>  × g for  <dig> min.
the supernatant was discarded, and the sediment was re-suspended in  <dig>  ml of sterile water.
the sediment was used to inoculate two löwestein-jensen with pyruvate solid medium.
lowëstein-jenssen slants were incubated at 37°c for  <dig> weeks and inspected weekly for growth.
when growth was detected, a smear was prepared to confirm the presence of acid-fast bacilli from suspect colonies by ziehl-neelsen staining.
identification
we identified m. bovis and mott to the species level and characterized m. bovis strains with spoligotyping and miru-vntr typing.
macroscopic morphology of the colonies and pigment production was recorded.
identification at species level was performed with the genotype®mtbc  for the mycobacterium complex strains that allows the differentiation of m. africanum i, m. bovis bcg, m. bovis ssp.
bovis, m. bovis ssp.
caprae and m. tuberculosis/m.
africanum ii/m.
canettii.
mott strains were identified by the genotype® mycobacterium cm and genotype® mycobacterium as mtbc .
the genotype assays were performed according to the manufacturer's instructions: dna extraction by the dna sss method  was followed by pcr amplification of a trait of the 23s rrna gene, as recommended.
reverse hybridization and detection were carried out on a shaking water bath .
the final identification was obtained by comparison of line probe patterns with the provided evaluation sheet  <cit> .
typing
the m. bovis isolates were further characterized by spoligotyping  <cit> .
the amplified product was detected by hybridization of the biotin-labelled pcr product onto spoligotyping membrane .
purified sterile water and chromosomal dna of m. tuberculosis h37rv and m. bovis bcg p <dig> were included as controls in each batch of tests.
the patterns were allocated a number in the m. bovis spoligotyping database.
the results were recorded in sb  code, followed by a field of  <dig> digits as defined on the m. bovis spoligotype database website .
all wildlife isolates  were also subjected to miru-vntr analysis .
extensive documentation  on the service and the genotyping methods is available at the miru-vntrplus website .
the analysis was carried out as described  <cit> , and primers previously reported were used to amplify the etr-a and etr-b loci   <cit> , miru  <dig>   <dig>   <dig>   <dig>   <dig>   <dig> and  <dig>  <cit>  using the recommended pcr conditions.
amplicon sizes were estimated by electrophoresis on a  <dig> % agarose gel at  <dig> v during  <dig> h, using 100-bp ladder .
figure  <dig> presents the spoligotyping patterns, vntr allelic profiles and typing pattern  codes defined for this study.
statistics
chi-square tests were used for between-pair comparisons of prevalences.
to test for the effect of host species vs site regarding the mycobacterial isolates, we used the czechanovsky similarity index  <cit> .
this index considers the list of mycobacterial species recorded in a given host type or in a given study area.
it is calculated by dividing two times the species shared between two lists, by the total number of species of both lists, as follows:
cz=2*/
considering the animals in which any mycobacterial infection was diagnosed, three generalized linear mixed models  were explored to test different explanatory variables that affect the presence of a mycobacterial type or group.
the most common mycobacterial groups were:  m. bovis  m bovis a <dig> and  m. scrofulaceum.
the presence or absence of infection in a mycobacterial group was considered as a binary variable.
the model was fitted using a logit link function.
the model considered social group as a random effect.
the model included host species , the study area and age  as categorical explanatory variables.
the distance to the water  was included as a continuous predictor.
to compare the spatial associations of infection by specific mycobacterial type and hosts, we included as explanatory continuous variable the ratio  between the nearest neighbor distance from host to a different host species with the same type of mycobacteria relative to the nearest distance to a con-specific host with the same type of mycobacteria .
a ratio > <dig> indicates that the nearest distance to a host with the same spoligotype is higher for a different host species.
all the aforementioned explanatory variables we also included in the models interacting with the host species.
due to over-parameterization of the models and zero inflated data, no interactions were included in the m. bovis a <dig> and m. scrofulaceum models.
p-value was set as ≤  <dig> .
we estimated exact confidence limits for prevalence  using sterne's exact method.
RESULTS
mycobacteria species and molecular types
we obtained a total of  <dig> mycobacterial isolates from dnp wildlife.
this included  <dig> m. bovis isolates belonging to  <dig> different typing patterns , and  <dig> isolates belonging to four mott .
m. bovis tps and mott species were isolated from wild boar , red deer , and fallow deer  .
wild boar and red deer had  <dig> m. bovis tps each, whereas fallow deer presented only  <dig> tps.
the number of different isolates per host  was  <dig> in wild boar,  <dig> in red deer and  <dig> in fallow deer .
regarding m. bovis, we identified  <dig> different spoligotyping patterns and  <dig> different vntr allelic profiles .
one spoligotyping pattern was new according to the m. bovis database, and was therefore introduced with code sb <dig>  co-infection of a single host by two m. bovis tps occurred in all three wild ungulate species.
one adult male red deer was infected with tps a <dig> and b <dig>  one adult male and one adult female fallow deer were co-infected with tps a <dig> and e <dig>  and two wild boar  were co-infected with tps a <dig> and b <dig>
mott species found in wildlife hosts included m. scrofulaceum  and m. intracellulare , both found in all host species, m. interjectum , and m. xenopi .
in four deer and four wild boar, m. bovis and mott were isolated concurrently .
in a single wild boar, both types of mycobacteria were simultaneously isolated from the two tissue samples collected and cultured, while in the remaining cases m. bovis was isolated from either lymph nodes or tonsils and the mott from the tissue where m. bovis was absent.
we recorded no cases of co-infection by different mott.
the proportion of wild ungulate m. bovis typing patterns  other than the dominant a <dig> and b <dig> did not differ statistically from 1998- <dig> to 2006- <dig> .
no spoligotyping patterns other than the two dominant ones  were detected among  <dig> cattle isolates in  <dig> and  <dig>
changes in mycobacterial typing patterns over time in dnp
all three m. bovis typing patterns recorded in dnp wildlife between  <dig> and  <dig>  were still evidenced in similar proportions in 2006- <dig> .
however, while only three different tps had been detected in dnp wildlife in the first period, up to  <dig> different ones were found in the second period .
two of these "new" tps  had already been recorded in cattle sampled in dnp between  <dig> and  <dig>  however,  <dig> other tps  had never before been reported in dnp.
in contrast with the situation in wildlife and to data from 1998- <dig>  when  <dig> out of  <dig> cattle spoligotyping patterns were different from a and b, no spoligotyping patterns other than the two dominant ones  were detected among  <dig> cattle isolates in  <dig> and  <dig> .
spatial structure
regarding the mott , m. interjectum was only found in wild boar from eb, in the central part of dnp.
in contrast, m. scrofulaceum was found in all three wildlife hosts  in cr , so , eb , and pu .
the only mott found in cattle  was isolated from a cow raised in pu.
m. intracellulare was often isolated from wild boar in pu and eb, and also from one fallow deer in eb and two red deer in so and ma, respectively.
tp identification  as in figure  <dig>
factors affecting the presence of m. bovis tps and motts
social groups, defined as animals sampled the same day at the same site, and with characteristics that were compatible with forming a stable  or seasonal  social group were characterized .
in total, 52% of  <dig> social groups were infected with m. bovis.
among red deer,  <dig> groups were identified, m. bovis infecting deer in one third of these groups .
in fallow deer, we identified  <dig> groups with  <dig>  including at least one m. bovis infected individual.
up to  <dig> wild boar groups were identified, and  <dig>  were found infected with m. bovis.
regarding the proportion of social groups with more than one tp, no between host species statistical differences were evidenced.
after eliminating one infected individual per group, red deer  and wild boar chi-square =  <dig> ,  <dig> d.f., p< <dig> ) from infected groups were more probably infected than con-specifics from non infected ones but this did not occur in fallow deer .
reference levels for 'area' and 'host species' are 'so ' and 'wild boar' respectively.
fd = fallow deer, rd = red deer.
cr = coto del rey, eb = estación biológica, ma = marismillas, pu = el puntal.
statistics concerning the glmms to test the factors affecting the presence of a given mycobacterial type or group are shown in table  <dig>  concerning the m. bovis vs mott glmm, the distance to water was statistically higher in mott infected individuals than in m. bovis ones .
the ratio of the minimum distances to similarly infected hosts  statistically interacted with the host.
the ratios  were similar for mott and m. bovis in both deer species , whereas they were higher for m. bovis than mott in the wild boar .
this would indicate that in wild boar the intraspecific spatial aggregation of m. bovis is higher than for mott.
when attending to specific mycobacterial types, there were statistical differences between zones for bovis tp a <dig>  so that it was dominant in wild ungulates from the north of dnp .
there were statistical differences in the probability of infection by m. scrofulaceum relative to other types among host species .
m. scrofulaceum presented a lower intraspecific spatial aggregation than the remaining mycobacterial types .
discussion
this study provided new insights into the ecology of m. bovis and environmental mycobacteria in complex host and pathogen communities, showing that mycobacteria are structured by host species and sampling site, even at very small spatial scales.
the study also showed that host species differences in spatial patterns may greatly depend on behavioral and/or specific host-pathogen-environment interactions, for which our molecular and ecological approach allowed obtaining valuable information on the involved risk factors.
mycobacterial species and typing patterns
contrary to most previous studies in wildlife, where single tps tend to dominate in each geographical region  we detected a high richness of both mott and m. bovis tps in dnp.
whereas single tps are indicative of single introduction events of m. bovis, in our case the high identified tp richness is probably a consequence of  historical cattle breeding and consequent exchanges with breeders from outside the park,  variable conditions provided by high environmental diversity, and  the diversity and abundance of suitable wildlife hosts.
multiple infection of a wildlife host with several m. bovis tps had recently been found in one wild boar from this study area  <cit> .
this observation is rare in wildlife m. bovis hosts  <cit> .
to the best of our knowledge, this is the first study reporting co-infection of red deer and fallow deer with several m. bovis tps.
moreover, the efficiency of isolating mycobacteria could have been improved with the inclusion of liquid media, suggesting that we detected only part of the true co-infections.
the relevance of these findings is that they demonstrate that m. bovis infected wildlife hosts may become infected more than once under natural conditions, at least in areas of high infection pressure such as dnp.
these results also suggest that cross-protection between different m. bovis strains is limited, further underlining the importance of genetic factors rather than immune responses in controlling mycobacterial infections in wildlife  <cit> .
additionally, the infection exclusion reported for closely related genotypes of other intracellular bacteria of the genus anaplasma  <cit>  did not appear to occur for m. bovis tps.
co-existence of members of the m. tuberculosis complex and mott, such as m. intracellulare, had already been reported in human patients  <cit> .
as previously discussed, the fact that we found several m. bovis - mott co-infections suggests that infection by one organism does not impede infection by the other in these wildlife host species.
however, in all three wildlife hosts, isolation of one group of mycobacteria occurred more frequently in individuals not infected by the other group, suggesting that either some competition between mycobacteria or some laboratory bias towards the first identifiable growth may exist.
this is also suggested by the finding that in all cases but in one, m. bovis was isolated from either lymph nodes or tonsils and the mott from the tissue where m. bovis was absent.
in humans, it has been suggested that bcg vaccination protects children against cervical lymph node infection by mott  <cit> .
several authors have reported infection of wild boar with m. scrofulaceum, m. interjectum, m. xenopi and m. intracellulare  <cit> .
all four mott recorded in this study had also already been reported in other wildlife species  <cit> .
however, this is the first report of m. xenopi in deer.
changes over time in dnp
apparently, the community of m. bovis in domestic cattle lost strain richness from time one  to time two , which may result from the application of the official test and slaughter program.
however, the alternative hypothesis of some rare strains going undetected at any sampling period cannot be completely excluded.
part of the new tps isolated from wildlife had been reported in cattle in the earlier survey .
this suggests cases of spill-over from cattle to wild ungulates, and subsequent maintenance of these tps in wildlife reservoir hosts.
other tps had been detected neither in dnp cattle nor in wildlife, but are widespread in spain .
this would suggest a recent introduction, possibly via infected cattle.
however, tp e <dig> is of particular interest.
this tp had never been detected, but is similar to the dominant tp a <dig> except for one spacer.
more sampling and long term studies are needed in order to test whether pathogen evolution resulted in higher tp richness in wildlife species when compared to cattle  <cit> .
spatial structure
our finding that different wildlife species were infected with the same types at a very local scale suggests that transmission is likely to occur between the species.
fallow deer differed from red deer and wild boar in showing more homogeneity in their mycobacterial isolates, regardless of the sampling area.
this may be due to a higher rate of movement of fallow deer between areas and therefore relates to specific territorial and aggregation behaviors as commented above.
this in turn would be relevant for disease control, suggesting a higher capacity of this host for spreading pathogens over long distances.
the different distribution patterns of m. bovis tps may be due to historical introduction of different tps, presumably by infected cattle, in different parts of dnp or, alternatively, if environmental survival of mycobacteria plays a role, to a better adaptation of certain tps to the varying habitat characteristics of northern and southern dnp.
factors affecting the presence of m. bovis tps and motts
in a previous paper we found that infection risk in wild boar was dependent on wild boar m. bovis prevalence in the buffer area containing interacting individuals.
however, this was not evidenced for deer  <cit> .
in this study, using the identification of social groups as an alternative tool, we found that red deer and wild boar were more prone to be infected with m. bovis if they were part of a group with at least one infected deer, while as commented above, this was not the case for fallow deer.
high intraspecific transmission rates at early ages within wild boar social groups have been suggested in wild boar from spain  <cit> , and this probably relates to close interaction when foraging or routing.
animal behavior is an important aspect of disease/host dynamics that as yet has not been well documented but may play an important role in the transmission in free-ranging wildlife populations  <cit> .
owing to higher contact rates and common environmental risk factors, btb transmission should occur more frequently within certain social groups.
recently,  <cit>  used host population genetics to show that contact within family groups probably was a significant mechanism of m. bovis transmission among white-tailed deer  in michigan .
in dnp, modelling suggested that wild boar infection probability depends on wild boar btb prevalence in a buffer zone of interacting individuals, while no such effect was observed in deer  <cit> .
the fallow deer was the only species whose mycobacterial community showed more intra-species similarity throughout dnp than site similarity.
although fallow deer displayed the lowest prevalence , its highly gregarious behavior and subsequent increased transmission risk  may cause mycobacterial strains to be shared by many social groups after social disruption.
this is consistent with the finding that fallow deer displayed the lowest m. bovis prevalence, but a disproportionally high social group prevalence  as compared to that of red deer.
that is, the findings that fallow deer belonging to groups with infected individuals were only rarely infected, or that most infected fallow deer groups had only one infected animal, strongly suggest that either the intra-specific intra-group transmission rate or the susceptibility of fallow deer to btb is lower than in red deer.
however, the alternative explanation that culturing from head lymphoid tissue only missed to detect infection disproportionally more in fallow deer than in red deer or wild boar cannot be excluded.
confirming the above discussed, a spatial structuring in the mycobacterial isolates was evidenced for m. bovis a <dig> type, so that it was dominant in wild ungulates from the north of dnp while b <dig> was dominant in the south .
when we assessed the spatial associations  of m. bovis tps, mott, and m. scrofulaceum, our findings were consistent with spatial aggregation of the host species with the same types.
the spatial distribution of m. bovis versus mott was probably linked with the water, since wildlife hosts infected with m. bovis were most often sampled closer to the marshland than mott.
environmental water sources could act not only as environmental sources of mycobacteria but also by favoring closer contact between the species  <cit> , and this could promote more the transmission of m. bovis by close contact than indirect transmission of mott, which one would expect to be more dependent on external factors.
there were statistical differences in the probability of infection by m. scrofulaceum relative to other types among host species.
m scrofulaceum is a slow-growing atypical mycobacteria that is found in environmental water sources.
nonetheless, no association was evidenced with distance to marshland.
we speculate that the rooting behavior of wild boar may relate to increased exposure to this mycobacteria than other hosts.
nonetheless, our study does not discard that advanced host species-pathogen interactions may also result in different relative occurrences of mycobacterial types across the studied host species.
CONCLUSIONS
the diversity of mycobacteria described herein is indicative of multiple introduction events and a complex multi-host and multi-pathogen epidemiology in dnp.
fine-tuning the epidemiology of mycobacterial infections allowed us to answer a number of relevant questions: first, co-infection of a single host by two m. bovis tps occurred in all three wild ungulate species, confirming that one host can get infected twice.
second, significant changes in the mycobacterial isolate community may have taken place, even in a short time period .
third, we confirmed that red deer and wild boar, but not fallow deer from infected social groups were more probably infected than those from non infected groups.
hence, we agree with the views of several authors suggesting that aspects of host social organization should be taken into account in wildlife epidemiology  <cit> .
fourth, we got insights of spatial structure in mycobacteria distribution, and discussed both habitat-related and host-related explanations for the observed differences.
finally, we conclude that wildlife in dnp is frequently exposed to different species of non-tuberculous, environmental mycobacteria, which could interact with the immune response to pathogenic mycobacteria, although the effects are unknown  <cit> .
in the present study we found evidence of mixed infection, i.e., co-infection of a single host by two m. bovis tps in all three wild ungulate species, and also four deer and four wild boar concurrently presented m. bovis and mott.
the possibility of cross contamination at laboratory or dna level was ruled out.
nonetheless the sensitivity of bacterial culture and dna fingerprinting for the identification of more than one mycobacteria species or m. tuberculosis complex strain may be limited when the strains are not present in the particular cultured organ/tissue.
in the context of the endemicity of tuberculosis and the possibility of repeated exposure to newer infections, one can expect more cases with multiple strain infection  <cit> .
therefore more research concerning whether infection with one strain would protect against infection with another strain is needed.
molecular typing did not allow inferring the direction of transmission  <cit> .
however, findings of rare tps such as e <dig> among both fallow deer and wild boar strongly suggest that interspecies transmission and/or common sources of infection do occur among wild ungulates.
conversely, the lack of isolation of rare m. bovis spoligotype patterns from cattle of the 2006- <dig> sample suggests that spill-back from the wildlife reservoir to livestock may not be a very usual event.
the results highlight the suitability of molecular typing for surveys at small spatial and temporal scales.
however, increased surveillance along with a better understanding of the transmission routes, environmental persistence, and associated risk factors  are needed if we are to effectively control bovine tb in dnp.
one remaining question relates to the influence of the genotype of mycobacteria on the virulence  <cit> , which may be mediated by secondary infections, which should be addressed by future research.
competing interests
the authors declare that they have no competing interests.
authors' contributions
conceived and designed the study: cg, mt, jn, jf, jv.
participated in sampling and field work: cg, mt, jn, jv.
carried out the laboratory work: mt, ja analyzed the data: cg, jn, pa, jv.
draft the manuscript: cg, mt, jn, pa, jf, jv.
all authors read and approved the final manuscript.

 identifying mechanisms of pathogen transmission, including potential environmental sources, is critical to control disease. molecular epidemiology integrates conventional epidemiological approaches with molecular techniques to track specific strains of pathogens in order to understand the distribution of pathogens in populations and environments. this can be used to elucidate inter- and intra-specific transmission pathways and environmental risk factors, from individual to population, and from local to broader spatial scales. the genus mycobacterium comprises over species and several subspecies. over of these can cause disease in livestock, wildlife and humans, occurring worldwide. mycobacterial diseases such as bovine tuberculosis have become a major sanitary and conservation problem even in relatively unmanaged natural areas across the world. similar to other shared diseases, the existence of wildlife reservoirs is limiting the effectiveness of eradication schemes in livestock. in btb, known risk factors for wild ungulates include age, gender, density, spatial aggregation, intra and inter-specific contact, fencing and other habitat features as well as genetic factors. however, most data derive from large scale studies, while detailed information at small spatial scales is still very scarce, and usually, fine associations with spatial and environmental factors are not addressed. even less information is available regarding the effects on wildlife and livestock of mycobacteria other than tuberculosis. environmental mycobacteria or mott include a large number of species that can cause serious illnesses in humans, particularly in immunocompromised patients. for example, mycobacterium interjectum has been identified as a causative agent of cervical lymphadenitis in children, and of cutaneous infections in immunosuppressed patients. m. xenopi may cause pulmonary disease in humans, and m. scrofulaceum may cause cutaneous infections and lymphadenitis. in humans, risk factors for mott infections include immunosuppression, contaminated water and aerosol exposure, and short or old age. mott are widely distributed in the environment, particularly in wet soil, marshland, streams, rivers and estuaries, but each species shows different preferences. because of its habitat characteristics, extension and their sizeable wild and domestic animal populations, doñana national park in southern spain has been proposed as a good natural laboratory for studying wildlife mycobacteriosis. molecular typing of m. bovis isolates for the period 1998- showed that wildlife species in dnp were infected only with those m. bovis typing patterns that were more prevalent in local cattle. furthermore, the results were suggestive of micro-evolutionary events in the local m. bovis population. in the same period, m. bovis infection prevalence in dnp was 33% in european wild boar, 21% in red deer, and 26% in fallow deer. in a more recent study, we confirmed infection with m. bovis in 52% wild boar, 27% red deer and 18% fallow deer from dnp in 2006- and evidenced that m. bovis prevalence decreased from north to south in wild boar and red deer, whereas no clear spatial pattern was observed for fallow deer. three wild ungulates coexist in dnp, wild boar, fallow deer and red deer, along with domestic cattle subjected to btb eradication programs. we included the wild species as our study models as all are highly susceptible to btb and are known to show high prevalence in the area. in addition, their different ecology and behavior peculiarities can play a role in the epidemiology of mycobacteria, for example, variations in sociability or gregariousness, and scavenging habits. in addition, different habitats could provide variable environmental suitability for m. bovis persistence. in this sense, scrublands and woodlands are preferably used by red deer and wild boar compared with fallow deer. in this study we used molecular epidemiological techniques to establish the extent of m. bovis strain richness and other environmental mycobacterial species in isolates collected in wildlife and cattle from the dnp, so as the association with social, spatial and environmental factors in this multi-host and multi-pathogen situation. we hypothesized that infection by mycobacteria would differ among hosts and sites depending on host, pathogen and local environment ecology. fine-tuning mycobacterial epidemiology in dnp allowed rising a number of relevant questions: do hosts get infected twice by m. bovis and mott, and can this interfere in m. bovis infection or vice versa? have new m. bovis types appeared or have any changes in type composition taken place in recent years? is there an effect of the social group on infection risk? is there a spatial structure in mycobacteria distribution? are there species-specific variants of mycobacteria that could be attributed to species-specific behavior patterns and/or to advanced host species-pathogen interactions? methods study area the study was carried out in dnp, located in south-western spain and covering ha. this is a flat region of sandy soils bordering the atlantic ocean, with a maximum elevation of m. the climate is mediterranean sub-humid with marked seasons. in the wet season, most of the marshlands are flooded and wildlife and cattle tend to graze in the more elevated scrublands. in summer, the wetter and more productive ecotone between the scrublands and the marshes supports aggregations of wild and domestic ungulates. human access is restricted and management is carried out by park authorities. limited traditional exploitation of some natural resources, such as logging, and cattle and horse rising are allowed. after when btb in wildlife was first diagnosed in dnp a government-sponsored program was initiated to eradicate btb-positive cattle. ungulate populations have been culled by shooting. animal sampling from april to april european wild boar, red deer, and fallow deer were sampled within the park by shooting. the culling of wild ungulates was approved by the research commission of doñana national park in accordance with management rules established by the autonomous government of andalucía. for each animal we recorded the exact position with gps. sex and age, based on tooth eruption patterns, were recorded in the field. a necropsy was performed on site and the presence of tuberculosis-like lesions recorded by macroscopic inspection of lymph nodes and abdominal and thoracic organs. this protocol included the examination of the lungs for the presence of tb-compatible macroscopic lesions during field inspection and a sample was collected. a tonsil and a head lymph node sample from each individual were collected for culture. in wild boar, one piece of the tonsils and one from both mandibular lymph nodes were submitted for microbiological studies. in red deer and fallow deer, one piece of the tonsils and head lymphnode samples, always containing at least half left and half right medial retropharyngeal lymph node, were submitted for culture. due to logistic and budget constraints, no thoracic or abdominal lymphoid tissues were cultured except when tb-compatible macroscopic lesions were evidenced. m. scr. = mycobacterium scrofulaceum; m. int. = mycobacterium interjectum, m. xen. = mycobacterium xenopi, m. int. = mycobacterium intracellulare social groups were defined as animals sampled the same day at the same site, and with characteristics that were compatible with forming a stable or seasonal group. only part of the individuals belonging to a given social group was sampled. sampling was performed according to european and spanish laws, and current guidelines for ethical use of animals in research and uclm animal experimentation committee. microbiological procedures culture: the specimens were liquefied and decontaminated with an equal volume of n-acetyl-l-cysteine and 2% naoh, mixed by centrifugal swirling, and then incubated for min. the reaction was neutralized by adding m phosphate-buffered saline, to a final volume of ml. the specimens were concentrated by centrifugation at × g for min. the supernatant was discarded, and the sediment was re-suspended in ml of sterile water. the sediment was used to inoculate two löwestein-jensen with pyruvate solid medium. lowëstein-jenssen slants were incubated at 37°c for weeks and inspected weekly for growth. when growth was detected, a smear was prepared to confirm the presence of acid-fast bacilli from suspect colonies by ziehl-neelsen staining. identification we identified m. bovis and mott to the species level and characterized m. bovis strains with spoligotyping and miru-vntr typing. macroscopic morphology of the colonies and pigment production was recorded. identification at species level was performed with the genotype®mtbc for the mycobacterium complex strains that allows the differentiation of m. africanum i, m. bovis bcg, m. bovis ssp. bovis, m. bovis ssp. caprae and m. tuberculosis/m. africanum ii/m. canettii. mott strains were identified by the genotype® mycobacterium cm and genotype® mycobacterium as mtbc. the genotype assays were performed according to the manufacturer's instructions: dna extraction by the dna sss method was followed by pcr amplification of a trait of the 23s rrna gene, as recommended. reverse hybridization and detection were carried out on a shaking water bath. the final identification was obtained by comparison of line probe patterns with the provided evaluation sheet. typing the m. bovis isolates were further characterized by spoligotyping. the amplified product was detected by hybridization of the biotin-labelled pcr product onto spoligotyping membrane. purified sterile water and chromosomal dna of m. tuberculosis h37rv and m. bovis bcg p were included as controls in each batch of tests. the patterns were allocated a number in the m. bovis spoligotyping database. the results were recorded in sb code, followed by a field of digits as defined on the m. bovis spoligotype database website. all wildlife isolates were also subjected to miru-vntr analysis. extensive documentation on the service and the genotyping methods is available at the miru-vntrplus website. the analysis was carried out as described, and primers previously reported were used to amplify the etr-a and etr-b loci, miru and using the recommended pcr conditions. amplicon sizes were estimated by electrophoresis on a % agarose gel at v during h, using 100-bp ladder. figure presents the spoligotyping patterns, vntr allelic profiles and typing pattern codes defined for this study. statistics chi-square tests were used for between-pair comparisons of prevalences. to test for the effect of host species vs site regarding the mycobacterial isolates, we used the czechanovsky similarity index. this index considers the list of mycobacterial species recorded in a given host type or in a given study area. it is calculated by dividing two times the species shared between two lists, by the total number of species of both lists, as follows: cz=2*/ considering the animals in which any mycobacterial infection was diagnosed, three generalized linear mixed models were explored to test different explanatory variables that affect the presence of a mycobacterial type or group. the most common mycobacterial groups were: m. bovis m bovis a and m. scrofulaceum. the presence or absence of infection in a mycobacterial group was considered as a binary variable. the model was fitted using a logit link function. the model considered social group as a random effect. the model included host species, the study area and age as categorical explanatory variables. the distance to the water was included as a continuous predictor. to compare the spatial associations of infection by specific mycobacterial type and hosts, we included as explanatory continuous variable the ratio between the nearest neighbor distance from host to a different host species with the same type of mycobacteria relative to the nearest distance to a con-specific host with the same type of mycobacteria. a ratio > indicates that the nearest distance to a host with the same spoligotype is higher for a different host species. all the aforementioned explanatory variables we also included in the models interacting with the host species. due to over-parameterization of the models and zero inflated data, no interactions were included in the m. bovis a and m. scrofulaceum models. p-value was set as ≤. we estimated exact confidence limits for prevalence using sterne's exact method.  mycobacteria species and molecular types we obtained a total of mycobacterial isolates from dnp wildlife. this included m. bovis isolates belonging to different typing patterns, and isolates belonging to four mott. m. bovis tps and mott species were isolated from wild boar, red deer, and fallow deer. wild boar and red deer had m. bovis tps each, whereas fallow deer presented only tps. the number of different isolates per host was in wild boar, in red deer and in fallow deer. regarding m. bovis, we identified different spoligotyping patterns and different vntr allelic profiles. one spoligotyping pattern was new according to the m. bovis database, and was therefore introduced with code sb co-infection of a single host by two m. bovis tps occurred in all three wild ungulate species. one adult male red deer was infected with tps a and b one adult male and one adult female fallow deer were co-infected with tps a and e and two wild boar were co-infected with tps a and b mott species found in wildlife hosts included m. scrofulaceum and m. intracellulare, both found in all host species, m. interjectum, and m. xenopi. in four deer and four wild boar, m. bovis and mott were isolated concurrently. in a single wild boar, both types of mycobacteria were simultaneously isolated from the two tissue samples collected and cultured, while in the remaining cases m. bovis was isolated from either lymph nodes or tonsils and the mott from the tissue where m. bovis was absent. we recorded no cases of co-infection by different mott. the proportion of wild ungulate m. bovis typing patterns other than the dominant a and b did not differ statistically from 1998- to 2006-. no spoligotyping patterns other than the two dominant ones were detected among cattle isolates in and changes in mycobacterial typing patterns over time in dnp all three m. bovis typing patterns recorded in dnp wildlife between and were still evidenced in similar proportions in 2006-. however, while only three different tps had been detected in dnp wildlife in the first period, up to different ones were found in the second period. two of these "new" tps had already been recorded in cattle sampled in dnp between and however, other tps had never before been reported in dnp. in contrast with the situation in wildlife and to data from 1998- when out of cattle spoligotyping patterns were different from a and b, no spoligotyping patterns other than the two dominant ones were detected among cattle isolates in and. spatial structure regarding the mott, m. interjectum was only found in wild boar from eb, in the central part of dnp. in contrast, m. scrofulaceum was found in all three wildlife hosts in cr, so, eb, and pu. the only mott found in cattle was isolated from a cow raised in pu. m. intracellulare was often isolated from wild boar in pu and eb, and also from one fallow deer in eb and two red deer in so and ma, respectively. tp identification as in figure factors affecting the presence of m. bovis tps and motts social groups, defined as animals sampled the same day at the same site, and with characteristics that were compatible with forming a stable or seasonal social group were characterized. in total, 52% of social groups were infected with m. bovis. among red deer, groups were identified, m. bovis infecting deer in one third of these groups. in fallow deer, we identified groups with including at least one m. bovis infected individual. up to wild boar groups were identified, and were found infected with m. bovis. regarding the proportion of social groups with more than one tp, no between host species statistical differences were evidenced. after eliminating one infected individual per group, red deer and wild boar chi-square =, d.f., p< ) from infected groups were more probably infected than con-specifics from non infected ones but this did not occur in fallow deer. reference levels for'area' and'host species' are'so' and'wild boar' respectively. fd = fallow deer, rd = red deer. cr = coto del rey, eb = estación biológica, ma = marismillas, pu = el puntal. statistics concerning the glmms to test the factors affecting the presence of a given mycobacterial type or group are shown in table concerning the m. bovis vs mott glmm, the distance to water was statistically higher in mott infected individuals than in m. bovis ones. the ratio of the minimum distances to similarly infected hosts statistically interacted with the host. the ratios were similar for mott and m. bovis in both deer species, whereas they were higher for m. bovis than mott in the wild boar. this would indicate that in wild boar the intraspecific spatial aggregation of m. bovis is higher than for mott. when attending to specific mycobacterial types, there were statistical differences between zones for bovis tp a so that it was dominant in wild ungulates from the north of dnp. there were statistical differences in the probability of infection by m. scrofulaceum relative to other types among host species. m. scrofulaceum presented a lower intraspecific spatial aggregation than the remaining mycobacterial types. discussion this study provided new insights into the ecology of m. bovis and environmental mycobacteria in complex host and pathogen communities, showing that mycobacteria are structured by host species and sampling site, even at very small spatial scales. the study also showed that host species differences in spatial patterns may greatly depend on behavioral and/or specific host-pathogen-environment interactions, for which our molecular and ecological approach allowed obtaining valuable information on the involved risk factors. mycobacterial species and typing patterns contrary to most previous studies in wildlife, where single tps tend to dominate in each geographical region we detected a high richness of both mott and m. bovis tps in dnp. whereas single tps are indicative of single introduction events of m. bovis, in our case the high identified tp richness is probably a consequence of historical cattle breeding and consequent exchanges with breeders from outside the park, variable conditions provided by high environmental diversity, and the diversity and abundance of suitable wildlife hosts. multiple infection of a wildlife host with several m. bovis tps had recently been found in one wild boar from this study area. this observation is rare in wildlife m. bovis hosts. to the best of our knowledge, this is the first study reporting co-infection of red deer and fallow deer with several m. bovis tps. moreover, the efficiency of isolating mycobacteria could have been improved with the inclusion of liquid media, suggesting that we detected only part of the true co-infections. the relevance of these findings is that they demonstrate that m. bovis infected wildlife hosts may become infected more than once under natural conditions, at least in areas of high infection pressure such as dnp. these results also suggest that cross-protection between different m. bovis strains is limited, further underlining the importance of genetic factors rather than immune responses in controlling mycobacterial infections in wildlife. additionally, the infection exclusion reported for closely related genotypes of other intracellular bacteria of the genus anaplasma did not appear to occur for m. bovis tps. co-existence of members of the m. tuberculosis complex and mott, such as m. intracellulare, had already been reported in human patients. as previously discussed, the fact that we found several m. bovis - mott co-infections suggests that infection by one organism does not impede infection by the other in these wildlife host species. however, in all three wildlife hosts, isolation of one group of mycobacteria occurred more frequently in individuals not infected by the other group, suggesting that either some competition between mycobacteria or some laboratory bias towards the first identifiable growth may exist. this is also suggested by the finding that in all cases but in one, m. bovis was isolated from either lymph nodes or tonsils and the mott from the tissue where m. bovis was absent. in humans, it has been suggested that bcg vaccination protects children against cervical lymph node infection by mott. several authors have reported infection of wild boar with m. scrofulaceum, m. interjectum, m. xenopi and m. intracellulare. all four mott recorded in this study had also already been reported in other wildlife species. however, this is the first report of m. xenopi in deer. changes over time in dnp apparently, the community of m. bovis in domestic cattle lost strain richness from time one to time two, which may result from the application of the official test and slaughter program. however, the alternative hypothesis of some rare strains going undetected at any sampling period cannot be completely excluded. part of the new tps isolated from wildlife had been reported in cattle in the earlier survey. this suggests cases of spill-over from cattle to wild ungulates, and subsequent maintenance of these tps in wildlife reservoir hosts. other tps had been detected neither in dnp cattle nor in wildlife, but are widespread in spain. this would suggest a recent introduction, possibly via infected cattle. however, tp e is of particular interest. this tp had never been detected, but is similar to the dominant tp a except for one spacer. more sampling and long term studies are needed in order to test whether pathogen evolution resulted in higher tp richness in wildlife species when compared to cattle. spatial structure our finding that different wildlife species were infected with the same types at a very local scale suggests that transmission is likely to occur between the species. fallow deer differed from red deer and wild boar in showing more homogeneity in their mycobacterial isolates, regardless of the sampling area. this may be due to a higher rate of movement of fallow deer between areas and therefore relates to specific territorial and aggregation behaviors as commented above. this in turn would be relevant for disease control, suggesting a higher capacity of this host for spreading pathogens over long distances. the different distribution patterns of m. bovis tps may be due to historical introduction of different tps, presumably by infected cattle, in different parts of dnp or, alternatively, if environmental survival of mycobacteria plays a role, to a better adaptation of certain tps to the varying habitat characteristics of northern and southern dnp. factors affecting the presence of m. bovis tps and motts in a previous paper we found that infection risk in wild boar was dependent on wild boar m. bovis prevalence in the buffer area containing interacting individuals. however, this was not evidenced for deer. in this study, using the identification of social groups as an alternative tool, we found that red deer and wild boar were more prone to be infected with m. bovis if they were part of a group with at least one infected deer, while as commented above, this was not the case for fallow deer. high intraspecific transmission rates at early ages within wild boar social groups have been suggested in wild boar from spain, and this probably relates to close interaction when foraging or routing. animal behavior is an important aspect of disease/host dynamics that as yet has not been well documented but may play an important role in the transmission in free-ranging wildlife populations. owing to higher contact rates and common environmental risk factors, btb transmission should occur more frequently within certain social groups. recently, used host population genetics to show that contact within family groups probably was a significant mechanism of m. bovis transmission among white-tailed deer in michigan. in dnp, modelling suggested that wild boar infection probability depends on wild boar btb prevalence in a buffer zone of interacting individuals, while no such effect was observed in deer. the fallow deer was the only species whose mycobacterial community showed more intra-species similarity throughout dnp than site similarity. although fallow deer displayed the lowest prevalence, its highly gregarious behavior and subsequent increased transmission risk may cause mycobacterial strains to be shared by many social groups after social disruption. this is consistent with the finding that fallow deer displayed the lowest m. bovis prevalence, but a disproportionally high social group prevalence as compared to that of red deer. that is, the findings that fallow deer belonging to groups with infected individuals were only rarely infected, or that most infected fallow deer groups had only one infected animal, strongly suggest that either the intra-specific intra-group transmission rate or the susceptibility of fallow deer to btb is lower than in red deer. however, the alternative explanation that culturing from head lymphoid tissue only missed to detect infection disproportionally more in fallow deer than in red deer or wild boar cannot be excluded. confirming the above discussed, a spatial structuring in the mycobacterial isolates was evidenced for m. bovis a type, so that it was dominant in wild ungulates from the north of dnp while b was dominant in the south. when we assessed the spatial associations of m. bovis tps, mott, and m. scrofulaceum, our findings were consistent with spatial aggregation of the host species with the same types. the spatial distribution of m. bovis versus mott was probably linked with the water, since wildlife hosts infected with m. bovis were most often sampled closer to the marshland than mott. environmental water sources could act not only as environmental sources of mycobacteria but also by favoring closer contact between the species, and this could promote more the transmission of m. bovis by close contact than indirect transmission of mott, which one would expect to be more dependent on external factors. there were statistical differences in the probability of infection by m. scrofulaceum relative to other types among host species. m scrofulaceum is a slow-growing atypical mycobacteria that is found in environmental water sources. nonetheless, no association was evidenced with distance to marshland. we speculate that the rooting behavior of wild boar may relate to increased exposure to this mycobacteria than other hosts. nonetheless, our study does not discard that advanced host species-pathogen interactions may also result in different relative occurrences of mycobacterial types across the studied host species.  the diversity of mycobacteria described herein is indicative of multiple introduction events and a complex multi-host and multi-pathogen epidemiology in dnp. fine-tuning the epidemiology of mycobacterial infections allowed us to answer a number of relevant questions: first, co-infection of a single host by two m. bovis tps occurred in all three wild ungulate species, confirming that one host can get infected twice. second, significant changes in the mycobacterial isolate community may have taken place, even in a short time period. third, we confirmed that red deer and wild boar, but not fallow deer from infected social groups were more probably infected than those from non infected groups. hence, we agree with the views of several authors suggesting that aspects of host social organization should be taken into account in wildlife epidemiology. fourth, we got insights of spatial structure in mycobacteria distribution, and discussed both habitat-related and host-related explanations for the observed differences. finally, we conclude that wildlife in dnp is frequently exposed to different species of non-tuberculous, environmental mycobacteria, which could interact with the immune response to pathogenic mycobacteria, although the effects are unknown. in the present study we found evidence of mixed infection, i.e., co-infection of a single host by two m. bovis tps in all three wild ungulate species, and also four deer and four wild boar concurrently presented m. bovis and mott. the possibility of cross contamination at laboratory or dna level was ruled out. nonetheless the sensitivity of bacterial culture and dna fingerprinting for the identification of more than one mycobacteria species or m. tuberculosis complex strain may be limited when the strains are not present in the particular cultured organ/tissue. in the context of the endemicity of tuberculosis and the possibility of repeated exposure to newer infections, one can expect more cases with multiple strain infection. therefore more research concerning whether infection with one strain would protect against infection with another strain is needed. molecular typing did not allow inferring the direction of transmission. however, findings of rare tps such as e among both fallow deer and wild boar strongly suggest that interspecies transmission and/or common sources of infection do occur among wild ungulates. conversely, the lack of isolation of rare m. bovis spoligotype patterns from cattle of the 2006- sample suggests that spill-back from the wildlife reservoir to livestock may not be a very usual event. the results highlight the suitability of molecular typing for surveys at small spatial and temporal scales. however, increased surveillance along with a better understanding of the transmission routes, environmental persistence, and associated risk factors are needed if we are to effectively control bovine tb in dnp. one remaining question relates to the influence of the genotype of mycobacteria on the virulence, which may be mediated by secondary infections, which should be addressed by future research. competing interests the authors declare that they have no competing interests. authors' contributions conceived and designed the study: cg, mt, jn, jf, jv. participated in sampling and field work: cg, mt, jn, jv. carried out the laboratory work: mt, ja analyzed the data: cg, jn, pa, jv. draft the manuscript: cg, mt, jn, pa, jf, jv. all authors read and approved the final manuscript. 
we describe the diversity of two kinds of mycobacteria isolates, environmental mycobacteria and mycobacterium bovis collected from wild boar, fallow deer, red deer and cattle in doñana national park, analyzing their association with temporal, spatial and environmental factors. high diversity of environmental mycobacteria species and m. bovis typing patterns were found. when assessing the factors underlying the presence of the most common types of both environmental mycobacteria and m. bovis tps in dnp, we evidenced host species differences in the occurrence, spatial structuration and differences in the degree of spatial association of specific types between host species. wildlife in dnp is frequently exposed to different species of non-tuberculous, environmental mycobacteria, which could interact with the immune response to pathogenic mycobacteria, although the effects are unknown. co-infection of a single host by two m. bovis tps occurred in all three wild ungulate species. while only three tps were detected in wildlife between and up to different ones were found during 2006- the opposite was observed in cattle. the diversity of mycobacteria described herein is indicative of multiple introduction events and a complex multi-host and multi-pathogen epidemiology in dnp. belonging to an m. bovis-infected social group was a significant risk factor for mycobacterial infection in red deer and wild boar, but not for fallow deer. m. bovis tps were usually found closer to water marshland than mott. in wild boar and red deer, isolation of one group of mycobacteria occurred more frequently in individuals not infected by the other group. this research highlights the suitability of molecular typing for surveys at small spatial and temporal scales. significant changes in the mycobacterial isolate community may have taken place, even in a short time period. aspects of host social organization should be taken into account in wildlife epidemiology.
 we describe the diversity of two kinds of mycobacteria isolates, environmental mycobacteria and mycobacterium bovis collected from wild boar, fallow deer, red deer and cattle in doñana national park, analyzing their association with temporal, spatial and environmental factors.  high diversity of environmental mycobacteria species and m. bovis typing patterns were found. when assessing the factors underlying the presence of the most common types of both environmental mycobacteria and m. bovis tps in dnp, we evidenced host species differences in the occurrence, spatial structuration and differences in the degree of spatial association of specific types between host species. co-infection of a single host by two m. bovis tps occurred in all three wild ungulate species. in wild boar and red deer, isolation of one group of mycobacteria occurred more frequently in individuals not infected by the other group. while only three tps were detected in wildlife between and up to different ones were found during 2006- the opposite was observed in cattle. belonging to an m. bovis-infected social group was a significant risk factor for mycobacterial infection in red deer and wild boar, but not for fallow deer. m. bovis tps were usually found closer to water marshland than mott.  the diversity of mycobacteria described herein is indicative of multiple introduction events and a complex multi-host and multi-pathogen epidemiology in dnp. significant changes in the mycobacterial isolate community may have taken place, even in a short time period. aspects of host social organization should be taken into account in wildlife epidemiology. wildlife in dnp is frequently exposed to different species of non-tuberculous, environmental mycobacteria, which could interact with the immune response to pathogenic mycobacteria, although the effects are unknown. this research highlights the suitability of molecular typing for surveys at small spatial and temporal scales. 