# New Cytoplasmic Virus-Like Elements (VLEs) in the Yeast Debaryomyces hansenii

Xymena Polomska, Cécile Neuvéglise, Joanna Zyzak, Barbara Żarowska, Serge Casaregola, Zbigniew Lazar

► To cite this version:

Xymena Polomska, Cécile Neuvéglise, Joanna Zyzak, Barbara Żarowska, Serge Casaregola, et al. New Cytoplasmic Virus-Like Elements (VLEs) in the Yeast Debaryomyces hansenii. Toxins, 2021, 13 (9), pp.615. 10.3390/toxins13090615. [https://doi.org/10.33431/teh.03343156](https://doi.org/10.33431/teh.03343156)

HAL Id: hal-03343156

https://hal.inrae.fr/hal-03343156

Submitted on 10 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# New Cytoplasmic Virus-Like Elements (VLEs) in the Yeast *Debaryomyces hansenii*

Xymena Połomska <fnref n="1" />,* ID, Cécile Neuvéglise <fnref n="2" />, Joanna Zyzak <fnref n="3" />, Barbara Żarowska <fnref n="1" />, Serge Casaregola <fnref n="4" /> and Zbigniew Lazar <fnref n="1" /> ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Połomska, X.; Neuvéglise, C.; Zyzak, J.; Żarowska, B.; Casaregola, S.; Lazar, Z. New Cytoplasmic Virus-Like Elements (VLEs) in the Yeast *Debaryomyces hansenii*. *Toxins* **2021**, *13*, 615. https://doi.org/10.3390/toxins13090615

Received: 20 July 2021
Accepted: 27 August 2021
Published: 1 September 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Yeasts can have additional genetic information in the form of cytoplasmic linear dsDNA molecules called virus-like elements (VLEs). Some of them encode killer toxins. The aim of this work was to investigate the prevalence of such elements in *D. hansenii* killer yeast deposited in culture collections as well as in strains freshly isolated from blue cheeses. Possible benefits to the host from harboring such VLEs were analyzed. VLEs occurred frequently among fresh *D. hansenii* isolates (15/60 strains), as opposed to strains obtained from culture collections (0/75 strains). Eight new different systems were identified: four composed of two elements and four of three elements. Full sequences of three new VLE systems obtained by NGS revealed extremely high conservation among the largest molecules in these systems except for one ORF, probably encoding a protein resembling immunity determinant to killer toxins of VLE origin in other yeast species. ORFs that could be potentially involved in killer activity due to similarity to genes encoding proteins with domains of chitin-binding/digesting and deoxyribonuclease NucA/NucB activity, could be distinguished in smaller molecules. However, the discovered VLEs were not involved in the biocontrol of *Yarrowia lipolytica* and *Penicillium roqueforti* present in blue cheeses.

**Keywords:** yeast; *Debaryomyces hansenii*; virus-like elements (VLEs); linear dsDNA plasmids; killer activity; osmotolerance; killer toxins

**Key Contribution:** Eight new systems of virus-like elements (VLEs) have been discovered in *Debaryomyces hansenii*; three of them have been sequenced and characterized.

## 1. Introduction

In addition to the genetic information in the form of DNA existing in the nucleus, yeast may have extra genetic elements located in the cytoplasm in the form of virus-like particles (VLPs) containing dsRNA as well as circular or linear dsDNA plasmids [1,2]. Linear dsDNA plasmids at present called virus-like elements (VLEs) due to their likely viral origin represent the predominant supplementary genetic material in yeasts. However, the number of discovered VLEs carrying strains is not very large. As a result of analysing over 1800 strains belonging to approximately 600 yeast species, molecules of this type have been found only in 40 strains [1,3–13]. Most of them belong to the Ascomycota phylum (34 strains) and the rest to the Basidiomycota phylum. In the latter, all discovered VLEs containing strains pertain only to one species, *Tausonia pullulans*, formerly

<footnote n="1">¹ Department of Biotechnology and Food Microbiology, Faculty of Biotechnology and Food Science, Wrocław University of Environmental and Life Sciences (WUELS), 50-375 Wrocław, Poland; barbara.zarowska@upwr.edu.pl (B.Z.); zbigniew.lazar@upwr.edu.pl (Z.L.)</footnote>

<footnote n="2">² SPO, INRAE, Montpellier SupAgro, Université de Montpellier, 34060 Montpellier, France; cecile.neuveglise@inrae.fr</footnote>

<footnote n="3">³ Department of Microbiology, Laboratory of Microbiome Immunobiology, Ludwik Hirszfeld Institute of Immunology and Experimental Therapy, Polish Academy of Sciences, 53-114 Wrocław, Poland; j.m.zyzak@gmail.com</footnote>

<footnote n="4">⁴ INRAE, AgroParisTech, Micalis Institute, CIRM-Levures, Université Paris-Saclay, 78350 Jouy-en-Josas, France; s.casaregola@gmail.com</footnote>

<footnote>* Correspondence: xymena.polomska@upwr.edu.pl; Tel.: +48-71-3207-791</footnote>
*Trichosporon pullulans* [1,14]. In Ascomycota, the largest group (17 strains) belongs to the Debaryomycetaceae family. According to the current taxonomy, these strains belong to the following species: *Babjeviella inositovora* (formerly *Pichia inositovora*), *Millerozyma (Pichia) acaciae*, *Schwanniomyces (Debaryomyces) etchellsii*, *Schwanniomyces (Debaryomyces) occidentalis*, *Schwanniomyces (Debaryomyces) polymorphus*, *Schwanniomyces vanrijiae* var. *yarrowii*, *Debaryomyces (Wingea) robertsiae*, and *Debaryomyces hansenii* [15–17]. Other strains are located within four families: Saccharomycetaceae, Pichiaceae, Phaffomycetaceae, Saccharomycopsidaceae, and Trichomonascaceae [1,3–5,9–12].

VLEs are usually found in pairs in yeast cells. There are also some cases of individual elements or systems consisting of three components [1,8,12]. Except for pPH1 from *Pichia heedii* and pPK1 and pPK2 from *Pichia kluyveri* discovered in mitochondria [1,12], all these elements are located in the cytoplasm [9,12,18,19]. Fragments of VLEs, so-called NUPAVs (nuclear sequences of plasmid and viral origin) can also be found in nuclear chromosomes [20]. However, the expression of nuclear genes that originated from such elements is abolished because of their high A+T content leading to the internal cleavage at poly(A) sites of mRNA and finally transcript fragmentation [21].

The largest elements in VLEs systems are considered autonomous, as they have genes encoding all enzymes necessary for their replication and transcription in cytoplasm regardless of the nuclear DNA [22–25]. Other molecules (non-autonomous) require the support of the enzymes originating from autonomous elements, e.g., RNA polymerase, helicase, capping enzyme, terminal recognition factor (TRF1), or single-strand DNA binding (SSB) proteins. However, some molecules such as pGKL1 from *Kluyveromyces lactis*, pPin1-3 from *B. inositovora*, pDHL1, and pDH1A from *D. hansenii*, may encode their own B-type DNA polymerase together with the terminal protein (TP), required for their replication. TP, along with SSB proteins and TRF1, is an essential part of the replication initiation complex of VLEs [26]. In some cases, e.g., pGKL1 or pPin1-3, it was also found that genes encoding killer toxins acting as anticodon nucleases (ACNases) and immunity determinants to these killer toxins may be located on non-autonomous elements [21,27–30].

One of the yeast species in which different VLE systems have been detected is *D. hansenii*, which is quite commonly found in low temperature or low water activity environments such as seawater, glaciers, refrigerated food, non-alcoholic beverages, wine, dairy, and meat products [31–34]. Although yeasts from this species are generally considered non-pathogenic receiving from EFSA the status of qualified presumption of safety (QPS) and the permission for use in food production [35], they have also been sporadically isolated from humans [36,37]. They occur in large numbers in ripening cheeses where they constitute the major yeast population at the end of the process, despite relatively low proteolytic and lipolytic activities [38]. This is due to their ability to assimilate lactose as well as lactic and citric acids, growing at low temperature and relatively high salt concentrations and, due to the formation of killer toxins at acidic pH [39,40].

It is known that these yeasts, depending on the strain, can secrete various killer toxin proteins. Their common feature is that they hydrolyze β-D-glucan present in the fungal cell walls. Two types of toxins have been identified acting on either (1→6) or (1→3)-β-D-glycosidic bonds. The enzyme of the first type may have different sizes, depending on the strain, ranging from 23 kDa [41] to even 232 kDa [42]. Additionally, one strain can secrete two killer proteins simultaneously [42]. The second type of toxin has been detected by Çorbacı and Uçar [43]. The molecular weight of the protein is 31.5 kDa. All mentioned *D. hansenii* toxins are encoded in the nuclear genome and although *D. hansenii* can also have additional genetic information in the form of VLEs in their cytoplasm, no connection between the killer activity of *D. hansenii* and the possession of VLEs has been found [8,9].

So far, four different plasmid systems have been discovered in representatives of *D. hansenii* species. The first of them, the three-elements system pDHL1/2/3 (8.4/9.2/15.0 kb), was isolated by Gunge et al. [8] from strain TK (CBS 7848) selected from Japanese pickles "Takuan". The second, detected in strain CBS 770, isolated from the throat of angina patient, consists of two plasmids pDH1A/B of 8.2 and 14.4 kb, respectively [9]. Fukuhara [1]
detected two additional strains having such molecules, CBS 1780 (pDH2A/B; 6.5/15 kb) and CBS 4890 (pDH3A/B; 4.7/13 kb). However, the function of all discovered VLEs in *D. hansenii* cells remains unknown.

Previously, we have been investigating the function of different yeast species including *D. hansenii* in the cheese ripening process and the possibility to utilize selected strains to shorten the cheese maturation [44–46]. We have found that *D. hansenii* strains originating from the cheese environment express a strong killer activity [39,47,48]. Furthermore, we have developed an industrial technology to produce preparations containing *D. hansenii* killer toxins for plant protection against fungal phytopathogens [49]. The technology has been implemented by Skotan S.A. (Poland) and the product named HanseFlora®is now available in Poland. In the current study, we have expanded our interests to the dissemination of VLEs in *D. hansenii* strains mainly occurring in cheese but also in other ecological niches. We tested strains deposited for years in culture collections as well as freshly isolated from the environment. By analyzing the sequence of three selected new VLE systems, we tried to find/predict a benefit for the host coming from having such elements in the cell.

**2. Results and Discussion**

**2.1. Frequency of VLEs in *D. hansenii* Yeast**

A set of 135 *D. hansenii* yeast strains was analyzed for the occurrence of VLEs; 75 of them were obtained from culture collections and 60 strains were freshly isolated from the environment (Tables 1 and 2). Culture collection strains originated mainly from dairy sources but also from ice and seawater, alcoholic and non-alcoholic drinks, meat, insects, and human skin. New isolates were gained from blue-veined cheeses, which are known to be abundant yeast microbiota [40]. To collect a large diversity of samples, eight different types of cheese produced in four European countries were selected. Isolation of yeast strains was performed on a specific medium named YORS and developed particularly for this work. This medium was based on YPD agar (Y) supplemented with oxytetracycline (O), rose bengal (R), and salt (S). Such composition allowed the selective isolation of halophilic yeast such as *D. hansenii* while inhibiting the growth of bacteria, filamentous fungi, and other yeast species present in blue-veined cheeses.

The type strain of *D. hansenii*, CLIB 197T, forms colonies of characteristic morphology on YORS agar: pink color, smooth surface, circular shape with pulvinate elevation, and entire margin (Figure S1). Randomly selected strains looking similar to the type strain were chosen for further experiments. All of them were identified as *D. hansenii* by colony PCR with a species-specific pair of primers—DhPADF and DhPADR developed by Wrent et al. [50], giving the product of 402 bp. In each case, a band of approximately 400 bp was obtained (Figure S2). The application of this method significantly reduced the time and cost of the analysis compared to amplification and sequencing of the D1/D2 domain of the large ribosomal subunit [51]. The relatively small size of the PCR product allowed us to carry out the PCR with biomass instead of purified DNA, thus, avoiding the time-consuming nucleic acid isolation required by the aforementioned method but also during the identification of *D. hansenii* by means of intergenic spacer rDNA amplification with subsequent restriction fingerprinting [52,53].

Among yeasts obtained from culture collections, no strain possessing VLEs was detected directly after reviving or after repeated cultivation in a medium containing sodium chloride, which can facilitate the proliferation of cells containing such elements [54]. The procedure of repeated cultures was used for *D. hansenii* CBS 770 and *D. hansenii* CBS 7848. It proved to be effective since no VLEs were observed in DNA preparations isolated from cultures directly after recovery from the freeze-dried preparation, while after six passages in YPD with 12% NaCl VLEs were successfully detected (Figure S3).
<table><caption>**Table 1.** *D. hansenii* strains obtained from culture collections used in the study.</caption><thead><tr><th>Culture Collection</th><th>Strain Designation</th><th>Source of Isolation</th></tr></thead><tbody><tr><td>DBFM <sup>1</sup></td><td>AI1a, AI4a, AI4b, AI4c, AI6c, AI12c, AI13b, AI14a, AI14b, AI14b-1, AI14c, BI6b, EI2b, EI3a, EI4a, EI12a, EI13c, MI1a, MI1a-1, MI1a-2, MI2a, MI4a, MI5b, MI6b, MI7a, OI1a, OI1b, OI5c, OI11c, OI13a, PI1a, PI4a, PI5b, PII3a<br>DRIP1b, DRIP2c, DRIIP1b, DRIIIP3c, DRIIW4c, DRIIW3b, DRIIW4c, DRIIIW6a, KI2a<br>DNIP4b, DNIIP2b</td><td>Dairy products<br>Rokpol cheese</td></tr><tr><td>CIRM-Levures <sup>2</sup></td><td>CLIB 236<br>CLIB 380, 594<br>CLIB 607<br>CLIB 608<br>CLIB 609, 684<br>CLIB 611<br>CLIB 613<br>CLIB 622<br>CLIB 920<br>CLIB 539, 1302, 1465, 1277, 1298, 1301<br>CLIB 543, 944, 1086, 1389<br>CLIB 545, 907<br>CLIB 542, 1295<br>CLIB 1142, 1143<br>CLIB 195, 1144<br>CLIB 1296</td><td>Dairy products and environment<br>Dorblu cheese<br>Danablu cheese<br>Roncal cheese<br>Goat's cheese<br>Camembert<br>Air<br>Milk<br>Dairy brine<br>Forage<br>Saint Nectaire<br>Cheese curd<br>Ice, glacial, and seawater<br>Beverages<br>Human<br>Meat products<br>Insects<br>Unknown<br>Salted pickles</td></tr></tbody></table>

<table><caption>**Table 2.** The origin of new *D. hansenii* isolates.</caption><thead><tr><th>Strain</th><th>Cheese Type</th><th>Country</th></tr></thead><tbody><tr><td>1 (a-d)</td><td>Fourme d'Ambert</td><td>France</td></tr><tr><td>2 (a-d)</td><td>Bleu d'Auvergne</td><td>France</td></tr><tr><td>3 (a-d)</td><td>Rokpol Lazur</td><td>Poland</td></tr><tr><td>4 (a-l)</td><td>Edelpilz</td><td>Germany</td></tr><tr><td>5 (a-f)</td><td>Rokpol KG</td><td>Poland</td></tr><tr><td>6 (a-j)</td><td>Dorblu</td><td>Germany</td></tr><tr><td>7 (a-j)</td><td>Bloose</td><td>Denmark</td></tr><tr><td>8 (a-j)</td><td>Turek niebieski</td><td>Poland</td></tr></tbody></table>

At present, it is difficult to find out whether the lack of the cytoplasmic DNA elements in the studied strains was due to their lack in the primary culture or have been caused by isolating these strains on standard yeast media or the following methods of preservation and storage. Strains of CIRM-Levures are stored in liquid nitrogen or in freeze-dried form, while cultures from DBFM are henceforth stored at -80 °C in a glycerol solution. However, for many years they have been stored under paraffin on YM agar slants at a temperature of 4 °C and have been regularly transferred once a year onto a fresh medium. Such a procedure could have led to the gradual loss of plasmids. No experiments about VLEs content have been conducted before culture conservation in both collections, but it should be implemented from now on. Consideration should, therefore, be given to the development of proper storage or revival conditions for strains deposited in culture collections, which could subsequently enhance chances for maintaining cytoplasmic DNA molecules in cells.

It appeared far more efficient to search for VLEs in fresh isolates of *D. hansenii* species obtained directly from cheeses with mould overgrowth. In this case, 25% of them contained such molecules (Table 3). There was high variability in the number of VLEs possessing strains between cheeses; VLEs were observed in half of the cheese samples tested. In most

<footnote n="1">¹ DBFM: Department of Biotechnology and Food Microbiology culture collection, Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland. ² CIRM-Levures: Centre International de Ressources Microbiennes-Levures, INRA, Jouy-en-Josas, France.</footnote>
cases, low detection rates in the rest of the trials could be correlated with the small number of analyzed strains (four isolates). In cheeses, where yeasts with VLEs were detected, their frequency ranged from 17% to as much as 50% of isolates, which suggested an important contribution of these elements to yeast cell functioning in such an ecosystem.

<table><caption>**Table 3.** Size of new VLEs systems discovered in *D. hansenii* yeast strains.</caption><thead><tr><td>System</td><td>Strains</td><td>VLE Size</td></tr></thead><tbody><tr><td>pDH4A/B/C</td><td>4a, 4c, 4d, 4e, 4f, 4i</td><td>4.8; 6.9; 15.1</td></tr><tr><td>pDH5A/B</td><td>5c</td><td>8.9; 15.1</td></tr><tr><td>pDH6A/B</td><td>7f, 7g, 7i</td><td>7.2; 15.1</td></tr><tr><td>pDH7A/B</td><td>7j</td><td>8.4; 15.1</td></tr><tr><td>pDH8A/B/C</td><td>8e</td><td>7.4; 8.0; 15.1</td></tr><tr><td>pDH9A/B</td><td>8g</td><td>7.6; 15.1</td></tr><tr><td>pDH10A/B/C</td><td>8h</td><td>7.6; 10.0; 15.1</td></tr><tr><td>pDH11A/B/C</td><td>8i</td><td>9.2; 10.0; 15.1</td></tr></tbody></table>

We also observed that in a given cheese both may occur, strains having the same VLE system, as in the case of pDH4A/B/C in the German Edelpilz cheese, as well as coexisting strains with different sets of such elements, as in strains isolated from the Polish Turek niebieski cheese (Table 3, Figure 1). In the latter, four VLEs possessing strains were detected, i.e., 8e, 8g, 8h, and 8i, and each of them had a different set of such molecules. Based on current data, it cannot be distinguished whether they arise from different ancient rearrangements of an ancestral VLEs system or whether there are undergoing current modifications.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** New VLEs systems discovered in *D. hansenii* strains freshly isolated from mould ripened cheeses. Names of yeast strains possessing particular systems are given in Table 3. * Genetic material isolated from strain without VLEs. M. Gene Ruler 1kb Plus DNA Ladder (Thermo Fisher Scientific, Warszawa, Poland).</figcaption></figure>

In total, eight different VLEs systems have been distinguished: four binary and four ternary (Figure 1, Table 3). The biggest molecule in each system had a consistent size of 15.1 kb for all tested strains, which is longer than the longest previously known element pSKL from *Lachancea kluyveri* of 14.2 kb [55]. The size variation was observed for the smallest and medium-sized elements. The length of the smallest one ranged from 4.8 to 9.2 kb and the medium one from 7.5 to 10.0 kb.

## 2.2. *The Characteristics of Selected VLEs Systems*

Three newly discovered sets of VLEs (pDH4A/B/C, pDH5A/B, and pDH6A/B) were selected for next-generation sequencing. They were isolated from strains acquired from cheeses produced in Germany, Poland, and Denmark (strains 4e, 5c, and 7g). The libraries were made from DNA fractions containing a mixture of linear DNA plasmids separated
from mitochondrial and nuclear DNA by ultracentrifugation in CsCl-bisbenzimide gradient (Figure S4). This approach was not previously used in the preparation of yeast VLEs for sequencing. So far, labour-intensive purification of individual molecules and their subsequent Sanger-based sequencing has been conducted [24,26,55–58].

Overlapping reads were assembled into contigs which were compared with GenBank nucleotide database using blastn [59]. It allowed identifying VLE sequences and contaminations with mitochondrial DNA. Such contamination was probably created at the fraction collection stage after separation in the CsCl-bisbenzimide gradient. Due to the small difference in buoyant density between mtDNA and VLEs, both fractions are faintly divided from one another and therefore easy to mix with each other. Nuclear DNA fragments have also been identified among the obtained sequences. Their presence indicates that despite careful operation chromosome fragmentation has not been avoided. Apart from the random fragments of different chromosomes, the sequences of nuclear origin belonged in many cases to the rRNA encoding genes, which are located in many copies in three loci in the genome of *D. hansenii* CBS 767T [60]. In addition, fragments of LTR-retrotransposon Tdh5 have also been identified [61]. Frequent detection of the latter may be associated with the presence of more than 20 copies of this element in the *D. hansenii* genome and subsequent generation of many extrachromosomal DNA copies of this element. While the length of Tdh5 retrotransposons is approximately 5.3 kb, their isolation together with the plasmid fraction was also possible.

Contigs assigned to VLEs constituted almost complete plasmid sequences. In the case of four of them (pDH4C, pDH5B, pDH6A, pDH6B), only the fragment of one terminal inverted repeat (TIR) sequence was missing, while for the other plasmids, sequences at both ends had to be filled by primer walking. Open reading frames were annotated using CLC Sequence Viewer v.8.0.0 (Qiagen Digital Insights). This software was selected based on previously analysed ORF detection accuracy on already known VLEs sequences: pGKL2, pPE1B, pPac1-1, and pSKL. The minimal ORF coding capacity was set to 70 aa as genes encoding 70 aa polypeptides were discovered on pGKL2 and pSKL and proved to be functional [62]. The annotated sequences of pDH4A/B/C, pDH5A/B, and pDH6A/B were deposited in GenBank under the following accession numbers: MF795091/MF795092/MF795093, KX904874/KX904875, and KX858805/KX904876, respectively.

All new VLEs turned out to be extremely A+T rich (73–76%) and compact in structure, exactly as for previously known plasmids. Their coding capacity was also very high; 75–83% and 88–90% for non-autonomous and autonomous molecules, respectively. At both ends of all analysed VLEs, terminal inverted repeats (TIR) strictly specific to the element type were found (Figures 2 and 3). Those with a length of 715–716 bp detected on autonomous elements turned out to be about 230 bp longer than the longest TIRs previously reported on pSKL. In contrast, TIRs from the smallest VLEs in the systems also proved to be the smallest known (91 bp) and were identical to the 91-first bases of pDH1A TIRs, which are 96 bp in length [58]. In the case of the middle non-autonomous elements, inverted repeats located at the ends were composed of TIR sequences from the small molecules extended for 33 bp.

Each of the newly sequenced VLEs had its own DNA polymerase-encoding gene that enables its replication in the cytoplasm independently of the nucleus using a previously disclosed protein-primed mechanism [reviewed by 30]. The sequence of this gene was closely related to the type of element, highly conserved for autonomous (99–100% identity) and non-autonomous (95–98% id.) molecules, while much more diverse between genes from these two types of VLEs (~ 29% id.).

The confirmation of distinctiveness of these two kinds of DNApol can also be found in the phylogenetic tree illustrating the relationships between various groups of proteins that are representatives of the type B superfamily of DNA polymerases (Figure 4). Similar phylogenetic analyses have been performed before but on a smaller number of sequences [57,58,63,64]. This time, the analysis was expanded due to an increase in the number of available data. Nevertheless, proteins encoded by autonomous and non-autonomous
plasmids formed separated groups within the clade of VLEs-origin DNApol s. It can be hypothesized that both groups derived from a common ancestor in which the DNApol gene was probably duplicated, then single copies were divided into two linear molecules of different functions and further differentiated. Moreover, our results indicated that VLEs-origin DNA polymerases are closely related to DNApol of adenoviruses and bacteriophages of both Gram-positive and -negative bacteria, which suggests a common origin. The structure of the phylogenetic tree obtained in this work showed that mitochondrial DNApol s differentiated earlier from the common progenitor than those of yeast cytoplasmic VLEs, Adenoviridae, and phages. Furthermore, the evolution of plasmid DNApol of mitochondria apparently followed more complicated paths. The basis for this presumption is that the DNApol of mitochondrial linear plasmids of fungi (Ascomycota and Basidiomycota) are arranged in three different clades distant from cytoplasmic and viral B-type polymerases. Therefore, our results confirm the conclusions already reached by Klassen and Meinhardt [64] from the analysis of a smaller pool of sequences.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** The schematic layout of open reading frames (ORFs) location on newly sequenced autonomous VLEs of *D. hansenii* (pDH4C, pDH5B, pDH6B) and previously known from other yeast species. The characteristics of individual ORFs is described in Tables 4 and 5. The arrows of light blue color correspond to ORFs encoding enzymes involved in replication, dark blue in transcription, and yellow of unknown function. The products of ORFs colored in purple show some similarities with immunity determinants to VLE killer toxins. The direction of arrows indicates the transcription direction; black triangles represent terminal inverted repeats (TIRs) at the ends of plasmids.</figcaption></figure>

<table><caption>**Table 4.** Characteristics of open reading frames (ORFs) detected on new autonomous VLEs from *D. hansenii* yeast (GenBank Acc. No.: pDH4C-MF795093, pDH5B-KX904875, pDH6B-KX904876).</caption><thead><tr><th>ORF</th><th>DNA Strand</th><th>ORF Length [bp]</th><th>ORF Identity [%]</th><th>Stop Codon</th><th>UCS</th><th>UCS Position <sup>1</sup> [bp]</th></tr></thead><tbody><tr><td>1</td><td>—</td><td>2997</td><td>100</td><td>TGA</td><td>ATATGA</td><td>-34</td></tr><tr><td>2</td><td>—</td><td>1698</td><td>100</td><td>TAA</td><td>ATATGA</td><td>-52</td></tr><tr><td>3</td><td>+</td><td>1761</td><td>99 <sup>2</sup></td><td>TAA</td><td>ATCTGA</td><td>-25</td></tr><tr><td>4</td><td>+</td><td>486</td><td>100</td><td>TGA</td><td>ATGTGA</td><td>-108</td></tr><tr><td>5</td><td>+</td><td>2964</td><td>99 <sup>3</sup></td><td>TAA</td><td>ATGTGA</td><td>-50</td></tr><tr><td>6</td><td>—</td><td>393</td><td>100</td><td>TAA</td><td>ATGTGA</td><td>-26</td></tr><tr><td>7</td><td>+</td><td>1392</td><td>100</td><td>TAA</td><td>ATTTGA</td><td>-43</td></tr><tr><td>8</td><td>+</td><td>318</td><td>100</td><td>TAA</td><td>ATTTGA</td><td>-8</td></tr><tr><td>9</td><td>+</td><td>861/786/1104</td><td>100 <sup>4</sup></td><td>TAA</td><td>ATATGA</td><td>-38</td></tr><tr><td>10</td><td>—</td><td>234</td><td>100</td><td>TAA</td><td>ATTTGA</td><td>-34</td></tr><tr><td>11</td><td>+</td><td>186</td><td>100</td><td>TAA</td><td>ATTTGA</td><td>-95</td></tr></tbody></table>

<footnote n="1">¹ The distance of UCS (upstream conserved sequence) first base from the start codon, ² 1760/1761 identities, ³ 2963/2964, ⁴ 100% identities on the first 785 bp of ORF9.</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** The scheme of open reading frames (ORFs) organization on newly sequenced non-autonomous VLEs of *D. hansenii* (pDH4A/B, pDH5A, pDH6A) and previously known from other yeast species. The characteristics of individual ORFs is described in Tables 6 and 7. The direction of arrows indicates the transcription direction; black triangles represent terminal inverted repeats (TIRs) at the ends of plasmids. Homologous ORFs or functional domains are marked with one color. DNA fragments encoding full domains are ended with vertical lines from both sides, while partial or with some similarities with arrows.</figcaption></figure>

It seems that genes of DNA polymerases have been pulled into VLEs that were sourced from a different origin. The research of Sýkora et al. [65] and Vopálenský et al. [66], who analyzed the transcriptional mechanism of pGKL from *K. lactis*, showed that RNA polymerase encoded on the autonomous molecules as well as the promoter sequences of all genes lying on both autonomous and non-autonomous molecules probably derived from the protoplast of poxviruses. Furthermore, the DNA polymerases found in viruses of the Poxviridae family do not show homology to the polymerases of VLEs origin.

Sequence analysis of pDH4C, pDH5B, and pDH6B showed their extraordinary conservation, with 99% identity. Therefore, in those molecules, all open reading frames were located precisely in the same orientation and order, and 10/11 of them had the same length, start, and stop codons (Figure 2, Table 4). Individual plasmids differed only in the length of the ORF9 located on the right end. In pDH6B, this ORF was the longest (1104 bp), while in pDH4C and pDH5B, two and one nucleotide insertions, respectively, resulted in the shift of the reading frame and the early appearance of a stop codon, which shortened the ORF length to 861 and 786 bp. A similar genetic organization of core 10 ORFs has already been reported for pGKL2, pPE1B, pPac1-1, and pSKL [22,24,25,55]. However, no homology to DhORF9 was found in these plasmids, although, the additional ORF (ORF1) was found on the left side of pSKL and pGKL2.

The predicted products of all extra ORFs located on VLEs of different species show some similarities to proteins providing immunity to plasmid-encoded killer toxins (Table 5). The putative DhORF9p appears to be most closely related to the immunity determinant of *W. robertsiae* encoded by pWR1A ORF5 (35% identical amino acids on 90% of protein length) and of *M. acaciae* pPac1-2 ORF4 (36% id. aa of 83% of p.l.), which are also closely related to each other (49% id. aa of 97% p.l.). No significant similarity of DhORF9p was found to the zymocin immunity protein encoded by pGKL1 ORF4. The deduced immunity-like proteins of the pSKL and pGKL2 probably belong to a different group of proteins resembling immunity as they have 48% common amino acids on the length of 80% of the protein sequence and they are homologous to the mentioned zymocin resistance protein (25% id./52% p.l.).
<figure><img src="image_8.png" /><figcaption>**Figure 4.** The phylogenetic relationship between B-type DNA polymerases of yeast cytoplasmic VLEs, viruses, and mitochondrial plasmids of different organisms. The phylogram was constructed with the neighbor-joining method using aligned full amino acid sequences of 51 proteins obtained from GenBank and from this study. GenBank accession numbers are given in brackets. The tree topology was verified with a bootstrap method with 1000 replications and the bootstrap percentage was presented on each node. The scale bar shows the number of changes per position.</figcaption></figure>
<table><caption>**Table 5.** Functional characteristics of autonomous VLEs from *D. hansenii* yeast (GenBank Acc. No.: pDH4C-MF795093, pDH5B-KX904875, pDH6B-KX904876).</caption><thead><tr><th>ORF</th><th>Protein Length<br>[aa]</th><th>Predicted Function/Similarity <sup>1</sup></th></tr></thead><tbody><tr><td>1</td><td>998</td><td>B-type DNA polymerase</td></tr><tr><td>2</td><td>565</td><td>Similarity to putative mRNA capping-enzymes from <i>M. acaciae</i> pPac1-1 and <i>S. etchellsii</i> pPE1B</td></tr><tr><td>3</td><td>586</td><td>Helicase-DEXdc and HELICc domain containing protein</td></tr><tr><td>4</td><td>161</td><td>Similarity to putative ssDNA binding protein from <i>M. acaciae</i> pPac1-1 and <i>S. etchellsii</i> pPE1B</td></tr><tr><td>5</td><td>987</td><td>RNA polymerase larger subunit</td></tr><tr><td>6</td><td>130</td><td>Similarity to putative RNA-polymerase subunits from <i>S. etchellsii</i> pPE1B (ORF7) and <i>M. acaciae</i> pPac1-1 (ORF7)</td></tr><tr><td>7</td><td>463</td><td>Unknown</td></tr><tr><td>8</td><td>105</td><td>DNA-binding protein TRF1 (terminal region recognition factor 1)</td></tr><tr><td>9</td><td>286/261/367</td><td>Similarity to immunity determinants to killer toxins encoded by <i>M. acaciae</i> pPac1-2 and <i>D. robertsiae</i> pWR1A, fragment of <i>M. farinosa</i> chromosome F (XP_004202381)</td></tr><tr><td>10</td><td>77</td><td>Unknown, high similarity to product of ORF8 from: <i>L. kluyveri</i> pSKL, <i>M. acaciae</i> pPac1-1, <i>K. lactis</i> pGKL2, <i>S. etchellsii</i> pPE1B</td></tr><tr><td>11</td><td>62</td><td>Unknown, high similarity to products of: <i>S. etchellsii</i> pPE1B ORF3 <i>M. acaciae</i> pPac1-1 ORF11</td></tr></tbody></table>

The function of the immunity-like genes is as yet unknown. However, Schaffrath et al. [70] demonstrated using pGKL2ORF1p as an example that this protein does not guarantee resistance to zymocin. In addition, they also showed that this protein does not participate in the replication processes and maintenance of the VLE system.

*DehaORF9* appears to be intact with a putative upstream control region (UCR) acting as yeast plasmid gene's promoter [71-73]. This sequence contains the ATA TGA key fragment located 38 bp upstream from the start codon. Such sequences (consensus: ATNTGA), called upstream conserved sequences (UCS), usually precede the start codon about 100 bp for autonomous and around 30 bp for non-autonomous plasmids [25,30,70,74,75]. UCS in the promoter sequence is necessary for the correct binding of RNA polymerase of VLEs origin and handling of the transcription of plasmid genes in the cytoplasm.

The sequences of new non-autonomous VLEs were searched for genes, which could be involved in killer toxin production. In all new systems, ORFs containing either a full domain of chitinase activity (GH18) or their fragments were detected (Figure 3, Tables 6 and 7). Polypeptides with such function are an inseparable part of VLE-encoded killer toxins [27,29]. However, they are also present on cryptic elements which function is still unknown [9,58]. The most interesting ORF containing GH18 is located on pDH5A. The predicted protein encoded by ORF2 contains the full domain of chitinase along with the domain responsible for chitin binding (ChBP1). There is also an incomplete autolysin domain of 1,4-beta-N-acetylmuramyl-hydrolase activity. The autolysin contains elements of the LysM domain that is responsible for binding to peptidoglycan in bacterial cell walls as well as fungal chitin which is structurally similar [76]. LysM domains were also detected in proteins encoded on plasmids pDH1A, pGKL1, pPin1-3, pPac1-2, pWR1A, and pPE1A (Figure 3).

Additionally, in pDH5AORF2p, a fragment of the deoxyribonuclease domain (DNase NucA/NucB) was also found. Parts of DNase NucA/NucB were also present in other pDH systems, i.e., in ORF2 of pDH4B and in ORF3 of pDH6A. The prediction of protein localization with DeepLoc-1,0 [77] revealed that pDH5AORF2p and pDH6AORF3p can be extracellular proteins. However, only in pDH5AORF2p was the signal peptide of 18 aa detected confirming the possible secretion of this protein outside the cell. Interestingly, in each of the pDH systems, a protein that probably functions in the nucleoplasm (pDH4AORF3p, pDH5AORF4p, and pDHORF7p) was detected. These proteins are almost identical and highly homologous to the hypothetically cytoplasmic pDH4AORF2p, pDH5AORF3p, and

<footnote n="1">¹ The function of open reading frames was at first predicted based on blastx [67], then in case of no hits Conserved Domains Database [68], as well as pfam database [69], were searched against protein queries for conserved domains.</footnote>
pDH6AORF6p proteins. It could be assumed that these proteins may be an element of the host's genetic material defense against the killer activity of DNases.

<table><caption>**Table 6.** The characteristics of open reading frames (ORFs) detected on new non-autonomous VLEs of *D. hansenii* yeast: pDH4A (GenBank Acc. No. MF795091), pDH4B (MF795092), pDH5A (KX904874), and pDH6A (KX858805).</caption><thead><tr><th>VLE</th><th>ORF</th><th>DNA Strand</th><th>ORF Length [bp]</th><th>Stop Codon</th><th>UCS</th><th>UCS Position [bp] <sup>1</sup></th></tr></thead><tbody><tr><td rowspan="3">pDH4A</td><td>1</td><td>+</td><td>3003</td><td>TGA</td><td>ATG TGA</td><td>-25</td></tr><tr><td>2</td><td>+</td><td>363</td><td>TAA</td><td>ATA TGA</td><td>-27</td></tr><tr><td>3</td><td>+</td><td>222</td><td>TGA</td><td>ATG TGA</td><td>-26</td></tr><tr><td rowspan="3">pDH4B</td><td>1</td><td>+</td><td>3003</td><td>TGA</td><td>ATG TGA</td><td>-25</td></tr><tr><td>2</td><td>+</td><td>1647</td><td>TAA</td><td>ATA TGA</td><td>-254</td></tr><tr><td>3</td><td>-</td><td>585</td><td>TAG</td><td>ATG TGA</td><td>-28</td></tr><tr><td rowspan="4">pDH5A</td><td>1</td><td>+</td><td>3000</td><td>TGA</td><td>ATG TGA</td><td>-25</td></tr><tr><td>2</td><td>+</td><td>3681</td><td>TAA</td><td>ATA TGA</td><td>-27</td></tr><tr><td>3</td><td>+</td><td>363</td><td>TAA</td><td>ATG TGA</td><td>-26</td></tr><tr><td>4</td><td>+</td><td>369</td><td>TGA</td><td>ATG TGA</td><td>-26</td></tr><tr><td rowspan="7">pDH6A</td><td>1</td><td>+</td><td>3003</td><td>TGA</td><td>ATG TGA</td><td>-27</td></tr><tr><td>2</td><td>+</td><td>396</td><td>TAA</td><td>ATT TGA</td><td>-34</td></tr><tr><td>3</td><td>+</td><td>1092</td><td>TGA</td><td>TTA TGA</td><td>-66</td></tr><tr><td>4</td><td>+</td><td>360</td><td>TGA</td><td>ATG TGA</td><td>-24</td></tr><tr><td>5</td><td>+</td><td>264</td><td>TAA</td><td>ATA TGA</td><td>-40</td></tr><tr><td>6</td><td>+</td><td>363</td><td>TAA</td><td>ATG TGA</td><td>-26</td></tr><tr><td>7</td><td>+</td><td>369</td><td>TAA</td><td>ATG TGA</td><td>-26</td></tr></tbody></table>

<table><caption>**Table 7.** The analysis of predicted proteins based on open reading frame (ORF) sequences detected on new non-autonomous VLEs of *D. hansenii* yeast: pDH4A (GenBank Acc. No. MF795091), pDH4B (MF795092), pDH5A (KX904874), and pDH6A (KX858805).</caption><thead><tr><th>VLE</th><th>ORF</th><th>Protein Length [aa]</th><th>Subcellular Localization Prediction <sup>1</sup></th><th>Predicted Function or Functional Domains/Similarity <sup>2</sup></th></tr></thead><tbody><tr><td rowspan="3">pDH4A</td><td>1</td><td>1000</td><td>Cytoplasm</td><td>B-type DNA polymerase</td></tr><tr><td>2</td><td>120</td><td>Cytoplasm</td><td>Unknown, identical to pDH5A ORF3p and pDH6A ORF6p, high homology to pDH4A ORF3p, pDH5A ORF4p, pDH6A ORF7p and B. inositovora XP_018982205</td></tr><tr><td>3</td><td>73</td><td>Nucleoplasm</td><td>Unknown, identical to the part of pDH5A ORF4p and pDH6A ORF7p</td></tr><tr><td rowspan="3">pDH4B</td><td>1</td><td>1000</td><td>Cytoplasm</td><td>B-type DNA polymerase</td></tr><tr><td>2</td><td>548</td><td>Peroxisome membrane</td><td>Similarity to GH18 chitinase-like superfamily proteins, fragment of DNase NucA/NucB</td></tr><tr><td>3</td><td>194</td><td>Cytoplasm</td><td>Unknown, some similarities to ORF1p from autonomous plasmids pSKL and pGKL1</td></tr><tr><td rowspan="4">pDH5A</td><td>1</td><td>999</td><td>Cytoplasm</td><td>B-type DNA polymerase</td></tr><tr><td>2</td><td>1226 <sup>3</sup></td><td>Extracellular</td><td>Fragment of autolysin domain, full chitin binding 1 domain (ChtBD1), full GH18 chitinase-like domain, fragment of deoxyribonuclease NucA/NucB domain</td></tr><tr><td>3</td><td>120</td><td>Cytoplasm</td><td>as described for pDH4A ORF2p</td></tr><tr><td>4</td><td>122</td><td>Nucleoplasm</td><td>Unknown, identical to pDH6A ORF7p, high homology to pDH4A ORF2p, pDH5A ORF3p and pDH6A ORF6p</td></tr><tr><td rowspan="7">pDH6A</td><td>1</td><td>1000</td><td>Cytoplasm</td><td>B-type DNA polymerase</td></tr><tr><td>2</td><td>131</td><td>Mitochondrion matrix</td><td>Contains fragment of GH18 chitinase domain</td></tr><tr><td>3</td><td>363</td><td>Extracellular</td><td>Contains fragment of DNase NucA/NucB</td></tr><tr><td>4</td><td>119</td><td>Cytoplasm</td><td>Unknown, homology to proteins of unknown function of D. hansenii</td></tr><tr><td>5</td><td>87</td><td>Extracellular</td><td>CBS 767<sup>T</sup>, pDH1A ORF4p and to the part of pPE1A ORF3p</td></tr><tr><td>6</td><td>120</td><td>Cytoplasm</td><td>Unknown, homology to the part of pDH1A ORF3p</td></tr><tr><td>7</td><td>122</td><td>Nucleoplasm</td><td>As described for pDH4A ORF2p</td></tr></tbody></table>

<footnote n="1">¹ described below Table 4.</footnote>

<footnote n="1">¹ Subcellular localization of hypothetical proteins was predicted with DeepLoc-1 [77]. ² The function of open reading frames was predicted as described under Table 3. ³ Protein with detected signal peptide.</footnote>
However, these assumptions need to be confirmed in laboratory experiments. DNase activity has not yet been identified in any protein derived from linear plasmids, including killer ones. However, PaT and zymocin killer toxins indirectly affect DNA integrity in *M. acaciae* and *K. lactis*. These toxins cleave selected tRNAs, which leads to damages in the DNA repair systems and accumulation of untidy mutations [29,78–80].

The pDH5AORF2p protein may undergo modifications with KEX2 protease because five motifs recognizable by this enzyme can be found in its sequence: one in the autolysin domain, two in the chitinolytic domain, and two before the fragment of the putative DNase. The post-translational modifications by KEX1/KEX2 proteases are necessary for the secretion of active killer toxins encoded both chromosomally and by VLEs (for a review, see [81]). Further, ORF2p secondary structure analysis with Phobius [82] showed that probably the fragment involved in chitin-binding and its cleavage is localized on the outer side of the cell membrane of the target yeast or mould cell, then, there is one transmembrane region followed by the element of the DNA digesting activity located on the cytoplasmic side.

In the case of the remaining proteins potentially associated with killer activity (pDH4BORF2p, pDH6AORF2-5p), none of them have a signal peptide, the presence of which would suggest the possibility of extracellular secretion. Only pDH4BORF2p has a 31-amino acid trans-membrane region located between the fragments identified by the Phobius software [82] as cytoplasmic and non-cytoplasmic. The transmembrane region (beta subunit) is present, e.g., in zymocin, and is probably responsible for the transport of the active gamma subunit to the cytoplasm of the attacked cell [30].

Some ORFs on the non-autonomous VLEs are not conserved in the other *D. hansenii* systems. An example is ORF3 from pDH4B, which has not been detected in other *D. hansenii* VLEs. However, ORF3p showed homology to the immunity-like protein of ORF1 from the autonomous plasmids pSKL and pGKL2. Their function is still unknown but blastp search showed some similarities with alpha-mannosyltransferases, which may suggest that they could be involved in protein glycosylation processes.

## 2.3. *The Relationship between New VLEs and Killer Activity*

Considering the above data, we decided to check the relationship between new VLEs and killer activity in practice. For that purpose, selected yeasts possessing VLEs were cultured for six weeks in media without selection pressure (without NaCl) and at an elevated temperature (30 °C). Similar conditions caused the loss of the pDHL1/2/3 system in cells of *D. hansenii* CBS 7848 [8]. The lack of VLEs was first confirmed by electrophoretic separation of total DNA isolated from yeast cells, and then, by polymerase chain reaction (Figure S5). Gunge et al. [8] have shown that cells can gradually lose their VLEs from a system, eventually losing their autonomous molecules. Therefore, primers enabling the amplification of a 963 bp DNA polymerase fragment from the autonomous plasmids were designed for PCR reactions. As it is shown in Figure S5, the applied procedure of treatment enabled the complete removal of VLEs from *D. hansenii* cells.

After confirmation of the VLEs' lack, a cross test was performed with *Y. lipolytica* PII6a yeast and *P. roqueforti* PR1 mould used as sensitive organisms (Figures S6–S9). Both species coexist with *D. hansenii* in blue-veined cheeses and could be the target of their VLE-encoded toxins due to the presence of chitin in their cell wall [83,84]. No killer activity was observed at neutral pH, which is optimal for VLE-origin killer toxins, zymocin, and PaT [30,85]. The test was performed at 14 °C, which is used for blue-veined cheese ripening as well as at 28 °C which is the standard temperature for yeast cultivation. In addition, the loss of immunity of *D. hansenii* strains to its own toxin due to the loss of linear plasmids was also not observed, neither at pH 4.6 nor 7 (Figure S10). The killer effect of the studied strains on *Y. lipolytica* and *P. roqueforti* was observed at pH 4.6, which is optimal for killer toxin activity of *D. hansenii* strains isolated from cheese [39,42,86]. However, no difference in lethal activity between VLEs possessing and cured of them isolates could be noted (Figures S11 and S12). Thus, the antifungal activity of the yeast *D. hansenii* is not encoded
within VLEs. It is in accordance with results obtained previously by Gunge et al. [8] for the
pDHL1/2/3 system.

2.4. *The Relationship between Plasmids and Halotolerance*

Bearing in mind that the function of linear cytoplasmic genomic elements of *D. hansenii*
remains cryptic, it was also examined if their presence increases the yeast resistance to a
high concentration of sodium chloride. Although sustaining VLEs in *D. hansenii* cells may
be dependent on the presence of salt in the environment [54], there was no change in growth
profiles of yeast deprived of such elements in comparison with strains possessing them
in YPD medium containing NaCl in the range of 0 to 24% (for an example see Figure S13).
Another possibility for *D. hansenii* to maintain VLEs is its putative toxicity against lactic
acid bacteria present in cheese. This hypothesis remains to be tested.

**3. Conclusions**

Although it is known that VLEs are the most widespread type of additional genetic
information in yeast, we still do not know everything about them. This is due, among
other things, to the fact that until now sequences of only three complete VLE systems
and four individual linear dsDNA molecules have been deposited and analyzed. In this
work, the available data were expanded to include sequences of the full three new systems
from one yeast species—*D. hansenii*, including, for the first time, a system consisting of
three elements. These plasmids, apart from the genes associated with their own replication
and transcription, have ORFs resembling genes of killer activity and resistance. However,
our experimental studies showed no relationship between killer activity and the presence
of VLEs in host cells and, therefore, the function of VLE-encoded proteins is still not
deciphered. Further analyses are needed to check the conditions under which VLE-derived
genes are expressed and correlate them with organisms that may be targeted for potential
VLE-encoded toxins. It would be also interesting to isolate *Debaryomyces* strains harboring
VLEs from different environments than dairy products. The linear plasmid purification
procedure followed by NGS used in this work can be utilized for rapid sequencing of similar
systems in other organisms and can contribute to a better understanding of their functions.

**4. Materials and Methods**

4.1. *Microorganisms*

The presence of VLEs was studied in the pool of culture collection *D. hansenii* strains
originated from Department of Biotechnology and Food Microbiology (DBFM) deposit,
WUELS, Poland and from CIRM-Levures, INRA, France (http://www6.inra.fr/cirm/
Levures) as well as of new isolates from different blue-veined cheeses produced in several
European countries (Tables 1 and 2). The following strains were used as references for
their identification: *D. hansenii* CLIB 197T (=CBS 767T), *D. hansenii* CLIB 907 (formerly
*C. famata* var. *famata* CBS 1795T), *D. hansenii* CBS 7848, *D. hansenii* CBS 770, *D. fabryi* CLIB
422T (=CBS 789T), *D. subglobosus* CLIB 908T (=CBS 1796T). *Yarrowia lipolytica* PII6a obtained
from DBFM, and *Penicillium roqueforti* PR1 (CHR Hansen) served as sensitive cultures for
antifungal activity tests. All yeasts used in this study were stored at 4 °C on YPD agar
slants (yeast extract 1%, peptone 2%, glucose 2%, agar 1.5%) supplemented with 3% NaCl,
except *Y. lipolytica* and *P. roqueforti*, which were kept on YPD agar without salt under the
same conditions.

4.2. *New Yeast Strain Isolation*

Yeast isolation was performed on YORS agar (yeast extract 10 g L⁻¹, peptone 20 g L⁻¹, glucose 20 g L⁻¹, oxytetracycline 0.1 g L⁻¹, rose bengal 0.033 g L⁻¹, NaCl 120 g L⁻¹, agar 15 g L⁻¹). The medium without oxytetracycline was sterilized at 121 °C for 20 min. After cooling to approximately 90 °C, sterilely weighed oxytetracycline was added to the medium and quickly mixed until dissolved. Cheese samples of 5 g were homogenized in 45 mL of physiological salt solution (Stomacher Lab Blender, Seward, Worthing, UK) and
diluted in the same liquid. Decimal dilutions were spread onto the surface of YORS agar. Colonies were taken from plates after 5 days incubation at 25 °C and streaked onto YPD agar plates with 3% NaCl.

### 4.3. *Yeast Isolate Identification*

As we wanted to analyze only *D. hansenii* strains, new yeast isolates were initially examined by colony PCR with *D. hansenii* species-specific pair of primers DhPadF/DhPadR developed by Wrent et al. [50] under their conditions except for the initial denaturation step, which was prolonged from 5 to 10 min. The reaction was performed with a small portion of 24 hour-old biomass taken from a plate with a 10 µL tip and mixed with 10 µL of 2 × Taq Nova Red Master Mix (Blirt-DNA-Gdansk, Gdansk, Poland), 0.4 µL of 10 µM each primer, and 9.2 µL of distilled water. The product was visualized with ethidium bromide after electrophoresis in 1.2% agarose gel at 120 V for 40 min.

### 4.4. *Plasmid Detection*

The VLEs presence in cells of new yeast strains was studied immediately after isolation, while for culture collection strains, it was carried out twice: first directly after reviving and then after four propagations in YPD medium with gradually increased NaCl concentration in the range of 0–9%. Each culture was shaken at 25 °C at 160 rpm (G10 Gyrotory Shaker, New Brunswick Scientific Co., Edison, NJ, USA) and after 7 days 1 mL was transferred to a new medium. The last culture was made in YPD containing 12% of NaCl and took 24 h. Then, the biomass was centrifuged at 6000 rpm (Sigma 16–32K), washed with distilled water, and used for DNA isolation.

The total DNA was obtained with GeneMATRIX Bacterial and Yeast Genomic DNA Purification Kit (Eur<sub>x</sub>, Gdansk, Poland) according to the manufacturer's protocol. Plasmids were separated from nuclear and mitochondrial DNA in 0.8% agarose gel with 1×TAE buffer at 90 V for 1.5 h and visualized with ethidium bromide.

VLEs for NGS sequencing were isolated according to an adapted protocol of Querol et al. [87]. Briefly, cells were harvested, lysed with zymolyase 100T (ICN Biochemicals, Aurora Ohio USA), and treated with SDS. Total DNA was precipitated within ethanol and resuspended in TE buffer. CsCl (1 g mL<sup>-1</sup>) and 10 µL mL<sup>-1</sup> of 10 mg mL<sup>-1</sup> stock solution of bisbenzimide (Hoechst; Thermo Fisher Scientific Poland, Warszawa, Poland) were added to total DNA. Then, the plasmid fraction along with mitochondrial DNA was separated by ultracentrifugation (Beckman L7-55 Ultracentrifuge, Beckman Coulter, Indianapolis, IND, USA) after 60 h at 50,000 rpm at 20 °C. The fraction with plasmid DNA was purified four times with isoamyl alcohol buffer made of 100 mL of isoamyl alcohol, 1 mL Tris-HCl 1 M pH 8, 20 g CsCl, and 20 mL water. DNA was precipitated with Na-acetate and isopropanol and finally resuspended in 50 µL of TE buffer.

### 4.5. *Sequencing, Assembly and Annotation*

VLE containing fractions were used for the construction of libraries, which were sequenced with Illumina MiSeq at Genome S.A. (Warszawa, Poland). The plasmid libraries were run 250 cycles giving paired-end reads of 2 × 250 bp. Raw Illumina reads were trimmed and assembled using CLC Genomic Workbench v. 8.0.0 (QIAGEN Digital Insights). The obtained contigs were blasted against the *D. hansenii* nuclear and mitochondrial genome to exclude contigs not belonging to VLEs [59]. Missing fragments were filled with the primer walking method. Obtained full sequences of linear plasmids were deposited in GenBank under the following accession numbers: pDH4A/B/C- MF795091/MF795092/MF795093, pDH5A/B-KX904874/KX904875, and pDH6A/B-KX858805/ KX904876.

Gene annotation on VLEs was performed with CLC Sequence Viewer (QIAGEN Digital Insights). The function of open reading frames was at first predicted by sequence similarity on the basis of blastx [67]. Then, in case of no hits, the Conserved Domains Database [68], as well as the Pfam database [69], were used to search for conserved domains.
Protein secondary structure and subcellular localization were analyzed with Phobius [82] and DeepLoc-1.0 [77].

### 4.6. Phylogenetic Analysis

Amino acid sequences of DNA polymerases were aligned in CLC sequence viewer v.8.0.0 (QIAGEN Digital Insights). The same software was used for the construction of the phylogenetic tree based on the neighbor-joining method with Jukes-Cantor protein distance measure after 1000 bootstrap replications.

### 4.7. VLE Curing

To force the yeast cells to lose linear cytoplasmic DNA, strains were cultivated on YPD medium without NaCl at 30 °C, with shaking at 160 rpm for 7 days and transferring 1 mL into a fresh medium. Such procedure was repeated 8 times. Then, the total DNA was isolated and inspected on agarose electrophoresis as described previously. Furthermore, the isolated DNA was used as a template for PCR reaction with primers designed with Primer3Plus [88] specific for the DNA polymerase gene from autonomous *D. hansenii* VLEs: pCpolDNAF1 (5' TTG GTG CGA TAT AGA TGG AAA A 3') and pCpolDNAR1 (5' TGC TTT TTA AAC CAG CTC CA 3'). The reaction mixture was: 0.4 µL DNA mixed with 10 µL of 2 × Taq Nova Red Master Mix (Blirt-DNA-Gdansk), 0.4 µL of 10 µM each primer, and 8.8 µL of distilled water. The PCR reaction was performed in a Personal Thermocycler (Biometra, Göttingen, Germany) under the following conditions: initial denaturation at 95 °C for 5 min; 35 cycles of denaturation at 94 °C for 1 min, primers annealing at 49 °C for 30s, extension at 72 °C for 1 min; and final elongation at 72 °C for 10 min. The product was visualized with ethidium bromide after electrophoresis in 1.2% agarose gel at 120 V for 40 min.

### 4.8. Antifungal Activity Test

The connection between antifungal activity and the presence of VLEs was tested according to the method described by Woods and Bevan [89]. Strains with such elements and cured of them were streaked in a thick line on a lawn of a sensitive *Y. lipolytica* PII6a strain (5 × 10⁵ cells mL⁻¹ of YPD-MB agar pH4.6 or pH 7.0 containing 2.5% NaCl) or *P. roqueforti* PR1 (1 × 10⁴ cells per plate of PDA + 2.5% NaCl pH4.6 and 7.0) and cultivated at 14 °C and 28 °C for 72 h. An analogous experiment was performed with *D. hansenii* strains possessing VLEs streaked on a lawn of *D. hansenii* cured of them. The experiment was performed on YMB pH 7.0 with 2.5% NaCl. The streaked strain was considered a killer yeast if the line of the biomass was surrounded by a clear growth inhibition zone of the sensitive strain. In the case of *P. roqueforti*, changes in sporulation were also inspected. All experiments were run in four replicates.

### 4.9. Halotolerance Test

Strains possessing virus-like elements and cured of them were cultivated on YPD medium with a concentration of NaCl in the range of 0–24%. Cultures were performed in an Automated Microbiology Workstation Bioscreen C (Labsystems Oy, Helsinki, Finland) in a total volume of 350 L in five replications at 28 °C for 96 h with continuous agitation. The BioLink software (Labsystems Oy, Helsinki, Finland) was applied to plot growth curves from optical density measurements performed with a wide band filter (420–580 nm).

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/toxins13090615/s1, Figure S1: The morphology of *D. hansenii* CLIB 197 (CBS 767) type strain colonies on YORS agar medium after 7 days of incubation at 25 °C. Figure S2: The products of colony PCR with *D. hansenii* species-specific pair of primers DhPadF and DhPadR: 1. GeneRuler 100bp Plus (Thermo Fisher Scientific, USA), 2. reagent control, 3. *D. hansenii* CLIB 197T (CBS 767T), 4. *D. hansenii* CLIB 907 (CBS 1795), 5. *D. hansenii* 5c, 6. *D. hansenii* 7g, 7. *D. fabryi* CLIB 422T (CBS 789T), 8. *D. subglobosus* CLIB 908T (CBS 1796T). Figure S3: Comparison of VLEs isolation from yeast biomass obtained after 24h culture in YPD broth inoculated with freeze-dried cells from the collection (2,6) and
from biomass gained after six passages in YPD broth with 12% NaCl (4,8). 2,4. *D. hansenii* CBS 7848; 6,8. *D. hansenii* CBS 770. 1,3,5,7. GeneRuler 1kb DNA Ladder (Thermo Fisher Scientific). Figure S4: The separation of virus like elements (VLEs) containing fraction from nuclear and mitochondrial DNA by centrifugation in CsCl-bisbenzimide gradient. Figure S5: Confirmation of *D. hansenii* cell curing of VLEs using PCR with primers specific for autonomous cytoplasmic DNA elements. Lines 1 and 25: GeneRuler 100 bp Plus (ThermoFisher Scientific). Lines 2-10: PCR with DNA isolated from strains before curing, 2: 4e, 3: 5c, 4: 7g, 5: 7j, 6: 8e, 7: 8g, 8: 8h, 9: 8i, 10: CBS 767, 11-12: CBS 7848 before and after curing, 13-14: CBS 770 before and after curing. Lines 15-22: PCR with DNA from VLE-cured strains, 15: 4e, 16: 5c, 17: 7g, 18: 7j, 19: 8e, 20: 8g, 21: 8h, 22: 8i, 23: PCR with purified pDH4C, 24: reagent control. Figure S6: Killer activity test of *D. hansenii* strains before (horizontal streaks) and after (vertical streaks) curing of VLEs against sensitive *Y. lipolytica* PII6a yeast strain (in lawn) on YPG agar pH 7.0 with 2.5% NaCl at 14°C. (A) 4e, (B) 5c, (C) 7g, (D) 8h, (E) CBS 7848, (F) CBS 770. Figure S7: Killer activity test of *D. hansenii* strains before (horizontal streaks) and after (vertical streaks) curing of VLEs against sensitive *Y. lipolytica* PII6a yeast strain (in lawn) on YPG agar with 2.5% NaCl pH 7.0 at 28 °C. (A) 4e, (B) 5c, (C) 7g, (D) 8h, (E) CBS 7848, (F) CBS 770. Figure S8: Killer activity test of *D. hansenii* strains before (horizontal streaks) and after (vertical streaks) curing of VLEs against sensitive *P. roqueforti* PR1 (in lawn) on PDA with 2.5% NaCl pH 7.0 at 14 °C. (A) 4e, (B) 5c, (C) 7g, (D) 8h, (E) CBS 7848, (F) CBS 770. Figure S9: Killer activity test of *D. hansenii* strains before (horizontal streaks) and after (vertical streaks) curing of VLEs against sensitive *P. roqueforti* PR1 (in lawn) on PDA with 2.5% NaCl pH 7.0 at 28 °C. (A) 4e, (B) 5c, (C) 7g, (D) 8h, (E) CBS 7848, (F) CBS 770. Figure S10: Killer activity test of *D. hansenii* strains before (streaks) and after (in lawn) curing of VLEs on YMB pH 7.0 with 2.5% NaCl at 14 °C. (A) 4e, (B) 5c, (C) 7g, (D) 8h. The same results were obtained on YMB pH 4.6 with 2.5% NaCl at 14 °C. Figure S11: Killer activity of *D. hansenii* strains before (horizontal streaks) and after (vertical streaks) curing of VLEs against sensitive *Y. lipolytica* PII6a yeast strain (in lawn, dead cells colored in blue) on YMB pH 4.6 with 2.5% NaCl at 14 °C for 72h. (A) 4e, (B) 5c, (C) 7g, (D) 8h. Figure S12: Antifungal activity of *D. hansenii* strains before (vertical streaks) and after (horizontal streaks) curing of linear plasmids against sensitive *P. roqueforti* PR1 (in lawn) on PDA + 2.5% NaCl pH 4.6 at 14 °C. (A) 4e left and up, 5c down and right, (B) 7g left and down, CBS 7848 up and right, (C) 8h left and up, CBS 770, (D) control plate. Figure S13: The growth of *D. hansenii* 5c strain before (red line) and after (yellow line) VLEs curing in YPD medium with different NaCl concentration. Cultures performed in BioscreenC plotting OD value using wide band filter of range 420–580 nm.

**Author Contributions:** Conceptualization, X.P.; methodology, X.P. and C.N.; investigation, X.P., C.N., J.Z. and B.Ż.; resources, X.P., C.N. and S.C.; data curation, X.P.; writing—original draft preparation, X.P.; writing—review and editing, X.P., C.N., S.C. and Z.L.; project administration, X.P. and C.N.; funding acquisition, X.P. and C.N. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Polish Ministry of Science and Higher Education [Project N N 312 258138] as well as the French "ministère de l'Europe et des Affaires étrangères (MEAE)" and "ministère de l'Enseignement supérieur, de la Recherche et de l'Innovation (MESRI)" for international cooperation in the framework of PHC program Polonium.

**Data Availability Statement:** The sequences of pDH4A/B/C, pDH5A/B and pDH6A/B are available in GenBank (https://www.ncbi.nlm.nih.gov/) under the following accession numbers: MF795091/MF795092/MF795093, KX904874/KX904875 and KX858805/KX904876, respectively.

**Acknowledgments:** The authors would like to thank Noémie Jacques for providing strains from the CIRM-Levures collection.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
## References

1. Fukuhara, H. Linear DNA Plasmids of Yeasts. *FEMS Microbiol. Lett.* **1995**, *131*. [CrossRef]

2. Meinhardt, F.; Klassen, R. *Microbial Linear Plasmids*; Meinhardt, F., Klassen, R., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; Volume 7, ISBN 978-3-540-72024-9.

3. Gunge, N.; Sakaguchi, K. Intergeneric Transfer of Deoxyribonucleic Acid Killer Plasmids, pGK11 and pGK12, from *Kluyveromyces lactis* into *Saccharomyces cerevisiae* by Cell Fusion. *J. Bacteriol.* **1981**, *147*, 155-160. [CrossRef] [PubMed]

4. Kitada, K.; Hishinuma, F. A New Linear DNA Plasmid Isolated from the Yeast *Saccharomyces kluyveri*. *MGG* **1987**, *206*, 377-381. [CrossRef]

5. Shepherd, H.S.; Ligon, J.M.; Bolen, P.L.; Kurtzman, C.P. Cryptic DNA Plasmids of the Heterothallic Yeast *Saccharomycopsis crataegensis*. *Curr. Genet.* **1987**, *12*, 297-304. [CrossRef]

6. Ligon, J.M.; Bolen, P.L.; Hill, D.S.; Bothast, R.J.; Kurtzman, C.P. Physical and Biological Characterization of Linear DNA Plasmids of the Yeast *Pichia inositovora*. *Plasmid* **1989**, *21*, 185-194. [CrossRef]

7. Worsham, P.L.; Bolen, P.L. Killer Toxin Production in Pichia Acaciae Is Associated with Linear DNA Plasmids. *Curr. Genet.* **1990**, *18*, 77-80. [CrossRef]

8. Gunge, N.; Fukuda, K.; Morikawa, S.; Murakami, K.; Takeda, M.; Miwa, A. Osmophilic Linear Plasmids from the Salt-Tolerant Yeast *Debaryomyces hansenii*. *Curr. Genet.* **1993**, *23*, 443-449. [CrossRef]

9. Cong, Y.-S.; Yarrow, D.; Li, Y.-Y.; Fukuhara, H. Linear DNA Plasmids from *Pichia etchellsii*, *Debaryomyces hansenii* and *Wingea robertsiae*. *Microbiology* **1994**, *140*, 1327-1335. [CrossRef]

10. Hagenson, M.J.; Barr, K.A.; Stroman, D.W.; Harpold, M.M.; Klein, R.D.; Gaertner, F.H. *Pichia pastoris* Linear Plasmids and DNA Fragments Thereof. European Patent Office. European Patent Application: 92115908.3, March 1993.

11. Banerjee, H.; Kopvak, C.; Curley, D. Identification of Linear DNA Plasmids of the Yeast *Pichia pastoris*. *Plasmid* **1998**, *40*, 58-60. [CrossRef]

12. Blaisonneau, J.; Nosek, J.; Fukuhara, H. Linear DNA plasmid pPK2 of *Pichia kluyveri*: Distinction between cytoplasmic and mitochondrial linear plasmids in yeasts. *Yeast* **1999**, *15*. [CrossRef]

13. Chen, W.-B.; Han, Y.-F.; Jong, S.-C.; Chang, S.-C. Isolation, Purification, and Characterization of a Killer Protein from *Schwanniomyces occidentalis*. *Appl. Environ. Microbiol.* **2000**, *66*, 5348-5352. [CrossRef] [PubMed]

14. Liu, X.-Z.; Wang, Q.-M.; Göker, M.; Groenewald, M.; Kachalkin, A.V.; Lumbsch, H.T.; Millanes, A.M.; Wedin, M.; Yurkov, A.M.; Boekhout, T.; et al. Towards an Integrated Phylogenetic Classification of the Tremellomycetes. *Stud. Mycol.* **2015**, *81*, 85-147. [CrossRef] [PubMed]

15. Kurtzman, C.P.; Robnett, C.J. Synonymy of the Yeast *Genera* *Wingea* and *Debaryomyces*. *Antonie Leeuwenhoek* **1994**, *66*, 337-342. [CrossRef]

16. Groenewald, M.; Daniel, H.-M.; Robert, V.; Poot, G.A.; Smith, M.T. Polyphasic Re-Examination of *Debaryomyces Hansenii* Strains and Reinstatement of *D. hansenii*, *D. fabryi* and *D. subglobosus*. *Pers. Mol. Phylogeny Evol. Fungi* **2008**, *21*, 17-27. [CrossRef] [PubMed]

17. Kurtzman, C.P.; Suzuki, M. Phylogenetic Analysis of Ascomycete Yeasts That Form Coenzyme Q-9 and the Proposal of the New Genera *Babjeviella*, *Meyerozyma*, *Millerozyma*, *Priceomyces*, and *Scheffersomyces*. *Mycoscience* **2010**, *51*, 2-14. [CrossRef]

18. Gunge, N. Linear DNA Killer Plasmids from the Yeast *Kluyveromyces*. *Yeast* **1986**, *2*, 153-162. [CrossRef]

19. Stam, J.C.; Kwakman, J.; Meijer, M.; Stuitje, A.R. Efficient Isolation of the Linear DNA Killer Plasmid of *Kluyveromyces Lactis*: Evidence for Location and Expression in the Cytoplasm and Characterization of Their Terminally Bound Proteins. *Nucleic. Acids Res.* **1986**, *14*, 6871-6884. [CrossRef]

20. Frank, A.C.; Wolfe, K.H. Evolutionary Capture of Viral and Plasmid DNA by Yeast Nuclear Chromosomes. *Eukaryot. Cell* **2009**, *8*, 1521-1531. [CrossRef]

21. Kast, A.; Voges, R.; Schroth, M.; Schaffrath, R.; Klassen, R.; Meinhardt, F. Autoselection of Cytoplasmic Yeast Virus Like Elements Encoding Toxin/Antitoxin Systems Involves a Nuclear Barrier for Immunity Gene Expression. *PLoS Genet.* **2015**, *11*. [CrossRef]

22. Tommasino, M.; Ricci, S.; Galeotti, C.L. Genome Organization of the Killer Plasmid PGK12 from Kluyveromyces Lactis. *Nucleic Acids Res.* **1988**, *16*, 5863-5878. [CrossRef]

23. Wilson, D.W.; Meacock, P.A. Extranuclear Gene Expression in Yeast: Evidence for a Plasmidencoded RNA Polymerase of Unique Structure. *Nucleic Acid Res.* **1988**, *16*, 8097-8112.

24. Klassen, R.; Jablonowski, D.; Schaffrath, R.; Meinhardt, F. Genome Organization of the Linear Pichia Etchellsii Plasmid PPE1A: Evidence for Expression of an Extracellular Chitin-Binding Protein Homologous to the α-Subunit of the Kluyveromyces Lactis Killer Toxin. *Plasmid* **2002**, *47*, 224-233. [CrossRef]

25. Jeske, S.; Meinhardt, F. Autonomous Cytoplasmic Linear Plasmid pPac1-1 of *Pichia acaciae*: Molecular Structure and Expression Studies. *Yeast* **2006**, *23*, 479-486. [CrossRef]

26. Schaffrath, R.; Meacock, P.A. An SSB Encoded by and Operating on Linear Killer Plasmids From *Kluyveromyces lactis*. *Yeast* **2001**, *18*, 1239-1247. [CrossRef]

27. Stark, M.J.R.; Boyd, A.; Pritchard, R.H. The Killer Toxin of *Kluyveromyces laclis*: Characterization of the Toxin Subunits and Identification of the Genes Which Encode Them. *EMBO J.* **1986**, *5*, 1995-2002. [CrossRef]

28. Klassen, R.; Meinhardt, F. Linear Plasmids pWR1A and pWR1B of the Yeast *Wingea robertsiae* Are Associated with a Killer Phenotype. *Plasmid* **2002**, *48*, 142-148. [CrossRef]
29. Klassen, R.; Teichert, S.; Meinhardt, F. Novel Yeast Killer Toxins Provoke S-Phase Arrest and DNA Damage Checkpoint Activation. *Mol. Microbiol.* **2004**, *53*, 263–273. [CrossRef]

30. Stark, M.J.R.; Boyd, A.; Mileham, A.J.; Ramonos, M.A. The Plasmid--encoded Killer System of *Kluyveromyces lactis*: A Review. *Yeast* **1990**, *6*, 1–29. [CrossRef] [PubMed]

31. Gardini, F.; Suzzi, G.; Lombardi, A.; Galgano, F.; Crudele, M.A.; Andrighetto, C.; Schirone, M.; Tofalo, R. A Survey of Yeasts in Traditional Sausages of Southern Italy. *FEMS Yeast Res.* **2001**, *1*. [CrossRef] [PubMed]

32. Simoncini, N.; Rotelli, D.; Virgili, R.; Quintavalla, S. Dynamics and Characterization of Yeasts during Ripening of Typical Italian Dry-Cured Ham. *Food Microbiol.* **2007**, *24*, 577–584. [CrossRef] [PubMed]

33. Brežná, B.; Ženišová, K.; Chovanová, K.; Chebeňová, V.; Kraková, L.; Kuchta, T.; Pangallo, D. Evaluation of Fungal and Yeast Diversity in Slovakian Wine-Related Microbial Communities. *Antonie Leeuwenhoek* **2010**, *98*, 519–529. [CrossRef]

34. Jacques, N.; Zenouche, A.; Gunde-Cimerman, N.; Casaregola, S. Increased Diversity in the Genus *Debaryomyces* from Arctic Glacier Samples. *Antonie Leeuwenhoek* **2015**, *107*, 487–501. [CrossRef] [PubMed]

35. Ricci, A.; Allende, A.; Bolton, D.; Chemaly, M.; Davies, R.; Girones, R.; Koutsoumanis, K.; Lindqvist, R.; Nørrung, B.; Robertson, L.; et al. Update of the List of QPS-recommended Biological Agents Intentionally Added to Food or Feed as Notified to EFSA 6: Suitability of Taxonomic Units Notified to EFSA until March 2017. *EFSA J. 2017*, *15*, e04884. [CrossRef]

36. Wagner, D.; Sander, A.; Bertz, H.; Finke, J.; Kern, W.V. Breakthrough Invasive Infection Due to *Debaryomyces hansenii* (Teleomorph *Candida famata*) and Scopulariopsis Brevicaulis in a Stem Cell Transplant Patient Receiving Liposomal Amphotericin B and Caspofungin for Suspected Aspergillosis. *Infection* **2005**, *33*, 397–400. [CrossRef] [PubMed]

37. Chan, G.F.; Sinniah, S.; Tengku Idr, T.I.N.; Ahmad Puad, M.S.; Abd Rahman, A.Z. Multiple Rare Opportunistic and Pathogenic Fungi in Persistent Foot Skin Infection. *Pak. J. Biol. Sci.* **2013**, *16*, 208–218. [CrossRef] [PubMed]

38. Van den Tempel, T.; Jakobsen, M. The Technological Characteristics of *Debaryomyces hansenii* and *Yarrowia lipolytica* and Their Potential as Starter Cultures for Production of Danablu. *Int. Dairy J.* **2000**, *10*, 263–270. [CrossRef]

39. Żarowska, B.; Wojtatowicz, M.; Połomska, X.; Juszczyk, P.; Chrzanowska, J. Factors Affecting Killer Activity of Some Yeast Species Occurring in Rokpol Cheese. *Folia Microbiol.* **2004**, *49*, 713–717. [CrossRef]

40. Fröhlich-Wyder, M.-T.; Arias-Roth, E.; Jakob, E. Cheese Yeasts. *Yeast* **2019**, *36*, 129–141. [CrossRef] [PubMed]

41. Santos, A.; Marquina, D.; Barroso, J.; Peinado, J.M. (1-6)-Beta-D-Glucan as the Cell Wall Binding Site for *Debaryomyces hansenii* Killer Toxin. *Lett. Appl. Microbiol.* **2002**, *34*, 95–99. [CrossRef] [PubMed]

42. Żarowska, B. *Biosynteza i Charakterystyka Toksyn Killerowych Drożdży Debaryomyces Hansenii*; Wydawnictwo Uniwersytetu Przyrodniczego we Wrocławiu: Wrocław, Poland, 2012; pp. 34–49.

43. Çorbacı, C.; Uçar, F.B. Purification, Characterization and in Vivo Biocontrol Efficiency of Killer Toxins from *Debaryomyces hansenii* Strains. *Int. J. Biol. Macromol.* **2018**, *119*, 1077–1082. [CrossRef]

44. Wojtatowicz, M.; Chrzanowska, J.; Juszczyk, P.; Skiba, A.; Gdula, A. Identification and Biochemical Characteristics of Yeast Microflora of Rokpol Cheese. *Int. J. Food Microbiol.* **2001**, *69*, 135–140. [CrossRef]

45. Połomska, X.; Wojtatowicz, M.; Zarowska, B.; Szołtysik, M.; Chrzanowska, J. Freeze-Drying Preservation of Yeast Adjunct Cultures for Cheese Production. *Pol. J. Food Nutr. Sci.* **2012**, *62*, 143–150. [CrossRef]

46. Szołtysik, M.; Dąbrowska, A.; Babij, K.; Pokora, M.; Zambrowicz, A.; Połomska, X.; Wojtatowicz, M.; Chrzanowska, J. Biochemical and Microbiological Changes in Cheese Inoculated with *Yarrowia lipolytica* Yeast. *Zywnosc. Nauka. Technol. Jakosc/Food. Sci. Technol. Qual.* **2013**, *89*, 49–64. [CrossRef]

47. Grzegorczyk, M.; Żarowska, B.; Restuccia, C.; Cirvilleri, G. Postharvest Biocontrol Ability of Killer Yeasts against *Monilinia fructigena* and *Monilinia fructicola* on Stone Fruit. *Food Microbiol.* **2017**, *61*. [CrossRef]

48. Czarnecka, M.; Żarowska, B.; Połomska, X.; Restuccia, C.; Cirvilleri, G. Role of Biocontrol Yeasts *Debaryomyces hansenii* and *Wickerhamomyces anomalus* in Plants’ Defence Mechanisms against *Monilinia fructicola* in Apple Fruits. *Food Microbiol.* **2019**, *83*, 1–8. [CrossRef]

49. Żarowska, B.; Połomska, X.; Grzegorczyk, M.; Regiec, P.; Gudarowska, E.; Czaplicka-Pędzich, M.; Sosna, I.; Figiel, A.; Pasławska, M.; Serowik, M.; et al. Method for obtaining a preparation containing killer toxins. Patent of Republic of Poland Pat.229949, September 2018.

50. Wrent, P.; Rivas, E.-M.; Gil de Prado, E.; Peinado, J.M.; de Silóniz, M.-I. Development of Species-Specific Primers for Rapid Identification of *Debaryomyces hansenii*. *Int. J. Food Microbiol.* **2015**, *193*, 109–113. [CrossRef] [PubMed]

51. O'Donnel, K. Fusarium and Its near Relatives. In *The Fungal Holomorph: Mitotic, Meiotic and Pleomorphic Speciation in Fungal Systematics*; Reynolds, D.R., Taylor, J.W., Eds.; CAB International: Wallingford, UK, 1993; pp. 225–233.

52. Romero, P.; Patino, B.; Quiros, M.; Gonzalezjaen, M.; Valderrama, M.; Desiloniz, M.; Penaido, J. Differential Detection of Isolated from Intermediate-Moisture Foods by PCR-RFLP of the IGS Region of rDNA. *FEMS Yeast Res.* **2005**, *5*, 455–461. [CrossRef] [PubMed]

53. Nguyen, H.-V.; Gaillardin, C.; Neuveglise, C. Differentiation of *Debaryomyces hansenii* and *Candida famata* by rRNA Gene Intergenic Spacer Fingerprinting and Reassessment of Phylogenetic Relationships among *D. hansenii*, *C. famata*, *D. fabryi*, *C. flareri* (*D. subglobosus*) and *D. prosopidis*: Description of *D. vietnamensis* sp. Nov. Closely Related to *D. nepalensis*. *FEMS Yeast Res.* **2009**, *9*, 641–662. [CrossRef]

54. Fukuda, K.; Jinshan, C.; Kawano, M.; Sudo, K.; Gunge, N. Stress Responses of Linear Plasmids From. *FEMS Microbiol. Lett.* **2004**, *237*, 243–248. [CrossRef]
55. Hishinuma, F.; Hirai, K. Genome Organization of the Linear Plasmid, pSKL, Isolated from *Saccharomyces kluyveri*. *MGG* **1991**, *226*, 97–106. [CrossRef]

56. Tommasino, M. Killer system of *Kluyveromyces lactis*: The open reading frame 10 of the pGKL2 plasmid encodes a putative DNA binding protein. *Yeast* **1991**, *7*, 245–252. [CrossRef]

57. Klassen, R.; Meinhardt, F. Structural and Functional Analysis of the Killer Element PPin1-3 from *Pichia inositovora*. *MGG* **2003**, *270*, 190–199. [CrossRef]

58. Satwika, D.; Klassen, R.; Meinhardt, F. Repeated Capture of a Cytoplasmic Linear Plasmid by the Host Nucleus in *Debaryomyces hansenii*. *Yeast* **2012**, *29*, 145–154. [CrossRef]

59. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

60. Dujon, B.; Sherman, D.; Fischer, G.; Durrens, P.; Casaregola, S.; Lafontaine, I.; de Montigny, J.; Marck, C.; Neuvéglise, C.; Talla, E.; et al. Genome Evolution in Yeasts. *Nature* **2004**, *430*, 35–44. [CrossRef]

61. Neuveglise, C.; Feldmann, H.; Bon, E.; Gaillardin, C.; Casaregola, S. Genomic Evolution of the Long Terminal Repeat Retrotransposons in Hemiascomycetous Yeasts. *Genome Res.* **2002**, *12*, 930–943. [CrossRef]

62. Larsen, M.; Meinhardt, F. *Kluyveromyces lactis* Killer System: Identification of a New Gene Encoded by pGKL2. *Curr. Genet.* **2000**, *38*, 271–275. [CrossRef]

63. Kempken, F.; Hermanns, J.; Osiewacz, H.D. Evolution of Linear Plasmids. *J. Mol. Evol.* **1992**, *35*, 502–513. [CrossRef] [PubMed]

64. Meinhardt, F.; Klassen, R. Linear protein-primed replicating plasmids in eukaryotic microbes. In *Microbial Linear Plasmids*; Meinhardt, F., Klassen, R., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; Volume 7, pp. 188–227. ISBN 978-3-540-72024-9.

65. Sýkora, M.; Pospíšek, M.; Novák, J.; Mrvová, S.; Krásný, L.; Vopálenský, V. Transcription Apparatus of the Yeast Virus-like Elements: Architecture, Function, and Evolutionary Origin. *PLoS Pathog.* **2018**, *14*, e1007377. [CrossRef] [PubMed]

66. Vopálenský, V.; Sýkora, M.; Mašek, T.; Pospíšek, M. Messenger RNAs of Yeast Virus-Like Elements Contain Non-Templated 5' Poly(A) Leaders, and Their Expression Is Independent of EIF4E and Pab1. *Front Microbiol.* **2019**, *10*, 2366. [CrossRef]

67. Gish, W.; States, D.J. Identification of Protein Coding Regions by Database Similarity Search. *Nat. Genet.* **1993**, *3*, 266–272. [CrossRef]

68. Marchler-Bauer, A.; Derbyshire, M.K.; Gonzales, N.R.; Lu, S.; Chitsaz, F.; Geer, L.Y.; Geer, R.C.; He, J.; Gwadz, M.; Hurwitz, D.I.; et al. CDD: NCBI's Conserved Domain Database. *Nucleic Acids Res.* **2015**, *43*, D222–D226. [CrossRef]

69. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam Protein Families Database: Towards a More Sustainable Future. *Nucleic Acids Res.* **2016**, *44*, D279–D285. [CrossRef]

70. Schaffrath, R.; Meacock, P.A.; Meinhardt, F. Yeast Killer Plasmid pGKL2: Molecular Analysis of UCS5, a Cytoplasmic Promoter Element Essential for ORF5 Gene Function. *MGG* **1996**, *250*, 286–294. [CrossRef]

71. Meinhardt, F.; Wodara, C.; Larsen, M.; Schickel, J. A Novel Approach to Express a Heterologous Gene on Kluyveromyces Lactis Linear Killer Plasmids: Expression of the Bacterial Aph Gene from a Cytoplasmic Promoter Fragment without In-Phase Fusion to the Plasmid Open Reading Frame. *Plasmid* **1994**, *32*, 318–327. [CrossRef]

72. Schrüberl, J.; Meinhardt, F. Extrachromosomal Inheritance: Yeast Linear Killer Plasmids as a Tool in Genetic Engineering. In *Progress in Botany*; Springer: Berlin/Heidelberg, Germany, 1995.

73. Schründer, J.; Gunge, N.; Meinhardt, F. Extranuclear Expression of the Bacterial Xylose Isomerase (Xyl A) and the UDP-Glucose Dehydrogenase (Has B) Genes in Yeast with *Kluyveromyces lactis* Linear Killer Plasmids as Vectors. *Curr. Microbiol.* **1996**, *33*, 323–330. [CrossRef]

74. Schickel, J.; Helmig, C.; Meinhardt, F. Kluyveromyces Lactis Killer System: Analysis of Cytoplasmic Promoters of the Linear Plasmids. *Nucleic Acids Res.* **1996**, *24*, 1879–1886. [CrossRef] [PubMed]

75. Gunge, N.; Tokunaga, M. Linear DNA Plasmids and Killer System of *Kluyveromyces lactis*. In *Genetics and Biotechnology*; Springer: Berlin/Heidelberg, Germany, 2004.

76. Mesnage, S.; Dellarole, M.; Baxter, N.J.; Rouget, J.-B.; Dimitrov, J.D.; Wang, N.; Fujimoto, Y.; Hounslow, A.M.; Lacroix-Desmazes, S.; Fukase, K.; et al. Molecular Basis for Bacterial Peptidoglycan Recognition by LysM Domains. *Nat. Commun.* **2014**, *5*, 4269. [CrossRef] [PubMed]

77. Almagro Armenteros, J.J.; Sønderby, C.K.; Sønderby, S.K.; Nielsen, H.; Winther, O. DeepLoc: Prediction of Protein Subcellular Localization Using Deep Learning. *Bioinformatics* **2017**, *33*, 3387–3395. [CrossRef]

78. Klassen, R.; Wemhoff, S.; Krause, J.; Meinhardt, F. DNA Repair Defects Sensitize Cells to Anticodon Nuclease Yeast Killer Toxins. *MGG* **2011**, *285*, 185–195. [CrossRef] [PubMed]

79. Wemhoff, S.; Klassen, R.; Meinhardt, F. DNA Damage Induced by the Anticodon Nuclease from a *Pichia acaciae* Killer Strain Is Linked to Ribonucleotide Reductase Depletion. *Cell. Microbiol.* **2016**, *18*, 211–222. [CrossRef] [PubMed]

80. Wemhoff, S.; Klassen, R.; Beetz, A.; Meinhardt, F. DNA Damage Responses Are Induced by tRNA Anticodon Nucleases and Hygromycin, B. *PLoS ONE* **2016**, *11*, e0157611. [CrossRef]

81. Schaffrath, R.; Meinhardt, F.; Klassen, R. Yeast Killer Toxins: Fundamentals and Applications. In *The Mycota. Physiology and Genetics: Selected Basic and Applied Aspects*; Anke, T., Schüffler, A., Eds.; Springer International Publishing: New York, NY, USA, 2018; pp. 87–118.
82. Kall, L.; Krogh, A.; Sonnhammer, E.L.L. Advantages of Combined Transmembrane Topology and Signal Peptide Prediction—the Phobius Web Server. *Nucleic. Acids Res* **2007**, *35*, W429–W432. [CrossRef]

83. Vega, R.; Dominguez, A. Cell Wall Composition of the Yeast and Mycelial Forms of *Yarrowia lipolytica*. *Arch. Microbiol.* **1986**, *144*, 124–130. [CrossRef]

84. Andriyanova, D.A.; Smirnova, G.P.; Shashkov, A.S.; Chizhov, A.O.; Galanina, L.A.; Feofilova, E.P.; Usov, A.I. Polysaccharide Composition of Mycelium and Cell Walls of the Fungus *Penicillium roqueforti*. *Russ. J. Bioorganic. Chem.* **2011**, *37*, 356–363. [CrossRef]

85. McCracken, D.A.; Martin, V.J.; Stark, M.J.R.; Bolen, P.L. The Linear-Plasmid-Encoded Toxin Produced the Yeast *Pichia acaciae*: Characterization an Comparison with the Toxin of *Kluyveromyces lactis*. *Microbiology* **1994**, *140*, 425–431. [CrossRef]

86. Çorbaci, C.; Uçar, F.B. Production and Optimization of Killer Toxin in *Debaryomyces hansenii* Strains. *Braz. Arch. Biol. Technol.* **2017**, *60*. [CrossRef]

87. Untergasser, A.; Nijveen, H.; Rao, X.; Bisseling, T.; Geurts, R.; Leunissen, J.A.M. Primer3Plus, an Enhanced Web Interface to Primer3. *Nucleic. Acids Res* **2007**, *35*, W71–W74. [CrossRef]

88. Querol, A.; Barrio, E.; Ramón, D. A Comparative Study of Different Methods of Yeast Strain Characterization. *Syst. Appl. Microbiol.* **1992**, *15*, 439–446. [CrossRef]

89. Woods, D.R.; Bevan, E.A. Studies on the Nature of the Killer Factor Produced by *Saccharomyces cerevisiae*. *J. Gen. Microbiol.* **1968**, *51*, 115–126. [CrossRef]