Cluster analysis of behavioural and event-related potentials during a contingent negative variation paradigm in remitting-relapsing and benign forms of multiple sclerosis

Background: Event-related potentials (ERPs) may be used as a highly sensitive way of detecting subtle degrees of cognitive dysfunction. On the other hand, impairment of cognitive skills is increasingly recognised as a hallmark of patients suffering from multiple sclerosis (MS). We sought to determine the psychophysiological pattern of information processing among MS patients with the relapsing-remitting form of the disease and low physical disability considered as two subtypes: 'typical relapsing-remitting' (RRMS) and 'benign MS' (BMS). Furthermore, we subjected our data to a cluster analysis to determine whether MS patients and healthy controls could be differentiated in terms of their psychophysiological profile.Methods: We investigated MS patients with RRMS and BMS subtypes using event-related potentials (ERPs) acquired in the context of a Posner visual-spatial cueing paradigm. Specifically, our study aimed to assess ERP brain activity in response preparation (contingent negative variation -CNV) and stimuli processing in MS patients. Latency and amplitude of different ERP components (P1, eN1, N1, P2, N2, P3 and late negativity -LN) as well as behavioural responses (reaction time -RT; correct responses -CRs; and number of errors) were analyzed and then subjected to cluster analysis. Results: Both MS groups showed delayed behavioural responses and enhanced latency for long-latency ERP components (P2, N2, P3) as well as relatively preserved ERP amplitude, but BMS patients obtained more important performance deficits (lower CRs and higher RTs) and abnormalities related to the latency (N1, P3) and amplitude of ERPs (eCNV, eN1, LN). However, RRMS patients also demonstrated abnormally high amplitudes related to the preparation performance period of CNV (cCNV) and post-processing phase (LN). Cluster analyses revealed that RRMS patients appear to make up a relatively homogeneous group with moderate deficits mainly related to ERP latencies, whereas BMS patients appear to make up a rather more heterogeneous group with more severe information processing and attentional deficits. Conclusions: Our findings are suggestive of a slowing of information processing for MS patients that may be a consequence of demyelination and axonal degeneration, which also seems to occur in MS patients that show little or no progression in the physical severity of the disease over time.

Genetic variability of the mTOR pathway and prostate cancer risk in the European Prospective Investigation on Cancer (EPIC)

The mTOR (mammalian target of rapamycin) signal transduction pathway integrates various signals, regulating ribosome biogenesis and protein synthesis as a function of available energy and amino acids, and assuring an appropriate coupling of cellular proliferation with increases in cell size. In addition, recent evidence has pointed to an interplay between the mTOR and p53 pathways. We investigated the genetic variability of 67 key genes in the mTOR pathway and in genes of the p53 pathway which interact with mTOR. We tested the association of 1,084 tagging SNPs with prostate cancer risk in a study of 815 prostate cancer cases and 1,266 controls nested within the European Prospective Investigation into Cancer and Nutrition (EPIC). We chose the SNPs (n = 11) with the strongest association with risk (p<0.01) and sought to replicate their association in an additional series of 838 prostate cancer cases and 943 controls from EPIC. In the joint analysis of first and second phase two SNPs of the PRKCI gene showed an association with risk of prostate cancer (ORallele = 0.85, 95% CI 0.78–0.94, p = 1.3×10−3 for rs546950 and ORallele = 0.84, 95% CI 0.76–0.93, p = 5.6×10−4 for rs4955720). We confirmed this in a meta-analysis using as replication set the data from the second phase of our study jointly with the first phase of the Cancer Genetic Markers of Susceptibility (CGEMS) project. In conclusion, we found an association with prostate cancer risk for two SNPs belonging to PRKCI, a gene which is frequently overexpressed in various neoplasms, including prostate cancer.

Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000-2010.

Most of the non-B HIV-1 subtypes are predominant in Sub-Saharan Africa and India although they have been found worldwide. In the last decade, immigration from these areas has increased considerably in Spain. The objective of this study was to evaluate the prevalence of non-B subtypes circulating in a cohort of HIV-1-infected immigrants in Seville, Southern Spain and to identify drug resistance-associated mutations. METHODS: Complete protease and first 220 codons of the reverse transcriptase coding regions were amplified and sequenced by population sequencing. HIV-1 subtypes were determined using Stanford University Drug Resistance Database, and phylogenetic analysis was performed comparing multiple reported sequences. Drug resistance mutations were defined according to the International AIDS Society-USA. RESULTS: From 2000 to 2010 a total of 1,089 newly diagnosed HIV-1-infected patients were enrolled in our cohort. Of these, 121 were immigrants, of which 98 had ethical approval and informed consent to include in our study. Twenty-nine immigrants (29/98, 29.6%) were infected with non-B subtypes, of which 15/29 (51.7%) were CRF02-AG, mostly from Sub-Saharan Africa, and 2/29 (6.9%) were CRF01-AE from Eastern Europe. A, C, F, J and G subtypes from Eastern Europe, Central-South America and Sub-Saharan Africa were also present. Some others harboured recombinant forms CRF02-AG/CRF01-AE, CRF2-AG/G and F/B, B/C, and K/G, in PR and RT-coding regions. Patients infected with non-B subtypes showed a high frequency of minor protease inhibitor resistance mutations, M36I, L63P, and K20R/I. Only one patient, CRF02_AG, showed major resistance mutation L90M. Major RT inhibitor resistance mutations K70R and A98G were present in one patient with subtype G, L100I in one patient with CRF01_AE, and K103N in another patient with CRF01_AE. Three patients had other mutations such as V118I, E138A and V90I. CONCLUSIONS: The circulation of non-B subtypes has significantly increased in Southern Spain during the last decade, with 29.6% prevalence, in association with demographic changes among immigrants. This could be an issue in the treatment and management of these patients. Resistance mutations have been detected in these patients with a prevalence of 7% among treatment-naïve patients compared with the 21% detected among patients under HAART or during treatment interruption.

ELOVL6 Genetic Variation Is Related to Insulin Sensitivity: A New Candidate Gene in Energy Metabolism

BACKGROUND: The elongase of long chain fatty acids family 6 (ELOVL6) is an enzyme that specifically catalyzes the elongation of saturated and monounsaturated fatty acids with 12, 14 and 16 carbons. ELOVL6 is expressed in lipogenic tissues and it is regulated by sterol regulatory element binding protein 1 (SREBP-1). OBJECTIVE: We investigated whether ELOVL6 genetic variation is associated with insulin sensitivity in a population from southern Spain. DESIGN: We undertook a prospective, population-based study collecting phenotypic, metabolic, nutritional and genetic information. Measurements were made of weight and height and the body mass index (BMI) was calculated. Insulin resistance was measured by homeostasis model assessment. The type of dietary fat was assessed from samples of cooking oil taken from the participants' kitchens and analyzed by gas chromatography. Five SNPs of the ELOVL6 gene were analyzed by SNPlex. RESULTS: Carriers of the minor alleles of the SNPs rs9997926 and rs6824447 had a lower risk of having high HOMA_IR, whereas carriers of the minor allele rs17041272 had a higher risk of being insulin resistant. An interaction was detected between the rs6824447 polymorphism and the intake of oil in relation with insulin resistance, such that carriers of this minor allele who consumed sunflower oil had lower HOMA_IR than those who did not have this allele (P = 0.001). CONCLUSIONS: Genetic variations in the ELOVL6 gene were associated with insulin sensitivity in this population-based study.

An antibody-based ELISA for quantification of Ani s 1, a major allergen from Anisakis simplex.

Anisakis simplex is a nematode parasite that can infect humans who have eaten raw or undercooked seafood. Larvae invading the gastrointestinal mucosa excrete/secrete proteins that are implicated in the pathogenesis of anisakiasis and can induce IgE-mediated symptoms. Since Ani s 1 is a potent secreted allergen with important clinical relevance, its measurement could assess the quality of allergenic products used in diagnosis/immunotherapy of Anisakis allergy and track the presence of A. simplex parasites in fish foodstuffs. An antibody-based ELISA for quantification of Ani s 1 has been developed based on monoclonal antibody 4F2 as capture antibody and biotin-labelled polyclonal antibodies against Ani s 1 as detection reagent. The dose-response standard curves, obtained with natural and recombinant antigens, ranged from 4 to 2000 ng/ml and were identical and parallel to that of the A. simplex extract. The linear portion of the dose-response curve with nAni s 1 was between 15 and 250 ng/ml with inter-assay and intra-assays coefficients of variation less than 20% and 10%, respectively. The assay was specific since there was no cross-reaction with other extracts (except Ascaris extracts) and was highly sensitive (detection limit of 1·8 ng/ml), being able to detect Ani s 1 in fish extracts from codfish and monkfish.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

CD209 in inflammatory bowel disease: a case-control study in the Spanish population

BACKGROUND The etiology of Ulcerative Colitis (UC) and Crohn's Disease (CD), considered together as Inflammatory Bowel Diseases (IBD), involves environmental and genetic factors. Although some genes are already known, the genetics underlying these diseases is complex and new candidates are continuously emerging. The CD209 gene is located in a region linked previously to IBD and a CD209 functional polymorphism (rs4804803) has been associated to other inflammatory conditions. Our aim was to study the potential involvement of this CD209 variant in IBD susceptibility. METHODS We performed a case-control study with 515 CD patients, 497 UC patients and 731 healthy controls, all of them white Spaniards. Samples were typed for the CD209 single nucleotide polymorphism (SNP) rs4804803 by TaqMan technology. Frequency comparisons were performed using chi2 tests. RESULTS No association between CD209 and UC or CD was observed initially. However, stratification of UC patients by HLA-DR3 status, a strong protective allele, showed that carriage of the CD209_G allele could increase susceptibility in the subgroup of HLA-DR3-positive individuals (p = 0.03 OR = 1.77 95% CI 1.04-3.02, vs. controls). CONCLUSION A functional variant in the CD209 gene, rs4804803, does not seem to be influencing Crohn's disease susceptibility. However, it could be involved in the etiology or pathology of Ulcerative Colitis in HLA-DR3-positive individuals but further studies are necessary.

Genetic polymorphisms of RANTES, IL1-A, MCP-1 and TNF-A genes in patients with prostate cancer

BACKGROUND Inflammation has been implicated as an etiological factor in several human cancers, including prostate cancer. Allelic variants of the genes involved in inflammatory pathways are logical candidates as genetic determinants of prostate cancer risk. The purpose of this study was to investigate whether single nucleotide polymorphisms of genes that lead to increased levels of pro-inflammatory cytokines and chemokines are associated with an increased prostate cancer risk. METHODS A case-control study design was used to test the association between prostate cancer risk and the polymorphisms TNF-A-308 A/G (rs 1800629), RANTES-403 G/A (rs 2107538), IL1-A-889 C/T (rs 1800587) and MCP-1 2518 G/A (rs 1024611) in 296 patients diagnosed with prostate cancer and in 311 healthy controls from the same area. RESULTS Diagnosis of prostate cancer was significantly associated with TNF-A GA + AA genotype (OR, 1.61; 95% CI, 1.09-2.64) and RANTES GA + AA genotype (OR, 1.44; 95% CI, 1.09-2.38). A alleles in TNF-A and RANTES influenced prostate cancer susceptibility and acted independently of each other in these subjects. No epistatic effect was found for the combination of different polymorphisms studied. Finally, no overall association was found between prostate cancer risk and IL1-A or MCP-1 polymorphisms. CONCLUSION Our results and previously published findings on genes associated with innate immunity support the hypothesis that polymorphisms in proinflammatory genes may be important in prostate cancer development.

Systematic survey of clonal complexity in tuberculosis at a populational level and detailed characterization of the isolates involved.

Clonally complex infections by Mycobacterium tuberculosis are progressively more accepted. Studies of their dimension in epidemiological scenarios where the infective pressure is not high are scarce. Our study systematically searched for clonally complex infections (mixed infections by more than one strain and simultaneous presence of clonal variants) by applying mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) analysis to M. tuberculosis isolates from two population-based samples of respiratory (703 cases) and respiratory-extrapulmonary (R+E) tuberculosis (TB) cases (71 cases) in a context of moderate TB incidence. Clonally complex infections were found in 11 (1.6%) of the respiratory TB cases and in 10 (14.1%) of those with R+E TB. Among the 21 cases with clonally complex TB, 9 were infected by 2 independent strains and the remaining 12 showed the simultaneous presence of 2 to 3 clonal variants. For the 10 R+E TB cases with clonally complex infections, compartmentalization (different compositions of strains/clonal variants in independent infected sites) was found in 9 of them. All the strains/clonal variants were also genotyped by IS6110-based restriction fragment length polymorphism analysis, which split two MIRU-defined clonal variants, although in general, it showed a lower discriminatory power to identify the clonal heterogeneity revealed by MIRU-VNTR analysis. The comparative analysis of IS6110 insertion sites between coinfecting clonal variants showed differences in the genes coding for a cutinase, a PPE family protein, and two conserved hypothetical proteins. Diagnostic delay, existence of previous TB, risk for overexposure, and clustered/orphan status of the involved strains were analyzed to propose possible explanations for the cases with clonally complex infections. Our study characterizes in detail all the clonally complex infections by M. tuberculosis found in a systematic survey and contributes to the characterization that these phenomena can be found to an extent higher than expected, even in an unselected population-based sample lacking high infective pressure.

Prospective universal application of mycobacterial interspersed repetitive-unit-variable-number tandem-repeat genotyping to characterize Mycobacterium tuberculosis isolates for fast identification of clustered and orphan cases.

The use of molecular tools for genotyping Mycobacterium tuberculosis isolates in epidemiological surveys in order to identify clustered and orphan strains requires faster response times than those offered by the reference method, IS6110 restriction fragment length polymorphism (RFLP) genotyping. A method based on PCR, the mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) genotyping technique, is an option for fast fingerprinting of M. tuberculosis, although precise evaluations of correlation between MIRU-VNTR and RFLP findings in population-based studies in different contexts are required before the methods are switched. In this study, we evaluated MIRU-VNTR genotyping (with a set of 15 loci [MIRU-15]) in parallel to RFLP genotyping in a 39-month universal population-based study in a challenging setting with a high proportion of immigrants. For 81.9% (281/343) of the M. tuberculosis isolates, both RFLP and MIRU-VNTR types were obtained. The percentages of clustered cases were 39.9% (112/281) and 43.1% (121/281) for RFLP and MIRU-15 analyses, and the numbers of clusters identified were 42 and 45, respectively. For 85.4% of the cases, the RFLP and MIRU-15 results were concordant, identifying the same cases as clustered and orphan (kappa, 0.7). However, for the remaining 14.6% of the cases, discrepancies were observed: 16 of the cases clustered by RFLP analysis were identified as orphan by MIRU-15 analysis, and 25 cases identified as orphan by RFLP analysis were clustered by MIRU-15 analysis. When discrepant cases showing subtle genotypic differences were tolerated, the discrepancies fell from 14.6% to 8.6%. Epidemiological links were found for 83.8% of the cases clustered by both RFLP and MIRU-15 analyses, whereas for the cases clustered by RFLP or MIRU-VNTR analysis alone, links were identified for only 30.8% or 38.9% of the cases, respectively. The latter group of cases mainly comprised isolates that could also have been clustered, if subtle genotypic differences had been tolerated. MIRU-15 genotyping seems to be a good alternative to RFLP genotyping for real-time interventional schemes. The correlation between MIRU-15 and IS6110 RFLP findings was reasonable, although some uncertainties as to the assignation of clusters by MIRU-15 analysis were identified.

Dynamic arterial elastance to predict arterial pressure response to volume loading in preload-dependent patients

INTRODUCTION Hemodynamic resuscitation should be aimed at achieving not only adequate cardiac output but also sufficient mean arterial pressure (MAP) to guarantee adequate tissue perfusion pressure. Since the arterial pressure response to volume expansion (VE) depends on arterial tone, knowing whether a patient is preload-dependent provides only a partial solution to the problem. The objective of this study was to assess the ability of a functional evaluation of arterial tone by dynamic arterial elastance (Ea(dyn)), defined as the pulse pressure variation (PPV) to stroke volume variation (SVV) ratio, to predict the hemodynamic response in MAP to fluid administration in hypotensive, preload-dependent patients with acute circulatory failure. METHODS We performed a prospective clinical study in an adult medical/surgical intensive care unit in a tertiary care teaching hospital, including 25 patients with controlled mechanical ventilation who were monitored with the Vigileo(®) monitor, for whom the decision to give fluids was made because of the presence of acute circulatory failure, including arterial hypotension (MAP ≤65 mmHg or systolic arterial pressure <90 mmHg) and preserved preload responsiveness condition, defined as a SVV value ≥10%. RESULTS Before fluid infusion, Ea(dyn) was significantly different between MAP responders (MAP increase ≥15% after VE) and MAP nonresponders. VE-induced increases in MAP were strongly correlated with baseline Ea(dyn) (r(2) = 0.83; P < 0.0001). The only predictor of MAP increase was Ea(dyn) (area under the curve, 0.986 ± 0.02; 95% confidence interval (CI), 0.84-1). A baseline Ea(dyn) value >0.89 predicted a MAP increase after fluid administration with a sensitivity of 93.75% (95% CI, 69.8%-99.8%) and a specificity of 100% (95% CI, 66.4%-100%). CONCLUSIONS Functional assessment of arterial tone by Ea(dyn), measured as the PVV to SVV ratio, predicted arterial pressure response after volume loading in hypotensive, preload-dependent patients under controlled mechanical ventilation.

Brachial artery peak velocity variation to predict fluid responsiveness in mechanically ventilated patients

INTRODUCTION Although several parameters have been proposed to predict the hemodynamic response to fluid expansion in critically ill patients, most of them are invasive or require the use of special monitoring devices. The aim of this study is to determine whether noninvasive evaluation of respiratory variation of brachial artery peak velocity flow measured using Doppler ultrasound could predict fluid responsiveness in mechanically ventilated patients. METHODS We conducted a prospective clinical research in a 17-bed multidisciplinary ICU and included 38 mechanically ventilated patients for whom fluid administration was planned due to the presence of acute circulatory failure. Volume expansion (VE) was performed with 500 mL of a synthetic colloid. Patients were classified as responders if stroke volume index (SVi) increased >or= 15% after VE. The respiratory variation in Vpeakbrach (DeltaVpeakbrach) was calculated as the difference between maximum and minimum values of Vpeakbrach over a single respiratory cycle, divided by the mean of the two values and expressed as a percentage. Radial arterial pressure variation (DeltaPPrad) and stroke volume variation measured using the FloTrac/Vigileo system (DeltaSVVigileo), were also calculated. RESULTS VE increased SVi by >or= 15% in 19 patients (responders). At baseline, DeltaVpeakbrach, DeltaPPrad and DeltaSVVigileo were significantly higher in responder than nonresponder patients [14 vs 8%; 18 vs. 5%; 13 vs 8%; P < 0.0001, respectively). A DeltaVpeakbrach value >10% predicted fluid responsiveness with a sensitivity of 74% and a specificity of 95%. A DeltaPPrad value >10% and a DeltaSVVigileo >11% predicted volume responsiveness with a sensitivity of 95% and 79%, and a specificity of 95% and 89%, respectively. CONCLUSIONS Respiratory variations in brachial artery peak velocity could be a feasible tool for the noninvasive assessment of fluid responsiveness in patients with mechanical ventilatory support and acute circulatory failure. TRIAL REGISTRATION ClinicalTrials.gov ID: NCT00890071.

Exposure to Trihalomethanes through Different Water Uses and Birth Weight, Small for Gestational Age, and Preterm Delivery in Spain

BACKGROUND Evidence associating exposure to water disinfection by-products with reduced birth weight and altered duration of gestation remains inconclusive. OBJECTIVE We assessed exposure to trihalomethanes (THMs) during pregnancy through different water uses and evaluated the association with birth weight, small for gestational age (SGA), low birth weight (LBW), and preterm delivery. METHODS Mother-child cohorts set up in five Spanish areas during the years 2000-2008 contributed data on water ingestion, showering, bathing, and swimming in pools. We ascertained residential THM levels during pregnancy periods through ad hoc sampling campaigns (828 measurements) and regulatory data (264 measurements), which were modeled and combined with personal water use and uptake factors to estimate personal uptake. We defined outcomes following standard definitions and included 2,158 newborns in the analysis. RESULTS Median residential THM ranged from 5.9 μg/L (Valencia) to 114.7 μg/L (Sabadell), and speciation differed across areas. We estimated that 89% of residential chloroform and 96% of brominated THM uptakes were from showering/bathing. The estimated change of birth weight for a 10% increase in residential uptake was -0.45 g (95% confidence interval: -1.36, 0.45 g) for chloroform and 0.16 g (-1.38, 1.70 g) for brominated THMs. Overall, THMs were not associated with SGA, LBW, or preterm delivery. CONCLUSIONS Despite the high THM levels in some areas and the extensive exposure assessment, results suggest that residential THM exposure during pregnancy driven by inhalation and dermal contact routes is not associated with birth weight, SGA, LBW, or preterm delivery in Spain.

Residential Exposure to Outdoor Air Pollution during Pregnancy and Anthropometric Measures at Birth in a Multicenter Cohort in Spain

BACKGROUND. A growing body of research suggests that prenatal exposure to air pollution may be harmful to fetal development. We assessed the association between exposure to air pollution during pregnancy and anthropometric measures at birth in four areas within the Spanish Children's Health and Environment (INMA) mother and child cohort study. METHODS. Exposure to ambient nitrogen dioxide (NO2) and benzene was estimated for the residence of each woman (n = 2,337) for each trimester and for the entire pregnancy. Outcomes included birth weight, length, and head circumference. The association between residential outdoor air pollution exposure and birth outcomes was assessed with linear regression models controlled for potential confounders. We also performed sensitivity analyses for the subset of women who spent more time at home during pregnancy. Finally, we performed a combined analysis with meta-analysis techniques. RESULTS. In the combined analysis, an increase of 10 µg/m3 in NO2 exposure during pregnancy was associated with a decrease in birth length of -0.9 mm [95% confidence interval (CI), -1.8 to -0.1 mm]. For the subset of women who spent ≥ 15 hr/day at home, the association was stronger (-0.16 mm; 95% CI, -0.27 to -0.04). For this same subset of women, a reduction of 22 g in birth weight was associated with each 10-µg/m3 increase in NO2 exposure in the second trimester (95% CI, -45.3 to 1.9). We observed no significant relationship between benzene levels and birth outcomes. CONCLUSIONS. NO2 exposure was associated with reductions in both length and weight at birth. This association was clearer for the subset of women who spent more time at home.

Bovine tuberculosis in Doñana Biosphere Reserve: the role of wild ungulates as disease reservoirs in the last Iberian lynx strongholds.

Doñana National Park (DNP) in southern Spain is a UNESCO Biosphere Reserve where commercial hunting and wildlife artificial feeding do not take place and traditional cattle husbandry still exists. Herein, we hypothesized that Mycobacterium bovis infection prevalence in wild ungulates will depend on host ecology and that variation in prevalence will reflect variation in the interaction between hosts and environmental risk factors. Cattle bTB reactor rates increased in DNP despite compulsory testing and culling of infected animals. In this study, 124 European wild boar, 95 red deer, and 97 fallow deer were sampled from April 2006 to April 2007 and analyzed for M. bovis infection. Modelling and GIS were used to identify risk factors and intra and inter-species relationships. Infection with M. bovis was confirmed in 65 (52.4%) wild boar, 26 (27.4%) red deer and 18 (18.5%) fallow deer. In the absence of cattle, wild boar M. bovis prevalence reached 92.3% in the northern third of DNP. Wild boar showed more than twice prevalence than that in deer (p<0.001). Modelling revealed that M. bovis prevalence decreased from North to South in wild boar (p<0.001) and red deer (p<0.01), whereas no spatial pattern was evidenced for fallow deer. Infection risk in wild boar was dependent on wild boar M. bovis prevalence in the buffer area containing interacting individuals (p<0.01). The prevalence recorded in this study is among the highest reported in wildlife. Remarkably, this high prevalence occurs in the absence of wildlife artificial feeding, suggesting that a feeding ban alone would have a limited effect on wildlife M. bovis prevalence. In DNP, M. bovis transmission may occur predominantly at the intra-species level due to ecological, behavioural and epidemiological factors. The results of this study allow inferring conclusions on epidemiological bTB risk factors in Mediterranean habitats that are not managed for hunting purposes. Our results support the need to consider wildlife species for the control of bTB in cattle and strongly suggest that bTB may affect animal welfare and conservation.