Skip to main content
  • Research article
  • Open access
  • Published:

Starvation-induced cell fusion and heterokaryosis frequently escape imperfect allorecognition systems in an asexual fungal pathogen



Asexual fungi include important pathogens of plants and other organisms, and their effective management requires understanding of their evolutionary dynamics. Genetic recombination is critical for adaptability and could be achieved via heterokaryosis — the co-existence of genetically different nuclei in a cell resulting from fusion of non-self spores or hyphae — and the parasexual cycle in the absence of sexual reproduction. Fusion between different strains and establishment of viable heterokaryons are believed to be rare due to non-self recognition systems. Here, we investigate the extent and mechanisms of cell fusion and heterokaryosis in the important asexual plant pathogen Verticillium dahliae.


We used live-cell imaging and genetic complementation assays of tagged V. dahliae strains to analyze the extent of non-self vegetative fusion, heterokaryotic cell fate, and nuclear behavior. An efficient CRISPR/Cas9-mediated system was developed to investigate the involvement of autophagy in heterokaryosis. Under starvation, non-self fusion of germinating spores occurs frequently regardless of the previously assessed vegetative compatibility of the partners. Supposedly “incompatible” fusions often establish viable heterokaryotic cells and mosaic mycelia, where nuclei can engage in fusion or transfer of genetic material. The molecular machinery of autophagy has a protective function against the destruction of “incompatible” heterokaryons.


We demonstrate an imperfect function of somatic incompatibility systems in V. dahliae. These systems frequently tolerate the establishment of heterokaryons and potentially the initiation of the parasexual cycle even between strains that were previously regarded as “incompatible.”


Generation of genetic and phenotypic diversity is of utmost importance for the survival and adaptability of species over evolutionary timescales [1, 2]. Genetic reshuffling in meiosis provides a constant supply of genetic substrates for natural selection and enables efficient purging of deleterious mutations from sexual populations [3, 4]. A complete lack of recombination in organisms such as obligate asexual species can result in the accumulation of problematic alleles, a phenomenon known as “Muller’s Ratchet” [5]. Some source of recombination is, therefore, usually deemed essential for the long-term survival of species.

In earlier days of fungal genetics, an alternative mechanism for genetic recombination, termed the parasexual cycle, was described in fungi [6]. This refers to a combination of non-sexual processes that lead to genetic diversification through the somatic cell and nuclear fusion, followed by mitotic recombination during ploidy reduction (i.e., haploidization) [6, 7]. Meiotic sexual cycles may have evolved from ancestral parasexual processes during the early eukaryotic evolution [8]. Variations of the parasexual cycle are present in several and diverse extant eukaryotes, i.e., many filamentous fungi and yeasts [9,10,11,12], protists [13, 14], and mammalian cells, where it can contribute to the phenotypic plasticity of tumor cell populations [15].

The parasexual cycle is normally initiated in heterokaryons, i.e., cells that host genetically distinct nuclei in a common cytoplasm. Heterokaryosis in fungi often arises from somatic non-self fusion (or anastomosis) between hyphae or spores of different strains [16, 17]. However, fusion between different individuals and establishment of viable heterokaryons are believed to be largely prevented in nature by non-self recognition systems, a phenomenon described as vegetative (or heterokaryon) incompatibility [18,19,20]. Incompatibility barriers can block cell communication to prevent fusion, or they can trigger post-fusion rejection mechanisms that result in the induction of programmed cell death in the incompatible heterokaryotic cell [19,20,21]. The versatile allorecognition mechanisms have been analyzed mostly in the sexual fungi Neurospora crassa, Podospora anserina, and Cryphonectria parasitica, where they are believed to function as defense mechanisms against the rapid spread of infectious elements and somatic parasitism [18, 19]. Fungal strains unable to form viable vegetative heterokaryons because of these barriers are referred to as vegetatively incompatible, whereas compatible strains that can fuse and establish stable heterokaryons are classified into the same Vegetative Compatibility Group (VCG) [22]. Despite the large perceived extent of incompatibility in natural populations [22, 23], many fungi seem competent in parasexual interactions in the laboratory [9,10,11], and the possibly heterokaryosis-mediated horizontal gene transfer has gained increasing attention in the postgenomic era as a driver of fungal evolution [24, 25].

The soil-borne asexual ascomycete Verticillium dahliae is the causal agent of Verticillium wilt, a disease that affects a wide range of plants and accounts for significant annual economic losses [26]. This haploid species is considered strictly asexual and propagates predominantly by clonal expansion [27]. Despite the apparent lack of sexual activity, V. dahliae populations exhibit signs of recombination, which, although not frequent, has the potential to create novel clonal lineages [28, 29]. The parasexual cycle has been described in Verticillium species and presents interesting features [30,31,32]. These include high instability of transient diploid nuclei, which are formed within heterokaryons following cell fusion, and frequent mitotic recombination during haploidization, often leading to recombined progeny [30, 33, 34]. Furthermore, overcoming cell wall-related barriers by protoplast fusion or nuclear microinjection resulted in interspecific heterokaryons between V. dahliae and V. albo-atrum at frequencies similar to intraspecific combinations [35]. Our previous finding that heterokaryosis in V. dahliae is occasionally possible even between strains of different VCGs (thus being considered incompatible) [36] raises the question whether parasexuality in asexual fungi may be more important than previously recognized. We recently established fusion between spores via conidial anastomosis tubes (CATs) as a convenient system for addressing this question [37].

In this work, our analyses of the extent and mechanisms of heterokaryosis in the strictly asexual fungus V. dahliae challenge the general assumption of absolute vegetative incompatibility barriers in asexual fungi. Frequent heterokaryon formation and transfer of genetic material even between supposedly “incompatible” partners, under non-selective conditions, suggest that parasexuality could play an important role in genomic and phenotypic diversification. Overall, our study paves the way for the elucidation of non-sexual genetic interactions in asexual fungi and the evolutionary significance of these processes in adaptation and survival.


Conidial fusion occurs frequently between “incompatible” V. dahliae strains under starvation

Starvation induces significantly CAT-mediated fusion of V. dahliae conidia/germlings [37], and this motivated us to investigate the possibility of conidial fusion between strains of different VCGs under such conditions. Seven representative strains (Additional file 1: Table S1 [36,37,38]; Fig. 1a) were labeled with sGFP- or mCherry-tagged histone H1 (nuclear localization) for their microscopic identification in pairings and, therefore, the discrimination of fusion between identical or unlike cells. In 14 pairings of the labeled strains, inter-strain fusion (Fig. 1b) was observed at significant levels (5.2-33.8% of fusions) not only in “compatible” pairings (i.e., same VCG), but also in all “incompatible” combinations (i.e., different VCGs) (no significant difference between “compatible” and “incompatible” pairings; Student’s t test, p value > 0.05; Fig. 1c). These results demonstrate that spontaneous cell fusion occurs frequently under starvation conditions between conidia or germlings of different strains, regardless of their previous VCG classification.

Fig. 1.
figure 1

Spontaneous non-self fusion via CATs occurs frequently between “incompatible” strains of V. dahliae under starvation. a Wild-type strains used in this experiment and their VCG classification. Single-letter strain codes (in brackets) denote derivative strains expressing sGFP- or mCherry-tagged histone H1 (details in Additional file 1: Table S1). The structure of the dendrogram reflects the phylogenetic relationships of VCGs as these were previously determined [39]. b Examples of CAT-mediated fusion between conidia/germlings of strains c-d (originating from the same wild-type strain), c–e (members of VCG 2B) and c–g (“incompatible” members of VCGs 2B and 4A, respectively). Bars = 5 μm. c Frequencies of intra- and inter-strain fusion in pairings of V. dahliae strains a–i. Control combinations indicated with an asterisk involved strains originating from the same wild-type strain. Each pairing was tested in triplicate, and 250-300 fusion events were analyzed per replicate. Bars = standard deviation (SD)

Heterokaryotic cells arisen from CAT-mediated fusion of “incompatible” strains often escape incompatibility-triggered cell death

We next asked whether fusion of “incompatible” conidia/germlings could lead to the establishment of viable heterokaryotic cells. To investigate this, we paired labeled strains (expressing fluorescently tagged histone H1 and/or cytoplasmic sGFP; Additional file 1: Table S1) and used live-cell imaging to monitor the cell fate and nuclear behavior in self, “compatible”, and “incompatible” fusions over time, up to 24 h after fusion. In total, we collected time-course data from 908 individual fusions, from three self (304 events), three “compatible” (80 events), and four “incompatible” (524 events) pairings (Additional file 2: Table S2). Overall, in this time window we identified four distinct types of cellular behavior following CAT-mediated fusion (Fig. 2a, b):

  1. I.

    Fusion led to a viable anastomosed cell. No nuclear migration occurred through the CAT, and septae were frequently formed within the CAT (70%). This behavior was frequently observed in self and “compatible” pairings (62.3%), but also in “incompatible” combinations (40.5%). In “compatible” pairings, cytoplasmic intermixing between the fused cells was always observed, whereas that was the case in 63.6% of “incompatible” cases (Additional file 3: Movie S1).

  2. II.

    Fusion was followed by a nuclear division in one of the interacting cells (20–150 min after fusion). One of the daughter nuclei then migrated through the CAT from the donor cell to its fusion partner (migration was completed in 10–20 min). Following a dikaryotic phase in the recipient cell (10 min–3 h), its resident (i.e., original) nucleus was degraded. This type was only detected in self-pairings (35.5%).

  3. III.

    Similar to (II), but degradation of the resident nucleus happened before migration of the other nucleus. This was often observed in “compatible” pairings (35.0%), less frequently in “incompatible” (7.4%), and only rarely in self-interactions. Nuclear migration through a CAT (types II and III) was often (59.3%) followed by formation of a septum in the CAT, regardless of the “compatibility” classification of the paired strains.

  4. IV.

    Fusion was followed by a cell death reaction observed in half of the “incompatible” pairings tested, and rarely in “compatible” interactions. The catastrophic reaction started 10–120 min after fusion, with or without obvious cytoplasmic mixing (Additional file 4: Fig. S1a, b). It was generally restricted to the fused cellular compartments without affecting the adjacent cells, and it typically involved thickening of the surrounding septae and cell walls, degradation of both nuclei and gradual cell shrinkage (Additional file 5: Movie S2; Additional file 4: Fig. S1c).

Fig. 2.
figure 2

“Incompatible” fused cells (via CATs) often escape incompatibility-triggered cell death. a Examples of the four observed types of cellular behavior following CAT-mediated fusion of “compatible” and “incompatible” cells. Strains shown: Ls.17 H1-mCherry paired with Cf.38 H1-sGFP (I), Ls.17 H1-sGFP (II), PH H1-sGFP (III), and BB H1-sGFP (IV). Arrows: nuclei undergoing degradation; arrowheads: migrating nuclei; asterisks: cell shrinkage. Bars = 5 μm. b Frequencies of the four types of cellular behavior in self, “compatible”, and “incompatible” pairings. The numbers of detected events in each case are provided in brackets (detailed results in Additional file 2: Table S2). c Frequency of nuclear migration through CATs in viable heterokaryons resulting from self, “compatible”, and “incompatible” pairings. Bars = SD. In b, c, statistical significance of differences between the compared groups was tested with one-way ANOVA followed by Tukey’s post hoc test; in b, groups marked by the same superscript character did not differ significantly; in c, ** p value ≤ 0.01; ns: non-significant (p value > 0.05)

Fatal incompatibility reactions following CAT-mediated fusion were essentially limited to “incompatible” fusions but affected only half of them. The considerable remaining fraction (47.9%) apparently escaped cell death and formed stable heterokaryotic cells without any signs of cellular degeneration in the analyzed time window, in all tested combinations of strains (Additional file 2: Table S2). Cell viability was confirmed by the addition at the end of each experiment of methylene blue, which generally fails to stain viable heterokaryons, in contrast to dead cells (Additional file 4: Fig. S1d). Nuclear translocation through the CAT was detected at significant levels in all types of pairings, including “incompatible” interactions (Fig. 2c). These results demonstrate that fusion of “incompatible” conidia/germlings can lead to the establishment of viable heterokaryons with the possibility of nuclear migration through their anastomosis bridge, at significant frequencies. Nuclear translocation is invariably associated with selective degradation of the nucleus of the recipient cell, similarly to previous observations in Fusarium oxysporum [40], both in “compatible” and in “incompatible” fusions (Fig. 2b).

“Incompatible” fusions that escape cell death support heterokaryotic growth in colonies

Since CAT-mediated fusion of “incompatible” strains can often lead to viable heterokaryotic cells, we hypothesized that these could further support hyphal growth to the formation of heterokaryotic colonies. We tested this using a collection of 19 nitrate non-utilizing (nit) strains of all V. dahliae VCGs, as well as two strains of the related species V. nonalfalfae (Additional file 1: Table S1). When hyphae of “incompatible” complementary nit strains are confronted on minimal medium (MM) that selects for prototrophy, typically no heterokaryotic growth occurs due to incompatibility barriers (Fig. 3a). We performed 110 pairings using this hyphal-based traditional assay (Fig. 3a, b), in comparison to an assay that permitted CAT-mediated cell fusion of the same pairs of strains (Fig. 3c, d). The results were drastically different (Fig. 3b, d; Additional file 6: Table S3), with the CAT-based method yielding 5.5 times more positive results than the hyphal-based assay, in combinations of all VCGs. In total, 49.0% of the “incompatible” pairings tested positive only in the CAT-based assay. The latter yielded several interspecific responses between V. dahliae and V. nonalfalfae, as well as interactions involving a heterokaryon self-incompatible (HSI) strain (Fig. 3d).

Fig. 3.
figure 3

Conidial/germling fusion (via CATs) supports stable growth of “incompatible” heterokaryotic colonies. The traditional (hyphal-based) VCG (a) and the new CAT-based (c) complementation assays for the detection of heterokaryotic growth are described schematically. The results from the two assays are summarized in b and d, respectively (see Additional file 6: Table S3 for details). A single-letter code has been assigned to each strain (details of strains are provided in Additional file 1: Table S1). HSI: heterokaryon self-incompatible. Vna: V. nonalfalfae. Each pairing (in both assays) was performed in at least three replicates. Pairings that produced very limited prototrophic growth in only one of the replicates were regarded as negative. Examples of strong and weak heterokaryons (21 days old) are shown for each method (right), as well as their morphological characteristics after subculturing them on fresh MM (arrows, last row) in d. Bars = 2 cm

Most “incompatible” interactions were weak, exhibiting sparse aerial mycelium, tufts of hyphae, or localized groups of microsclerotia in at least two out of the three replicates, as opposed to the strong interactions with denser prototrophic growth in all three replicates (Fig. 3b, d). We sought, therefore, to validate the heterokaryotic nature of such colonies by excluding the possibilities of reversion to prototrophy or cross-feeding of non-fused cells. When 30 randomly selected weak “incompatible” heterokaryons were analyzed, all retained the ability to grow slowly over successive reculturings on fresh MM, with denser growth in their central regions, irregular colony shape with frequent morphological sectoring, and uneven pigment distribution, which are all typical features of heterokaryons (Fig. 3d). We then plated approx. 104 uninucleate conidia of each heterokaryon on MM and we never detected prototrophic growth, but only nit-characteristic sparse growth. These results confirmed the heterokaryotic nature of the hyphal masses. Our experiments suggest that large numbers of fusions occur even between “incompatible” strains and that sufficiently many of them escape incompatibility-triggered cell death to support prototrophic growth of the resulting colonies.

Heterokaryosis arises from mixed infections in planta

With the exception of its dormant microsclerotia in the soil, V. dahliae completes its life cycle in its plant host [26]. To examine whether the parasexual cycle in nature could be initiated by heterokaryon formation within plants, we subjected eggplant seedlings to mixed artificial infections (Fig. 4a) with “compatible” (Ls.17 × Cf.38) and “incompatible” (Ls.17 × PH; Ls.17 × BB; Ls.17 × Ca.146) combinations of strains, which were marked with complementary antibiotic-resistance markers (Additional file 1: Table S1). Fungal material was re-isolated from stem sections of the infected plants 40 days after inoculation. Presumably heterokaryotic growth was detected on double-selective PDA medium (amended with both antibiotics) in 16 out of 425 tested xylem chips (3.8%), from all “compatible” and “incompatible” combinations, whereas no growth was observed in any of our single-infection or mock-infection controls (Fig. 4b, c). All 16 strains retained their ability to grow slowly over successive reculturings on double-selective medium, but no growth was detected when approx. 104 uninucleate conidia of each sample were plated on that medium. We also used resistance cassette-specific primers for PCR screening of 14 of these heterokaryons, which confirmed their heterokaryotic nature. Given the complete suppression of conidial fusion on PDA and the previously demonstrated possibility of conidial fusion in a xylem sap-simulating medium [37], our results suggest that the isolated “compatible” and “incompatible” heterokaryons arose from fusion between strains in planta.

Fig. 4.
figure 4

Heterokaryon formation in planta following mixed infection. a, b Eggplant seedlings were artificially infected with mixed inocula of “compatible” and “incompatible” combinations of V. dahliae strains. The pathogen was re-isolated from stem sections (from eight xylem chips along the stem of each plant) of 17 infected plants per strain combination, 40 days after inoculation. From each treated plant, three xylem chips were transferred to PDA (controls of re-isolation), while the remaining five were transferred to double-selective PDA supplemented with hygB and G418, to select for heterokaryotic growth. Mock-infected plants (with water) and single-strain infections were used as controls. c Similarly to b, for seedlings infected with mixed inocula of strains. d Summary of tested pairings and results

Optimization of an efficient CRISPR/Cas9-mediated gene targeting system for the analysis of autophagy in V. dahliae heterokaryons

The characteristics of incompatibility-triggered cell death and selective nuclear degradation following cell fusion (Fig. 2; Additional file 4: Fig. S1) indicated that the autophagic machinery of the cells may be an important component of the mechanisms that regulate the heterokaryotic cell fate decision and nuclear interaction. To investigate this, we first searched the V. dahliae genome for homologs of the core autophagy genes and discovered that V. dahliae possesses a full autophagy gene complement (Additional file 7: Table S4). We selected the hallmark autophagy genes atg1 (VDAG_05745) and atg8 (VDAG_01225), which are essential for the induction of autophagy and formation of autophagosomes, respectively [41], to investigate the possible roles of autophagy in the post-fusion processes.

To overcome practical limitations that complicate reverse genetic investigations in V. dahliae, we used the AMA1 replicator sequence from Aspergillus nidulans [42] to ensure autonomous plasmid maintenance in a CRISPR/Cas9-based system [43] for efficient gene targeting in V. dahliae. For plasmid delivery to the fungal cells we developed a lithium acetate/heat-shock-mediated method for conidial transformation. This technique yielded lower transformation efficiencies than our standard protoplast-based method (up to 17 and 106 transformants/μg of plasmid, respectively), but it proved much faster and more convenient. Using this method, we transformed V. dahliae Ls.17 with a plasmid containing the hygromycin B-resistance (hph) selection marker, and we studied plasmid stability over 14 successive reculturings of 25 randomly selected transformants. While full stability was observed on selective medium, a gradual decrease in plasmid stability occurred in the absence of selection (Additional file 4: Fig. S2a). This pattern proved very convenient, as it offered clear-cut selection of transformants and maintenance of strains under selective conditions, but also straightforward elimination of the plasmid containing the Spcas9 endonuclease gene and the selection marker, after the targeted genetic modification had been achieved.

We then constructed an AMA1-containing plasmid for the CRISPR/Cas9-mediated targeting of the N-terminal region of V. dahliae atg8 (Fig. 5a). In this plasmid, the constitutive PgpdA promoter was used for the pol-II-mediated transcription of the sgRNA cassette. The sgRNA was flanked by the ribozymes HH and HDV, which ensure the post-transcriptional release of the sgRNA without modifications (Fig. 5a). The plasmid was used to co-transform V. dahliae Ls.17 together with a linear repair substrate consisting of either the sgfp or the mCherry coding sequence, flanked by homologous arms corresponding to the genomic regions that surround the start codon of atg8 (Fig. 5b). Using this method, we achieved in situ seamless (i.e., without introducing any additional sequence) marker-free tagging of atg8 (Fig. 5c), with efficiencies that ranged between 12 and 66% of transformants for 100-1000-bp-long homologous flanks, respectively (Additional file 4: Fig. S2b). Properly tagged strains were validated using PCR and microscopic analyses that confirmed the induction of autophagy by starvation and treatment with rapamycin, which is a known inducer of autophagy [44] (Additional file 4: Fig. S2c), as well as the previously described involvement of autophagy in conidial germination [45] (Additional file 4: Fig. S2d).

Fig. 5.
figure 5

Schematic representation of CRISPR/Cas9-mediated gene targeting in V. dahliae. a Both SpCas9 and sgRNA cassettes were introduced into the same, AMA1-containing plasmid. The sgRNA was flanked by the hammerhead (HH) and hepatitis delta virus (HDV) ribozymes, which enable the release of the mature sgRNA after their self-catalysis (scissors indicate cleavage points). b Homologous repair substrates consist of either a fluorescent protein gene (for gene tagging) or an antibiotic-resistance cassette (for gene disruption), flanked by homologous arms of variable length. c Co-transformation of the AMA1-containing plasmid with the homologous repair substrate for gene tagging or disruption. The Cas9/sgRNA complex induces in the genomic locus targeted DSBs (scissors), which can be repaired using the homologous repair substrate to result in seamless gene tagging or gene disruption. In the latter case, stop codons (asterisk) were introduced in the homologous arm downstream of the antibiotic-resistance cassette, to ensure the knockout of the gene (example shown: V. dahliae atg8)

Next, we used the same system for inactivation of V. dahliae atg8 by gene disruption. For this, the G418-resistance (neoR) cassette, flanked by 100-bp-long homologous arms (Fig. 5b), was used for co-transformation of V. dahliae Ls.17 with the plasmid containing the sgRNA cassette (Fig. 5a). To ensure inactivation of the gene following homologous recombination, stop codons were introduced in the arm downstream of the neoR cassette (Fig. 5c). This approach for gene knockout proved faster than the commonly used Agrobacterium tumefaciens-based system (ATMT) (Additional file 4: Fig. S3a), which we used to delete atg1. In both cases, knockout mutants were verified by PCR and Southern hybridization analyses (Additional file 4: Fig. S3b, c). The CRISPR/Cas9-based system yielded correctly disrupted mutants at a relatively low frequency (1.1% of transformants) using 100-bp-long arms, which conveniently overcomes the need for multistep vector construction.

Autophagy is involved in post-fusion selective nuclear degradation but is not required for cell fusion and acts against incompatibility-triggered cell death

In order to gain further insight into the roles of autophagy in heterokaryosis, we first constructed a V. dahliae H1-mCherry sGFP-Atg8 strain (Ls.17 background), which we used for time-lapse imaging of CAT-mediated self-fusion. Following anastomosis, we invariably observed co-localization of nuclei undergoing selective degradation (Figs. 2a and 6a, b; Additional file 4: Fig. S1c, S4; Additional file 8: Movie S3) with Atg8-containing organelles. These were ranging from small autophagosomes associated with parts of the nuclear chromatin (Additional file 4: Fig. S4) to larger and sometimes ring-like structures that appeared to engulf the whole nucleus before its degradation (Fig. 6a; Additional file 8: Movie S3). During the Atg8-nuclear interaction, the Atg8-containing vesicles were often dynamically undergoing fusion and fission (Additional file 4: Fig. S4).

Fig. 6.
figure 6

Autophagy is involved in selective nuclear degradation, but not in cell fusion or incompatibility-triggered death. a Live-cell imaging of a self-fusion (Ls.17 H1-mCherry sGFP-Atg8). A ring-like structure with accumulated sGFP-Atg8 (arrows) surrounds the sequestered nucleus (arrowheads). b Live-cell imaging of an “incompatible” fusion (Ls.17 H1-mCherry sGFP-Atg8 × BB H1-sGFP). Localized accumulation of sGFP-Atg8 near the fusion point (arrows) shortly before the incompatibility reaction. Arrowheads: nucleus undergoing degradation. c Frequency of active conidia and their fraction involved in CAT-mediated fusion (left; n = 300 conidia per replicate), non-apical hyphal compartments and fused cells with more than one nuclei (middle; n = 300 hyphal cells or 150 anastomoses per replicate), and inviable fusions determined by staining with methylene blue (right; n = 150 anastomoses per replicate). Wild-type (wt) pairing: Ls.17 H1-mCherry × PH H1-sGFP; pairing of autophagy-deficient mutants (Δatg1): Ls.17 Δatg1 × PH Δatg1. Each strain/pairing was tested in triplicate. Statistical significance of differences was tested with one-way ANOVA followed by Tukey’s post hoc test (left and middle; bars with the same letter do not differ significantly, p value > 0.05) or Student’s t-test (right; * p ≤ 0.05). Error bars: SD. d, e Multinucleate cells (arrowheads) arise frequently in Δatg1 and Δatg8 autophagy-deficient mutants from cell fusion (d) or sub-apical nuclear division (e), in contrast to the wild type that exhibits strictly uninucleate organization. Calcofluor white was used for cell wall staining. Scale bars = 5 μm

In our autophagy-deficient mutants Δatg1 and Δatg8, conidial germination and CAT-mediated fusion occurred at wild-type levels (Fig. 6c), in contrast to F. oxysporum where deletion of atg8 was previously reported to impair its fusion frequency [46]. Both V. dahliae mutants exhibited significantly more bi- and multi-nucleate fusions than the wild type (Fig. 6c, d; Additional file 9: Movie S4), whereas gene complementation restored the wild-type phenotype in both cases (Additional file 4: Fig. S5). Similarly, bi-/multi-nucleate cells arose frequently from nuclear division in Δatg1 and Δatg8 hyphal compartments (Fig. 6c, e; Additional file 10: Movie S5), in contrast to the strictly uninucleate organization of the wild-type mycelium (Fig. 6c, e), which is ensured by degradation of one of the daughter nuclei of each non-apical nuclear division (Additional file 11: Movie S6). Overall, our results suggest that the nuclear number is strictly regulated in V. dahliae to one per cell, and this process is mediated by the autophagic machinery.

We then tested whether autophagy is involved in the post-fusion incompatibility-triggered reaction, which displays morphological features reminiscent of autophagic programmed cell death (Fig. 2a; Additional file 4: Fig. S1). Time-lapse imaging of “incompatible” fusions between strains Ls.17 H1-mCherry sGFP-Atg8 and BB H1-sGFP revealed in all cases localized accumulation of Atg8 in the vicinity of the fusion point, prior to the onset of the catastrophic reaction (Fig. 6b). Co-localization of nuclei undergoing degradation with Atg8-containing vacuoles was also observed in all “incompatible” fusions (Fig. 6b). To investigate whether autophagy mediates the incompatibility reaction, we also deleted atg1 from strain PH and analyzed a PH Δatg1 × Ls.17 Δatg1 pairing. Staining of inviable fusions with methylene blue (Additional file 4: Fig. S1d) revealed a significantly higher frequency of cell death in the autophagy-deficient pair than in the wild-type pair (Fig. 6c), while typical morphological characteristics of incompatibility were observed in affected heterokaryons in both cases. Therefore, although autophagy is induced before the manifestation of the incompatibility reaction, it is not required for cell death but actually protects some “incompatible” fusions from destruction. This is in line with previous studies of P. anserina that proposed a possibly protective role of autophagy against incompatibility-triggered cell death [47, 48].

Nuclear interactions in “compatible” and “incompatible” heterokaryons

Co-existence of genetically distinct nuclei in heterokaryons is a requirement for the parasexual cycle and may enable horizontal gene transfer, but the mechanisms involved in these nuclear interactions are not well understood. Despite our many attempts, we were unable to capture direct fusion of nuclei in our heterokaryons, which is explained by the very low frequencies of diploidization (10−6–10−8) in V. dahliae heterokaryons [33]. However, in a few cases we observed gradual accumulation of sGFP signal in mCherry-labeled nuclei or vice versa, both in “compatible” and in “incompatible” pairings (Fig. 7a). In these cases, transfer was strictly unidirectional and not always correlated temporally with nuclear division. Although the underlying mechanism is unknown, such interactions provide evidence for genetic transfer between nuclei. In addition, we analyzed CAT-mediated self-fusions of Ls.17 Δatg1 cells to examine whether their higher frequencies of binucleate cells facilitate the microscopic study of nuclear fusion (Fig. 6c). In rare cases, temporal interaction of nuclei possibly indicating their fusion was detected in this strain (Fig. 7b).

Fig. 7.
figure 7

Nuclear interactions in “compatible” and “incompatible” V. dahliae heterokaryons. a Gradual accumulation of sGFP signal in mCherry-labeled nuclei in heterokaryotic cells between strains Ls.17 H1-mCherry × Cf.38 H1-sGFP (“compatible”), and Ls.17 H1-mCherry × PH H1-sGFP (“incompatible”). Arrowheads: recipient nuclei. b Putative nuclear fusion (arrowhead) in strain Ls.17 Δatg1. Bars in a, b = 5 μm. c Schematic representation of the possible fates of a non-self cell fusion in V. dahliae. Fused cells can remain viable or be subjected to an incompatibility-triggered reaction. Division of rare hypothetical fused nuclei in viable heterokaryotic cells could generate distinct nuclear lineages in mosaic colonies. d Nuclei of “incompatible” (Ls.17 H1-mCherry × BB sGFP, Ls.17 H1-mCherry × PH sGFP) and self-pairing control (Ls.17 H1-mCherry × Ls.17 sGFP) heterokaryons following growth under non-selective conditions. Bars = 10 μm

As long as they are able to divide, putatively fused nuclei could form distinct nuclear lineages in chimeric heterokaryotic colonies (Fig. 7c). To investigate whether this is possible in “incompatible” interactions, we scanned microscopically heterokaryotic colonies formed via CATs, under non-selective conditions, between strains with sGFP- or mCherry-labeled nuclei. Both our self-pairing control (Ls.17) and the “incompatible” fusions (Ls.17 × BB, Ls.17 × PH) gave rise to mosaic colonies with mixed populations of nuclei with sGFP, mCherry, or both signals (Fig. 7d). These colonies mainly consisted of regions with parental properties, interspersed with sectors or individual hyphae (and uninucleate conidia) with mixed fluorescent signal, which is in agreement with the mosaic nature of V. dahliae heterokaryons described earlier [31]. Our results are consistent with the possibility of parasexual interactions in “incompatible” heterokaryons of V. dahliae.


Despite the evolutionary benefits and persistence of sexual reproduction across the eukaryotic tree of life [4, 5, 49], loss of sexual competence is observed in members of several eukaryotic lineages, including approx. one fifth of described fungi (i.e., the mitosporic fungi) that either appear to be strictly asexual or are suspected to engage only rarely in some cryptic form of sex [50, 51]. These organisms are very important to plant pathology and their effective management depends on our ability to predict successfully the modes and pace of their evolution. Central questions in this regard are what are the relative contributions, mechanisms, and spatiotemporal patterns of clonal reproduction and recombination in these organisms [50]. Here, we used V. dahliae, an asexual pathogen, to gain insights into the extent and mechanisms of heterokaryosis, an essential component of parasexuality, which enables genomic recombination in the absence of sex and meiosis [7, 9, 23, 52]. We found that starvation induces fusion between conidia or germlings of different strains, regardless of their traditionally assessed vegetative compatibility [22]. Cell fusion can often lead to the establishment of viable heterokaryons, where nuclei can interact to initiate the parasexual cycle.

The VCGs of asexual fungi, assessed by forcing hyphae of complementary auxotrophic mutants to anastomose on minimal medium, have been widely regarded as incompatible groups [22, 53]. However, cross (or “bridging”) inter-VCG reactions have been reported in some species [54,55,56,57,58], and their genetic analysis in V. dahliae previously revealed that such interactions can be due to heterokaryosis [36]. In addition, studies of CAT-mediated fusion of Colletotrichum lindemuthianum [59] and F. oxysporum [60] suggested that incompatibility may be suppressed at the conidial/germling stage of development. However, our previous investigation of conidial pairings of V. dahliae strains in minimal medium showed excellent congruence with the traditional hyphal pairings on agar medium [36], providing no evidence of incompatibility suppression during conidial germination in V. dahliae.

We recently discovered that starvation of cells induces significantly CAT-mediated fusion of V. dahliae conidia/germlings [37]. Here, we show that this happens regardless of their VCG classification, even between species, and that inter-VCG fusions often form viable heterokaryotic cells and colonies. No selection pressure for heterokaryons was exercised in our experiments, in sharp contrast to the traditional VCG assays that force heterokaryosis between auxotrophic strains under non-physiological selection conditions. In addition, we demonstrated that heterokaryosis can result from mixed infections in planta, which constitute the most likely opportunity for encounter of non-dormant V. dahliae strains in nature [61]. Overall, our study demonstrates that the traditionally assessed incompatibility systems do not constitute absolute barriers to heterokaryon formation in V. dahliae, where a fraction of starved conidia can escape allorecognition systems, anastomose, and lead to the formation of viable heterokaryons, regardless of their VCG classification. Our preliminary experiments revealed that starvation similarly induces anastomosis between mature hyphae of strains on water agar. A presumed extension of the imperfect nature of incompatibility systems, observed here in fusion of conidia, to hyphal fusion could explain the frequently observed cross-VCG hyphal interactions in this species [54, 58, 62,63,64]. It would also suggest that it is the environmental conditions (e.g., starvation) rather than the developmental stage that induces anastomosis and hypothetically parasexual phenomena even between previously considered “incompatible” strains. Consistently, mycelial compatibility assessed by the analysis of barrage zone formation between V. dahliae strains was previously shown to differ significantly from their VCG classification [65], providing further support to the hypothesis that somatic compatibility in this species depends on the environmental conditions and the level of selective pressure.

Multiple checkpoint systems have evolved in the sexual ascomycete N. crassa to essentially prevent heterokaryosis between genetically distinct individuals [20, 21]. Our study indicates that partner recognition and chemotropic attraction preceding cell fusion in V. dahliae are not subject to the same genetic control as in N. crassa, where an allorecognition system prevents fusion between incompatible individuals [66]. Nevertheless, we found evidence for the existence of two downstream allorecognition checkpoint systems following cell fusion. First, a fraction of “incompatible” interactions failed to achieve cytoplasmic mixing, which may have been prevented by a system similar to that of N. crassa, where allelic differences between individuals at the cwr locus can cause cell fusion arrest by blocking dissolution of the cell wall [67]. Secondly, in a fraction of “incompatible” interactions, fusion triggered a cell death reaction characterized by hyphal compartmentalization, thickening of the cell wall due to increased deposition of chitin [68], nuclear degradation, and cell shrinkage, therefore resembling the incompatibility reactions described in other fungi [69,70,71]. Such reactions are mediated by allelic interactions between sec-9/plp-1 [72] and rcd-1 [73] in N. crassa, as well as different repertoires of multiple het loci among fungi [18, 19, 70]. However, a fundamental difference in V. dahliae is that a significant fraction of “incompatible” conidial fusions escapes allorecognition control under starvation, ensuring the formation of viable heterokaryons, in agreement with similar observations in the other asexual fungi F. oxysporum [60] and C. lindemuthianum [59].

Overall, our data favor the hypothesis of a trade-off between two mutually exclusive processes that both provide advantages to the adaptive evolution of asexual fungi, such as V. dahliae. On the one hand, genetic intermixing between individuals is subject to allorecognition control, which triggers cell death in the majority of “incompatible” fusions, possibly to protect these syncytial organisms from the spread of infectious or parasitic elements [19, 20]. On the other hand, limited availability of nutrients induces cell fusion [37], and a fraction of the resulting heterokaryons overcome the imperfect incompatibility barriers hypothetically to grant populations with access to parasexuality and, therefore, genetic diversification that can be advantageous under stressful conditions. Even low levels of parasexual genomic reshuffling could have a marked effect on the genetic structure of fungal populations [74] and could explain the signs of recombination observed in V. dahliae populations [28, 29], as well as the imperfect clonality of its VCGs [39]. Parasexuality could also be critical for interspecific hybridization [75] and horizontal gene transfer [24, 25], which may have significant evolutionary roles. This is supported by our evidence of interspecific interactions between V. dahliae and V. nonalfalfae, similarly to those recently reported between Colletotrichum species [76], and the amphidiploid nature of the related hybrid species V. longisporum, presumed to have arisen parasexually between V. dahliae and other species [77, 78].

Our study provides evidence that the molecular machinery of autophagy might be involved in the proposed trade-off as a link between starvation and the cell fate decision in V. dahliae heterokaryons. We found that Atg1- and Atg8-mediated autophagy participates in the incompatibility reaction, but not as a mediator of cell death. Instead, autophagy seems to act against destruction of “incompatible” heterokaryons, which is in line with previous studies of P. anserina that revealed accelerated cell death in autophagy-deficient mutants [48]. The autophagic machinery could function to protect “incompatible” fusions, e.g., by degrading incompatibility effectors, and survival of each heterokaryon may depend on the success of this process. Since cell fusion is coupled with starvation, which also induces autophagy [41], the majority of fused cells must have high levels of autophagic activity and this could function in favor of their survival.

Previous studies of the parasexual cycle [6, 30] and the hypotheses of parasexual horizontal gene transfer [24, 25] predict exchange of genetic material between distinct nuclei within heterokaryons. Consistently, we observed frequent nuclear migration through the anastomosis bridge both in “compatible” and in “incompatible” heterokaryons, commonly generating hyphae with genetically different nuclei. The resulting dikaryons, however, were eventually subjected to a nuclear division-triggered selective degradation process targeting the host nucleus of the affected cellular compartment, as previously also observed in F. oxysporum [40]. This process involves Atg1- and Atg8-mediated macroautophagy of the whole nucleus (i.e., macronucleophagy), which acts as a homeostatic mechanism ensuring a uninucleate status, both in hyphae and in fused cells of V. dahliae. This mechanism appears conserved between uninucleate fungi, as suggested by similar findings in F. oxysporum [46].

Despite the eventual degradation of the resident nucleus, however, the preceding phase of nuclear co-existence in the heterokaryon could provide a sufficient window of opportunity for genetic exchanges. One possibility is nuclear fusion, a defining component of the sensu stricto parasexual cycle, originally described in A. nidulans [6] and later in several other fungi [10] including Verticillium species [30, 33, 34]. Those studies inferred nuclear fusion within heterokaryons to the formation of heterozygous diploids using segregation analyses of auxotrophic and drug-resistant genetic markers, as well as analyses of spore sizes as a proxy for ploidy levels [79]. Single-spore heterozygous diploid cultures were isolated and shown to generate frequently aneuploid and haploid sectors, presumably due to chromosomal non-disjunction, often leading to progeny with recombined markers [6, 7, 10, 30, 33]. However, the estimated frequency of parasexual nuclear fusion in those studies (in the range of 10−6–10−8 in the case of Verticillium species [30, 33, 34]) is prohibitively low for capturing the process in real time by microscopy, which also prevented us in the present study from doing so. Apart from direct nuclear fusion, other mechanisms of nuclear interaction might also facilitate smaller-scale so-called horizontal gene or chromosome transfers (HGT or HCT) in heterokaryons, for example uptake of particular chromosomes or chromatin remnants of the degraded nucleus by the surviving one. Such a mechanism could possibly be involved in the transfer of dispensable or supernumerary chromosomes reported in C. gloeosporioides [80] and F. oxysporum [81, 82]. We speculate that the previously observed fractionated nuclei (“micronuclei”) [60, 83], often observed in highly vacuolated fungal cells, could mediate such transfers, in a process hypothetically analogous to the lateral transfer of oncogenes in mammalian cells by apoptotic bodies [84].


This report contributes to a better understanding of the evolution of asexual fungi. Our findings have important implications for the interpretation of the traditional VCG data, which are commonly discussed under the assumption of genetic isolation between VCGs in asexual species. We propose the evolution of imperfect incompatibility systems in V. dahliae to preserve the ability of the species to engage in parasexuality. This conclusion suggests important roles of heterokaryosis and parasexuality in fungal evolution, and it welcomes future research aiming at a better understanding of the extent and modes of parasexual recombination achieved through this process. We propose V. dahliae as a particularly suitable organism for such analyses, and we expect future experimentation to be facilitated by the methodological toolkit for CRISPR/Cas9-mediated gene targeting that we optimized in this study to overcome experimental limitations in this and many other fungi.


Strains and culture conditions

All V. dahliae strains and plasmids constructed and used in this study are listed in Additional file 1: Table S1 and Additional file 12: Table S5 [43, 59, 68, 85,86,87,88,89], respectively. Their background wild-type isolates were previously classified into Vegetative Compatibility Groups (VCGs) based on traditional complementation assays between nitrate non-utilizing (nit) mutants [36,37,38]. Preparation, handling, and maintenance of monoconidial strains, as well as culture media and conditions have been previously described [37, 58].

Labeling of strains with fluorescent markers

Nuclei of V. dahliae strains were labeled with either sGFP- or mCherry-tagged histone H1. Carboxy-terminal tagging of V. dahliae histone H1 with the red fluorescent protein mCherry was performed using either a fusion PCR strategy that we have previously described [37] or a binary vector (pVV19) suitable for Agrobacterium tumefaciens-mediated transformation (ATMT).

The first method involved the generation of two PCR amplicons using the genomic DNA of V. dahliae isolate Ls.17 as template, containing the endogenous promoter and the coding region of the histone H1 gene (VDAG_09854), and its transcription termination sequence, respectively. The mCherry coding sequence and the phleomycin-resistance (bleR) cassette were amplified from plasmid pAN8.1-mCherry [40]. The two genomic amplicons and the mCherry-containing cassette were mixed in a PCR reaction without primers, which was followed by a second reaction that included the first PCR product as template and the primer pair VdH1FusF/VdH1FusR (Additional file 13: Table S6). The fusion PCR product was then used to transform protoplasts of V. dahliae strains, according to our previously reported protocol [68].

Alternatively, we constructed plasmid pVV19 for mCherry-tagging of histone H1 by amplifying the following fragments with the Herculase II Fusion DNA polymerase (Agilent, Santa Clara, CA, USA) (all PCR primers are listed in the Additional file 13: Table S6): (i) a 1972-bp-long fragment from genomic DNA of V. dahliae strain Ls.17, containing the coding region of the histone H1 gene (VDAG_09854) excluding the stop codon, and a 1000-bp-long upstream genomic fragment presumably including the endogenous promoter; (ii) a 717-bp-long fragment from plasmid pMaM330 [86] containing the mCherry gene; (iii) a 489-bp-long fragment from plasmid pFC332 [43] containing the transcription termination sequence of the A. nidulans tef1 gene; and (iv) a 1648-bp-long fragment from plasmid pBS-genR [68] containing the G418-resistance (neoR) cassette. The four PCR amplicons were subjected to an in vitro DNA assembly reaction in the backbone of plasmid pOSCAR [88] using the NEBuilder HiFi DNA Assembly Master Mix kit (New England Biolabs, Ipswich, MA, USA), according to the manufacturer’s recommendations. The resulting vector was used for ATMT of V. dahliae strains, according to our recently reported protocol [68]. Nuclear labeling of V. dahliae strains with heterologous sGFP-tagged histone H1 was performed using plasmid pMF357 [59]. This contains a fusion construct of the sgfp gene to the Neurospora crassa histone H1 gene and the hph cassette that confers resistance to hygromycin B. This plasmid was used for transformation of V. dahliae protoplasts [68].

To study cytoplasmic mixing following CAT-mediated cell fusion, V. dahliae strains were constructed to express sGFP with cytoplasmic localization. For this, we used plasmid pIGPAPA [89], which carries the sgfp gene under the control of the N. crassa ICL promoter and the hph cassette, to transform protoplasts of the desired V. dahliae strains [68].

Analysis of CAT-mediated cell fusion

Our methods for analysis and quantification of CAT-mediated fusion of V. dahliae conidia/germlings have recently been described [37]. Briefly, fresh conidial suspensions (in water) were always prepared from 7-day-old PDA cultures and diluted in CAT medium (0.75 g/l β-glycerophosphate disodium salt · 5H2O) to a final concentration of 1.0 × 106 conidia per ml. Suspensions were then mixed (1:1), and 100 μl of their mixture were transferred to a well of a 96-well glass-bottom plate (MGB096-1-2-LG-L from Matrical Bioscience, Spokane, WA, USA). The plate was incubated at 24 °C (in the dark) for 60 h before imaging. Each pairing was performed in triplicate, and 250-300 fusion events were analyzed per replicate. Staining with 0.005% methylene blue (Sigma-Aldrich, St. Louis, MO, USA) (incubation at 24 °C for 5 min) was used for quantifying cell viability of fused cells.


Image acquisition was performed using a Nikon (Tokyo, Japan) Ti-E epifluorescence microscope equipped with an autofocus system (Perfect Focus System, Nikon), a LED light engine (Spectra X from Lumencor, Beaverton, OR, USA), filter sets 390/18 and 435/48, 469/35 and 525/50, and 542/27 and 600/52 (excitation and emission, respectively; all from Semrock, Rochester, NY, USA, except for 525/50, which was from Chroma Technology, Bellows Falls, VT, USA), and a sCMOS camera (Flash4.0 from Hamamatsu, Honshu, Japan). At least 20 non-overlapping fields of view were imaged per well. Calcofluor white M2R (Sigma-Aldrich, St. Louis, MO, USA) was used at a final concentration of 10 μg/ml (incubation at 24 °C for 5 min) for staining of chitin in cell walls and septae. In time-lapse experiments, the imaging plates were scanned for image acquisition at 10 min-intervals for up to 24 h (exposure time: 20 ms for the blue channel, 50 ms for the green and red channels) and were incubated undisturbed at 24 °C (in the dark between image acquisitions). In some experiments we performed z-stack (step size = 0.5 μm) time-lapse imaging at each position, for up to 8 h. Image processing was performed using ImageJ [90]. Images were adjusted to a uniform contrast across all time points for each experiment, and the maximum intensity projection method was used for processing z-stacks.

Hyphal- and CAT-based complementation assays

We used two complementation methods for assessing the degree of vegetative genetic isolation between V. dahliae strains, namely the traditional VCG (hyphal-based) assay and a newly developed assay that involves CAT-mediated fusion of conidia/germlings of the interacting strains. For both methods, we used the same set of previously characterized nitrate non-utilizing nit mutants (nit1 and nitM; with reversion rates lower than 1.0 × 10−7) of a collection of V. dahliae isolates of all VCGs (Additional file 1: Table S1) [36,37,38].

In the VCG assay, sparse mycelial growth of complementary nit mutants on co-inoculated minimal medium (MM) plates enables mycelial confrontation of the paired isolates on the agar. When this interaction leads to heterokaryotic prototrophic growth along the confrontation zone (Fig. 3a), the corresponding strains are considered “compatible” and assigned to the same Vegetative Compatibility Group (VCG) [54]. We performed all VCG assays on Czapek-Dox MM agar [37] and scored plates for prototrophic growth four weeks after inoculation, according to standard procedures and previously described criteria [36, 58]. Each pairing was performed in triplicate and all used nit mutants were also inoculated separately on MM agar (in triplicate) to serve as controls.

We also developed and assessed in this study a CAT-based complementation assay of nit mutants. This method is based on pairing conidia of the desired V. dahliae strains under conditions that are favorable for CAT formation, and then transferring the samples to MM to select for heterokaryotic prototrophic growth (Fig. 3c). For this, we prepared a fresh conidial suspension (in water) from a 7-day-old PDA culture of each strain, and we diluted each suspension in CAT medium (0.75 g/l β-glycerophosphate disodium salt · 5H2O) to a final concentration of 1.0 × 106 conidia/ml. The desired suspensions were then mixed (1:1) for each pairing, and 100 μl of each mixture were transferred to a well of a polystyrene 96-well flat-bottom plate (2620/S from Kartell, Noviglio, Italy), with a cellophane disk placed at the bottom of the well. The plates were incubated at 24 °C (in the dark) for 72 h before the cellophane disks were aseptically transferred to the surface of Czapek-Dox MM agar. The plates were scored for prototrophic growth on the cellophane disks over four weeks after inoculation, while they were incubated at 24 °C (in the dark). Each pairing was performed in triplicate and all tested nit mutants were also subjected to the same procedure individually (in triplicate) as controls.

Presumed heterokaryons generated by both methods were randomly selected for further characterization to validate their heterokaryotic nature and exclude the possibilities of prototrophic growth because of marker reversion or cross-feeding of non-fused individuals [32]. For this, they were sub-cultured on selective medium (MM), conidial suspensions were prepared from these cultures, and approx. 1.0 × 104 conidia were spread on fresh MM. Conidia of V. dahliae are always unicellular and predominantly uninucleate (99.7% in the wild-type isolate Ls.17 according to our observations), while the rare bi- and multi-nucleate conidia are considered strictly homozygous due to their origin from the single nucleus of a phialide [34, 91]. Therefore, a heterokaryon is expected to grow on selective medium due to genetic complementation between its genetically distinct nuclei, whereas its conidia would be of one of the parental types, and would thus remain unable to grow under selection for prototrophy.

Mixed artificial infection of eggplant seedlings

Eggplant seedlings at the one true leaf-stage were drenched with 20 ml of mixed conidial suspensions (5.0 × 106 conidia per ml) or sterile water (mock-infection controls) [92]. Plants were incubated at 25 °C with a 12-h light-dark cycle. Fungal re-isolation was performed from sections along the stem of each treated plant (from eight xylem chips per plant) 40 days after inoculation, according to described procedures [92]. Three xylem chips from each plant were transferred to acidified PDA as controls, while the remaining five chips were transferred to acidified double-selective PDA (10 μg/ml hygB and 100 μg/ml G418) to select for heterokaryotic growth.

Identification of autophagy genes

The protein sequences of 41 genes involved in autophagy (Additional file 7: Table S4) were retrieved from the Saccharomyces Genome Database ( and the Ensembl Fungi database ( and used as queries in tBLASTn searches of the V. dahliae Ls.17 genome (NCBI Genome assembly GCA_000952015.1) to detect homologs of core autophagy genes, applying previously described criteria [93]. For reasons of comparison, we included in the Additional file 7: Table S4 data for various classes of the fungal phylogeny, adopted from [93].

Lithium acetate (LiAc)/heat shock-mediated transformation of V. dahliae conidia

At least 108 conidia of each V. dahliae strain were harvested from PDA plates, pelleted, resuspended in 30 ml of a phosphate-KCl buffer (0.6 M KCl, 1/15 M KH2PO4, and 1/15 M Na2HPO4, pH 5.8; supplemented with 25 mM DTT) and incubated for 2 h at 30 °C (50 rpm). Conidia were then washed once with water and once with a 0.1 M LiAc solution, before being resuspended in 0.5 ml of 0.1 M LiAc and incubated at room temperature for 20 min. Aliquots (100 μl) of this suspension were mixed with 5–10 μg of plasmid, 240 μl of 50% PEG-4000 and 35 μl of 1 M LiAc, and they were incubated on ice for 30 min before being transferred to 34 °C for 10 min (heat-shock). Each sample was then washed once and plated on a PDA medium overlaid with a cellophane sheet. Following incubation at 25 °C for 24 h, the sheet was transferred onto PDA supplemented with hygromycin B (15 μg/ml) for selection of transformants.

CRISPR/Cas9-mediated gene tagging and knockout

The components and the strategy of our CRISPR/Cas9-based method, adopted from [43], are presented in Fig. 5. An AMA1-containing plasmid (pVV27) with the SpCas9 cassette, an atg8-specific sgRNA cassette (for targeting V. dahliae atg8, VDAG_01225) and the hph selection marker was generated by the insertion of two PCR fragments into plasmid pFC332 [43]. For this, two PCR fragments were amplified from plasmid pFC334 using the Herculase II Fusion DNA Polymerase with primers containing the protospacer sequence of atg8 (Additional file 13: Table S6). The backbone plasmid pFC332 was linearized with PacI digestion, and the three fragments were assembled into plasmid pVV27 using the NEBuilder HiFi DNA Assembly Master Mix kit. Repair substrates with 1.0-kb-long homologous arms (Fig. 5b, c) were constructed by cloning the sgfp or mCherry gene (amplified from plasmid pIGPAPA or pMaM330, respectively; Additional file 12: Table S5) and the desired amplified flanks from genomic DNA of V. dahliae Ls.17 into plasmids pVV25 and pVV26, respectively (Additional file 12: Table S5), using the NEBuilder HiFi DNA Assembly Master Mix kit. Repair substrates with 0.1 kb-long homologous arms (Fig. 5b, c) were obtained directly by PCR amplification from plasmid pIGPAPA, pMaM330 or pSD1 (for the sgfp, mCherry or neoR gene, respectively; Additional file 12: Table S5) using primers containing the desired homologous arms (Additional file 13: Table S6). Co-transformation of V. dahliae strains with plasmid pVV27 and a repair substrate was performed using a protoplast-mediated method [68] or the LiAc/heat shock-mediated method. Transformants were screened by PCR and the fluorescently tagged strains were confirmed by microscopy. To induce autophagy in these experiments, we used rapamycin (Sigma-Aldrich, St. Louis, MO, USA) at a final concentration of 500 nM.

Gene deletion using Agrobacterium tumefaciens-mediated transformation (ATMT) and complementation of knockout strains

The atg1 homolog of V. dahliae (VDAG_05745) was deleted by the ATMT-based method (Additional file 4: Fig. S3a) that we have previously described [68]. Knockout strains were complemented by re-introducing the corresponding wild-type genes [68].

Statistical analyses

Statistical significance of differences between multiple compared groups was tested using the two-sided analysis of variance (ANOVA) test, following the evaluation of homogeneity of variance across samples using the Levene’s test (F-test, p ≤ 0.05). The Tukey’s honestly significant difference (HSD) test was performed to assess group mean separation. Mean comparison between two groups was performed using the two-sided Student’s t-test. All results are presented as mean value ± SD. Sample sizes and number of replicates are provided in the legends of the corresponding figures.

Availability of data and materials

The data supporting the findings of this study are presented within the manuscript and the additional supporting files of this article. Raw microscopy data are available from the corresponding author on reasonable request.


  1. Booy G, Hendriks RJJ, Smulders MJM, Groenendael JMV, Vosman B. Genetic Diversity and the Survival of Populations. Plant Biol. 2000;2(4):379–95.

    Article  Google Scholar 

  2. Frankham R. Genetics and extinction. Biol Conserv. 2005;126(2):131–40.

    Article  Google Scholar 

  3. Fisher RA. The genetical theory of natural selection. Oxford, England: Clarendon Press; 1930.

    Book  Google Scholar 

  4. Muller HJ. Some Genetic Aspects of Sex. Am Nat. 1932;66(703):118–38.

    Article  Google Scholar 

  5. Felsenstein J. The evolutionary advantage of recombination. Genetics. 1974;78(2):737–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Pontecorvo G, Roper JA. Genetic analysis without sexual reproduction by means of polyploidy in Aspergillus nidulans. J Gen Microbiol. 1952;6:7.

    Google Scholar 

  7. Pontecorvo G. The Parasexual Cycle in Fungi. Annu Rev Microbiol. 1956;10(1):393–400.

    Article  CAS  PubMed  Google Scholar 

  8. Goodenough U, Heitman J. Origins of Eukaryotic Sexual Reproduction. Cold Spring Harb Perspect Biol. 2014;6(3):a016154.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Strom NB, Bushley KE. Two genomes are better than one: history, genetics, and biotechnological applications of fungal heterokaryons. Fungal Biol Biotechnol. 2016;3(1):4.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Tinline RD, MacNeill BH. Parasexuality in plant pathogenic fungi. Annu Rev Phytopathol. 1969;7(1):147–68.

    Article  Google Scholar 

  11. Anderson MZ, Thomson GJ, Hirakawa MP, Bennett RJ. A ‘parameiosis’ drives depolyploidization and homologous recombination in Candida albicans. Nat Commun. 2019;10(1):4388.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zhao Y, Wang Y, Upadhyay S, Xue C, Lin X. Activation of Meiotic Genes Mediates Ploidy Reduction during Cryptococcal Infection. Curr Biol. 2020;30:1387–1396.e5.

    Article  CAS  Google Scholar 

  13. Gaunt MW, Yeo M, Frame IA, Stothard JR, Carrasco HJ, Taylor MC, et al. Mechanism of genetic exchange in American trypanosomes. Nature. 2003;421(6926):936–9.

    Article  CAS  PubMed  Google Scholar 

  14. Carpenter ML, Assaf ZJ, Gourguechon S, Cande WZ. Nuclear inheritance and genetic exchange without meiosis in the binucleate parasite Giardia intestinalis. J Cell Sci. 2012;125(Pt 10):2523–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Miroshnychenko D, Baratchart E, Ferrall-Fairbanks MC, Velde RV, Laurie MA, Bui MM, et al. Spontaneous cell fusions as a mechanism of parasexual recombination in tumour cell populations. Nat Ecol Evol. 2021;5(3):379–91.

    Article  PubMed  Google Scholar 

  16. Mela AP, Rico-Ramírez AM, Glass NL. Syncytia in Fungi. Cells. 2020;9(10):2255.

    Article  CAS  PubMed Central  Google Scholar 

  17. Roca MG, Read ND, Wheals AE. Conidial anastomosis tubes in filamentous fungi. FEMS Microbiol Lett. 2005;249(2):191–8.

    Article  CAS  Google Scholar 

  18. Saupe SJ. Molecular genetics of heterokaryon incompatibility in filamentous ascomycetes. Microbiol Mol Biol Rev. 2000;64(3):489–502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Paoletti M. Vegetative incompatibility in fungi: From recognition to cell death, whatever does the trick. Fungal Biol Rev. 2016;30(4):152–62.

    Article  Google Scholar 

  20. Daskalov A, Heller J, Herzog S, Fleißner A, Glass NL. Molecular Mechanisms Regulating Cell Fusion and Heterokaryon Formation in Filamentous Fungi. Microbiol Spectr. 2017;5(2).

  21. Gonçalves AP, Glass NL. Fungal social barriers: to fuse, or not to fuse, that is the question. Commun Integr Biol. 2020;13(1):39–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Leslie JF. Fungal Vegetative Compatibility. Annu Rev Phytopathol. 1993;31(1):127–50.

    Article  CAS  PubMed  Google Scholar 

  23. Clutterbuck AJ. Parasexual recombination in fungi. J Genet. 1996;75(3):281–6.

    Article  Google Scholar 

  24. Fitzpatrick DA. Horizontal gene transfer in fungi. FEMS Microbiol Lett. 2012;329(1):1–8.

    Article  CAS  PubMed  Google Scholar 

  25. Soanes D, Richards TA. Horizontal Gene Transfer in Eukaryotic Plant Pathogens. Annu Rev Phytopathol. 2014;52(1):583–614.

    Article  CAS  PubMed  Google Scholar 

  26. Pegg GF, Brady BL. Verticillium wilts. Wallingford, Oxon, UK: CABI Pub; 2002. Accessed 1 Nov 2020.

    Book  Google Scholar 

  27. Gurung S, Short DPG, Atallah ZK, Subbarao KV. Clonal Expansion of Verticillium dahliae in Lettuce. Phytopathology®. 2014;104:641–9.

    Article  CAS  Google Scholar 

  28. Atallah ZK, Maruthachalam K, du Toit L, Koike ST, Michael Davis R, Klosterman SJ, et al. Population analyses of the vascular plant pathogen Verticillium dahliae detect recombination and transcontinental gene flow. Fungal Genet Biol. 2010;47(5):416–22.

    Article  CAS  PubMed  Google Scholar 

  29. Milgroom MG, del Mar Jiménez-Gasco M, García CO, Drott MT, Jiménez-Díaz RM. Recombination between Clonal Lineages of the Asexual Fungus Verticillium dahliae Detected by Genotyping by Sequencing. PLoS ONE. 2014;9:e106740.

    Article  Google Scholar 

  30. Hastie AC. The parasexual cycle in Verticillium albo-atrum. Genet Res. 1964;5(2):305–15.

    Article  Google Scholar 

  31. Puhalla JE, Mayfield JE. The Mechanism of Heterokaryotic Growth in VERTICILLIUM DAHLIAE. Genetics. 1974;76(3):411–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Typas MA, Heale JB. Heterokaryosis and the role of cytoplasmic inheritance in dark resting structure formation in Verticillium spp. Mol Gen Genet MGG. 1976;146(1):17–26.

    Article  Google Scholar 

  33. Typas MA, Heale JB. Heterozygous diploid analyses via the parasexual cycle and a cytoplasmic pattern of inheritance in Verticillium spp. Genet Res. 1978;31(2):131–44.

    Article  Google Scholar 

  34. Hastie AC. Phialide analysis of mitotic recombination inVerticillium. Mol Gen Genet MGG. 1968;102(3):232–40.

    Article  CAS  PubMed  Google Scholar 

  35. Typas MA. Heterokaryon Incompatibility and Interspecific Hybridization between Verticillium albo-atrum and Verticillium dahliae Following Protoplast Fusion and Microinjection. Microbiology. 1983;129(10):3043–56.

    Article  Google Scholar 

  36. Papaioannou IA, Typas MA. High-Throughput Assessment and Genetic Investigation of Vegetative Compatibility in Verticillium dahliae. J Phytopathol. 2015;163(6):475–85.

    Article  CAS  Google Scholar 

  37. Vangalis V, Knop M, Typas MA, Papaioannou IA. Establishment of conidial fusion in the asexual fungus Verticillium dahliae as a useful system for the study of non-sexual genetic interactions. Curr Genet. 2021;67(3):471–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Papaioannou IA, Dimopoulou CD, Typas MA. “Cryptic” group-I introns in the nuclear SSU-rRNA gene of Verticillium dahliae. Curr Genet. 2014;60(3):135–48.

    Article  CAS  PubMed  Google Scholar 

  39. Papaioannou IA, Dimopoulou CD, Typas MA. Structural and phylogenetic analysis of the rDNA intergenic spacer region of Verticillium dahliae. FEMS Microbiol Lett. 2013;347(1):23–32.

    Article  CAS  PubMed  Google Scholar 

  40. Ruiz-Roldán MC, Köhli M, Roncero MIG, Philippsen P, Di Pietro A, Espeso EA. Nuclear dynamics during germination, conidiation, and hyphal fusion of Fusarium oxysporum. Eukaryot Cell. 2010;9(8):1216–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Feng Y, He D, Yao Z, Klionsky DJ. The machinery of macroautophagy. Cell Res. 2014;24(1):24–41.

    Article  CAS  PubMed  Google Scholar 

  42. Gems D, Johnstone IL, Clutterbuck AJ. An autonomously replicating plasmid transforms Aspergillus nidulans at high frequency. Gene. 1991;98(1):61–7.

    Article  CAS  PubMed  Google Scholar 

  43. Nødvig CS, Nielsen JB, Kogle ME, Mortensen UH. A CRISPR-Cas9 System for Genetic Engineering of Filamentous Fungi. PLoS ONE. 2015;10(7):e0133085.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Mizushima N. Autophagy: process and function. Genes Dev. 2007;21(22):2861–73.

    Article  CAS  PubMed  Google Scholar 

  45. Kikuma T, Ohneda M, Arioka M, Kitamoto K. Functional Analysis of the ATG8 Homologue Aoatg8 and Role of Autophagy in Differentiation and Germination in Aspergillus oryzae. Eukaryot Cell. 2006;5(8):1328–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Corral-Ramos C, Roca MG, Di Pietro A, Roncero MIG, Ruiz-Roldán C. Autophagy contributes to regulation of nuclear dynamics during vegetative growth and hyphal fusion in Fusarium oxysporum. Autophagy. 2015;11(1):131–44.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Pinan-Lucarré B, Paoletti M, Dementhon K, Coulary-Salin B, Clavé C. Autophagy is induced during cell death by incompatibility and is essential for differentiation in the filamentous fungus Podospora anserina: Autophagy, differentiation and cell death in Podospora. Mol Microbiol. 2003;47(2):321–33.

    Article  PubMed  Google Scholar 

  48. Pinan-Lucarré B, Balguerie A, Clavé C. Accelerated Cell Death in Podospora Autophagy Mutants. Eukaryot Cell. 2005;4(11):1765–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Otto SP, Lenormand T. Resolving the paradox of sex and recombination. Nat Rev Genet. 2002;3(4):252–61.

    Article  CAS  PubMed  Google Scholar 

  50. Taylor JW, Jacobson DJ, Fisher MC. THE EVOLUTION OF ASEXUAL FUNGI: Reproduction, Speciation and Classification. Annu Rev Phytopathol. 1999;37(1):197–246.

    Article  CAS  PubMed  Google Scholar 

  51. Dyer PS, Kück U. Sex and the Imperfect Fungi. Microbiol Spectr. 2017;5(3).

  52. Bennett RJ. The parasexual lifestyle of Candida albicans. Curr Opin Microbiol. 2015;28:10–7.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Puhalla JE. Classification of strains of Fusarium oxysporum on the basis of vegetative compatibility. Can J Bot. 1985;63(2):179–83.

    Article  Google Scholar 

  54. Joaquim TR, Rowe RC. Reassessment of Vegetative Compatibility Relationships Among Strains of Verticillium dahliae Using Nitrate-Nonutilizing Mutants. Phytopathology. 1990;80(11):1160–6.

    Article  Google Scholar 

  55. Molnár A, Sulyok L, Hornok L. Parasexual recombination between vegetatively incompatible strains in Fusarium oxysporum. Mycol Res. 1990;94(3):393–8.

    Article  Google Scholar 

  56. Kerényi Z, Táborhegyi É, Pomázi A, Hornok L. Variability amongst strains of Fusarium poae assessed by vegetative compatibility and RAPD polymorphism. Plant Pathol. 1997;46(6):882–9.

    Article  Google Scholar 

  57. Viji G, Uddin W, O’Neill NR, Mischke S, Saunders JA. Genetic Diversity of Sclerotinia homoeocarpa Isolates from Turfgrasses from Various Regions in North America. Plant Dis. 2004;88(11):1269–76.

    Article  CAS  PubMed  Google Scholar 

  58. Papaioannou IA, Ligoxigakis EK, Vakalounakis DJ, Markakis EA, Typas MA. Phytopathogenic, morphological, genetic and molecular characterization of a Verticillium dahliae population from Crete, Greece. Eur J Plant Pathol. 2013;136(3):577–96.

    Article  Google Scholar 

  59. Ishikawa FH, Souza EA, Shoji J, Connolly L, Freitag M, Read ND, et al. Heterokaryon Incompatibility Is Suppressed Following Conidial Anastomosis Tube Fusion in a Fungal Plant Pathogen. PLoS ONE. 2012;7(2):e31175.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Shahi S, Beerens B, Bosch M, Linmans J, Rep M. Nuclear dynamics and genetic rearrangement in heterokaryotic colonies of Fusarium oxysporum. Fungal Genet Biol. 2016;91:20–31.

    Article  CAS  PubMed  Google Scholar 

  61. Wheeler DL, Johnson DA. Does Coinoculation with Different Verticillium dahliae Genotypes Affect the Host or Fungus? Phytopathology. 2019;109:780–6.

    Article  CAS  Google Scholar 

  62. El-Bebany AF, Alkher H, Adam LR, Daayf F. Vegetative compatibility of Verticillium dahliae isolates from potato and sunflower using nitrate non-utilizing (nit) mutants and PCR-based approaches. Can J Plant Pathol. 2013;35(1):1–9.

    Article  CAS  Google Scholar 

  63. Hiemstra JA, Rataj-Guranowska M. Vegetative Compatibility Groups in Verticillium dahliae Isolates from the Netherlands as Compared to VCG Diversity in Europe and in the USA. Eur J Plant Pathol. 2003;109(8):827–39.

    Article  CAS  Google Scholar 

  64. Elena P. Vegetative compatibility groups within Verticillium dahliae isolates from different hosts in Greece. Plant Pathol. 1998;47(5):635–40.

    Article  Google Scholar 

  65. Papaioannou IA, Typas MA. Barrage formation is independent from heterokaryon incompatibility in Verticillium dahliae. Eur J Plant Pathol. 2015;141(1):71–82.

    Article  Google Scholar 

  66. Heller J, Zhao J, Rosenfield G, Kowbel DJ, Gladieux P, Glass NL. Characterization of Greenbeard Genes Involved in Long-Distance Kind Discrimination in a Microbial Eukaryote. PLoS Biol. 2016;14(4):e1002431.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Gonçalves AP, Heller J, Span EA, Rosenfield G, Do HP, Palma-Guerrero J, et al. Allorecognition upon Fungal Cell-Cell Contact Determines Social Cooperation and Impacts the Acquisition of Multicellularity. Curr Biol. 2019;29:3006–3017.e3.

    Article  Google Scholar 

  68. Vangalis V, Papaioannou IA, Markakis EA, Knop M, Typas MA. Hex1, the Major Component of Woronin Bodies, Is Required for Normal Development, Pathogenicity, and Stress Response in the Plant Pathogenic Fungus Verticillium dahliae. J Fungi. 2020;6(4):344.

    Article  CAS  Google Scholar 

  69. Dementhon K, Paoletti M, Pinan-Lucarré B, Loubradou-Bourges N, Sabourin M, Saupe SJ, et al. Rapamycin Mimics the Incompatibility Reaction in the Fungus Podospora anserina. Eukaryot Cell. 2003;2(2):238–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Glass NL, Dementhon K. Non-self recognition and programmed cell death in filamentous fungi. Curr Opin Microbiol. 2006;9(6):553–8.

    Article  CAS  PubMed  Google Scholar 

  71. Gonçalves AP, Heller J, Daskalov A, Videira A, Glass NL. Regulated Forms of Cell Death in Fungi. Front Microbiol. 2017;8.

  72. Heller J, Clavé C, Gladieux P, Saupe SJ, Glass NL. NLR surveillance of essential SEC-9 SNARE proteins induces programmed cell death upon allorecognition in filamentous fungi. Proc Natl Acad Sci. 2018;115(10):E2292–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Daskalov A, Gladieux P, Heller J, Glass NL. Programmed Cell Death in Neurospora crassa Is Controlled by the Allorecognition Determinant rcd-1. Genetics. 2019;213(4):1387–400.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Hoekstra RF. Population genetics of filamentous fungi. Antonie Van Leeuwenhoek. 1994;65(3):199–204.

    Article  CAS  PubMed  Google Scholar 

  75. Steensels J, Gallone B, Verstrepen KJ. Interspecific hybridization as a driver of fungal evolution and adaptation. Nat Rev Microbiol. 2021:1–16.

  76. Mehta N, Baghela A. Quorum sensing-mediated inter-specific conidial anastomosis tube fusion between Colletotrichum gloeosporioides and C. siamense. IMA Fungus. 2021;12:7.

    Article  Google Scholar 

  77. Karapapa VK, Bainbridge BW, Heale JB. Morphological and molecular characterization of Verticillium longisporum comb, nov., pathogenic to oilseed rape. Mycol Res. 1997;101(11):1281–94.

    Article  Google Scholar 

  78. Inderbitzin P, Davis RM, Bostock RM, Subbarao KV. The Ascomycete Verticillium longisporum Is a Hybrid and a Plant Pathogen with an Expanded Host Range. PLoS ONE. 2011;6(3):e18260.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Typas MA, Heale JB. Analysis of Ploidy Levels in Strains of Verticillium using a Coulter Counter. Microbiology. 1977;101:177–80.

    Google Scholar 

  80. He C, Rusu AG, Poplawski AM, Irwin JA, Manners JM. Transfer of a supernumerary chromosome between vegetatively incompatible biotypes of the fungus Colletotrichum gloeosporioides. Genetics. 1998;150(4):1459–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Ma L-J, van der Does HC, Borkovich KA, Coleman JJ, Daboussi M-J, Di Pietro A, et al. Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature. 2010;464(7287):367–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Vlaardingerbroek I, Beerens B, Rose L, Fokkens L, Cornelissen BJC, Rep M. Exchange of core chromosomes and horizontal transfer of lineage-specific chromosomes in Fusarium oxysporum. Environ Microbiol. 2016;18(11):3702–13.

    Article  CAS  PubMed  Google Scholar 

  83. Roca MG, Davide LC, Mendes-Costa MC, Wheals A. Conidial anastomosis tubes in Colletotrichum. Fungal Genet Biol. 2003;40(2):138–45.

    Article  PubMed  Google Scholar 

  84. Bergsmedh A, Szeles A, Henriksson M, Bratt A, Folkman MJ, Spetz A-L, et al. Horizontal transfer of oncogenes by uptake of apoptotic bodies. Proc Natl Acad Sci. 2001;98(11):6407–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Lu S, Lyngholm L, Yang G, Bronson C, Yoder OC, Turgeon BG. Tagged mutations at the Tox1 locus of Cochliobolus heterostrophus by restriction enzyme-mediated integration. Proc Natl Acad Sci U S A. 1994;91(26):12649–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Khmelinskii A, Meurer M, Ho C-T, Besenbeck B, Füller J, Lemberg MK, et al. Incomplete proteasomal degradation of green fluorescent proteins in the context of tandem fluorescent protein timers. Mol Biol Cell. 2016;27(2):360–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Nguyen QB, Kadotani N, Kasahara S, Tosa Y, Mayama S, Nakayashiki H. Systematic functional analysis of calcium-signalling proteins in the genome of the rice-blast fungus, Magnaporthe oryzae, using a high-throughput RNA-silencing system. Mol Microbiol. 2008;68(6):1348–65.

    Article  CAS  PubMed  Google Scholar 

  88. Paz Z, García-Pedrajas MD, Andrews DL, Klosterman SJ, Baeza-Montañez L, Gold SE. One Step Construction of Agrobacterium-Recombination-ready-plasmids (OSCAR), an efficient and robust tool for ATMT based gene deletion construction in fungi. Fungal Genet Biol. 2011;48(7):677–84.

    Article  CAS  PubMed  Google Scholar 

  89. Horwitz BA, Sharon A, Lu S-W, Ritter V, Sandrock TM, Yoder OC, et al. A G Protein Alpha Subunit fromCochliobolus heterostrophusInvolved in Mating and Appressorium Formation. Fungal Genet Biol. 1999;26(1):19–32.

    Article  CAS  PubMed  Google Scholar 

  90. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–82.

    Article  CAS  PubMed  Google Scholar 

  91. Hastie AC. Mitotic Recombination in Conidiophores of Verticillium albo-atrum. Nature. 1967;214(5085):249–52.

    Article  Google Scholar 

  92. Markakis EA, Ligoxigakis EK, Avramidou EV, Tzanidakis N. Survival, Persistence, and Infection Efficiency of Verticillium dahliae Passed Through the Digestive System of Sheep. Plant Dis. 2014;98(9):1235–40.

    Article  CAS  PubMed  Google Scholar 

  93. Wang Q, Liu H, Xu H, Hei R, Zhang S, Jiang C, et al. Independent losses and duplications of autophagy-related genes in fungal tree of life: ATG genes in fungal kingdom. Environ Microbiol. 2019;21(1):226–43.

    Article  CAS  PubMed  Google Scholar 

Download references


Plasmids pFC332 and pFC334 (U.H. Mortensen, Technical University of Denmark) were obtained from Addgene.


This research is co-financed by Greece and the European Union (European Social Fund-ESF) through the Operational Programme «Human Resources Development, Education and Lifelong Learning» in the context of the project “Strengthening Human Resources Research Potential via Doctorate Research” (MIS-5000432), implemented by the State Scholarships Foundation (IKY). Open Access funding enabled and organized by Projekt DEAL.

Author information

Authors and Affiliations



Conceptualization: IAP and MAT; Methodology: IAP, MAT, MK, and VV; Investigation: VV, IAP, and IL; Formal analysis: VV and IAP; Supervision: IAP and MAT; Visualization: IAP and VV; Writing – Original Draft: IAP and VV; Writing – review and editing: IAP, MAT, and VV; Funding acquisition: MK and MAT. All authors read and approved the final manuscript.

Authors’ information

Twitter handles:

Vasileios Vangalis @VVangalis

Ioannis A. Papaioannou @ioannis_ap

Corresponding author

Correspondence to Ioannis A. Papaioannou.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Table S1. List of strains used in this study.

Additional file 2.

Table S2. Results of time-lapse imaging experiments.

Additional file 3.

Movie S1. Establishment of cytoplasmic continuity following CAT-mediated fusion in an “incompatible” pairing (PH sGFP × BB H1-sGFP).

Additional file 4

Supplementary Figures. Figure S1: Typical features of incompatibility-triggered cell death in V. dahliae. Figure S2: Fluorescent tagging of V. dahliae Atg8 by a CRISPR/Cas9-based system. Figure S3: Construction and validation of V. dahliae atg1 and atg8 knockout mutants. Figure S4: Atg8-mediated nuclear degradation. Figure S5: Frequencies of multinucleate cells in complemented Δatg1-c and Δatg8-c strains.

Additional file 5.

Movie S2. Time-lapse imaging of a fusion event between the “incompatible” strains PH sGFP × BB H1-sGFP (cell wall staining: calcofluor white).

Additional file 6.

Table S3. Results of complementation assays.

Additional file 7

Table S4. Verticillium dahliae homologs of autophagy genes.

Additional file 8.

Movie S3. Time-lapse imaging of a self-fusion of strain Ls.17 H1-mCherry sGFP-atg8. The host nucleus is surrounded by a ring-like Atg8-containing structure and eventually degraded.

Additional file 9.

Movie S4. Time-lapse imaging of a self-fusion of strain Ls.17 H1-mCherry Δatg1. Both nuclei of the dikaryon remain viable following nuclear migration.

Additional file 10.

Movie S5. Nuclear division in a hyphal compartment of strain Ls.17 H1-mCherry Δatg8 results in the formation of a binucleate cell.

Additional file 11.

Movie S6. Nuclear division in a non-apical hyphal compartment of strain Ls.17 H1-mCherry is followed by degradation of one of the daughter nuclei. Only one nucleus remains viable in the compartment.

Additional file 12.

Table S5. List of plasmids.

Additional file 13.

Table S6. List of PCR primers.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vangalis, V., Likhotkin, I., Knop, M. et al. Starvation-induced cell fusion and heterokaryosis frequently escape imperfect allorecognition systems in an asexual fungal pathogen. BMC Biol 19, 169 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: