Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

A Novel Virus Alters Gene Expression and Vacuolar Morphology in Malassezia Cells and Induces a TLR3-Mediated Inflammatory Immune Response

Minji Park, Yong-Joon Cho, Donggyu Kim, Chul-Su Yang, Shi Mun Lee, Thomas L. Dawson Jr., Satoshi Nakamizo, Kenji Kabashima, Yang Won Lee, Won Hee Jung
Floyd L. Wormley Jr., Editor
Minji Park
aDepartment of Systems Biotechnology, Chung-Ang University, Anseong, South Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yong-Joon Cho
bThe Research Institute of Basic Sciences, Seoul National University, Seoul, South Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Donggyu Kim
cDepartment of Molecular and Life Science, Hanyang University, Ansan, South Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chul-Su Yang
cDepartment of Molecular and Life Science, Hanyang University, Ansan, South Korea
dDepartment of Bionano Technology, Hanyang University, Seoul, South Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shi Mun Lee
eSingapore Immunology Network (SIgN), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
jSkin Research Institute of Singapore (SRIS), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Thomas L. Dawson Jr.
eSingapore Immunology Network (SIgN), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
fCenter for Cell Death, Injury & Regeneration, Medical University of South Carolina, Charleston, South Carolina, USA
jSkin Research Institute of Singapore (SRIS), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
kDepartment of Drug Discovery & Biomedical Sciences, Medical University of South Carolina, Charleston, South Carolina, USA
lDepartment of Biochemistry & Molecular Biology, Medical University of South Carolina, Charleston, South Carolina, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Satoshi Nakamizo
eSingapore Immunology Network (SIgN), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
jSkin Research Institute of Singapore (SRIS), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kenji Kabashima
eSingapore Immunology Network (SIgN), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
gDepartment of Dermatology, Kyoto University Graduate School of Medicine, Kyoto, Japan
jSkin Research Institute of Singapore (SRIS), Agency for Science, Technology and Research (A*STAR), Biopolis, Singapore
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yang Won Lee
hDepartment of Dermatology, School of Medicine, Konkuk University, Seoul, South Korea
iResearch Institute of Medicine, Konkuk University, Seoul, South Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Won Hee Jung
aDepartment of Systems Biotechnology, Chung-Ang University, Anseong, South Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Won Hee Jung
Floyd L. Wormley Jr.
Texas Christian University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.01521-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Most fungal viruses have been identified in plant pathogens, whereas the presence of viral particles in human-pathogenic fungi is less well studied. In the present study, we observed extrachromosomal double-stranded RNA (dsRNA) segments in various clinical isolates of Malassezia species. Malassezia is the most dominant fungal genus on the human skin surface, and species in this group are considered etiological factors of various skin diseases including dandruff, seborrheic dermatitis, and atopic dermatitis. We identified novel dsRNA segments, and our sequencing results revealed that the virus, named MrV40, belongs to the Totiviridae family and contains an additional satellite dsRNA segment encoding a novel protein. The transcriptome of virus-infected Malassezia restricta cells was compared to that of virus-cured cells, and the results showed that transcripts involved in ribosomal biosynthesis were downregulated and those involved in energy production and programmed cell death were upregulated. Moreover, transmission electron microscopy revealed significantly larger vacuoles in virus-infected M. restricta cells, indicating that MrV40 infection dramatically altered M. restricta physiology. Our analysis also revealed that viral nucleic acid from MrV40 induced a TLR3 (Toll-like receptor 3)-mediated inflammatory immune response in bone marrow-derived dendritic cells, suggesting that a viral element contributes to the pathogenicity of Malassezia.

IMPORTANCE Malassezia is the most dominant fungal genus on the human skin surface and is associated with various skin diseases including dandruff and seborrheic dermatitis. Among Malassezia species, Malassezia restricta is the most widely observed species on the human skin. In the current study, we identified a novel dsRNA virus, named MrV40, in M. restricta and characterized the sequence and structure of the viral genome along with an independent satellite dsRNA viral segment. Moreover, expression of genes involved in ribosomal synthesis and programmed cell death was altered, indicating that virus infection affected the physiology of the fungal host cells. Our data also showed that the viral nucleic acid from MrV40 induces a TLR3-mediated inflammatory immune response in bone marrow-derived dendritic cells, indicating that a viral element likely contributes to the pathogenicity of Malassezia. This is the first study to identify and characterize a novel mycovirus in Malassezia.

INTRODUCTION

Viruses have been observed in many fungal species since their first identification in mushrooms (1). Fungal viruses, also known as mycoviruses, possess different forms of viral genomes including double-stranded RNA (dsRNA), single-stranded RNA (ssRNA), and single-stranded DNA (ssDNA). It is estimated that 30 to 80% of all fungal species, mainly endophytic fungi, are infected with viruses. Unlike viruses that infect other organisms, the transmission of fungal virus occurs vertically by cell division or horizontally via mating or hyphal anastomosis, with no extracellular phase of the virus life cycle. dsRNA segments have predominantly been found in fungal viruses, and taxonomically, the fungal dsRNA viruses are classified into seven families: Chrysoviridae, Endornaviridae, Megabirnaviridae, Quadriviridae, Partitiviridae, Reoviridae, and Totiviridae (2).

The model fungus Saccharomyces cerevisiae also carries a dsRNA virus that belongs to the Totiviridae family and is known as the L-A virus. A unique feature of fungal viruses of the Totiviridae family is their capability to produce a killer toxin that lyses susceptible neighbor strains, whereas the virus-containing strain (also known as the killer strain) is immune to the toxin. Studies of how the virus produces killer toxins in S. cerevisiae showed that the killer toxins are encoded by a satellite dsRNA segment, known as the M satellite, within the L-A virus. To date, four different killer toxins, K1, K2, K28, and Klus, have been described (3–6). The S. cerevisiae L-A virus forms icosahedral particles with a diameter of approximately 39 nm (7). The virus possesses a nonsegmented 4.6-kb dsRNA genome consisting of two open reading frames (ORFs), gag and pol, which overlap by 130 bp (8). A major 76-kDa capsid protein is encoded by gag, and a 180-kDa minor protein species is encoded as a Gag-Pol fusion protein by a −1 ribosomal frameshift event (9, 10). The C terminus of the Gag-Pol fusion protein possesses viral RNA-dependent RNA polymerase (RDRP) activity (8). The ribosomal frameshift is an interesting feature in a compact viral genome and has been commonly found in various viral genomes as a mechanism to allow viruses to express overlapping ORFs. Studies of S. cerevisiae L-A virus revealed that the mechanism of frameshifting is based on the sequence structures including a canonical slippery heptamer, 5′-X XXY YYZ-3′ (X = A, U, or G; Y = A or U; Z = A, U, or C) and RNA pseudoknot (11).

Fungal viruses have also been considered biocontrol agents in the field of agriculture. For example, a virus causes hypovirulence in the chestnut blight fungus Cryphonectria parasitica (12, 13), and a virus mediates the biocontrol of other phytopathogenic fungi such as Helminthosporium victoriae, Sclerotinia sclerotiorum, and Botrytis cinerea (14).

Although viral infections in fungal cells are widespread, the interactions between the fungal virus and its host are not well understood. One of the best-studied host defense mechanisms is RNA silencing. Several studies have shown that RNA silencing functions as an antiviral defense mechanism against C. parasitica in fungi. Disruption of the dicer pathway in C. parasitica increases the susceptibility to virus infections (15), and p29 was identified as a suppressor that inhibits expression of the genes required for RNA silencing-mediated viral defense in the fungus (16). Similarly, conserved RNA silencing-mediated antiviral defense systems have been identified in Aspergillus nidulans, Rosellinia necatrix, and Fusarium graminearum (17–19).

Malassezia is the most dominant fungal genus on the human skin surface and is considered an etiological factor in various skin diseases including dandruff, seborrheic dermatitis, and atopic dermatitis (20–23). Eighteen Malassezia species have been identified; among them, Malassezia restricta is the most abundant on the human skin (20). Recent studies showed an increased burden of M. restricta on the scalp of patients with dandruff, indicating an association between dandruff and the fungus, although its role as a pathogenic organism is still unclear, and host susceptibility should be taken into consideration (21, 24–27). Most fungal viruses are found in plant-pathogenic fungi, whereas few examples of viral particles have been identified in human-pathogenic fungi such as Candida albicans (28). In the present study, we observed extrachromosomal dsRNA segments in various M. restricta clinical isolates which represented a novel viral genome and its satellite. Sequence analysis revealed that the virus belongs to the Totiviridae family and that the additional satellite dsRNA segment encodes a novel protein. The interactions between the viral elements and the fungal host and the impact of the virus on fungal interactions with immune cells were evaluated.

RESULTS

Identification of extrachromosomal dsRNA segments in Malassezia.Extrachromosomal nucleic acid bands were observed in total nucleic acid extracts of the M. restricta strains isolated in our recent study (29). Among the strains, M. restricta KCTC 27540 was used to identify extrachromosomal segments. Total nucleic acids were extracted from the strain and digested with DNase I, RNase A, and RNase T1. The extrachromosomal segments and rRNA were resistant to DNase I, indicating that they were neither ssDNA nor dsDNA. RNase A degraded all nucleic acids except for genomic DNA, whereas RNase T1, which catalyzes the degradation of ssRNA, removed rRNA only (30). These results suggested that the extrachromosomal segments correspond to dsRNA and that M. restricta KCTC 27540 possesses two separate extrachromosomal segments as estimated by agarose gel electrophoresis to be approximately 4.5 and 1.5 kb (Fig. 1A).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Extrachromosomal dsRNA segments in Malassezia. (A) Nucleic acids from Malassezia restricta KCTC 27540 were separated on a 0.7% agarose gel. Lane 1, total nucleic acids; lane 2, total nucleic acids treated with DNase I; lane 3, total nucleic acids treated with RNase A; lane 4, total nucleic acids treated with RNase T1. (B) Presence of dsRNA in M. restricta strains. Total nucleic acids were extracted and treated with DNase I and RNase T1 and were separated on a 0.7% agarose gel. (C) Total nucleic acids were extracted from different Malassezia species, treated with DNase I and RNase T1, and separated on a 0.7% agarose gel.

To confirm whether the extrachromosomal segments observed in other M. restricta clinical isolates were also dsRNA, total nucleic acid extracts from strains other than M. restricta KCTC 27540 were treated with DNase I and RNase T1. The extrachromosomal segments remained unaffected after enzyme treatment, indicating that they are also dsRNA, as in M. restricta KCTC 27540 (Fig. 1B). Moreover, other Malassezia species including M. globosa, M. pachydermatis, and M. sympodialis showed similar extrachromosomal dsRNA segments, suggesting that these segments are common characteristics of Malassezia species (Fig. 1C). Agarose gel electrophoresis revealed extrachromosomal segments composed of at least two separate dsRNA fragments except for M. restricta KCTC 27543, which showed a single dsRNA fragment. Additionally, the large fragments of dsRNA showed similar sizes (∼5.0 kb), whereas the small dsRNA segments varied in size in different strains (Fig. 1B and C). We hypothesized that the dsRNA segments from Malassezia strains represent the dsRNA elements of mycoviruses, which are prevalent in all major fungal taxa (2). Indeed, a dsRNA mycovirus-infected soil fungus, Trichoderma atroviride NFCF028, displayed similar extrachromosomal segments (see Fig. S1 in the supplemental material) (31).

FIG S1

Extrachromosomal dsRNA segments in Malassezia restricta KCTC 27540 and T. atroviride NFCF028. Nucleic acids from the strains were separated on a 0.7% agarose gel. Lane 1, total nucleic acids; lane 2, total nucleic acids treated with DNase I; lane 3, total nucleic acids treated with RNase A; lane 4, total nucleic acids treated with RNase T1. Download FIG S1, PDF file, 0.3 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Sucrose gradient ultracentrifugation was conducted to purify virus particles to confirm that the dsRNA segments in the Malassezia strains were indeed viral elements. The separated nucleic acids and proteins in each fraction were analyzed. Two dsRNA fragments (∼5.0 and ∼1.7 kb) were clearly visible in fractions 1 to 6 following agarose gel electrophoresis (Fig. 2A). Moreover, the results of sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) showed that fractions 3 to 6 contained protein bands with an estimated molecular weight of ∼77 kDa (Fig. 2B). This molecular weight is similar to that of the known capsid protein of the S. cerevisiae mycovirus (9, 32, 33). Fractions 3 to 6 were subsequently evaluated by microscopy to visualize mycovirus particles in M. restricta. Transmission electron microscopy (TEM) images showed virus-like particles with an isometric shape and a diameter of 43 nm (Fig. 2C). These results support the hypothesis that the extrachromosomal dsRNA segments form the genome of mycovirus in M. restricta KCTC27540. We named the viral particle MrV40 (M. restricta KCTC 27540 mycovirus). The large and small dsRNA viral fragments were named MrV40L and MrV40S, respectively. In addition to evaluating images of the purified virus particles, we examined the morphology of the virus-infected and the virus-cured M. restricta KCTC 27540 strain pairs by TEM. The virus-cured M. restricta KCTC 27540 strain was prepared by sterilizing the virus-infected M. restricta KCTC 27540 through numerous serial passages (Text S1 and Fig. S2). In general, there was no morphological difference between the virus-infected and the virus-cured KCTC 27540 strains. However, we noticed that the size of vacuoles in the virus-infected strain was significantly larger than that in the virus-cured strain, suggesting that the virus influences vacuole size in M. restricta (Fig. 3).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Isolation of virus particles from Malassezia restricta KCTC 27540 and their TEM images. The collected fractions after sucrose gradient ultracentrifugation were analyzed for their nucleic acids and proteins on an agarose gel (A) and SDS-PAGE gel (B), respectively. Images of the viral particles (red arrowheads) were obtained using a transmission electron microscope (C).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Morphology of vacuoles in virus-infected and virus-cured Malassezia restricta cells. Ultrathin sections of virus-infected (A) and virus-cured (B) M. restricta KCTC 27540 strains were observed by TEM. Cytoplasmic vacuole formation in a cell, expressed as the percentage of the area occupied by vacuoles (C) and the number of vacuoles (D) per cell, was estimated (n = 100 cells for each strain). The vacuole area was estimated using ImageJ (https://imagej.nih.gov/ij/), and the plot was created with BioVinci software (BioTuring Inc., San Diego, CA, USA). Statistical analysis for differences between the strains was performed with unpaired t tests. **, P < 0.001.

TEXT S1

Supplemental materials and methods. Cell fixation, transmission electron microscopy (TEM), cDNA synthesis, RNA isolation and transcriptome analysis by RNA sequencing, quantitative RT-PCR (Q-RT-PCR), annexin V-FITC/propidium iodide (PI) staining and flow cytometric analysis, heterologous expression and purification of the capsid protein, heterologous expression and purification of the MrV40S satellite proteins, and preparation of fungal cell lysates. Download Text S1, PDF file, 0.1 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S2

Confirmation of the virus-cured Malassezia restricta strain. Total nucleic acids from the virus-infected and the virus-cured M. restricta KCTC 27540 and KCTC 27524 strains were extracted and treated with RNase T1. The dsRNA viral segments were not detected in the virus-cured M. restricta KCTC 27540 and KCTC 27524 strains (lane 2 and lane 4, respectively). Download FIG S2, PDF file, 0.3 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Determination of dsRNA sequence of MrV40L.The complete sequence of MrV40L was determined by a combination of the Illumina MiSeq technique and the Sanger sequencing method using purified viral dsRNA. The length of the complete assembled sequence of MrV40L was 4,606 bp, and two overlapping open reading frames (ORFs), designated ORF1 and ORF2, were identified (Fig. 4A). ORF1 corresponds to the region from nucleotides (nt) 28 to 2097 and encodes a polypeptide of 689 amino acids with a molecular weight of 77 kDa. ORF2 corresponds to the region from nt 1949 to 4538 and encodes a polypeptide of 862 amino acids with a molecular weight of 98 kDa. The results of BLAST analysis showed that the protein sequences of ORF1 and ORF2 were highly similar to that of the capsid protein (the Pfam families of LA-virus_coat, PF09220) and viral RNA-directed RNA polymerase (RDRP; the Pfam family of RDRP_4, PF02123) of Scheffersomyces segobiensis virus L belonging to the genus Totivirus (Totiviridae family) with 53% (YP_009507830.1) and 52% (YP_009507831.1) identities, respectively (34). Eight conserved motifs, which are commonly found in totiviruses, were found within ORF2, supporting the classification of MrV40 as a totivirus (Fig. 4B) (35). Moreover, phylogenetic analysis with known amino acid sequences of RDRPs from dsRNA mycoviruses demonstrated that the RDRP encoded by the MrV40L genome was clustered with totiviruses (Fig. 4C). Thus, MrV40L is a dsRNA viral genome encoding a capsid protein and an RDRP, and we propose that MrV40 belongs to the Totivirus genus.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Genome and phylogenetic analysis of MrV40L. (A) Comparison of the genomic organization of MrV40L with that of other totiviruses. (B) Eight conserved motifs within the viral RDRP of MrV40L. The red boxes indicate conserved motifs. The amino acid sequences of ScVL-A, ScVL-BC, SsVL, and MrV40L were aligned using Jalview. (C) Multiple alignment of 18 amino acid sequences of RDRPs from dsRNA viruses was analyzed using the neighbor-joining method with bootstrap test (1,000 replicates) using MEGA X (81, 83, 84). The evolutionary distances are in units of the number of amino acid substitutions per site. All ambiguous positions were removed for each sequence pair. CHV1, Cryphonectria hypovirus 1; CHV2, Cryphonectria hypovirus 2; AsV, Armigeres subalbatus virus; DTV, Drosophila melanogaster totivirus; OMRV, Omono River virus; ToV, Tianjin totivirus; GSTV, Golden shiner totivirus; ScVL-A, S. cerevisiae virus L-A; ScVL-BC, S. cerevisiae virus L-BC; SsVL, S. segobiensis virus L; GaRV-L1, Gremmeniella abietina virus L1; Hv190SV, Helminthosporium victoriae virus-190S; MoV1, Magnaporthe oryzae virus 1; HmV-17, Helicobasidium mompa totivirus 1-17; PcV, Penicillium chrysogenum virus; Hv145SV, Helminthosporium victoriae 145S virus; AbV1, Agaricus bisporus virus 1 (8, 33, 34, 85–97).

Next, we analyzed the genomic sequences of dsRNA viruses in other clinical M. restricta isolates to confirm that they were similar to those of MrV40L and thus belong to the Totivirus genus. Reverse transcription-PCR (RT-PCR) was performed using four primer sets corresponding to the sequence of the conserved regions of the capsid protein and RDRP in the MrV40L genome (Fig. 5A). The results showed that six of 11 viruses generated RT-PCR products, suggesting the presence of viral genomes similar to MrV40Ls, while the remaining samples had dissimilar viral genomes (Fig. 5B).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Phylogenetic classification of viruses found in Malassezia restricta. (A and B) The cDNAs of dsRNAs from M. restricta strains were amplified using primers for genes encoding CP and RDRP (ORF1 and ORF2, respectively) homologs in the viral genome. The following primers were used: lane 1, MrV40L_CP_F1 and MrV40L_CP_R1; lane 2, MrV40L_CP_F2 and MrV40L_CP_R2; lane 3, MrV40L_RDRP_F1 and MrV40L_RDRP_R1; lane 4, MrV40L_RDRP_F2 and MrV40L_RDRP_R2 (see Table S2 in the supplemental material). (C) MrV-Ls were clustered into three clades (clades I, II, and III). Multiple alignment of nucleotide sequences of combined parts of gag and pol from 11 MrV-Ls was analyzed by the neighbor-joining method with bootstrap test (1,000 replicates) in MEGA X (81, 83, 84). The evolutionary distances are in units of the number of base substitutions per site.

Furthermore, we performed phylogenetic classification of the same viruses found in other clinical M. restricta isolates by multilocus sequence typing of the 1,075-bp region, 638 bp of ORF1 and 437 bp of ORF2, corresponding to gag and pol, respectively, within the viral genomes. The results revealed that the viruses were classified into three clades: clade I (MrV12L, MrV16L, MrV40L, MrV79L, MrV80L, and MrV82L), clade II (MrV18L, MrV43L, MrV83L, and MrV50L), and clade III (MrV24L) (Fig. 5C). Additionally, the sequences of the 1,075-bp region of MrV40L, MrV79L, and MrV82L in clade I were 100% identical, suggesting that they originated from the same lineage.

Determination of dsRNA sequence of MrV40S.To determine the sequence of MrV40S, the dsRNA segments of MrV40S were extracted from agarose gels and then subjected to cDNA cloning and sequencing (see Materials and Methods). Using the partial sequences obtained from cDNA clones and the sequences obtained from repeated 5′ rapid amplification of cDNA ends (RACE), we successfully determined the complete sequence of MrV40S. The sequence of MrV40S was 1,355 bp, and a single ORF was identified in the 3′ region (from nt 773 to 1097), which encoded a polypeptide of 124 amino acids with a molecular weight of 15.6 kDa. The ORF was designated ORF3. Although we obtained the complete sequence of MrV40S ORF3, no homologous protein sequence was identified by BLAST analysis using all currently available databases.

Mycoviruses belonging to Totivirus, particularly the S. cerevisiae L-A virus, often possess a satellite dsRNA segment known as M dsRNA, which is responsible for producing a killer toxin that excludes neighboring yeast cells. Because MrV40L resembles M dsRNA, we predicted that MrV40S produces a protein that inhibits other Malassezia strains and/or other fungal and/or bacterial cells residing on the skin. To test the toxin-like activity of the protein produced from MrV40S, ORF3 was cloned and the protein was heterologously expressed in Escherichia coli and purified (Fig. 6A). The activity of the purified protein was evaluated against several pathogenic fungi and bacteria including M. restricta, C. albicans, Cryptococcus neoformans, E. coli, and Staphylococcus aureus. Unexpectedly, the purified protein displayed no growth-inhibitory effect on the microbial cells tested (data not shown). Based on these results, we concluded that the novel protein produced from ORF3 likely has no toxin-like activity against the microorganisms tested, and further functional studies are required to characterize its function.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Characterization and similarity analysis of MrV40S. (A) The protein encoded by ORF3 was heterologously expressed in E. coli and purified using a His tag column. (B) RT-PCR results. The primers specific to MrV40S were used to amplify satellite dsRNA from other M. restricta strains containing small dsRNA segments. The following primers were used: lane 1, MrV40S_F5 and V40S_SP1; lane 2, MrV40S_ORF_F1 and MrV40S_R1 (see Table S2 in the supplemental material).

M dsRNA in S. cerevisiae produces several variants of killer toxins, namely, M1, M2, M3, and Mlus, and they showed limited homology in their protein sequences (4–6). Based on this information, we compared the sequences of satellite dsRNA genomes between MrV40S and those in other M. restricta strains containing the small dsRNA segments. RT-PCR was conducted using a series of primers specific to MrV40S. The results suggested that among the nine strains tested, only the satellite dsRNA in M. restricta KCTC 27582 possessed a similar dsRNA segment. These data indicate that the satellite dsRNA genome sequences are highly variable (Fig. 6B).

MrV40 altered transcriptome profiles in the host M. restricta strain.To understand the influence of MrV40 on the host physiology, we compared the transcriptomes of virus-infected and virus-cured M. restricta KCTC 27540 strain pairs by RNA sequencing. Transcriptome analysis showed that 258 and 809 genes were up- and downregulated by more than 2-fold, respectively, in the virus-infected strain compared to those in the virus-cured strain (Table 1 and Table S1). These results indicate that the presence of the virus may impact a number of physiological processes in M. restricta.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Number of differentially expressed genes

TABLE S1

Expression of genes calculated from transcriptome sequencing (RNA-seq) data. Download Table S1, PDF file, 1.8 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

In S. cerevisiae, numerous genes are known to be involved in maintaining its dsRNA virus (36) and are categorized into two groups: MAK (maintenance of killer) and SKI (superkiller). At least 25 MAK genes have been reported (37); among them, MAK3, MAK10, and PET18 were shown to be required for the maintenance of both S. cerevisiae L-A virus and its M satellite, whereas all other MAK genes are responsible only for the M satellite (38–43). We could not find homologs of MAK10 or PET18 in the genome of M. restricta, suggesting that a different or modified virus maintenance mechanism may be present in the fungus compared to S. cerevisiae (Table 2). The results also showed that among the MAK homologs, MAK1, MAK5, MAK11, MAK16, MAK21, MAK7, and MAK8 were downregulated by more than 2-fold. In S. cerevisiae, these genes, except for MAK1, were involved in 60S ribosomal subunit biosynthesis and M satellite maintenance. Thus, we predicted that MrV40 reduced ribosomal biogenesis within the host. Additionally, M. restricta contained a series of SKI homologs. However, no SKI homologs were differentially expressed in the virus-infected M. restricta KCTC 27540 compared to those in the virus-cured strain, indicating that the function of SKI genes was not critical for maintaining the virus. In addition to the genes involved in virus maintenance, we observed upregulation of numerous genes required for the tricarboxylic acid (TCA) cycle and the electron transport chain, including MRET_1104 [NADH dehydrogenase (ubiquinone) Fe-S protein 7], MRET_1378 [succinate dehydrogenase (ubiquinone) iron-sulfur subunit], MRET_1953 [NADH dehydrogenase (ubiquinone) Fe-S protein 1], MRET_2042 (fumarate hydratase, class II), MRET_2097 [succinate dehydrogenase (ubiquinone)], MRET_2956 (2-oxoglutarate dehydrogenase E1 component), MRET_3173 (dihydrolipoamide dehydrogenase), and MRET_4117 (aconitate hydratase) (Table 3). These results suggest that virus maintenance and propagation may require higher energy production in the host cell.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Differential expression of genes involved in maintaining dsRNA virus

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 3

Differential expression of genes involved in TCA cycle and electron transport chain

The overall dysregulation of primary metabolism may disturb the normal cell physiology in the M. restricta strain infected with MrV40. Indeed, we observed upregulation of genes involved in programmed cell death in the fungal host cells. For example, the expression of MRET_3200 (p38 mitogen-activated kinase [MAP] kinase), MRET_1134 (programmed cell death 6-interacting protein), and MRET_2499 (autophagy-related protein, ATG101), which are associated with programmed cell death, was upregulated by 4.43-, 3.14-, and 2.83-fold in the virus-infected fungal cells, respectively. Moreover, MRET_0131 (Fas-associated factor), which is involved in Fas-induced apoptosis, was found to be strongly upregulated (8.80-fold) in the presence of the virus within the fungal host (44). It is well known that programmed cell death is triggered during virus infection (45, 46), and our findings agree with this observation. To confirm that there might be an increase in programmed cell death, we stained cells with annexin V and propidium iodide, which is commonly used to quantitate the apoptosis-like phenomenon (47, 48). Using flow cytometric analysis, we found a higher percentage of positive population in the virus-infected KCTC 27540 strain than in the virus-cured strain, implying that programmed cell death might be increased in the virus-infected M. restricta strain (Fig. S4). We also tested whether increased programmed cell death influences other physiological characteristics of the virus-infected cells such as susceptibility to an antifungal drug. However, antifungal susceptibility of virus-infected KCTC 27540 was similar to that of the virus-cured strain (data not shown), suggesting that the viral infection may have a limited effect on the physiological characteristics of the host cells.

FIG S3

Validation of the differential expression in transcriptome analysis. Differential expression of selected genes was confirmed by Q-RT-PCR. MRET_0131 (Fas-associated factor); MRET_2499 (autophagy-related protein, ATG101); MRET_3200 (p38 MAP kinase); MRET_1953 (NADH dehydrogenase Fe-S protein 1); MRET_2956 (2-oxoglutarated dehydrogenase E1 component); MRET_0230 (large subunit ribosomal protein L8e); MRET_1468 (large subunit ribosomal protein L22e). Download FIG S3, PDF file, 0.2 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG S4

Evaluation of apoptosis in M. restricta strains. The virus-infected and the virus-cured M. restricta KCTC 27540 cells were stained with annexin V and propidium iodide (PI) to evaluate apoptosis. Annexin V binds to phosphatidylserine (PS), which translocates from the inner leaflet of the plasma membrane to the outer leaflet upon apoptosis while the plasma membrane integrity is still maintained (1). PI is permeable when cell membrane was ruptured by damage or upon cell death (2). Overall apoptosis was quantified by annexin V (Ann), and early and late apoptosis was quantified by annexin V (Ann)/PI and by PI, respectively. Statistical analysis for differences between the strains was performed with unpaired t tests. *, P < 0.01. Download FIG S4, PDF file, 0.10 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

MrV40 induced the TLR3-mediated immune system.Since it was first reported that fungal viral dsRNA induces cytokine production in rabbits (49–51), several studies have demonstrated that TLR3 plays a central role in viral dsRNA recognition and production of inflammatory cytokines in innate immune cells (52, 53). Additionally, a recent study showed that S. cerevisiae viral dsRNA stimulates the immune system through TLR3 in a human embryonic kidney cell line (54). We therefore investigated whether MrV40 itself or the virus-containing Malassezia cells alter the expression patterns of TLRs and cytokine production in mammalian cells. To this end, bone marrow-derived dendritic cells (BMDCs) obtained from C57BL/6 mice, which have been used as a model to study the interactions between fungal cells including Malassezia and the innate immune system in a mammalian host, were used in our study (55, 56).

Two independent strain pairs, virus-infected and virus-cured M. restricta KCTC 27540 strains as well as virus-infected and virus-cured M. restricta KCTC 27524 strains, were used for the study. Live fungal cells, fungal cell lysates, purified capsid protein, purified satellite protein (produced from MrV40S ORF3), and purified MrV40 dsRNA were prepared and coincubated with BMDCs, and the expression of TLRs and cytokines was analyzed. We found that TLR3 expression was significantly upregulated by live fungal cells and fungal cell lysates of the virus-infected M. restricta strains, while live cells and their lysates of the virus-cured strains did not enhance the expression of TLR3. Considering that the expression of other TLRs was unchanged upon addition of live cells and cell lysates, we concluded that host response upon treatment with virus-infected M. restricta cells was mainly mediated by TLR3 in murine BMDCs. Moreover, upregulation of TLR3 expression was observed in BMDCs treated with purified MrV40 dsRNA but not with purified capsid and satellite proteins, suggesting that MrV40 dsRNA is the causative agent that enhances the expression of TLR3 (Fig. 7A). Moreover, our data suggested that purified capsid and satellite proteins from MrV40 did not influence the expression of TLRs.

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

Evaluation of TLR and cytokine levels upon treatment with the viral elements. TLR expression (A) and cytokine production (B) in BMDCs coincubated with live fungal cells, total fungal cell lysates, purified MrV40 capsid protein, purified satellite protein, and purified MrV40 dsRNA. PBS and poly(I·C) served as negative and positive controls, respectively. WT, wild-type mice (C57BL/6); tlr3−/−, TLR3 knockout mice (B6;129S1-Tlr3tm1Flv/J, 005217); FC, live fungal cells; LY, cell lysate; 40+, virus-infected KCTC 27540; 40-, virus-cured KCTC 27540; 24+, virus-infected KCTC 27524; 24-, virus-cured KCTC 27524; CP, purified capsid protein; p15, purified MrV40S satellite protein; dsRNA, purified MrV40 dsRNA. Values represent the average from quadruplicates with standard deviations. Statistical analysis for differences between PBS and each sample was performed with unpaired t tests. *, P < 0.001; **, P < 0.0001.

The expression of several inflammatory cytokines including tumor necrosis factor alpha (TNF-α), interleukin-6 (IL-6), IL-10, alpha interferon (IFN-α), and IFN-γ, in response to live fungal cells, fungal cell lysates, purified capsid protein, purified satellite protein, and purified MrV40 dsRNA, was evaluated. Particularly, IFN-α and IFN-γ were included because of their involvement in the antiviral response in mammals (57). The results of the cytokine analysis indicated that virus-infected live M. restricta cells and their lysates induced the production of TNF-α, IL-6, IL-10, IFN-α, and IFN-γ, while no change in cytokine expression was observed in BMDCs treated with virus-cured live M. restricta cells and their cell lysates. Similarly, capsid protein and MrV40 dsRNA also enhanced the expression of all cytokines tested (Fig. 7B). Cytokine profiles were also measured in homozygous TLR3 knockout mice, as our observations suggested a TLR3-mediated host response in BMDCs treated with virus-infected fungal cells. As shown in Fig. 7B, cytokine levels in BMDCs isolated from TLR3 knockout mice were not significantly altered except in those treated with the capsid and satellite proteins. These data suggest that MrV40 increased the production of the cytokines TNF-α, IL-6, IL-10, IFN-α, and IFN-γ in a TLR3-dependent manner in BMDCs. However, as suggested above, the cytokine production induced by the capsid protein seemed to be TLR3 independent. Moreover, the purified satellite protein also caused a TLR3-independent increase in cytokine production except of TNF-α. These results suggest that the capsid and the satellite proteins also contributed to upregulation of cytokine expression in a TLR3-independent manner. Overall, our data suggest that dsRNA in MrV40 triggers an increase in the production of cytokines involved in inflammation and that TLR3 plays a central role in the host response.

DISCUSSION

In the current study, we detected dsRNA virus in several clinical isolates of Malassezia species. Among them, MrV40 identified in M. restricta KCTC 27540 was selected, and its genome structure and effects on host gene expression and on the mammalian immune response were evaluated. Our data showed that MrV40 consists of two RNA segments, which we named MrV40L and MrV40S. The results of the genome sequence analysis suggested that these segments were 4,606 and 1,355 bp in length, respectively, and belong to the genus Totivirus. Typically, the genomes of the viruses belonging to the genus Totivirus consist of nonsegmented dsRNA with sizes between 4.6 and 7.0 kb and contain two ORFs, gag and pol. Studies have specifically examined the genome structure of Totivirus because of the overlapping nature of the two ORFs, where a −1 frameshift occurs, resulting in translation of the fusion protein (2). The overlapping ORFs and frameshift were frequently observed in a compact viral genome to translate proteins and were found in several dsRNA and ssRNA viruses; these ORFs allow ribosomes to translate CP and RDRP continuously with a missing CP termination codon (11, 58, 59). The mechanism of frameshifting in the viral genome is based on the pseudoknot structure of the mRNA for efficient slipping via a slippery site (60). For example, in the genome of the S. cerevisiae L-A virus, the 5′ ORF (gag) encodes a 76-kDa CP and the 3′ ORF (pol) encodes an RDRP, which is expressed as a 180-kDa CP-RDRP fusion protein generated by a −1 ribosomal frameshift (8, 11). In our study, MrV40L (the major dsRNA segment in MrV40) contained two overlapping ORFs, ORF1 and ORF2. We identified a putative slippery site heptamer, 5′-GGGTTTT-3′, at the region from nt 1968 to 1974 for the −1 ribosomal frameshift, which may be associated with production of the fusion protein of 170 kDa in MrV40L. A previous study suggested that in the S. cerevisiae L-A virus, the rate of ribosomal frameshifts was approximately 1.8%, giving 120 CP and 2 CP-RDRP fusion protein molecules per virus particle (11). Considering the low efficiency of producing the fusion proteins by ribosomal frameshifting, we expected to observe a significantly lower translation rate of the fusion protein; indeed, the putative 170-kDa band was not detected by SDS-PAGE. In addition to MrV40L, we determined the sequence of the satellite dsRNA segment MrV40S and found that it consists of 1,355 nt containing a single ORF, ORF3, producing a novel 15.6-kDa protein. As observed for other totiviruses, the possible toxin-like activity of the protein was investigated in our study, but no growth-inhibitory activity against several bacteria and fungi was observed.

It has been estimated that 30 to 80% of fungal species in nature are infected with viruses, and a fungal host normally shows no specific symptoms upon infection (14). However, several genes were shown to be required for maintaining and propagating viruses in the host fungal cells. In S. cerevisiae, numerous chromosomal genes are known to be involved in viral propagation and expression of the viral killer toxin (36). Furthermore, several genes are known to be responsible for maintaining the L-A virus and M dsRNA in S. cerevisiae. The MAK genes are required for the propagation and maintenance of the L-A virus and M dsRNA in S. cerevisiae (61). Among the MAK genes, MAK3, MAK10, and PET18 are required for the maintenance of both L-A virus and M dsRNA, whereas all other MAK genes are responsible only for M dsRNA (38–43). Particularly, MAK3 encodes an N-acetyltransferase and is required for N-acetylation of the coat protein (38, 62). A previous study showed that the coat proteins without acetylation failed to self-assemble, resulting in the loss of all dsRNA viruses (39). MAK10 and PET18 (MAK30+MAK31) encode a noncatalytic subunit of N-terminal acetyltransferase and a protein of unknown function, respectively. Mutant strains lacking each gene contained unstable viral particles, indicating that the genes are involved in the structural stability of L-A virus and M dsRNA (40). MAK1 (MAK17, TOP1) encodes DNA topoisomerase I, and other MAK genes including MAK2, MAK5, MAK11, MAK16, MAK21, MAK7 (RPL4A), MAK8 (TCM1), and MAK18 (RPL41B) are related to 60S ribosomal subunit assembly (43, 63, 64). All mutant strains lacking the above genes showed decreased levels of free 60S ribosomal subunits and the inability to maintain M dsRNA, suggesting that stable propagation of the satellite dsRNA depends on 60S ribosome synthesis (42, 43). In addition to the MAK genes, the SKI gene family has been shown to be involved in the maintenance and propagation of virus in S. cerevisiae. SKI1 (XRN1) encoding a 5′-to-3′ exonuclease is involved in the degradation of uncapped mRNA including viral mRNA, and SKI2, SKI3, SKI6, SKI7, and SKI8 block the translation of viral mRNAs (65–67).

In the current study, homologs of most MAK and SKI genes were identified in M. restricta. The results of the transcriptome analysis suggested that most MAK genes were downregulated, which may in turn reduce ribosome synthesis in the M. restricta strain containing MrV40. In contrast, no SKI homolog showed significantly altered transcript levels in the M. restricta strain harboring MrV40. Moreover, we found enhanced expression of genes involved in energy metabolism and programmed cell death in the M. restricta strain containing the virus. S. cerevisiae also displayed relatively small changes in fungal host gene expression upon virus infection, possibly because of coadaptation of the virus within the fungal host (68). Maintenance and propagation of virus within the fungal host may be involved in the posttranscriptional mechanism and may contribute to the minimal changes in host gene expression. Notably, the possibility that an RNA silencing pathway in M. restricta cells influences virus maintenance was excluded because of the absence of homologous genes required for the pathway in the genome of the fungus.

In addition to transcriptome analysis, we directly investigated whether the virus influences the cellular morphology of M. restricta and the structures of its intracellular organelles by TEM. We observed significantly larger vacuoles in virus-infected M. restricta cells. An increased vacuole size upon virus infection has been reported previously. The phytopathogenic fungus Botrytis porri infected with dsRNA virus 1 showed the formation of abundant vacuoles (69), and turnip mosaic virus induced the formation of a large central vacuole in Nicotiana benthamiana plant cells (70). Particularly, in N. benthamiana, turnip mosaic virus particles were shown to accumulate in vacuoles and be protected by the organelle membranes against the harsh host environment, particularly during xylem vessel differentiation. Furthermore, virus accumulation in vacuoles in plant-to-plant virus transmission has been suggested (70). Although the expression of genes involved in vacuole biogenesis was not altered in transcriptome analysis in virus-infected M. restricta cells, there may be a relationship between vacuole functions and MrV40 in M. restricta. Therefore, additional physiological studies are needed.

TLR3 is well conserved in most vertebrates, localizes on the endosomal transmembrane, and plays a role in immune and nonimmune cells. It has been suggested that TLR3 senses viral dsRNA taken up through endocytosis and contributes to defending the host against viral infection by regulating the expression of a range of cytokines (52, 71). A previous study demonstrated a significant decrease in cytokines (IFN-α/β, IFN-γ, and IL-12p40) and increase in viral PFU in the spleen of Tlr3−/− mice infected with mouse cytomegalovirus (72). We found that TLR3 expression was significantly increased by the live cells of two independent virus-infected M. restricta strains and their cell lysates, while no change was found by the live cells and cell lysates of the virus-cured strains. The causative agent was very likely to be dsRNA of MrV40 in the virus-infected M. restricta strains since purified dsRNA from the MrV40 particles enhanced TLR3 expression. The virus-infected live M. restricta cells, their cell lysates, and purified dsRNA from the MrV40 particles also showed increased production of cytokines involved in inflammation. However, in Tlr3−/− BMDCs, the production of all cytokines in the cells treated with the virus-infected live M. restricta cells, their cell lysates, and purified dsRNA from the MrV40 particles was attenuated, suggesting that TLR3 plays a central role in the host response against dsRNA of MrV40 and mediates the production of the cytokines.

Similarly, the Leishmania parasite L. guyanensis infected with the dsRNA virus belonging to Totivirus induced a TLR3-mediated inflammatory immune response within the vertebrate host, indicating that the Leishmania RNA virus (LRV1) functions as an innate immunogen. Moreover, in Leishmania, LRV1 was suggested to stimulate the inflammatory immune response and increase the severity and persistence of the parasite (73). However, whether MrV40 contributes to hyper- or hypopathogenicity of M. restricta remains unclear. Therefore, further studies to characterize the potential function of MrV40 as an innate immunogen are needed. In addition, since our data suggest contribution of the MrV40 capsid and the satellite proteins to host cytokine expression in a TLR3-independent manner, additional studies are needed to define a role of the capsid and the satellite proteins in interactions between these viral components and the mammalian host cells.

Overall, our study demonstrated the existence of a dsRNA virus within M. restricta. We also determined the sequences and structure of the viral genome along with the independent RNA segment. Moreover, we identified viruses not only from different strains of M. restricta but also from other Malassezia species, although variations were observed in the viral genomes. Evidence that the viral nucleic acid from MrV40 induces a TLR3-mediated inflammatory immune response was obtained, suggesting that a viral element contributes to the pathogenicity of Malassezia.

MATERIALS AND METHODS

Strains and growth conditions.M. restricta KCTC 27512, KCTC 27516, KCTC 27518, KCTC 27524, KCTC 27540, KCTC 27543, KCTC 27550, KCTC 27879, KCTC 27880, KCTC 27882, and KCTC 27883, M. globosa CBS 7966, M. pachydermatis KCTC 27587, and M. sympodialis KCTC 27817 were obtained as previously described and cultured on Leeming and Notman agar (LNA) medium at 34°C for 3 days (29, 74–76). Among these strains, M. restricta KCTC 27540 was used to identify dsRNA viruses. E. coli DH10B and BL21 were grown in Luria-Bertani broth at 37°C.

Purification of virus particles.The virus particles were purified as previously described with some modifications (69). Briefly, approximately 10 g of M. restricta cells was harvested and suspended in 10 ml of extraction buffer (0.1 M sodium phosphate buffer containing 3% [vol/vol] Triton X-100, pH 7.0). The suspended cells were vortexed with glass beads and centrifuged at 10,000 × g at 4°C for 20 min to remove cell debris. The supernatant containing crude extracts was ultracentrifuged at 119,000 × g at 4°C for 2 h to precipitate all particles. The resulting pellet was suspended in 1 ml of 0.1 M sodium phosphate buffer (pH 7.0), and the suspension was centrifuged at 16,000 × g at 4°C for 30 min to remove large particles. The supernatant was overlaid on sucrose solution with a gradient concentration ranging from 10% to 40% (wt/vol) and centrifuged at 70,000 × g at 4°C for 2 h. Fractions (700 μl/each fraction) were collected from the top of the ultracentrifuged sucrose gradient solution. Each fraction was analyzed by 1% agarose gel electrophoresis and 8% SDS-PAGE to detect dsRNA segments and protein, respectively. Methods for cell fixation and TEM are described in Text S1 in the supplemental material.

Determination of MrV40L and MrV40S sequences.The reference genome sequence was required to obtain the complete sequence of MrV40L by Illumina MiSeq. We first obtained the partial fragment of MrV40L and determined its nucleotide sequence. Briefly, 1 μg of cDNA obtained as described above was amplified using 1 μM tagged oligo_1 primer to bind the tag and 1.25 units of Taq DNA polymerase (Bioneer, Daejeon, South Korea). The amplicons were cloned into the pGEM-T vector (Promega, Madison, WI, USA), transformed into E. coli DH10B, and sequenced with the universal M13 primers. The resulting partial sequences of MrV40L showed high similarity with the genome of Scheffersomyces segobiensis virus L isolate NRRL Y-11571, which was used to assemble the reference genome in our study. To sequence the entire genome of the virus from M. restricta KCTC 27540, dsRNA was purified using cellulose A (Advantec, Taipei, Taiwan) from total cell extract as described previously (77). After purification, purified dsRNA was treated with Ambion DNase I (Thermo Fisher Scientific, Waltham, MA, USA) to remove remaining genomic DNA. The sequencing library was constructed using the TruSeq stranded total RNA sample prep kit (Illumina, San Diego, CA, USA) following the manufacturer’s instructions and excluding poly(A) selection. The resulting library was sequenced on an Illumina MiSeq instrument according to the manufacturer’s instructions. The generated raw sequences, 250-bp paired-end reads, were assembled by CLC Genomics Workbench v7.5 (Qiagen, Hilden, Germany), and the resulting contigs were identified by BLASTN search of the NCBI nucleotide database. Based on the reference genome, two contigs of MrV40L were assembled. The gap between the contigs was sequenced by gap-filling RT-PCR using primers Linkage1_V40 and Linkage2_V40) (see Table S2), and both termini were sequenced by RACE PCR using a 5′/3′ RACE kit, second generation (Roche, Basel, Switzerland), according to the 5′ RACE protocol described by the manufacturer.

TABLE S2

Primers used for sequencing and heterologous expression of viral proteins. Download Table S2, PDF file, 0.02 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

To determine the MrV40S genome, the RACE PCR method was used because we failed to find any contig of MrV40S. Like MrV40L, the partial sequence of cDNA from MrV40S was amplified using the tagged oligo_1 primer. The amplicons were cloned into the pGEM-T vector (Promega) and sequenced. Using the partial sequence obtained, 5′ RACE PCR was performed and repeated until no additional sequence was found.

In silico analysis.To identify the ORFs in MrV40L and MrV40S, the nucleotide sequences of the genomes were analyzed by ORFfinder (https://www.ncbi.nlm.nih.gov/orffinder/). To predict ribosomal frameshift signals, a putative slippery site was identified using FSFinder2 (78). The amino acid sequences of CP and RDRP of MrV40L were analyzed using the Pfam database (http://pfam.xfam.org/) and BLASTP (https://blast.ncbi.nlm.nih.gov/). For sequence alignment, Jalview2.0 was used (79). The phylogenetic tree showing the relationship of MrV40L with RDRP of other dsRNA fungal viruses was constructed using MEGA X, and phylogenetic reconstruction analysis was performed by the neighbor-joining method with 1,000 bootstrap replications (80, 81).

Mice and cell culture.Wild-type C57BL/6 mice were purchased from Orient Bio Co. (Gyeonggi-do, South Korea). TLR3−/− (B6;129S1-Tlr3tm1Flv/J, 005217) mice were obtained from Jackson Laboratory (Bar Harbor, ME, USA). All animal experimental procedures were reviewed and approved by the Institutional Animal Care and Use Committee of Hanyang University (protocol 2018-0085) and performed in accordance with Korean Food and Drug Administration guidelines. All animals were maintained in a specific-pathogen-free environment. Primary BMDCs were isolated from C57BL/6 mice and cultured in Dulbecco’s modified Eagle’s medium for 3 to 5 days in the presence of 20 ng/ml recombinant granulocyte-macrophage colony-stimulating factor (R&D Systems, Minneapolis, MN, USA) as described previously (56). Cell cultures were stained to detect dendritic cells with CD11c-fluorescein isothiocyanate (FITC) (eBiosciences, San Diego, CA, USA).

Cytokine measurement.Mouse cytokines in the culture supernatants were measured with a BD OptEIA enzyme-linked immunosorbent assay (ELISA) kit (BD Biosciences, Franklin Lakes, NJ, USA) as described previously (82). All assays were performed according to the manufacturer’s instructions. Two independent strain pairs, virus-infected and virus-cured M. restricta KCTC 27540 and virus-infected and virus-cured M. restricta KCTC 27540, were used. To treat BMDCs with live fungal cells, a 20:1 ratio (2.0 × 106 and 1.0 × 105 cells/ml of M. restricta and BMDCs, respectively) was used. A final concentration of 100 μg/ml (total protein/volume) of fungal cell lysates was used to treat BMDCs. A final concentration of 1 μM purified MrV40 capsid and satellite proteins and 500 ng/ml of purified MrV40 dsRNA from virus-infected M. restricta KCTC 27540 strain were incubated with BMDCs. Phosphate-buffered saline (PBS) and poly(I)-poly(C) [poly(I·C), 25 μg/ml] served as negative and positive controls, respectively. Cytokines were measured 6 h after treatment. All experiments were repeated four times.

Other methods.For descriptions of cDNA synthesis, RNA isolation, transcriptome analysis by RNA sequencing, quantitative real-time PCR, heterologous expression and purification of MrV40 proteins, cell fixation and transmission electron microscopy, and flow cytometric analysis, see Text S1 in the supplemental material.

Data availability.The complete genome sequences of MrV40L and MrV40S were deposited into GenBank under accession numbers MN603497 and MN603498, respectively. The transcriptome data were deposited in Gene Expression Omnibus under the accession number GSE138985.

FIG S5

Purification of the MrV40 capsid protein. The capsid protein encoded by ORF1 of MrV40 was heterologously expressed in E. coli, purified using a His tag column, and used to analyze the expression of TLRs and cytokines in BMDCs. Download FIG S5, PDF file, 0.02 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S3

Primers used for Q-RT-PCR. Download Table S3, PDF file, 0.09 MB.
Copyright © 2020 Park et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

We thank James W. Kronstad for critical reading of the manuscript and Taegun Seo for helpful discussions. There is a companion paper from the Joe Heitman group at Duke University, and we thank Joe Heitman and Giuseppe Ianiri for communication.

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT), 2019R1F1A1061930 (W.H.J.), 2019R1A4A1024764 (W.H.J.), 2019R1I1A2A01064237 (C.-S.Y.), and Agency for Science, Technology and Research IAF-PP (HBMS) grant Asian Skin Microbiome Project H18/01a0/016 (T.L.D.).

FOOTNOTES

    • Received 9 June 2020
    • Accepted 30 July 2020
    • Published 1 September 2020
  • Copyright © 2020 Park et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Hollings M
    . 1962. Viruses associated with a die-back disease of cultivated mushroom. Nature 196:962–965. doi:10.1038/196962a0.
    OpenUrlCrossRefWeb of Science
  2. 2.↵
    1. Ghabrial SA,
    2. Caston JR,
    3. Jiang D,
    4. Nibert ML,
    5. Suzuki N
    . 2015. 50-plus years of fungal viruses. Virology 479-480:356–368. doi:10.1016/j.virol.2015.02.034.
    OpenUrlCrossRef
  3. 3.↵
    1. Bussey H
    . 1991. K1 killer toxin, a pore-forming protein from yeast. Mol Microbiol 5:2339–2343. doi:10.1111/j.1365-2958.1991.tb02079.x.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    1. Hannig EM,
    2. Leibowitz MJ
    . 1985. Structure and expression of the M2 genomic segment of a type 2 killer virus of yeast. Nucleic Acids Res 13:4379–4400. doi:10.1093/nar/13.12.4379.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Schmitt MJ,
    2. Tipper DJ
    . 1990. K28, a unique double-stranded RNA killer virus of Saccharomyces cerevisiae. Mol Cell Biol 10:4807–4815. doi:10.1128/mcb.10.9.4807.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Rodríguez-Cousiño N,
    2. Maqueda M,
    3. Ambrona J,
    4. Zamora E,
    5. Esteban R,
    6. Ramírez M
    . 2011. A new wine Saccharomyces cerevisiae killer toxin (Klus), encoded by a double-stranded RNA virus, with broad antifungal activity is evolutionarily related to a chromosomal host gene. Appl Environ Microbiol 77:1822–1832. doi:10.1128/AEM.02501-10.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Cheng RH,
    2. Caston JR,
    3. Wang G-J,
    4. Gu F,
    5. Smith TJ,
    6. Baker TS,
    7. Bozarth RF,
    8. Trus BL,
    9. Cheng N,
    10. Wickner RB
    . 1994. Fungal virus capsids, cytoplasmic compartments for the replication of double-stranded RNA, formed as icosahedral shells of asymmetric Gag dimers. J Mol Biol 244:255–258. doi:10.1006/jmbi.1994.1726.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Icho T,
    2. Wickner RB
    . 1989. The double-stranded RNA genome of yeast virus L-A encodes its own putative RNA polymerase by fusing two open reading frames. J Biol Chem 264:6716–6723.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Hopper JE,
    2. Bostian K,
    3. Rowe L,
    4. Tipper D
    . 1977. Translation of the L-species dsRNA genome of the killer-associated virus-like particles of Saccharomyces cerevisiae. J Biol Chem 252:9010–9017.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Fujimura T,
    2. Wickner RB
    . 1988. Gene overlap results in a viral protein having an RNA binding domain and a major coat protein domain. Cell 55:663–671. doi:10.1016/0092-8674(88)90225-5.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Dinman JD,
    2. Icho T,
    3. Wickner RB
    . 1991. A -1 ribosomal frameshift in a double-stranded RNA virus of yeast forms a gag-pol fusion protein. Proc Natl Acad Sci U S A 88:174–178. doi:10.1073/pnas.88.1.174.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Van Alfen NK,
    2. Jaynes RA,
    3. Anagnostakis SL,
    4. Day PR
    . 1975. Chestnut blight: biological control by transmissible hypovirulence in Endothia parasitica. Science 189:890–891. doi:10.1126/science.189.4206.890.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Choi GH,
    2. Nuss DL
    . 1992. Hypovirulence of chestnut blight fungus conferred by an infectious viral cDNA. Science 257:800–803. doi:10.1126/science.1496400.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Ghabrial SA,
    2. Suzuki N
    . 2009. Viruses of plant pathogenic fungi. Annu Rev Phytopathol 47:353–384. doi:10.1146/annurev-phyto-080508-081932.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Segers GC,
    2. Zhang X,
    3. Deng F,
    4. Sun Q,
    5. Nuss DL
    . 2007. Evidence that RNA silencing functions as an antiviral defense mechanism in fungi. Proc Natl Acad Sci U S A 104:12902–12906. doi:10.1073/pnas.0702500104.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Zhang X,
    2. Segers GC,
    3. Sun Q,
    4. Deng F,
    5. Nuss DL
    . 2008. Characterization of hypovirus-derived small RNAs generated in the chestnut blight fungus by an inducible DCL-2-dependent pathway. J Virol 82:2613–2619. doi:10.1128/JVI.02324-07.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Hammond TM,
    2. Andrewski MD,
    3. Roossinck MJ,
    4. Keller NP
    . 2008. Aspergillus mycoviruses are targets and suppressors of RNA silencing. Eukaryot Cell 7:350–357. doi:10.1128/EC.00356-07.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Yaegashi H,
    2. Shimizu T,
    3. Ito T,
    4. Kanematsu S
    . 2016. Differential inductions of RNA Silencing among encapsidated double-stranded RNA mycoviruses in the white root rot fungus Rosellinia necatrix. J Virol 90:5677–5692. doi:10.1128/JVI.02951-15.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Yu J,
    2. Lee KM,
    3. Cho WK,
    4. Park JY,
    5. Kim KH
    . 2018. Differential contribution of RNA interference components in response to distinct Fusarium graminearum virus infections. J Virol 92:e01756-17. doi:10.1128/JVI.01756-17.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Findley K,
    2. Oh J,
    3. Yang J,
    4. Conlan S,
    5. Deming C,
    6. Meyer JA,
    7. Schoenfeld D,
    8. Nomicos E,
    9. Park M, NIH Intramural Sequencing Center Comparative Sequencing Program,
    10. Kong HH,
    11. Segre JA
    . 2013. Topographic diversity of fungal and bacterial communities in human skin. Nature 498:367–370. doi:10.1038/nature12171.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    1. Clavaud C,
    2. Jourdain R,
    3. Bar-Hen A,
    4. Tichit M,
    5. Bouchier C,
    6. Pouradier F,
    7. El Rawadi C,
    8. Guillot J,
    9. Ménard-Szczebara F,
    10. Breton L,
    11. Latgé J-P,
    12. Mouyna I
    . 2013. Dandruff is associated with disequilibrium in the proportion of the major bacterial and fungal populations colonizing the scalp. PLoS One 8:e58203. doi:10.1371/journal.pone.0058203.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Patiño-Uzcátegui A,
    2. Amado Y,
    3. Cepero de García M,
    4. Chaves D,
    5. Tabima J,
    6. Motta A,
    7. Cárdenas M,
    8. Bernal A,
    9. Restrepo S,
    10. Celis A
    . 2011. Virulence gene expression in Malassezia spp from individuals with seborrheic dermatitis. J Invest Dermatol 131:2134–2136. doi:10.1038/jid.2011.178.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    1. Zhang E,
    2. Tanaka T,
    3. Tajima M,
    4. Tsuboi R,
    5. Nishikawa A,
    6. Sugita T
    . 2011. Characterization of the skin fungal microbiota in patients with atopic dermatitis and in healthy subjects. Microbiol Immunol 55:625–632. doi:10.1111/j.1348-0421.2011.00364.x.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Park T,
    2. Kim HJ,
    3. Myeong NR,
    4. Lee HG,
    5. Kwack I,
    6. Lee J,
    7. Kim BJ,
    8. Sul WJ,
    9. An S
    . 2017. Collapse of human scalp microbiome network in dandruff and seborrhoeic dermatitis. Exp Dermatol 26:835–838. doi:10.1111/exd.13293.
    OpenUrlCrossRef
  25. 25.↵
    1. Wang L,
    2. Clavaud C,
    3. Bar-Hen A,
    4. Cui M,
    5. Gao J,
    6. Liu Y,
    7. Liu C,
    8. Shibagaki N,
    9. Guéniche A,
    10. Jourdain R,
    11. Lan K,
    12. Zhang C,
    13. Altmeyer R,
    14. Breton L
    . 2015. Characterization of the major bacterial–fungal populations colonizing dandruff scalps in Shanghai, China, shows microbial disequilibrium. Exp Dermatol 24:398–400. doi:10.1111/exd.12684.
    OpenUrlCrossRef
  26. 26.↵
    1. Xu Z,
    2. Wang Z,
    3. Yuan C,
    4. Liu X,
    5. Yang F,
    6. Wang T,
    7. Wang J,
    8. Manabe K,
    9. Qin O,
    10. Wang X,
    11. Zhang Y,
    12. Zhang M
    . 2016. Dandruff is associated with the conjoined interactions between host and microorganisms. Sci Rep 6:srep24877. doi:10.1038/srep24877.
    OpenUrlCrossRef
  27. 27.↵
    1. Grice EA,
    2. Dawson TL, Jr..
    2017. Host-microbe interactions: Malassezia and human skin. Curr Opin Microbiol 40:81–87. doi:10.1016/j.mib.2017.10.024.
    OpenUrlCrossRef
  28. 28.↵
    1. Sharma S,
    2. Gupta S,
    3. Shrivastava JN
    . 2011. Presence of virus like particles in human pathogenic fungi: Chrysosporium sps and Candida albicans. Indian J Virol 22:104–110. doi:10.1007/s13337-011-0048-x.
    OpenUrlCrossRef
  29. 29.↵
    1. Park M,
    2. Cho YJ,
    3. Lee YW,
    4. Jung WH
    . 2017. Whole genome sequencing analysis of the cutaneous pathogenic yeast Malassezia restricta and identification of the major lipase expressed on the scalp of patients with dandruff. Mycoses 60:188–197. doi:10.1111/myc.12586.
    OpenUrlCrossRef
  30. 30.↵
    1. Edy VG,
    2. Szekely M,
    3. Loviny T,
    4. Dreyer C
    . 1976. Action of nucleases on double‐stranded RNA. Eur J Biochem 61:563–572. doi:10.1111/j.1432-1033.1976.tb10051.x.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Lee SH,
    2. Yun S-H,
    3. Chun J,
    4. Kim D-H
    . 2017. Characterization of a novel dsRNA mycovirus of Trichoderma atroviride NFCF028. Arch Virol 162:1073–1077. doi:10.1007/s00705-016-3214-z.
    OpenUrlCrossRef
  32. 32.↵
    1. Sommer SS,
    2. Wickner RB
    . 1982. Yeast L dsRNA consists of at least three distinct RNAs; evidence that the non-Mendelian genes [HOK], [NEX] and [EXL] are on one of these dsRNAs. Cell 31:429–441. doi:10.1016/0092-8674(82)90136-2.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Park CM,
    2. Lopinski JD,
    3. Masuda J,
    4. Tzeng TH,
    5. Bruenn JA
    . 1996. A second double-stranded RNA virus from yeast. Virology 216:451–454. doi:10.1006/viro.1996.0083.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Taylor DJ,
    2. Ballinger MJ,
    3. Bowman SM,
    4. Bruenn JA
    . 2013. Virus-host co-evolution under a modified nuclear genetic code. PeerJ 1:e50. doi:10.7717/peerj.50.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Bruenn JA
    . 1993. A closely related group of RNA-dependent RNA polymerases from double-stranded RNA viruses. Nucleic Acids Res 21:5667–5669. doi:10.1093/nar/21.24.5667.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Schmitt MJ,
    2. Breinig F
    . 2002. The viral killer system in yeast: from molecular biology to application. FEMS Microbiol Rev 26:257–276. doi:10.1111/j.1574-6976.2002.tb00614.x.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Wickner RB,
    2. Leibowitz MJ
    . 1979. Mak mutants of yeast: mapping and characterization. J Bacteriol 140:154–160. doi:10.1128/JB.140.1.154-160.1979.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Tercero J,
    2. Dinman J,
    3. Wickner R
    . 1993. Yeast MAK3 N-acetyltransferase recognizes the N-terminal four amino acids of the major coat protein (gag) of the LA double-stranded RNA virus. J Bacteriol 175:3192–3194. doi:10.1128/jb.175.10.3192-3194.1993.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Tercero J,
    2. Wickner RB
    . 1992. MAK3 encodes an N-acetyltransferase whose modification of the LA gag NH2 terminus is necessary for virus particle assembly. J Biol Chem 267:20277–20281.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Fujimura T,
    2. Wickner RB
    . 1987. LA double-stranded RNA viruslike particle replication cycle in Saccharomyces cerevisiae: particle maturation in vitro and effects of mak10 and pet18 mutations. Mol Cell Biol 7:420–426. doi:10.1128/mcb.7.1.420.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Lee Y-J,
    2. Wickner R
    . 1992. MAK10, a glucose-repressible gene necessary for replication of a dsRNA virus of Saccharomyces cerevisiae, has T cell receptor alpha-subunit motifs. Genetics 132:87–96.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Wickner RB,
    2. Ridley SP,
    3. Fried HM,
    4. Ball SG
    . 1982. Ribosomal protein L3 is involved in replication or maintenance of the killer double-stranded RNA genome of Saccharomyces cerevisiae. Proc Natl Acad Sci U S A 79:4706–4708. doi:10.1073/pnas.79.15.4706.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Ohtake Y,
    2. Wickner RB
    . 1995. Yeast virus propagation depends critically on free 60S ribosomal subunit concentration. Mol Cell Biol 15:2772–2781. doi:10.1128/mcb.15.5.2772.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Ryu S-W,
    2. Lee S-J,
    3. Park M-Y,
    4. Jun J-I,
    5. Jung Y-K,
    6. Kim E
    . 2003. Fas-associated factor 1, FAF1, is a member of Fas death-inducing signaling complex. J Biol Chem 278:24003–24010. doi:10.1074/jbc.M302200200.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Fujikura D,
    2. Miyazaki T
    . 2018. Programmed cell death in the pathogenesis of influenza. Int J Mol Sci 19:2065. doi:10.3390/ijms19072065.
    OpenUrlCrossRef
  46. 46.↵
    1. Jorgensen I,
    2. Rayamajhi M,
    3. Miao EA
    . 2017. Programmed cell death as a defence against infection. Nat Rev Immunol 17:151–164. doi:10.1038/nri.2016.147.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Carmona-Gutierrez D,
    2. Bauer MA,
    3. Zimmermann A,
    4. Aguilera A,
    5. Austriaco N,
    6. Ayscough K,
    7. Balzan R,
    8. Bar-Nun S,
    9. Barrientos A,
    10. Belenky P,
    11. Blondel M,
    12. Braun RJ,
    13. Breitenbach M,
    14. Burhans WC,
    15. Büttner S,
    16. Cavalieri D,
    17. Chang M,
    18. Cooper KF,
    19. Côrte-Real M,
    20. Costa V,
    21. Cullin C,
    22. Dawes I,
    23. Dengjel J,
    24. Dickman MB,
    25. Eisenberg T,
    26. Fahrenkrog B,
    27. Fasel N,
    28. Fröhlich K-U,
    29. Gargouri A,
    30. Giannattasio S,
    31. Goffrini P,
    32. Gourlay CW,
    33. Grant CM,
    34. Greenwood MT,
    35. Guaragnella N,
    36. Heger T,
    37. Heinisch J,
    38. Herker E,
    39. Herrmann JM,
    40. Hofer S,
    41. Jiménez-Ruiz A,
    42. Jungwirth H,
    43. Kainz K,
    44. Kontoyiannis DP,
    45. Ludovico P,
    46. Manon S,
    47. Martegani E,
    48. Mazzoni C,
    49. Megeney LA,
    50. Meisinger C,
    51. Nielsen J,
    52. Nyström T,
    53. Osiewacz HD,
    54. Outeiro TF,
    55. Park H-O,
    56. Pendl T,
    57. Petranovic D,
    58. Picot S,
    59. Polčic P,
    60. Powers T,
    61. Ramsdale M,
    62. Rinnerthaler M,
    63. Rockenfeller P,
    64. Ruckenstuhl C,
    65. Schaffrath R,
    66. Segovia M,
    67. Severin FF,
    68. Sharon A,
    69. Sigrist SJ,
    70. Sommer-Ruck C,
    71. Sousa MJ,
    72. Thevelein JM,
    73. Thevissen K,
    74. Titorenko V,
    75. Toledano MB,
    76. Tuite M,
    77. Vögtle F-N,
    78. Westermann B,
    79. Winderickx J,
    80. Wissing S,
    81. Wölfl S,
    82. Zhang ZJ,
    83. Zhao RY,
    84. Zhou B,
    85. Galluzzi L,
    86. Kroemer G,
    87. Madeo F
    . 2018. Guidelines and recommendations on yeast cell death nomenclature. Microb Cell 5:4–31. doi:10.15698/mic2018.01.607.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Newbold A,
    2. Martin BP,
    3. Cullinane C,
    4. Bots M
    . 2014. Detection of apoptotic cells using propidium iodide staining. Cold Spring Harb Protoc 2014:1202–1206. doi:10.1101/pdb.prot082545.
    OpenUrlCrossRef
  49. 49.↵
    1. Banks GT,
    2. Buck KW,
    3. Chain EB,
    4. Himmelweit F,
    5. Marks JE,
    6. Tyler JM,
    7. Hollings M,
    8. Last FT,
    9. Stone OM
    . 1968. Viruses in fungi and interferon stimulation. Nature 218:542–545. doi:10.1038/218542a0.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Kleinschmidt WJ,
    2. Cline JC,
    3. Murphy EB
    . 1964. Interferon production induced by statolon. Proc Natl Acad Sci U S A 52:741–744. doi:10.1073/pnas.52.3.741.
    OpenUrlFREE Full Text
  51. 51.↵
    1. Lampson GP,
    2. Tytell AA,
    3. Field AK,
    4. Nemes MM,
    5. Hilleman MR
    . 1967. Inducers of interferon and host resistance. I. Double-stranded RNA from extracts of Penicillium funiculosum. Proc Natl Acad Sci U S A 58:782–789. doi:10.1073/pnas.58.2.782.
    OpenUrlFREE Full Text
  52. 52.↵
    1. Perales-Linares R,
    2. Navas-Martin S
    . 2013. Toll-like receptor 3 in viral pathogenesis: friend or foe? Immunology 140:153–167. doi:10.1111/imm.12143.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Thompson MR,
    2. Kaminski JJ,
    3. Kurt-Jones EA,
    4. Fitzgerald KA
    . 2011. Pattern recognition receptors and the innate immune response to viral infection. Viruses 3:920–940. doi:10.3390/v3060920.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Claudepierre MC,
    2. Hortelano J,
    3. Schaedler E,
    4. Kleinpeter P,
    5. Geist M,
    6. Remy-Ziller C,
    7. Brandely R,
    8. Tosch C,
    9. Laruelle L,
    10. Jawhari A,
    11. Menguy T,
    12. Marchand JB,
    13. Romby P,
    14. Schultz P,
    15. Hartmann G,
    16. Rooke R,
    17. Bonnefoy JY,
    18. Preville X,
    19. Rittner K
    . 2014. Yeast virus-derived stimulator of the innate immune system augments the efficacy of virus vector-based immunotherapy. J Virol 88:5242–5255. doi:10.1128/JVI.03819-13.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Ishikawa T,
    2. Itoh F,
    3. Yoshida S,
    4. Saijo S,
    5. Matsuzawa T,
    6. Gonoi T,
    7. Saito T,
    8. Okawa Y,
    9. Shibata N,
    10. Miyamoto T,
    11. Yamasaki S
    . 2013. Identification of distinct ligands for the C-type lectin receptors Mincle and Dectin-2 in the pathogenic fungus Malassezia. Cell Host Microbe 13:477–488. doi:10.1016/j.chom.2013.03.008.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Vargas G,
    2. Rocha JD,
    3. Oliveira DL,
    4. Albuquerque PC,
    5. Frases S,
    6. Santos SS,
    7. Nosanchuk JD,
    8. Gomes AM,
    9. Medeiros LC,
    10. Miranda K,
    11. Sobreira TJ,
    12. Nakayasu ES,
    13. Arigi EA,
    14. Casadevall A,
    15. Guimaraes AJ,
    16. Rodrigues ML,
    17. Freire-de-Lima CG,
    18. Almeida IC,
    19. Nimrichter L
    . 2015. Compositional and immunobiological analyses of extracellular vesicles released by Candida albicans. Cell Microbiol 17:389–407. doi:10.1111/cmi.12374.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Samuel CE
    . 2001. Antiviral actions of interferons. Clin Microbiol Rev 14:778–809. doi:10.1128/CMR.14.4.778-809.2001.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    1. Brierley I,
    2. Boursnell ME,
    3. Binns MM,
    4. Bilimoria B,
    5. Blok VC,
    6. Brown TD,
    7. Inglis SC
    . 1987. An efficient ribosomal frame-shifting signal in the polymerase-encoding region of the coronavirus IBV. EMBO J 6:3779–3785. doi:10.1002/j.1460-2075.1987.tb02713.x.
    OpenUrlCrossRef
  59. 59.↵
    1. Jacks T,
    2. Varmus HE
    . 1985. Expression of the Rous sarcoma virus pol gene by ribosomal frameshifting. Science 230:1237–1242. doi:10.1126/science.2416054.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    1. Jacks T,
    2. Madhani HD,
    3. Masiarz FR,
    4. Varmus HE
    . 1988. Signals for ribosomal frameshifting in the Rous sarcoma virus gag-pol region. Cell 55:447–458. doi:10.1016/0092-8674(88)90031-1.
    OpenUrlCrossRefPubMedWeb of Science
  61. 61.↵
    1. Magliani W,
    2. Conti S,
    3. Gerloni M,
    4. Bertolotti D,
    5. Polonelli L
    . 1997. Yeast killer systems. Clin Microbiol Rev 10:369–400. doi:10.1128/CMR.10.3.369-400.1997.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Tercero JC,
    2. Riles LE,
    3. Wickner RB
    . 1992. Localized mutagenesis and evidence for post-transcriptional regulation of MAK3. A putative N-acetyltransferase required for double-stranded RNA virus propagation in Saccharomyces cerevisiae. J Biol Chem 267:20270–20276.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    1. Wickner RB
    . 1996. Double-stranded RNA viruses of Saccharomyces cerevisiae. Microbiol Rev 60:250–265. doi:10.1128/MMBR.60.1.250-265.1996.
    OpenUrlFREE Full Text
  64. 64.↵
    1. Thrash C,
    2. Bankier AT,
    3. Barrell BG,
    4. Sternglanz R
    . 1985. Cloning, characterization, and sequence of the yeast DNA topoisomerase I gene. Proc Natl Acad Sci U S A 82:4374–4378. doi:10.1073/pnas.82.13.4374.
    OpenUrlAbstract/FREE Full Text
  65. 65.↵
    1. Widner WR,
    2. Wickner RB
    . 1993. Evidence that the SKI antiviral system of Saccharomyces cerevisiae acts by blocking expression of viral mRNA. Mol Cell Biol 13:4331–4341. doi:10.1128/mcb.13.7.4331.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    1. Masison DC,
    2. Blanc A,
    3. Ribas JC,
    4. Carroll K,
    5. Sonenberg N,
    6. Wickner RB
    . 1995. Decoying the cap-mRNA degradation system by a double-stranded RNA virus and poly(A)-mRNA surveillance by a yeast antiviral system. Mol Cell Biol 15:2763–2771. doi:10.1128/mcb.15.5.2763.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    1. Benard L,
    2. Carroll K,
    3. Valle RC,
    4. Masison DC,
    5. Wickner RB
    . 1999. The ski7 antiviral protein is an EF1-α homolog that blocks expression of non-poly(A) mRNA in Saccharomyces cerevisiae. J Virol 73:2893–2900. doi:10.1128/JVI.73.4.2893-2900.1999.
    OpenUrlAbstract/FREE Full Text
  68. 68.↵
    1. McBride RC,
    2. Boucher N,
    3. Park DS,
    4. Turner PE,
    5. Townsend JP
    . 2013. Yeast response to LA virus indicates coadapted global gene expression during mycoviral infection. FEMS Yeast Res 13:162–179. doi:10.1111/1567-1364.12019.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Wu M,
    2. Jin F,
    3. Zhang J,
    4. Yang L,
    5. Jiang D,
    6. Li G
    . 2012. Characterization of a novel bipartite double-stranded RNA mycovirus conferring hypovirulence in the phytopathogenic fungus Botrytis porri. J Virol 86:6605–6619. doi:10.1128/JVI.00292-12.
    OpenUrlAbstract/FREE Full Text
  70. 70.↵
    1. Wan J,
    2. Basu K,
    3. Mui J,
    4. Vali H,
    5. Zheng H,
    6. Laliberte JF
    . 2015. Ultrastructural characterization of turnip mosaic virus-induced cellular rearrangements reveals membrane-bound viral particles accumulating in vacuoles. J Virol 89:12441–12456. doi:10.1128/JVI.02138-15.
    OpenUrlAbstract/FREE Full Text
  71. 71.↵
    1. Matsumoto M,
    2. Oshiumi H,
    3. Seya T
    . 2011. Antiviral responses induced by the TLR3 pathway. Rev Med Virol 21:67–77. doi:10.1002/rmv.680.
    OpenUrlCrossRefPubMed
  72. 72.↵
    1. Tabeta K,
    2. Georgel P,
    3. Janssen E,
    4. Du X,
    5. Hoebe K,
    6. Crozat K,
    7. Mudd S,
    8. Shamel L,
    9. Sovath S,
    10. Goode J,
    11. Alexopoulou L,
    12. Flavell RA,
    13. Beutler B
    . 2004. Toll-like receptors 9 and 3 as essential components of innate immune defense against mouse cytomegalovirus infection. Proc Natl Acad Sci U S A 101:3516–3521. doi:10.1073/pnas.0400525101.
    OpenUrlAbstract/FREE Full Text
  73. 73.↵
    1. Hartley MA,
    2. Ronet C,
    3. Zangger H,
    4. Beverley SM,
    5. Fasel N
    . 2012. Leishmania RNA virus: when the host pays the toll. Front Cell Infect Microbiol 2:99. doi:10.3389/fcimb.2012.00099.
    OpenUrlCrossRef
  74. 74.↵
    1. Kim M,
    2. Cho Y-J,
    3. Park M,
    4. Choi Y,
    5. Hwang SY,
    6. Jung WH
    . 2018. Genomic tandem quadruplication is associated with ketoconazole resistance in Malassezia pachydermatis. J Microbiol Biotechnol 28:1937–1945. doi:10.4014/jmb.1810.10019.
    OpenUrlCrossRef
  75. 75.↵
    1. Park M,
    2. Jung WH,
    3. Han SH,
    4. Lee YH,
    5. Lee YW
    . 2015. Characterisation and expression analysis of MrLip1, a class 3 family lipase of Malassezia restricta. Mycoses 58:671–678. doi:10.1111/myc.12412.
    OpenUrlCrossRef
  76. 76.↵
    1. Leeming JP,
    2. Notman FH
    . 1987. Improved methods for isolation and enumeration of Malassezia furfur from human skin. J Clin Microbiol 25:2017–2019. doi:10.1128/JCM.25.10.2017-2019.1987.
    OpenUrlAbstract/FREE Full Text
  77. 77.↵
    1. Okada R,
    2. Kiyota E,
    3. Moriyama H,
    4. Fukuhara T,
    5. Natsuaki T
    . 2015. A simple and rapid method to purify viral dsRNA from plant and fungal tissue. J Gen Plant Pathol 81:103–107. doi:10.1007/s10327-014-0575-6.
    OpenUrlCrossRef
  78. 78.↵
    1. Moon S,
    2. Byun Y,
    3. Kim HJ,
    4. Jeong S,
    5. Han K
    . 2004. Predicting genes expressed via -1 and +1 frameshifts. Nucleic Acids Res 32:4884–4892. doi:10.1093/nar/gkh829.
    OpenUrlCrossRefPubMedWeb of Science
  79. 79.↵
    1. Waterhouse AM,
    2. Procter JB,
    3. Martin DM,
    4. Clamp M,
    5. Barton GJ
    . 2009. Jalview version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics 25:1189–1191. doi:10.1093/bioinformatics/btp033.
    OpenUrlCrossRefPubMedWeb of Science
  80. 80.↵
    1. Kumar S,
    2. Nei M,
    3. Dudley J,
    4. Tamura K
    . 2008. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform 9:299–306. doi:10.1093/bib/bbn017.
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    1. Kumar S,
    2. Stecher G,
    3. Li M,
    4. Knyaz C,
    5. Tamura K
    . 2018. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol Biol Evol 35:1547–1549. doi:10.1093/molbev/msy096.
    OpenUrlCrossRefPubMed
  82. 82.↵
    1. Koh H-J,
    2. Kim Y-R,
    3. Kim J-S,
    4. Yun J-S,
    5. Kim S,
    6. Kim SY,
    7. Jang K,
    8. Yang C-S
    . 2018. CD82 hypomethylation is essential for tuberculosis pathogenesis via regulation of RUNX1-Rab5/22. Exp Mol Med 50:1–15. doi:10.1038/s12276-018-0091-4.
    OpenUrlCrossRefPubMed
  83. 83.↵
    1. Saitou N,
    2. Nei M
    . 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406–425. doi:10.1093/oxfordjournals.molbev.a040454.
    OpenUrlCrossRefPubMedWeb of Science
  84. 84.↵
    1. Felsenstein J
    . 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783–791. doi:10.1111/j.1558-5646.1985.tb00420.x.
    OpenUrlCrossRefPubMedWeb of Science
  85. 85.↵
    1. Sanderlin RS,
    2. Ghabrial SA
    . 1978. Physicochemical properties of two distinct types of virus-like particles from Helminthosporium victoriae. Virology 87:142–151. doi:10.1016/0042-6822(78)90166-6.
    OpenUrlCrossRefPubMed
  86. 86.↵
    1. Shapira R,
    2. Choi GH,
    3. Nuss DL
    . 1991. Virus-like genetic organization and expression strategy for a double-stranded RNA genetic element associated with biological control of chestnut blight. EMBO J 10:731–739. doi:10.1002/j.1460-2075.1991.tb08004.x.
    OpenUrlCrossRefPubMedWeb of Science
  87. 87.↵
    1. Hillman BI,
    2. Halpern BT,
    3. Brown MP
    . 1994. A viral dsRNA element of the chestnut blight fungus with a distinct genetic organization. Virology 201:241–250. doi:10.1006/viro.1994.1289.
    OpenUrlCrossRefPubMedWeb of Science
  88. 88.↵
    1. Van der Lende TR,
    2. Duitman EH,
    3. Gunnewijk MG,
    4. Yu L,
    5. Wessels JG
    . 1996. Functional analysis of dsRNAs (L1, L3, L5, and M2) associated with isometric 34-nm virions of Agaricus bisporus (white button mushroom). Virology 217:88–96. doi:10.1006/viro.1996.0096.
    OpenUrlCrossRefPubMed
  89. 89.↵
    1. Nomura K,
    2. Osaki H,
    3. Iwanami T,
    4. Matsumoto N,
    5. Ohtsu Y
    . 2003. Cloning and characterization of a totivirus double-stranded RNA from the plant pathogenic fungus, Helicobasidium mompa Tanaka. Virus Genes 26:219–226. doi:10.1023/a:1024453111809.
    OpenUrlCrossRefPubMedWeb of Science
  90. 90.↵
    1. Tuomivirta TT,
    2. Hantula J
    . 2003. Two unrelated double-stranded RNA molecule patterns in Gremmeniella abietina type A code for putative viruses of the families Totiviridae and Partitiviridae. Arch Virol 148:2293–2305. doi:10.1007/s00705-003-0194-6.
    OpenUrlCrossRefPubMed
  91. 91.↵
    1. Jiang D,
    2. Ghabrial SA
    . 2004. Molecular characterization of Penicillium chrysogenum virus: reconsideration of the taxonomy of the genus Chrysovirus. J Gen Virol 85:2111–2121. doi:10.1099/vir.0.79842-0.
    OpenUrlCrossRefPubMedWeb of Science
  92. 92.↵
    1. Urayama S,
    2. Ohta T,
    3. Onozuka N,
    4. Sakoda H,
    5. Fukuhara T,
    6. Arie T,
    7. Teraoka T,
    8. Moriyama H
    . 2012. Characterization of Magnaporthe oryzae chrysovirus 1 structural proteins and their expression in Saccharomyces cerevisiae. J Virol 86:8287–8295. doi:10.1128/JVI.00871-12.
    OpenUrlAbstract/FREE Full Text
  93. 93.↵
    1. Wu Q,
    2. Luo Y,
    3. Lu R,
    4. Lau N,
    5. Lai EC,
    6. Li WX,
    7. Ding SW
    . 2010. Virus discovery by deep sequencing and assembly of virus-derived small silencing RNAs. Proc Natl Acad Sci U S A 107:1606–1611. doi:10.1073/pnas.0911353107.
    OpenUrlAbstract/FREE Full Text
  94. 94.↵
    1. Isawa H,
    2. Kuwata R,
    3. Hoshino K,
    4. Tsuda Y,
    5. Sakai K,
    6. Watanabe S,
    7. Nishimura M,
    8. Satho T,
    9. Kataoka M,
    10. Nagata N,
    11. Hasegawa H,
    12. Bando H,
    13. Yano K,
    14. Sasaki T,
    15. Kobayashi M,
    16. Mizutani T,
    17. Sawabe K
    . 2011. Identification and molecular characterization of a new nonsegmented double-stranded RNA virus isolated from Culex mosquitoes in Japan. Virus Res 155:147–155. doi:10.1016/j.virusres.2010.09.013.
    OpenUrlCrossRefPubMedWeb of Science
  95. 95.↵
    1. Mor SK,
    2. Phelps NB
    . 2016. Molecular detection of a novel totivirus from golden shiner (Notemigonus crysoleucas) baitfish in the USA. Arch Virol 161:2227–2234. doi:10.1007/s00705-016-2906-8.
    OpenUrlCrossRef
  96. 96.↵
    1. Yang X,
    2. Zhang Y,
    3. Ge X,
    4. Yuan J,
    5. Shi Z
    . 2012. A novel totivirus-like virus isolated from bat guano. Arch Virol 157:1093–1099. doi:10.1007/s00705-012-1278-y.
    OpenUrlCrossRefPubMed
  97. 97.↵
    1. Zhai Y,
    2. Attoui H,
    3. Mohd Jaafar F,
    4. Wang HQ,
    5. Cao YX,
    6. Fan SP,
    7. Sun YX,
    8. Liu LD,
    9. Mertens PP,
    10. Meng WS,
    11. Wang D,
    12. Liang G
    . 2010. Isolation and full-length sequence analysis of Armigeres subalbatus totivirus, the first totivirus isolate from mosquitoes representing a proposed novel genus (Artivirus) of the family Totiviridae. J Gen Virol 91:2836–2845. doi:10.1099/vir.0.024794-0.
    OpenUrlCrossRefPubMedWeb of Science
  98. 98.↵
    1. Dihanich M,
    2. Van Tuinen E,
    3. Lambris J,
    4. Marshallsay B
    . 1989. Accumulation of viruslike particles in a yeast mutant lacking a mitochondrial pore protein. Mol Cell Biol 9:1100–1108. doi:10.1128/mcb.9.3.1100.
    OpenUrlAbstract/FREE Full Text
  99. 99.↵
    1. Leibowitz MJ,
    2. Wickner RB
    . 1978. pet 18: a chromosomal gene required for cell growth and for the maintenance of mitochondrial DNA and the killer plasmid of yeast. Mol Gen Genet 165:115–121. doi:10.1007/BF00269899.
    OpenUrlCrossRefPubMed
  100. 100.↵
    1. Toh-e A,
    2. Sahashi Y
    . 1985. The PET18 locus of Saccharomyces cerevisiae: a complex locus containing multiple genes. Yeast 1:159–171. doi:10.1002/yea.320010204.
    OpenUrlCrossRefPubMed
  101. 101.↵
    1. Thrash C,
    2. Voelkel K,
    3. DiNardo S,
    4. Sternglanz R
    . 1984. Identification of Saccharomyces cerevisiae mutants deficient in DNA topoisomerase I activity. J Biol Chem 259:1375–1377.
    OpenUrlAbstract/FREE Full Text
  102. 102.↵
    1. Icho T,
    2. Wickner RB
    . 1988. The MAK11 protein is essential for cell growth and replication of M double-stranded RNA and is apparently a membrane-associated protein. J Biol Chem 263:1467–1475.
    OpenUrlAbstract/FREE Full Text
  103. 103.↵
    1. Wickner RB
    . 1988. Host function of MAK16: G1 arrest by a mak16 mutant of Saccharomyces cerevisiae. Proc Natl Acad Sci U S A 85:6007–6011. doi:10.1073/pnas.85.16.6007.
    OpenUrlAbstract/FREE Full Text
  104. 104.↵
    1. Carroll K,
    2. Wickner RB
    . 1995. Translation and M1 double-stranded RNA propagation: MAK18= RPL41B and cycloheximide curing. J Bacteriol 177:2887–2891. doi:10.1128/jb.177.10.2887-2891.1995.
    OpenUrlAbstract/FREE Full Text
  105. 105.↵
    1. Toh-E A,
    2. Guerry P,
    3. Wickner RB
    . 1978. Chromosomal superkiller mutants of Saccharomyces cerevisiae. J Bacteriol 136:1002–1007. doi:10.1128/JB.136.3.1002-1007.1978.
    OpenUrlAbstract/FREE Full Text
  106. 106.↵
    1. Ball SG,
    2. Tirtiaux C,
    3. Wickner RB
    . 1984. Genetic control of LA and L-(BC) dsRNA copy number in killer systems of Saccharomyces cerevisiae. Genetics 107:199–217.
    OpenUrlAbstract/FREE Full Text
  107. 107.↵
    1. Ridley SP,
    2. Sommer SS,
    3. Wickner RB
    . 1984. Superkiller mutations in Saccharomyces cerevisiae suppress exclusion of M2 double-stranded RNA by LA-HN and confer cold sensitivity in the presence of M and LA-HN. Mol Cell Biol 4:761–770. doi:10.1128/mcb.4.4.761.
    OpenUrlAbstract/FREE Full Text
  108. 108.↵
    1. Dmochowska A,
    2. Dignard D,
    3. Henning D,
    4. Thomas DY,
    5. Bussey H
    . 1987. Yeast KEX1 gene encodes a putative protease with a carboxypeptidase B-like function involved in killer toxin and α-factor precursor processing. Cell 50:573–584. doi:10.1016/0092-8674(87)90030-4.
    OpenUrlCrossRefPubMedWeb of Science
  109. 109.↵
    1. Wickner RB
    . 1974. Chromosomal and nonchromosomal mutations affecting the “killer character” of Saccharomyces cerevisiae. Genetics 76:423–432.
    OpenUrlAbstract/FREE Full Text
  110. 110.↵
    1. Wickner RB,
    2. Leibowitz MJ
    . 1976. Two chromosomal genes required for killing expression in killer strains of Saccharomyces cerevisiae. Genetics 82:429–442.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
A Novel Virus Alters Gene Expression and Vacuolar Morphology in Malassezia Cells and Induces a TLR3-Mediated Inflammatory Immune Response
Minji Park, Yong-Joon Cho, Donggyu Kim, Chul-Su Yang, Shi Mun Lee, Thomas L. Dawson Jr., Satoshi Nakamizo, Kenji Kabashima, Yang Won Lee, Won Hee Jung
mBio Sep 2020, 11 (5) e01521-20; DOI: 10.1128/mBio.01521-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
A Novel Virus Alters Gene Expression and Vacuolar Morphology in Malassezia Cells and Induces a TLR3-Mediated Inflammatory Immune Response
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
A Novel Virus Alters Gene Expression and Vacuolar Morphology in Malassezia Cells and Induces a TLR3-Mediated Inflammatory Immune Response
Minji Park, Yong-Joon Cho, Donggyu Kim, Chul-Su Yang, Shi Mun Lee, Thomas L. Dawson Jr., Satoshi Nakamizo, Kenji Kabashima, Yang Won Lee, Won Hee Jung
mBio Sep 2020, 11 (5) e01521-20; DOI: 10.1128/mBio.01521-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Malassezia
Malassezia restricta
mycovirus
TLR3
cytokine
double-stranded RNA virus
virus
yeast virus

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511