Genomic Analyses of the Microsporidian Nosema ceranae,an Emergent Pathogen of Honey Bees

Recent steep declines in honey bee health have severely impacted the beekeeping industry, presenting new risks for agricultural commodities that depend on insect pollination. Honey bee declines could reflect increased pressures from parasites and pathogens. The incidence of the microsporidian pathogen Nosema ceranae has increased significantly in the past decade. Here we present a draft assembly (7.86 MB) of the N. ceranae genome derived from pyrosequence data, including initial gene models and genomic comparisons with other members of this highly derived fungal lineage. N. ceranae has a strongly AT-biased genome (74% A+T) and a diversity of repetitive elements, complicating the assembly. Of 2,614 predicted protein-coding sequences, we conservatively estimate that 1,366 have homologs in the microsporidian Encephalitozoon cuniculi, the most closely related published genome sequence. We identify genes conserved among microsporidia that lack clear homology outside this group, which are of special interest as potential virulence factors in this group of obligate parasites. A substantial fraction of the diminutive N. ceranae proteome consists of novel and transposable-element proteins. For a majority of well-supported gene models, a conserved sense-strand motif can be found within 15 bases upstream of the start codon; a previously uncharacterized version of this motif is also present in E. cuniculi. These comparisons provide insight into the architecture, regulation, and evolution of microsporidian genomes, and will drive investigations into honey bee–Nosema interactions.     

Lentinula edodes-Derived Polysaccharide Alters the Spatial Structure of Gut Microbiota in Mice

Lentinula edodes-derived polysaccharides possess many therapeutic characteristics, including anti-tumor and immuno-modulation. The gut microbes play a critical role in modulation of immune function. However, the impact of Lentinula edodes-derived polysaccharides on the gut microbes have not yet been explored. In this study, high-throughput pyrosequencing technique was employed to investigate the effects of a new heteropolysaccharide L2 from Lentinula edodes on microbiota diversity and composition of small intestine, cecum, colon and distal end of colon (feces) in mice. The results demonstrated that along mouse intestine the microbiota exhibit distinctly different space distribution. L2 treatment reduced the diversity and evenness of gut microbiota along the intestine, especially in the cecum and colon. In the fecal microbial communities, the decrease of Bacteroidetes by significantly increasing Proteobacteria were observed, which were characterized by the increased Helicobacteraceae and reduced S24-7 at family level. Some OTUs, corresponding to Bacteroides acidifaciens, Alistipes and Helicobacter suncus, were found to be significantly increased in L2 treated-mice. In particular, 4 phyla Chloroflexi, Gemmatimonadetes, Nitrospirae and Planctomycetes are exclusively present in L2-treated mice. This is helpful for further demonstrating healthy action mechanism of Lentinula edodes-derived polysaccharide L2.     

Production of the Catechol Type Siderophore Bacillibactin by the Honey Bee Pathogen Paenibacillus larvae

The Gram-positive bacterium Paenibacillus larvae is the etiological agent of American Foulbrood. This bacterial infection of honey bee brood is a notifiable epizootic posing a serious threat to global honey bee health because not only individual larvae but also entire colonies succumb to the disease. In the recent past considerable progress has been made in elucidating molecular aspects of host pathogen interactions during pathogenesis of P. larvae infections. Especially the sequencing and annotation of the complete genome of P. larvae was a major step forward and revealed the existence of several giant gene clusters coding for non-ribosomal peptide synthetases which might act as putative virulence factors. We here present the detailed analysis of one of these clusters which we demonstrated to be responsible for the biosynthesis of bacillibactin, a P. larvae siderophore. We first established culture conditions allowing the growth of P. larvae under iron-limited conditions and triggering siderophore production by P. larvae. Using a gene disruption strategy we linked siderophore production to the expression of an uninterrupted bacillibactin gene cluster. In silico analysis predicted the structure of a trimeric trithreonyl lactone (DHB-Gly-Thr)₃ similar to the structure of bacillibactin produced by several Bacillus species. Mass spectrometric analysis unambiguously confirmed that the siderophore produced by P. larvae is identical to bacillibactin. Exposure bioassays demonstrated that P. larvae bacillibactin is not required for full virulence of P. larvae in laboratory exposure bioassays. This observation is consistent with results obtained for bacillibactin in other pathogenic bacteria.     

Four Common Pesticides,Their Mixtures and a Formulation Solvent in the Hive Environment Have High Oral Toxicity to Honey Bee Larvae

Recently, the widespread distribution of pesticides detected in the hive has raised serious concerns about pesticide exposure on honey bee (Apis mellifera L.) health. A larval rearing method was adapted to assess the chronic oral toxicity to honey bee larvae of the four most common pesticides detected in pollen and wax - fluvalinate, coumaphos, chlorothalonil, and chloropyrifos - tested alone and in all combinations. All pesticides at hive-residue levels triggered a significant increase in larval mortality compared to untreated larvae by over two fold, with a strong increase after 3 days of exposure. Among these four pesticides, honey bee larvae were most sensitive to chlorothalonil compared to adults. Synergistic toxicity was observed in the binary mixture of chlorothalonil with fluvalinate at the concentrations of 34 mg/L and 3 mg/L, respectively; whereas, when diluted by 10 fold, the interaction switched to antagonism. Chlorothalonil at 34 mg/L was also found to synergize the miticide coumaphos at 8 mg/L. The addition of coumaphos significantly reduced the toxicity of the fluvalinate and chlorothalonil mixture, the only significant non-additive effect in all tested ternary mixtures. We also tested the common ‘inert’ ingredient N-methyl-2-pyrrolidone at seven concentrations, and documented its high toxicity to larval bees. We have shown that chronic dietary exposure to a fungicide, pesticide mixtures, and a formulation solvent have the potential to impact honey bee populations, and warrants further investigation. We suggest that pesticide mixtures in pollen be evaluated by adding their toxicities together, until complete data on interactions can be accumulated.     

Susceptibility of Various Species of Lepidopterous Pupae to the Entomogenous Nematode Neoaplectana carpocapsae

The susceptibility of certain species of lepidopterous pupae occurring in different ecological situations to the entomogenous nematode, Neoaplectana carpocapsae, was tested. Soil- or litter-pupating lepidopterous insects were not highly susceptible to N. carpocapsae. The most susceptible insect pupating in the soil was Spodoptera exigua with 63% pupal mortality, while Harrisinia brillians, which pupates in litter, had 55% mortality. Other soil- and litter-pupating insects had mortalities of less than 25%. Some insect species that pupate above ground were highly susceptible (> 84% mortality) to N. carpocapsae infection.     

Manipulation of the Xanthophyll Cycle Increases Plant Susceptibility to Sclerotinia sclerotiorum

The xanthophyll cycle is involved in dissipating excess light energy to protect the photosynthetic apparatus in a process commonly assessed from non-photochemical quenching (NPQ) of chlorophyll fluorescence. Here, it is shown that the xanthophyll cycle is modulated by the necrotrophic pathogen Sclerotinia sclerotiorum at the early stage of infection. Incubation of Sclerotinia led to a localized increase in NPQ even at low light intensity. Further studies showed that this abnormal change in NPQ was closely correlated with a decreased pH caused by Sclerotinia-secreted oxalate, which might decrease the ATP synthase activity and lead to a deepening of thylakoid lumen acidification under continuous illumination. Furthermore, suppression (with dithiothreitol) or a defect (in the npq1-2 mutant) of violaxanthin de-epoxidase (VDE) abolished the Sclerotinia-induced NPQ increase. HPLC analysis showed that the Sclerotinia-inoculated tissue accumulated substantial quantities of zeaxanthin at the expense of violaxanthin, with a corresponding decrease in neoxanthin content. Immunoassays revealed that the decrease in these xanthophyll precursors reduced de novo abscisic acid (ABA) biosynthesis and apparently weakened tissue defense responses, including ROS induction and callose deposition, resulting in enhanced plant susceptibility to Sclerotinia. We thus propose that Sclerotinia antagonizes ABA biosynthesis to suppress host defense by manipulating the xanthophyll cycle in early pathogenesis. These findings provide a model of how photoprotective metabolites integrate into the defense responses, and expand the current knowledge of early plant-Sclerotinia interactions at infection sites.     

Multiparameter Viability Assay for Stress Profiling Applied to the Food Pathogen Listeria monocytogenes F2365

A novel generic approach for stress profiling was applied to Listeria monocytogenes strain F2365. This food-borne pathogen was exposed to gradients of five different stresses of increasing intensity, typically ranging from moderate to lethal conditions. The stress factors included heat, acidic pH, a detergent disinfectant, an oxidant, and hyperosmotic conditions. In addition to CFU counts and lag time, five different molecular viability parameters were measured by fluorescence-based assays, including membrane integrity, membrane potential, esterase activity, redox activity, and intracellular pH stability. The last was measured by our recently invented real-time viability assay. Exposure to all stresses resulted in clear dose-response relationships for all viability parameters with the exception of hyperosmotic conditions. A statistical analysis showed strong correlations for (i) the growth parameters plate counts and lag times, (ii) the enzyme-associated functions redox and esterase activity, and (iii) the membrane-associated pH stability and membrane integrity. Results indicated a pronounced difference in the susceptibilities of the measured parameters depending on the stress factor applied. However, at relatively high stress intensities, all of the viability parameters became affected independent of the stress factor. Applications of the approach presented here include studies on the mechanism of action of unknown compounds with biocidal activity and a comparative analysis of the severities of the impact of stress conditions of interest. It appears that a meaningful evaluation of the impact of mild stress conditions can be obtained only through measurement of multiple viability parameters.     

Genome-Wide Transcriptional Response of Silkworm (Bombyx mori) to Infection by the Microsporidian Nosema bombycis

Microsporidia have attracted much attention because they infect a variety of species ranging from protists to mammals, including immunocompromised patients with AIDS or cancer. Aside from the study on Nosema ceranae, few works have focused on elucidating the mechanism in host response to microsporidia infection. Nosema bombycis is a pathogen of silkworm pébrine that causes great economic losses to the silkworm industry. Detailed understanding of the host (Bombyx mori) response to infection by N. bombycis is helpful for prevention of this disease. A genome-wide survey of the gene expression profile at 2, 4, 6 and 8 days post-infection by N. bombycis was performed and results showed that 64, 244, 1,328, 1,887 genes were induced, respectively. Up to 124 genes, which are involved in basal metabolism pathways, were modulated. Notably, B. mori genes that play a role in juvenile hormone synthesis and metabolism pathways were induced, suggesting that the host may accumulate JH as a response to infection. Interestingly, N. bombycis can inhibit the silkworm serine protease cascade melanization pathway in hemolymph, which may be due to the secretion of serpins in the microsporidia. N. bombycis also induced up-regulation of several cellular immune factors, in which CTL11 has been suggested to be involved in both spore recognition and immune signal transduction. Microarray and real-time PCR analysis indicated the activation of silkworm Toll and JAK/STAT pathways. The notable up-regulation of antimicrobial peptides, including gloverins, lebocins and moricins, strongly indicated that antimicrobial peptide defense mechanisms were triggered to resist the invasive microsporidia. An analysis of N. bombycis-specific response factors suggested their important roles in anti-microsporidia defense. Overall, this study primarily provides insight into the potential molecular mechanisms for the host-parasite interaction between B. mori and N. bombycis and may provide a foundation for further work on host-parasite interaction between insects and microsporidia.     

Genetic Resistance to Bacillus thuringiensis Alters Feeding Behaviour in the Cabbage Looper,Trichoplusia ni

Evolved resistance to xenobiotics and parasites is often associated with fitness costs when the selection pressure is absent. Resistance to the widely used microbial insecticide Bacillus thuringiensis (Bt) has evolved in several insect species through the modification of insect midgut binding sites for Bt toxins, and reports of costs associated with Bt resistance are common. Studies on the costs of Bt-resistance restrict the insect to a single artificial diet or host-plant. However, it is well documented that insects can self-select appropriate proportions of multiple nutritionally unbalanced foods to optimize life-history traits. Therefore, we examined whether Bt-resistant and susceptible cabbage loopers Trichoplusia ni differed in their nutrient intake and fitness costs when they were allowed to compose their own protein:carbohydrate diet. We found that Bt-resistant T. ni composed a higher ratio of protein to carbohydrate than susceptible T. ni. Bt-resistant males exhibited no fitness cost, while the fitness cost (reduced pupal weight) was present in resistant females. The absence of the fitness cost in resistant males was associated with increased carbohydrate consumption compared to females. We demonstrate a sex difference in a fitness cost and a new behavioural outcome associated with Bt resistance.     

XPF-673C>T Polymorphism Effect on the Susceptibility to Esophageal Cancer in Chinese Population

Purpose

Xeroderma pigmentsum group F (XPF) plays a pivotal role in DNA nucleotide excision repair and has been linked to the development of various cancers. This study aims to assess the association of XPF genetic variants with the susceptibility to esophageal squamous cell carcinoma (ESCC) in Chinese population.

Methods

This two-stage case-control study was conducted in a total of 1524 patients with ESCC and 1524 controls. Genotype of XPF -673C>T and 11985A>G variants were determined by polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP). Logistic regression analysis was performed to estimate odd ratios (ORs) and 95% confidence intervals (95% CI).

Results

Our case-control study showed that XPF -673TT genotype was associated with a decreased risk of ESCC compared with CC genotype in both case-control sets (Tangshan set: OR = 0.58; 95%CI = 0.34–0.99, P = 0.040; Beijing set: OR = 0.66; 95%CI = 0.46–0.95, P = 0.027). Stratified analyses revealed that a multiplicative interaction between -673C>T variant and age, sex or smoking status was evident (Gene-age: P_interaction = 0.002; Gene-sex: P_interaction = 0.002; Gene-smoking: P_interaction = 0.002). For XPF 11985A>G polymorphism, there was no significant difference of genotype distribution between ESCC cases and controls.

Conclusion

These findings indicated that genetic variants in XPF might contribute to the susceptibility to ESCC.     

A Comparative Chemogenomics Strategy to Predict Potential Drug Targets in the Metazoan Pathogen,Schistosoma mansoni

Schistosomiasis is a prevalent and chronic helmintic disease in tropical regions. Treatment and control relies on chemotherapy with just one drug, praziquantel and this reliance is of concern should clinically relevant drug resistance emerge and spread. Therefore, to identify potential target proteins for new avenues of drug discovery we have taken a comparative chemogenomics approach utilizing the putative proteome of Schistosoma mansoni compared to the proteomes of two model organisms, the nematode, Caenorhabditis elegans and the fruitfly, Drosophila melanogaster. Using the genome comparison software Genlight, two separate in silico workflows were implemented to derive a set of parasite proteins for which gene disruption of the orthologs in both the model organisms yielded deleterious phenotypes (e.g., lethal, impairment of motility), i.e., are essential genes/proteins. Of the 67 and 68 sequences generated for each workflow, 63 were identical in both sets, leading to a final set of 72 parasite proteins. All but one of these were expressed in the relevant developmental stages of the parasite infecting humans. Subsequent in depth manual curation of the combined workflow output revealed 57 candidate proteins. Scrutiny of these for ‘druggable’ protein homologs in the literature identified 35 S. mansoni sequences, 18 of which were homologous to proteins with 3D structures including co-crystallized ligands that will allow further structure-based drug design studies. The comparative chemogenomics strategy presented generates a tractable set of S. mansoni proteins for experimental validation as drug targets against this insidious human pathogen.     

A Dichotomy in Cortical Actin and Chemotactic Actin Activity between Human Memory and Naive T Cells Contributes to Their Differential Susceptibility to HIV-1 Infection

Human memory and naive CD4 T cells can mainly be identified by the reciprocal expression of the CD45RO or CD45RA isoforms. In HIV-1 infection, blood CD45RO memory CD4 T cells are preferentially infected and serve as a major viral reservoir. The molecular mechanism dictating this differential susceptibility to HIV-1 remains largely obscure. Here, we report that the different susceptibility of memory and naive T cells to HIV is not determined by restriction factors such as Apobec3G or BST2. However, we observed a phenotypic distinction between human CD45RO and CD45RA resting CD4 T cells in their cortical actin density and actin dynamics. CD45RO CD4 T cells possess a higher cortical actin density and can be distinguished as CD45RO⁺Actin^high. In contrast, CD45RA T cells are phenotypically CD45RA⁺Actin^low. In addition, the cortical actin in CD45RO memory CD4 T cells is more dynamic and can respond to low dosages of chemotactic induction by SDF-1, whereas that of naive cells cannot, despite a similar level of the chemokine receptor CXCR4 present on both cells. We further demonstrate that this difference in the cortical actin contributes to their differential susceptibility to HIV-1; resting memory but not naive T cells are highly responsive to HIV-mediated actin dynamics that promote higher levels of viral entry and early DNA synthesis in resting memory CD4 T cells. Furthermore, transient induction of actin dynamics in resting naive T cells rescues HIV latent infection following CD3/CD28 stimulation. These results suggest a key role of chemotactic actin activity in facilitating HIV-1 latent infection of these T cell subsets.     

Investigating Differences across Host Species and Scales to Explain the Distribution of the Amphibian Pathogen Batrachochytrium dendrobatidis

Many pathogens infect more than one host species, and clarifying how these different hosts contribute to pathogen dynamics can facilitate the management of pathogens and can lend insight into the functioning of pathogens in ecosystems. In this study, we investigated a suite of native and non-native amphibian hosts of the pathogen Batrachochytrium dendrobatidis (Bd) across multiple scales to identify potential mechanisms that may drive infection patterns in the Colorado study system. Specifically, we aimed to determine if: 1) amphibian populations vary in Bd infection across the landscape, 2) amphibian community composition predicts infection (e.g., does the presence or abundance of any particular species influence infection in others?), 3) amphibian species vary in their ability to produce infectious zoospores in a laboratory infection, 4) heterogeneity in host ability observed in the laboratory scales to predict patterns of Bd prevalence in the landscape. We found that non-native North American bullfrogs (Lithobates catesbeianus) are widespread and have the highest prevalence of Bd infection relative to the other native species in the landscape. Additionally, infection in some native species appears to be related to the density of sympatric L. catesbeianus populations. At the smaller host scale, we found that L. catesbeianus produces more of the infective zoospore stage relative to some native species, but that this zoospore output does not scale to predict infection in sympatric wild populations of native species. Rather, landscape level infection relates most strongly to density of hosts at a wetland as well as abiotic factors. While non-native L. catesbeianus have high levels of Bd infection in the Colorado Front Range system, we also identified Bd infection in a number of native amphibian populations allopatric with L. catesbeianus, suggesting that multiple host species are important contributors to the dynamics of the Bd pathogen in this landscape.     

Human Chondrocytes Respond Discordantly to the Protein Encoded by the Osteoarthritis Susceptibility Gene GDF5

A genetic deficit mediated by SNP rs143383 that leads to reduced expression of GDF5 is strongly associated with large-joint osteoarthritis. We speculated that this deficit could be attenuated by the application of exogenous GDF5 protein and as a first step we have assessed what effect such application has on primary osteoarthritis chondrocyte gene expression. Chondrocytes harvested from cartilage of osteoarthritic patients who had undergone joint replacement were cultured with wildtype recombinant mouse and human GDF5 protein. We also studied variants of GDF5, one that has a higher affinity for the receptor BMPR-IA and one that is insensitive to the GDF5 antagonist noggin. As a positive control, chondrocytes were treated with TGF-β1. Chondrocytes were cultured in monolayer and micromass and the expression of genes coding for catabolic and anabolic proteins of cartilage were measured by quantitative PCR. The expression of the GDF5 receptor genes and the presence of their protein products was confirmed and the ability of GDF5 signal to translocate to the nucleus was demonstrated by the activation of a luciferase reporter construct. The capacity of GDF5 to elicit an intracellular signal in chondrocytes was demonstrated by the phosphorylation of intracellular Smads. Chondrocytes cultured with TGF-β1 demonstrated a consistent down regulation of MMP1, MMP13 and a consistent upregulation of TIMP1 and COL2A1 with both culture techniques. In contrast, chondrocytes cultured with wildtype GDF5, or its variants, did not show any consistent response, irrespective of the culture technique used. Our results show that osteoarthritis chondrocytes do not respond in a predictable manner to culture with exogenous GDF5. This may be a cause or a consequence of the osteoarthritis disease process and will need to be surmounted if treatment with exogenous GDF5 is to be advanced as a potential means to overcome the genetic deficit conferring osteoarthritis susceptibility at this gene.     

Mycobacterium tuberculosis Genotypes Determined by Spoligotyping to Be Circulating in Colombia between 1999 and 2012 and Their Possible Associations with Transmission and Susceptibility to First-Line Drugs

Introduction

Tuberculosis (TB) remains a primary public health problem worldwide. The number of multidrug-resistant tuberculosis (MDR TB) cases has increased in recent years in Colombia. Knowledge of M. tuberculosis genotypes defined by spoligotyping can help determine the circulation of genotypes that must be controlled to prevent the spread of TB.

Objective

To describe the genotypes of M. tuberculosis using spoligotyping in resistant and drug-sensitive isolates and their possible associations with susceptibility to first-line drugs.

Methods

An analytical observational study was conducted that included 741 isolates of M. tuberculosis from patients. The isolates originated from 31 departments and were obtained by systematic surveillance between 1999 and 2012.

Results

In total 61.94% of the isolates were resistant to 1 or more drugs, and 147 isolates were MDR. In total, 170 genotypes were found in the population structure of Colombian M. tuberculosis isolates. The isolates were mainly represented by four families: LAM (39.9%), Haarlem (19%), Orphan (17%) and T (9%). The SIT42 (LAM 9) was the most common genotype and contained 24.7% of the isolates, followed by the genotypes SIT62 (Haarlem1), SIT53 (T1), and SIT50 (H3). A high clustering of isolates was evident with 79.8% of the isolates classified into 32 groups. The Beijing family was associated with resistant isolates, whereas the Haarlem and T families were associated with sensitive isolates. The Haarlem family was also associated with grouped isolates (p = 0.031).

Conclusions

A high proportion (approximately 80%) of isolates was found in clusters; these clusters were not associated with resistance to first-line drugs. The Beijing family was associated with drug resistance, whereas the T and Haarlem families were associated with susceptibility in the Colombian isolates studied.     

A New Pathogen Transmission Mechanism in the Ocean: The Case of Sea Otter Exposure to the Land-Parasite Toxoplasma gondii

Toxoplasma gondii is a land-derived parasite that infects humans and marine mammals. Infections are a significant cause of mortality for endangered southern sea otters (Enhydra lutris nereis), but the transmission mechanism is poorly understood. Otter exposure to T. gondii has been linked to the consumption of marine turban snails in kelp (Macrocystis pyrifera) forests. It is unknown how turban snails acquire oocysts, as snails scrape food particles attached to surfaces, whereas T. gondii oocysts enter kelp beds as suspended particles via runoff. We hypothesized that waterborne T. gondii oocysts attach to kelp surfaces when encountering exopolymer substances (EPS) forming the sticky matrix of biofilms on kelp, and thus become available to snails. Results of a dietary composition analysis of field-collected snails and of kelp biofilm indicate that snails graze the dense kelp-biofilm assemblage composed of pennate diatoms and bacteria inserted within the EPS gel-like matrix. To test whether oocysts attach to kelp blades via EPS, we designed a laboratory experiment simulating the kelp forest canopy in tanks spiked with T. gondii surrogate microspheres and controlled for EPS and transparent exopolymer particles (TEP - the particulate form of EPS). On average, 19% and 31% of surrogates were detected attached to kelp surfaces covered with EPS in unfiltered and filtered seawater treatments, respectively. The presence of TEP in the seawater did not increase surrogate attachment. These findings support a novel transport mechanism of T. gondii oocysts: as oocysts enter the kelp forest canopy, a portion adheres to the sticky kelp biofilms. Snails grazing this biofilm encounter oocysts as ‘bycatch’ and thereby deliver the parasite to sea otters that prey upon snails. This novel mechanism can have health implications beyond T. gondii and otters, as a similar route of pathogen transmission may be implicated with other waterborne pathogens to marine wildlife and humans consuming biofilm-feeding invertebrates.     

Adaptive Potential of Maritime Pine (Pinus pinaster) Populations to the Emerging Pitch Canker Pathogen,Fusarium circinatum

There is a concern on how emerging pests and diseases will affect the distribution range and adaptability of their host species, especially due to different conditions derived from climate change and growing globalization. Fusarium circinatum, which causes pitch canker disease in Pinus species, is an exotic pathogen of recent introduction in Spain that threatens its maritime pine (P. pinaster) stands. To predict the impact this disease will have on the species, we examine host resistance traits and their genetic architecture. Resistance phenotyping was done in a clonal provenance/progeny trial, using three-year-old cuttings artificially inoculated with the pathogen and maintained under controlled environmental conditions. A total number of 670 ramets were assessed, distributed in 10 populations, with a total of 47 families, 2 to 5 half-sibs per family, and 3–7 ramets per clone. High genetic variation was found at the three hierarchical levels studied: population, family and clone, being both additive and non-additive effects important. Narrow-sense and broad-sense heritability estimates were relatively high, with respective values of 0.43–0.58 and 0.51–0.8, depending on the resistance traits measured (lesion length, lesion length rate, time to wilting, and survival). These values suggest the species'' high capacity of evolutionary response to the F. circinatum pathogen. A population originated in Northern Spain was the most resistant, while another from Morocco was the most susceptible. The total number of plants that did not show lesion development or presented a small lesion (length<30 mm) was 224 out of 670, indicating a high proportion of resistant trees in the offspring within the analyzed populations. We found large differences among populations and considerable genetic variation within populations, which should allow, through natural or artificial selection, the successful adaptation of maritime pine to pitch canker disease.