Pitting the olive seed microbiome

BACKGROUND

The complex and co-evolved interplay between plants and their microbiota is crucial for the health and fitness of the plant holobiont. However, the microbiota of the seeds is still relatively unexplored and no studies have been conducted with olive trees so far. In this study, we aimed to characterize the bacterial, fungal and archaeal communities present in seeds of ten olive genotypes growing in the same orchard through amplicon sequencing to test whether the olive genotype is a major driver in shaping the seed microbial community, and to identify the origin of the latter. Therefore, we have developed a methodology for obtaining samples from the olive seed's endosphere under sterile conditions.

RESULTS

A diverse microbiota was uncovered in olive seeds, the plant genotype being an important factor influencing the structure and composition of the microbial communities. The most abundant bacterial phylum was Actinobacteria, accounting for an average relative abundance of 41%. At genus level, Streptomyces stood out because of its potential influence on community structure. Within the fungal community, Basidiomycota and Ascomycota were the most abundant phyla, including the genera Malassezia, Cladosporium, and Mycosphaerella. The shared microbiome was composed of four bacterial (Stenotrophomonas, Streptomyces, Promicromonospora and Acidipropionibacterium) and three fungal (Malassezia, Cladosporium and Mycosphaerella) genera. Furthermore, a comparison between findings obtained here and earlier results from the root endosphere of the same trees indicated that genera such as Streptomyces and Malassezia were present in both olive compartments.

CONCLUSIONS

This study provides the first insights into the composition of the olive seed microbiota. The highly abundant fungal genus Malassezia and the bacterial genus Streptomyces reflect a unique signature of the olive seed microbiota. The genotype clearly shaped the composition of the seed's microbial community, although a shared microbiome was found. We identified genera that may translocate from the roots to the seeds, as they were present in both organs of the same trees. These findings set the stage for future research into potential vertical transmission of olive endophytes and the role of specific microbial taxa in seed germination, development, and seedling survival.

Thriving beneath olive trees: The influence of organic farming on microbial communities

Soil health and root-associated microbiome are interconnected factors involved in plant health. The use of manure amendment on agricultural fields exerts a direct benefit on soil nutrient content and water retention, among others. However, little is known about the impact of manure amendment on the root-associated microbiome, particularly in woody species. In this study, we aimed to evaluate the effects of ovine manure on the microbial communities of the olive rhizosphere and root endosphere. Two adjacent orchards subjected to conventional (CM) and organic (OM) management were selected. We used metabarcoding sequencing to assess the bacterial and fungal communities. Our results point out a clear effect of manure amendment on the microbial community. Fungal richness and diversity were increased in the rhizosphere. The fungal biomass in the rhizosphere was more than doubled, ranging from 1.72 × 10⁶ ± 1.62 × 10⁵ (CM) to 4.54 × 10⁶ ± 8.07 × 10⁵ (OM) copies of the 18 S rRNA gene g^-1 soil. Soil nutrient content was also enhanced in the OM orchard. Specifically, oxidable organic matter, total nitrogen, nitrate, phosphorous, potassium and sulfate concentrations were significantly increased in the OM orchard. Moreover, we predicted a higher abundance of bacteria in OM with metabolic functions involved in pollutant degradation and defence against pathogens. Lastly, microbial co-occurrence network showed more positive interactions, complexity and shorter geodesic distance in the OM orchard. According to our results, manure amendment on olive orchards represents a promising tool for positively modulating the microbial community in direct contact with the plant.

Co-occurrence network analysis unveils the actual differential impact on the olive root microbiota by two Verticillium wilt biocontrol rhizobacteria

BACKGROUND

Verticillium wilt of olive (VWO), caused by Verticillium dahliae Kleb, is one of the most threatening diseases affecting olive cultivation. An integrated disease management strategy is recommended for the effective control of VWO. Within this framework, the use of biological control agents (BCAs) is a sustainable and environmentally friendly approach. No studies are available on the impact that the introduction of BCAs has on the resident microbiota of olive roots. Pseudomonas simiae PICF7 and Paenibacillus polymyxa PIC73 are two BCAs effective against VWO. We examined the effects of the introduction of these BCAs on the structure, composition and co-occurrence networks of the olive (cv. Picual) root-associated microbial communities. The consequences of the subsequent inoculation with V. dahliae on BCA-treated plants were also assessed.

RESULTS

Inoculation with any of the BCAs did not produce significant changes in the structure or the taxonomic composition of the 'Picual' root-associated microbiota. However, significant and distinctive alterations were observed in the topologies of the co-occurrence networks. The introduction of PIC73 provoked a diminution of positive interactions within the 'Picual' microbial community; instead, PICF7 inoculation increased the microbiota's compartmentalization. Upon pathogen inoculation, the network of PIC73-treated plants decreased the number of interactions and showed a switch of keystone species, including taxa belonging to minor abundant phyla (Chloroflexi and Planctomycetes). Conversely, the inoculation of V. dahliae in PICF7-treated plants significantly increased the complexity of the network and the number of links among their modules, suggestive of a more stable network. No changes in their keystone taxa were detected.

CONCLUSION

The absence of significant modifications on the structure and composition of the 'Picual' belowground microbiota due to the introduction of the tested BCAs underlines the low/null environmental impact of these rhizobacteria. These findings may have important practical consequences regarding future field applications of these BCAs. Furthermore, each BCA altered the interactions among the components of the olive belowground microbiota in idiosyncratic ways (i.e. PIC73 strongly modified the number of positive relations in the 'Picual' microbiota whereas PICF7 mostly affected the network stability). These modifications may provide clues on the biocontrol strategies used by these BCAs.

Unveiling Differences in Root Defense Mechanisms Between Tolerant and Susceptible Olive Cultivars to Verticillium dahliae

Verticillium wilt of olive (VWO), caused by the soil-borne vascular pathogen Verticillium dahliae, is one of the most devastating diseases affecting olive, the woody crop. One of the best VWO management measures is the use of tolerant cultivars. Yet, our knowledge about defense mechanisms that operate at the root level to explain tolerance to this disease is incomplete. Moreover, most of the approaches so far followed focus only on a specific mechanistic level (e.g., genetic, physiological, or biochemical) rather than on a holistic/multilevel perspective. In this study, eighteen root functional traits, the time-course expression of nine defense-related genes, the root lignin content, and the root membrane permeability were evaluated in six olive varieties differing in their level of tolerance/susceptibility to VWO. The aim was to find links between the level of tolerance to VWO and specific root defense mechanisms at the structural, genetic, biochemical, and physiological levels. Tolerant and susceptible cultivars showed substantial differences in the root system architecture and root lignin content. VWO-susceptible cultivars presented roots with higher specific length and area, but lower diameter and larger number of forks and tips compared to tolerant varieties that also showed less branched roots, higher root diameter, and larger basal content of lignin. Interestingly, VWO-tolerant varieties significantly increased their root lignin content and root membrane permeability after inoculation with V. dahliae. These results were seldom (or not at all) observed in the susceptible plants. At the genetic level, genes related to defense mechanisms, such as cell wall lignin biosynthesis (C4H and CO-MT), production of hydrolytic enzymes able to degrade the fungal cell wall (β-1.3-glucanase), and activation of innate immunity (BAK1 and WRKY5) increased their expression in tolerant cultivars from early moments after inoculation, in contrast to the susceptible ones. These results showed that differences in the root system architecture and lignin content may greatly determine the performance of olive against colonization and invasion by V. dahliae. Moreover, the increase in root membrane permeability in the presence of the pathogen was a typical response of tolerant cultivars. Finally, VWO-tolerant cultivars were able to mount a more intense and rapid defense-related genetic response to respond to the attack by V. dahliae.

Impacts of the Biocontrol Strain Pseudomonas simiae PICF7 on the Banana Holobiont: Alteration of Root Microbial Co-occurrence Networks and Effect on Host Defense Responses

The impact of the versatile biocontrol and plant-growth-promoting rhizobacteria Pseudomonas simiae PICF7 on the banana holobiont under controlled conditions was investigated. We examine the fate of this biological control agent (BCA) upon introduction in the soil, the effect on the banana root microbiota, and the influence on specific host genetic defense responses. While the presence of strain PICF7 significantly altered neither the composition nor the structure of the root microbiota, a significant shift in microbial community interactions through co-occurrence network analysis was observed. Despite the fact that PICF7 did not constitute a keystone, the topology of this network was significantly modified-the BCA being identified as a constituent of one of the main network modules in bacterized plants. Gene expression analysis showed the early suppression of several systemic acquired resistance and induced systemic resistance (ISR) markers. This outcome occurred at the time in which the highest relative abundance of PICF7 was detected. The absence of major and permanent changes on the banana holobiont upon PICF7 introduction poses advantages regarding the use of this beneficial rhizobacteria under field conditions. Indeed a BCA able to control the target pathogen while altering as little as possible the natural host-associated microbiome should be a requisite when developing effective bio-inoculants.

The Banana Root Endophytome: Differences between Mother Plants and Suckers and Evaluation of Selected Bacteria to Control Fusarium oxysporum f.sp. cubense

This study aimed to disentangle the structure, composition, and co-occurrence relationships of the banana (cv. Dwarf Cavendish) root endophytome comparing two phenological plant stages: mother plants and suckers. Moreover, a collection of culturable root endophytes (>1000) was also generated from Canary Islands. In vitro antagonism assays against Fusarium oxysporum f.sp. cubense (Foc) races STR4 and TR4 enabled the identification and characterization of potential biocontrol agents (BCA). Eventually, three of them were selected and evaluated against Fusarium wilt of banana (FWB) together with the well-known BCA Pseudomonas simiae PICF7 under controlled conditions. Culturable and non-culturable (high-throughput sequencing) approaches provided concordant information and showed low microbial diversity within the banana root endosphere. Pseudomonas appeared as the dominant genus and seemed to play an important role in the banana root endophytic microbiome according to co-occurrence networks. Fungal communities were dominated by the genera Ophioceras, Cyphellophora, Plecosphaerella, and Fusarium. Overall, significant differences were found between mother plants and suckers, suggesting that the phenological stage determines the recruitment and organization of the endophytic microbiome. While selected native banana endophytes showed clear antagonism against Foc strains, their biocontrol performance against FWB did not improve the outcome observed for a non-indigenous reference BCA (strain PICF7).

Assessing the Involvement of Selected Phenotypes of Pseudomonas simiae PICF7 in Olive Root Colonization and Biological Control of Verticillium dahliae

Pseudomonas simiae PICF7 is an indigenous inhabitant of the olive (Olea europaea L.) rhizosphere/root endosphere and an effective biocontrol agent against Verticillium wilt of olive (VWO), caused by the soil-borne fungus Verticillium dahliae. This study aimed to evaluate the potential involvement of selected phenotypes of strain PICF7 in root colonization ability and VWO biocontrol. Therefore, a random transposon-insertion mutant bank of P. simiae PICF7 was screened for the loss of phenotypes likely involved in rhizosphere/soil persistence (copper resistance), root colonization (biofilm formation) and plant growth promotion (phytase activity). Transposon insertions in genes putatively coding for the transcriptional regulator CusR or the chemotaxis protein CheV were found to affect copper resistance, whereas an insertion in fleQ gene putatively encoding a flagellar regulatory protein hampered the ability to form a biofilm. However, these mutants displayed the same antagonistic effect against V. dahliae as the parental strain. Remarkably, two mutants impaired in biofilm formation were never found inside olive roots, whereas their ability to colonize the root exterior and to control VWO remained unaffected. Endophytic colonization of olive roots was unaltered in mutants impaired in copper resistance and phytase production. Results demonstrated that the phenotypes studied were irrelevant for VWO biocontrol.

Comparative study of neighboring Holm oak and olive trees-belowground microbial communities subjected to different soil management

It is well-known that different plant species, and even plant varieties, promote different assemblages of the microbial communities associated with them. Here, we investigate how microbial communities (bacteria and fungi) undergo changes within the influence of woody plants (two olive cultivars, one tolerant and another susceptible to the soilborne fungal pathogen Verticillium dahliae, plus wild Holm oak) grown in the same soil but with different management (agricultural versus native). By the use of metabarcoding sequencing we determined that the native Holm oak trees rhizosphere bacterial communities were different from its bulk soil, with differences in some genera like Gp4, Gp6 and Solirubrobacter. Moreover, the agricultural management used in the olive orchard led to belowground microbiota differences with respect to the natural conditions both in bulk soils and rhizospheres. Indeed, Gemmatimonas and Fusarium were more abundant in olive orchard soils. However, agricultural management removed the differences in the microbial communities between the two olive cultivars, and these differences were minor respect to the olive bulk soil. According to our results, and at least under the agronomical conditions here examined, the composition and structure of the rhizospheric microbial communities do not seem to play a major role in olive tolerance to V. dahliae.

Epigenetic Regulation of Verticillium dahliae Virulence: Does DNA Methylation Level Play A Role?

Verticillium dahliae is the etiological agent of Verticillium wilt of olive. The virulence of Defoliating V. dahliae isolates usually displays differences and high plasticity. This work studied whether an epigenetic mechanism was involved in this plasticity. An inverse correlation between virulence and DNA methylation of protein-coding genes was found. A set of 831 genes was selected for their highly consistent inverse methylation profile and virulence in the five studied isolates. Of these genes, ATP-synthesis was highly represented, which indicates that the more virulent D isolates are, the more energy requirements they may have. Furthermore, there were numerous genes in the protein biosynthesis process: genes coding for the chromatin structure, which suggests that epigenetic changes may also affect chromatin condensation; many transmembrane transporter genes, which is consistent with denser compounds, traffic through membranes in more virulent isolates; a fucose-specific lectin that may play a role in the attachment to plant cell walls during the host infection process; and pathogenic cutinases that facilitate plant invasion and sporulation genes for rapid spreading alongside plants. Our findings support the notion that differences in the virulence of the Defoliating V. dahliae isolates may be controlled, at least to some extent, by an epigenetic mechanism.

Linking belowground microbial network changes to different tolerance level towards Verticillium wilt of olive

BACKGROUND

Verticillium wilt of olive (VWO) is caused by the soilborne fungal pathogen Verticillium dahliae. One of the best VWO management measures is the use of tolerant/resistant olive cultivars. Knowledge on the olive-associated microbiome and its potential relationship with tolerance to biotic constraints is almost null. The aims of this work are (1) to describe the structure, functionality, and co-occurrence interactions of the belowground (root endosphere and rhizosphere) microbial communities of two olive cultivars qualified as tolerant (Frantoio) and susceptible (Picual) to VWO, and (2) to assess whether these communities contribute to their differential disease susceptibility level.

RESULTS

Minor differences in alpha and beta diversities of root-associated microbiota were detected between olive cultivars regardless of whether they were inoculated or not with the defoliating pathotype of V. dahliae. Nevertheless, significant differences were found in taxonomic composition of non-inoculated plants' communities, "Frantoio" showing a higher abundance of beneficial genera in contrast to "Picual" that exhibited major abundance of potential deleterious genera. Upon inoculation with V. dahliae, significant changes at taxonomic level were found mostly in Picual plants. Relevant topological alterations were observed in microbial communities' co-occurrence interactions after inoculation, both at structural and functional level, and in the positive/negative edges ratio. In the root endosphere, Frantoio communities switched to highly connected and low modularized networks, while Picual communities showed a sharply different behavior. In the rhizosphere, V. dahliae only irrupted in the microbial networks of Picual plants.

CONCLUSIONS

The belowground microbial communities of the two olive cultivars are very similar and pathogen introduction did not provoke significant alterations in their structure and functionality. However, notable differences were found in their networks in response to the inoculation. This phenomenon was more evident in the root endosphere communities. Thus, a correlation between modifications in the microbial networks of this microhabitat and susceptibility/tolerance to a soilborne pathogen was found. Moreover, V. dahliae irruption in the Picual microbial networks suggests a stronger impact on the belowground microbial communities of this cultivar upon inoculation. Our results suggest that changes in the co-occurrence interactions may explain, at least partially, the differential VWO susceptibility of the tested olive cultivars. Video abstract.

Defining the root endosphere and rhizosphere microbiomes from the World Olive Germplasm Collection

The bacterial and fungal communities from the olive (Olea europaea L.) root systems have not yet been simultaneously studied. We show in this work that microbial communities from the olive root endosphere are less diverse than those from the rhizosphere. But more relevant was to unveil that olive belowground communities are mainly shaped by the genotype of the cultivar when growing under the same environmental, pedological and agronomic conditions. Furthermore, Actinophytocola, Streptomyces and Pseudonocardia are the most abundant bacterial genera in the olive root endosphere, Actinophytocola being the most prevalent genus by far. In contrast, Gp6, Gp4, Rhizobium and Sphingomonas are the main genera in the olive rhizosphere. Canalisporium, Aspergillus, Minimelanolocus and Macrophomina are the main fungal genera present in the olive root system. Interestingly enough, a large number of as yet unclassified fungal sequences (class level) were detected in the rhizosphere. From the belowground microbial profiles here reported, it can be concluded that the genus Actinophytocola may play an important role in olive adaptation to environmental stresses. Moreover, the huge unknown fungal diversity here uncovered suggests that fungi with important ecological function and biotechnological potential are yet to be identified.

Indigenous Pseudomonas spp. Strains from the Olive (Olea europaea L.) Rhizosphere as Effective Biocontrol Agents against Verticillium dahliae: From the Host Roots to the Bacterial Genomes

The use of biological control agents (BCA), alone or in combination with other management measures, has gained attention over the past decades, driven by the need to seek for sustainable and eco-friendly alternatives to confront plant pathogens. The rhizosphere of olive (Olea europaea L.) plants is a source of bacteria with potential as biocontrol tools against Verticillium wilt of olive (VWO) caused by Verticillium dahliae Kleb. A collection of bacterial isolates from healthy nursery-produced olive (cultivar Picual, susceptible to VWO) plants was generated based on morphological, biochemical and metabolic characteristics, chemical sensitivities, and on their in vitro antagonistic activity against several olive pathogens. Three strains (PIC25, PIC105, and PICF141) showing high in vitro inhibition ability of pathogens' growth, particularly against V. dahliae, were eventually selected. Their effectiveness against VWO caused by the defoliating pathotype of V. dahliae was also demonstrated, strain PICF141 being the rhizobacteria showing the best performance as BCA. Genotypic and phenotypic traits traditionally associated with plant growth promotion and/or biocontrol abilities were evaluated as well (e.g., phytase, xylanase, catalase, cellulase, chitinase, glucanase activities, and siderophore and HCN production). Multi-locus sequence analyses of conserved genes enabled the identification of these strains as Pseudomonas spp. Strain PICF141 was affiliated to the “Pseudomonas mandelii subgroup,” within the “Pseudomonas fluorescens group,” Pseudomonas lini being the closest species. Strains PIC25 and PIC105 were affiliated to the “Pseudomonas aeruginosa group,” Pseudomonas indica being the closest relative. Moreover, we identified P. indica (PIC105) for the first time as a BCA. Genome sequencing and in silico analyses allowed the identification of traits commonly associated with plant-bacteria interactions. Finally, the root colonization ability of these olive rhizobacteria was assessed, providing valuable information for the future development of formulations based on these strains. A set of actions, from rhizosphere isolation to genome analysis, is proposed and discussed for selecting indigenous rhizobacteria as effective BCAs.

Tolerance of olive (Olea europaea) cv Frantoio to Verticillium dahliae relies on both basal and pathogen-induced differential transcriptomic responses

Verticillium wilt of olive (VWO) is one of the most serious biotic constraints for this tree crop. Our knowledge of the genetics of the tolerance/resistance to this disease is very limited. Here we show that tolerance of the cv Frantoio relies on both basal and early pathogen-induced differential transcriptomic responses. A comparative transcriptomic analysis (RNA-seq) was conducted in root tissues of cvs Frantoio (VWO-tolerant) and Picual (VWO-susceptible). RNA samples originated from roots of inoculated olive plants during the early infection stages by Verticillium dahliae (highly virulent, defoliating pathotype). A huge number of differentially expressed genes (DEGs) were found between 'Frantoio' and 'Picual' (27 312 unigenes) in the absence of the pathogen. Upon infection with V. dahliae, 'Picual' and 'Frantoio' plants responded differently too. In the early infection stages, four clusters of DEGs could be identified according to their time-course expression patterns. Among others, a pathogenesis-related protein of the Bet v I family and a dirigent-like protein involved in lignification, and several BAK1, NHL1, reactive oxygen species stress response and BAM unigenes showed noticeable differences between cultivars. Tolerance of 'Frantoio' plants to VWO is a consequence of a complex and multifaceted process which involves many plant traits.

Systemic responses in a tolerant olive (Olea europaea L.) cultivar upon root colonization by the vascular pathogen Verticillium dahliae

Verticillium wilt of olive (VWO) is caused by the vascular pathogen Verticillium dahliae. One of the best VWO management measures is the use of tolerant cultivars; however, our knowledge on VWO tolerance/resistance genetics is very limited. A transcriptomic analysis was conducted to (i) identify systemic defense responses induced/repressed in aerial tissues of the tolerant cultivar Frantoio upon root colonization by V. dahliae, and (ii) determine the expression pattern of selected defense genes in olive cultivars showing differential susceptibility to VWO. Two suppression subtractive hybridization cDNA libraries, enriched in up-regulated (FU) and down-regulated (FD) genes respectively, were generated from "Frantoio" aerial tissues. Results showed that broad systemic transcriptomic changes are taking place during V. dahliae-"Frantoio" interaction. A total of 585 FU and 381 FD unigenes were identified, many of them involved in defense response to (a)biotic stresses. Selected genes were then used to validate libraries and evaluate their temporal expression pattern in "Frantoio." Four defense genes were analyzed in cultivars Changlot Real (tolerant) and Picual (susceptible). An association between GRAS1 and DRR2 gene expression patterns and susceptibility to VWO was observed, suggesting that these transcripts could be further evaluated as markers of the tolerance level of olive cultivars to V. dahliae.

Early and delayed long-term transcriptional changes and short-term transient responses during cold acclimation in olive leaves

Low temperature severely affects plant growth and development. To overcome this constraint, several plant species from regions having a cool season have evolved an adaptive response, called cold acclimation. We have studied this response in olive tree (Olea europaea L.) cv. Picual. Biochemical stress markers and cold-stress symptoms were detected after the first 24 h as sagging leaves. After 5 days, the plants were found to have completely recovered. Control and cold-stressed plants were sequenced by Illumina HiSeq 1000 paired-end technique. We also assembled a new olive transcriptome comprising 157,799 unigenes and found 6,309 unigenes differentially expressed in response to cold. Three types of response that led to cold acclimation were found: short-term transient response, early long-term response, and late long-term response. These subsets of unigenes were related to different biological processes. Early responses involved many cold-stress-responsive genes coding for, among many other things, C-repeat binding factor transcription factors, fatty acid desaturases, wax synthesis, and oligosaccharide metabolism. After long-term exposure to cold, a large proportion of gene down-regulation was found, including photosynthesis and plant growth genes. Up-regulated genes after long-term cold exposure were related to organelle fusion, nucleus organization, and DNA integration, including retrotransposons.

The biocontrol endophytic bacterium Pseudomonas fluorescens PICF7 induces systemic defense responses in aerial tissues upon colonization of olive roots

Pseudomonas fluorescens PICF7, a native olive root endophyte and effective biocontrol agent (BCA) against Verticillium wilt of olive, is able to trigger a broad range of defense responses in root tissues of this woody plant. In order to elucidate whether strain PICF7 also induces systemic defense responses in above-ground organs, aerial tissues of olive plants grown under non-gnotobiotic conditions were collected at different time points after root bacterization with this endophytic BCA. A suppression subtractive hybridization (SSH) cDNA library, enriched in up-regulated genes, was generated. This strategy enabled the identification of 376 ESTs (99 contigs and 277 singlets), many of them related to response to different stresses. Five ESTs, involved in defense responses, were selected to carry out time-course quantitative real-time PCR (qRT-PCR) experiments aiming to: (1) validate the induction of these genes, and (2) shed light on their expression pattern along time (from 1 to 15 days). Induction of olive genes potentially coding for lipoxygenase 2, catalase, 1-aminocyclopropane-1-carboxylate oxidase, and phenylananine ammonia-lyase was thus confirmed at some time points. Computational analysis also revealed that different transcription factors were up-regulated in olive aerial tissues (i.e., JERF, bHLH, WRKY), as previously reported for roots. Results confirmed that root colonization by this endophytic bacterium does not only trigger defense responses in this organ but also mounts a wide array of systemic defense responses in distant tissues (stems, leaves). This sheds light on how olive plants respond to the "non-hostile" colonization by a bacterial endophyte and how induced defense response can contribute to the biocontrol activity of strain PICF7.

Molecular Variability Within and Among Verticillium dahliae Vegetative Compatibility Groups Determined by Fluorescent Amplified Fragment Length Polymorphism and Polymerase Chain Reaction Markers

ABSTRACT A degree of genetic diversity may exist among Verticillium dahliae isolates within vegetative compatibility groups (VCGs) that bears phytopathological significance and is worth investigating using molecular tools of a higher resolution than VCG characterization. The molecular variability within and among V. dahliae VCGs was studied using 53 artichoke isolates from eastern-central Spain, 96 isolates from cotton, 7 from cotton soil, and 45 from olive trees in countries of the Mediterranean Basin. Isolates were selected to represent the widest available diversity in cotton- and olive-defoliating (D) and -nondefoliating (ND) pathotypes, as well as for VCG. The VCG of 96 cotton and olive isolates was determined in this present study. Molecular variability among V. dahliae isolates was assessed by fluorescent amplified fragment length polymorphism (AFLP) analysis and by polymerase chain reaction (PCR) assays for DNA fragments associated with the D (462 bp) and ND (824 bp) pathotypes, as well as a 334-bp amplicon associated with D pathotype isolates but also present in some VCG2B isolates. Isolates from cotton were in VCG1A, VCG1B, VCG2A, VCG2B, and VCG4B and those from olive trees were in VCG1A, VCG2A, and VCG4B. Artichoke isolates included representatives of VCG1A, VCG2A, VCG2B (including a newly identified VCG2Ba), and VCG4B. AFLP data were used to generate matrixes of genetic distance among isolates for cluster analysis using the neighbor-joining method and for analysis of molecular variance. Results demonstrated that V. dahliae isolates within a VCG subgroup are molecularly similar, to the extent that clustering of isolates correlated with VCG subgroups regardless of the host source and geographic origin. VCGs differed in molecular variability, with the variability being highest in VCG2B and VCG2A. For some AFLP/VCG subgroup clusterings, V. dahliae isolates from artichoke grouped in subclusters clearly distinct from those comprising isolates from cotton and olive trees. In addition, VCG2B isolates from artichoke formed two distinct clusters that correlated with PCR markers of 334 bp (VCG2B(334)) or 824 bp (VCG2B(824)). Artichoke isolates in the VCG2B(334)/2beta(334) cluster were molecularly similar to isolates of VCG1A. The molecular difference found among artichoke isolates in VCG2B correlates with virulence of isolates to artichoke and cotton cultivars demonstrated in a previous study.