{'fraction_non_alphanumeric': 0.03292548979539707, 'fraction_numerical': 0.019208664333699457, 'mean_word_length': 5.433853738701726, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 4, 'https://': 6, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Evolution of zygomycete secretomes and the origins of terrestrial fungal ecologies

Summary Fungi survive in diverse ecological niches by secreting proteins and other molecules into the environment to acquire food and interact with various biotic and abiotic stressors. Fungal secretome content is, therefore, believed to be tightly linked to fungal ecologies. We sampled 132 genomes from the early-diverging terrestrial fungal lineage zygomycetes (Mucoromycota and Zoopagomycota) and characterized their secretome composition. Our analyses revealed that phylogeny played an important role in shaping the secretome composition of zygomycete fungi with trophic mode contributing a smaller amount. Reconstruction of the evolution of secreted digestive enzymes revealed lineage-specific expansions, indicating that Mucoromycota and Zoopagomycota followed different trajectories early in their evolutionary history. We identified the presence of multiple pathogenicity-related proteins in the lineages known as saprotrophs, suggesting that either the ecologies of these fungi are incompletely known, and/or that these pathogenicity-related proteins have important functions associated with saprotrophic ecologies, both of which invite further investigation.


INTRODUCTION
Trophic modes and ecological guilds are some of the most defining traits of fungi. Ecologists typically classify fungi into trophic modes based on their association and outcome with hosts and substrates (Nguyen et al., 2016). Symbiotrophs associate with living hosts forming beneficial relationships with both host and fungus, such as mycorrhizae and lichens. Pathotrophs cause disease and ultimately harm their hosts, ranging from biotrophy to necrotrophy; whereas, saprotrophs degrade non-living organic matter, playing an essential role in nutrient cycling. Trophic modes are further divided into ecological guilds, such as mycorrhizae, plant pathogens, animal pathogens, and mycoparasites, which provide greater specificity with respect to host/substrate affiliation and nature of association. Although these categories may be imprecise in some situations, and some fungi may occupy more than one trophic mode or ecological guild, the ecology of a particular fungus is largely based on its ability to secrete biologically active molecules into the environment. These molecules define the heterotrophic nutritional mode of a fungus through its ability to interact with other organisms, decompose substrates comprising complex biopolymers, or obtain nutrients from living hosts through antagonistic, beneficial, or not-yet-defined strategies. As such, fungi are what they secrete.
The importance of fungal secretome has been analyzed in numerous ecological guilds of the fleshy, sporocarp-producing Dikarya. Mycorrhizae are characterized by the reduction of enzymes associated with the breakdown of plant cell walls, and the production of small secreted proteins and effectors that interact with specific plant hormone pathways (Plett and Martin, 2015;Romero et al., 2021). The combination of these features is hypothesized to reduce or inhibit the host defense response, allowing the fungus-plant symbiosis to establish and function. Wood decay fungi possess an arsenal of degradative enzymes acting on different constituents of plant cell walls, and the suite of enzymes employed by a given fungus will determine what plant cell walls biopolymers (e.g., cellulose, lignin, and so forth) will be degraded (Floudas et al., 2012). Mycorrhizae and wood decay represent relatively well-studied ecological guilds in extensively characterized systems of Dikarya. The majority of fungal species, however, are poorly studied and their function in nature is characterized by limited observations or extrapolations from related species. In this context, the analysis of fungal secretomes represents an opportunity to predict the potential ecologies of many fungi whose function(s) in ecosystems are difficult to observe and characterize.
Zygomycete fungi comprise the first diverging lineages of Kingdom Fungi to emerge after the loss of the flagellum and the origin of the filamentous fungal thallus (Kiss et al., 2019). The two phylum level clades, Mucoromycota and Zoopagomycota, have been resolved as a paraphyletic grade Spatafora et al., 2016) or sister groups , depending on taxon sampling and the molecular markers employed. Zygomycetes are ubiquitous in nature and with diverse ecologies. Zoopagomycota comprises pathogens of insects and herptile gut fungi of Entomophthoromycotina; symbionts of aquatic arthropods and mycoparasites of Kickxellomycotina; and mycoparasites and parasites of amoebae, rotifers, and nematodes of Zoopagomycotina. Based on our current understanding of their ecologies, Zoopagomycota is uniquely devoid of plant associations. However, some species can be maintained on simple media consisting of plant and fungal carbohydrates, consistent with a potential for saprotrophy. Conversely, the dominant ecologies of Mucoromycota are plant-associated and include arbuscular mycorrhizae of Glomeromycotina; root endophytes and chitinolytic soil saprotrophs of Mortierellomycotina; and saprotrophs of plant material, mycoparasites, ectomycorrhizae, plant pathogens and opportunistic animal pathogens of Mucoromycotina.
This diversification of ecologies occurred early in the evolution of nonflagellated, terrestrial fungi, and set these two fungal phyla on different evolutionary trajectories. Unfortunately, many of these fungi remain understudied owing to their cryptic nature, small number of economically important species, and a lack of genetic tools. Recently, however, genome sequencing efforts have significantly increased the phylogenetic diversity of zygomycete genomes. This windfall of data provides the unique opportunity to explore the ecological potentials of zygomycete fungi through analysis of their secretomes. Here we report the analyses of the secretome of Kingdom Fungi (Table S1) with a focus on digestive enzymes, putative small secreted proteins (SSPs), and other proteins hypothesized to function in the environment. Our sampling emphasizes the understudied, non-Dikarya phyla Mucoromycota and Zoopagomycota, and testing hypotheses about the evolution of fungal trophic modes and ecological guilds of lineages associated with the origins of terrestrial fungi.

Overview of fungal secretomes
On average, 5.7% of the total number predicted proteins (referred as proteome later in discussion) were secreted to the extracellular environment, with the lowest being 1.1% (Catenaria, Blastocladiomycota) and the highest being 12.0% (Acaulopage, Zoopagomycotina, Zoopagomycota) (Table S2). Overall, Zoopagomycota had a higher relative secretome abundance than other phyla (Table S2 and Figure S2), with an average of 8.5% (Table S2 and Figure S2). The average secretome:proteome ratio for Mucoromycota was 4.7% and varied across different subphyla (2.8% for Glomeromycotina, 4.7% for Mucoromycotina, and 6.9% for Mortierellomycotina). Secretome composition comprised on average 25% digestive enzymes, 46% SSPs, and 29% OTHER proteins ( Figure S2). Glomeromycotina stood out by having the highest SSP content (67%) and the lowest digestive enzyme content (11%). Among zygomycetes, there was a strong negative correlation between the relative abundance of digestive enzymes and the relative abundance of SSPs ( Figure 1).

Phylogenetic vs. ecological signals in secretome variation
Non-metric multidimensional scaling (NMDS) analyses revealed a strong influence of phylogenetic relatedness on the secretome content regardless of whether analyzed as total secretome (Figure 2) or by grouping proteins into digestive enzymes, SSPs, or OTHER ( Figure S3). Genomes from Mucoromycota and Zoopagomycota were separated from each other in the NMDS space. At the subphylum level, the Mucoromycota secretomes showed similar secretome composition within each subphylum. In contrast, the genomes from each Zoopagomycota subphylum had a more diverse secretome composition and were more loosely grouped in the NMDS plot ( Figure 2A). NMDS analyses also revealed patterns associated with taxon primary guild ( Figure 2A) among Zoopagomycota genomes. Within Entomophthoromycotina, the three herptile gut fungi (Basidiobolus) were clustered apart from the animal pathogenic members. Similarly, the animal symbionts and the saprobes from Kickxellomycotina formed two separate groups. Within Zoopagomycotina, mycoparasites were separated from the animal pathogens. The dual influence from phylogenetic relatedness and trophic ecologies was corroborated by results from Mantel tests. Phylogenetic relatedness was strongly correlated with total secretome content (r = 0.66-0.72, p < 0.01). Ecological guild showed significant but weaker correlation with total secretome content in partial and regular Mantel tests for both primary guild coding (r = 0.41 & 0.56; p < 0.01) and multiple guild coding (r = 0.44 & 0.47; p < 0.01). Similar patterns were recovered for individual protein groups of digestive enzymes, SSPs, and OTHER proteins ( Figure 2B).
In the NMDS analysis, when enzyme vectors were fit to the ordination plot, 41 enzymes were identified with strong r values (squared correlation coefficient, r > 0.2 & p < 0.001), including 33 cazymes, seven peptidases, and one lipase (Table S4, Figure 3C). Plant cell wall degrading enzymes (PCWDEs) were enriched in Mucoromycotina. In contrast, secretomes of Zoopagomycota, Glomeromycotina, and Mortierellomycotina were enriched in chitinases and chitin-binding modules, including AA11, CBM18, and GH19 (Zoopagomycota only). Chitinase family GH18 was widely distributed across the sampled zygomycete genomes but missing from Glomeromycotina. Chitin-binding CBM5 was commonly found in Mucoromycota but not in Zoopagomycota except for in two of the Basidiobolus genomes. The two peptidase families S01A and A01A had the highest NMDS r scores among all the digestive enzymes (r = 0.71 & p < 0.001 for both). Serine peptidase S01A was enriched in Kickxellomycotina and Entomophthoromycotina, with average copy numbers of 23 and 20, respectively, whereas this protein was missing or depauperate in the other fungi (average less than one copy per genome). In contrast, aspartic peptidase A01A was enriched in Mucoromycota (average of 17 copies) but found in low copy numbers in Zoopagomycota (average of two copies). Peptidase S10 had an even distribution across most of the genomes, with relatively low copy numbers (present in 122 out of 132 genomes, average of four copies per genome), but was absent from 70% of Zoopagomycotina genomes.

Small secreted proteins in zygomycete fungal genomes
A total of 47,494 SSPs were identified across the 132 zygomycete genomes. They formed 3,144 nonsingleton orthologous groups, including 35,134 SSPs (74%). 97.5% of the OGs were shared by two or more species (3,066). Among them, over two-thirds (2,222) did not show sequence similarity to known pfam proteins.
2 3 2 clustering analyses of the SSP non-singleton OGs identified the enrichment of SSPs shared by closely related taxa with a similar ecological guild ( Figure 4), but convergent ecologies across different subphyla did not share enriched orthologous OGs. No indicator OGs (i.e., the proteins OGs that are found mainly associated with particular taxa) were identified at the phylum level. At the subphylum level, Entomophthoromycotina, Glomeromycotina, and Mortierellomycotina had over 180 SSP-OGs identified as indicator OGs, while the other three subphyla had less than 70 ( Figure S5). The majority of indicator OGs were found with high specificity A values ( Figure S5) but lower fidelity scores (B values), indicating that most of the indicator SSP OGs were specific to, or enriched in a particular subphylum, but these OGs had limited distribution within this subphylum ( Figure S5). Figure 3. Digestive enzyme composition and evolution across sampled zygomycete genomes (A) Enzyme evolution was reconstructed using NOTUNG based on the RAxML species tree and best enzyme trees with bootstrap support. Bars above internodes are the total number of net gains (light grey to black) or net losses (pink to red) for cazymes, peptidases, and lipases. Colored bars by taxon names represent the primary guild of the taxon. (B) Heatmap of raw counts of PCWDEs, chitinases, and peptidases identified as strong vectors in NMDS analysis (r > 0.2). (C) Plot of NMDS analysis of digestive enzyme composition across sampled genomes. Each genome is color coded by its subphylum placement with the shape representing its primary guild. Arrows represent the strong vectors of PCWDEs, chitinases, and peptidases (r > 0.2).

OPEN ACCESS
iScience 25, 104840, August 19, 2022 5 iScience Article For Glomeromycotina, many of the indicator OGs with pfam hits were related to pathogenesis and/or defense (Table S5). Similarly, many of the indicator OGs for Entomophthoromycotina and Mortierellomycotina showed sequence similarity to proteins involved in defense (Table S5). In addition, some indicator proteins for Entomophthoromycotina and Mortierellomycotina were related to nutrient acquisition (Table S5). Though not identified as an indicator protein, LysM (PF01476) was found enriched in Mucoromycotina, Mortierellomycotina, and Basidiobolus (Entomophthoromycotina, Figure S6). When taxa were grouped by their primary guild, none or a small number of indicator OGs were identified for each category, except for the arbuscular mycorrhizal fungi which were present only in Glomeromycotina ( Figures S5C  and S5D).

OTHER proteins in zygomycete fungal genomes
A total of 22,712 OTHER proteins were identified across the 132 zygomycete genomes. They formed 1,246 non-singleton OGs, including 20,134 proteins (88.6%). Among the non-singleton OGs, 1,213 (97.35%) OGs were shared by two or more species. Approximately half of the OTHER-OGs (657) had hits to known protein families in the Pfam database.
In the IndicSpecies analysis, six OTHER-OGs were identified as indicator OGs for Zoopagomycota. Three of them showed sequence similarity to Pfam proteins (PF00264-a common central domain of tyrosinase, PF00188-a cysteine-rich secretory protein family, and PF03009-a diester phosphodiesterase family). In contrast, no indicator OTHER-OGs were identified for Mucoromycota. At the subphylum level, Mortierellomycotina had the highest number of indicator OGs (109), followed by Entomophthoromycotina (68), Glomeromycotina (39), Kickxellomycotina (35), Zoopagomycotina (23), and Mucoromycotina (12) (Figure S7). Many of these indicator OTHER-OGs were associated with interactions with other organisms or dealing with environmental stress (Table S6), including cotH proteins (PF08757) that are known to be involved in mucormycosis ( Figure S6)

DISCUSSION
Phylogeny is the major driver of secretome content Earlier studies of Dikarya found that secretome size was directly related to phylogeny but not to trophic modes (e.g., Krijger et al., 2014;Stolzer et al., 2012), whereas Lowe and Howlett (2012) reported a correlation between trophic ecology and secretome size for some fungi. Alfaro et al. (2014) argued that lifestyle shapes secretome composition, rather than secretome size. Here, we found that zygomycete secretome composition is primarily determined by phylogeny, with trophic mode and ecological guild playing an important secondary role in the diversification of secretome content (Figures 2, 5 and S3).
Mantel tests suggest that phylogenetic relatedness strongly affects the secretome composition and that trophic ecology played a minor role ( Figure 2). This is consistent with the adaptive potential of a fungus being largely determined by the genetic background inherited from its ancestors, which results in fewer shared homologies in secretome content among fungi with similar trophic ecologies but with different phylogenetic backgrounds. Although fungi with the same trophic ecology may utilize common nutrient resources, and encounter similar environmental stresses, they have developed independent, non-orthologous enzymatic and molecular mechanisms to do so.

Early diversification of nutritional modes with the emergence of nonflagellated, terrestrial fungi
Gene tree -species tree reconciliation analyses support early diversifications of digestive enzymes prior to the loss of the flagellum and origin of nonflagellated fungi (Figure 3). Reconstruction of digestive enzyme evolution concurrent with the emergence of nonflagellated, terrestrial fungi supports an early split in nutritional modes between Mucoromycota and Zoopagomycota. The majority of Mucoromycota taxa are associated with plants or plant-derived substrates, whereas Zoopagomycota taxa are mainly associated with animals or other fungi. PCWDEs are enriched in Mucoromycotina taxa, consistent with ecologies dependent on the breakdown of plant cell wall carbohydrates (Figure 2 and Table S7). In contrast, Zoopagomycota taxa possess more chitinases and peptidases, consistent with animal-and/or fungal-associated ecologies and the utilization of chitin and proteins as nutritional sources.
Rather than being phylum-level secretome synapomorphies, however, the emergence of these two ecological trajectories is characterized by ancient losses and differential retention of proteins in the MRCAs of Mucoromycota and Zoopagomycota, followed by more derived, lineage-specific secretome diversifications. For example, even though both lineages underwent losses of PCWDEs in their early evolutionary history, the MRCA of Zoopagomycota was estimated to possess only 18 copies of PCWDEs and few or no pectinases, while the MRCA of Mucoromycota is hypothesized to have retained more PCWDEs (75 copies), targeting a diverse array of plant-derived polysaccharides (Table S7). Subsequent diversifications of extant ecologies are the product of more derived expansions of enzyme families within subphyla, as evidenced in the increase in PCWDEs, chitinases, and peptidases within Mucorales, chitinases, and peptidases in Mortierellomycotina, and independent expansions of chitinases and peptidases within the subphyla Entomophthoromycotina and Kickxellomycotina (Figure 3).
These results are consistent with the hypothesis that Mucoromycota fungi were associated with streptophytes early in their evolutionary history. The earliest Zoopagomycota fungi, on the other hand, adapted to animal-based and fungal-based ecologies and utilized nutrient sources other than biopolymers derived from streptophyte cell walls. As such, early diversification events in Kingdom Fungi involved two different phylogenetic trajectories defined by host association and substrate utilization.

Evolution of fungal cell wall masking LysM domain
LysM effectors are small effector proteins with one or more Lysin motifs (LysM; PF01476). LysM motifs are typically between 44 and 65 amino acids long and bind to peptidoglycan, chitin, and their derivatives (Dubey et al., 2020 iScience Article the establishment of symbiosis or to defense response against pathogens (Hu et al., 2021). Fungi, on the other hand, also use LysM effectors for a diversity of functions. By binding to chitin in the fungal cell walls, LysM effector helps the fungi to avoid detection by the plant immune system and to protect fungi against mycoparasitism (Kombrink and Thomma, 2013). LysM effectors can also immobilize bacterial competitors through binding of bacterial peptidoglycan (Kombrink and Thomma, 2013). LysM effector domains are enriched in Mucoromycota but with poor representation in Zoopagomycota with the exception of Basidiobolus ( Figure S6). This striking difference in the distribution of LysM effector domains positively correlates with Mucoromycota ecologies being primarily associated primarily with plants as symbionts (e.g., mycorrhizae and root endophytes) and saprotrophs, while Zoopagomycota are defined by animal and fungi-associated ecologies.
A parallel pattern of enrichment is observed for the CE4/PF015522 domain, which functions in cell-wall remodeling and the deacetylation of chitin to chitosan, consistent with chitosan being an abundant cell wall carbohydrate of Mucoromycota ( Figure S6). But this domain is also involved in protecting fungi from attack by plant hydrolases by de-N-acetylating chitin, and avoidance of the plant immune system (Toro and Brachmann, 2016).
Zoopagomycota with animal-associated ecologies have higher secretome:proteome ratios Secretome:proteome ratios of the zygomycete genomes in this study are comparable to those of Dikarya fungi (e.g., Lo Presti et al., 2015; Toro and Brachmann, 2016; Figure S2). Glomeromycotina stands out in having lower secretome:proteome ratios than the other zygomycete subphyla, consistent with the findings of previous studies (Kamel et al., 2017;Miyauchi et al., 2020;Tisserant et al., 2013). Zoopagomycota, especially the subphyla Zoopagomycotina and Entomophthoromycotina, has a higher relative secretome content compared to other fungi. Many members of Zoopagomycota are associated with animals, as either pathogens or symbionts. The high relative secretome content in these taxa is different from what has been reported for animal pathogens from Dikarya, which generally have low secretome:proteome ratios (Krijger et al., 2014;Lowe and Howlett, 2012). It should be noted that the majority of the Dikarya animal pathogens studied so far are mostly from Eurotiomycetes and are opportunistic pathogens of mammals, which only represents a small portion of the animal pathogen diversity in Dikarya.
Entomophthoromycosis fungi are best known as insect pathogens, such as Entomophthora muscae, a pathogen of the common house fly (Gryganskyi et al., 2017). These insect pathogens infect their hosts through spores, grow vegetatively within the host, ultimately cause harm and death, and then sporulate on the insect exoskeleton to infect other hosts. Many species within the subphylum exhibit complex life histories, however, which involve cycling between complex animal associations characterized as symbiotrophs (e.g., insect dispersed, herptile gut symbionts of Basidiobolus) and pathotrophs (e.g., insect pathogens of Conidiobolus), and saprotrophic life history stages associated with soil and plant litter environments. Furthermore, both Basidiobolus and Conidiobolus can function as opportunistic pathogens of humans and other animals (Vilela and Mendoza, 2018). In contrast, the sampled Zoopagomycotina fungi are mycoparasites, amoeba parasites, and predators of amoeba or rotifers. Although little is known about the cellular interactions between Zoopagomycota fungi and their animal or fungal hosts, the high relative abundance of the Zoopagomycota secretome is consistent with ecologies characterized by complex fungal-host interactions and nutrient acquisition strategies.

Secretome content of GPI-anchored proteins involved in maintaining cell integrity and interacting with the environment
A portion of the predicted secreted proteins are not released as free extracellular proteins, but are bound to the cell wall or plasma membrane. A proportion of these proteins are attached to the cell walls or plasma membrane by a glycerophospholipid anchor (GPI anchor) and extend to the extracellular side of the fungal cell. These proteins play important roles in retaining the integrity of fungal cells, in protecting fungi from biotic and abiotic environmental stresses, and are also related to pathogenicity (Garcia-Rubio et al., 2020;Gow et al., 2017;Vogt et al., 2020).
An average of 10% of the secretome is predicted to be GPI-anchored. Pfam annotation of the GPIanchored proteins is consistent with the aforementioned functions of the fungal cell wall proteins. The most common GPI-anchored proteins identified in this study, Ser-Thr-rich GPI-anchored proteins (PF10342), are known to function in cell wall organization (Samalova et al., 2020 iScience Article chitin deacetylases from PF01522 are also involved in the remodeling of the fungal cell walls (Toro and Brachmann, 2016). All these proteins are also involved in pathogenesis, and/or protection of pathogenic fungi from attack by their host plants (Samalova et al., 2017). In addition, we recovered several other pathogenesis-related protein families, including fasciclin proteins (PF02469) and cotH proteins (PF08757). Both proteins were described as surface proteins with functions in cell-to-cell interactions. Fasciclin proteins mediate interactions between cell exterior and cell surface and they function in pathogenesis and/or symbiosis in fungi (Seifert, 2018;Wang et al., 2018). CotH proteins play key roles in the invasion of human endothelial cells during mucormycosis by adhering the fungal cell walls to the endothelial cell membrane (Baldin and Ibrahim, 2017). Both of these pathogenesis-related proteins are present in many non-pathogenic fungi sampled in this study ( Figure S6). Importantly, most of the predicted fasciclin proteins (89%) and cotH proteins (75%) are not GPI-anchored, consistent with these proteins likely having functions other than pathogenesis and are consistent with a co-opted function in opportunistic pathogens.

The trade-off between small secreted proteins and digestive enzymes
SSPs and digestive enzymes displayed a strong, negative logarithmic relationship with respect to secretome content (Figure 1). Taxa with abundant digestive enzymes possessed relatively few SSPs, and viceversa. The highest SSP:DE ratios were predictably found in symbiotrophs of arbuscular mycorrhizal Glomeromycotina, consistent with previous observations and SSP importance in the mycorrhizal symbiosis (Miyauchi et al., 2020). Compared to the Dikarya, zygomycete genomes generally have higher proportions of SSPs and relatively lower proportions of digestive enzymes in their secretomes ( Figure S2). In addition, a high proportion of zygomycete SSPs is predicted to be putative effector proteins ( Figure S8), regardless of their subphylum, trophic mode, or ecological guild. Effector proteins are best characterized in pathogenic and symbiotic fungi, and function in the interaction and colonization of host cells (Lo Presti et al., 2015). Saprotrophic zygomycetes maintain a wide array of effector-like proteins, suggesting that either effector proteins play important roles in their survival in the environment away from host cells, or that their ecologies are more complex in nature than currently appreciated, or both. Together with previous reports on effectors shared by pathogenic fungi and their saprotrophic relatives (Baldin and Ibrahim, 2017;Samalova et al., 2017;Seifert, 2018;Wang et al., 2018), our study lends support to the hypothesis that effectors play dual functions in both pathogen and symbiont fitness and survival of saprotrophs (Kim et al., 2016;Morris et al., 2009;Stergiopoulos et al., 2012).

Diversification of small secreted proteins and divergent and convergent ecologies
SSPs, including effectors, function in symbiotrophic and pathotrophic ecologies through regulating biochemical and cellular processes associated with host defense mechanisms (Feldman et al., 2020). Our analyses revealed abundant SSPs and putative effectors for zygomycete fungi as well as three main patterns of SSP distribution. First, we found that related species (either the members of the same genus or family) exhibited strong divergence and shared only a limited set of SSP orthologous clusters. This finding is consistent with a conserved core function among closely related taxa that possess common ecologies. Second, species classified within a subphylum, but with different trophic modes and ecological guilds, do not share SSP OGs. Instead, each ecological group was characterized by a unique set of SSPs. For example, animal pathogens and soil saprotrophs of Kickxellomycotina, and herptile gut symbionts and insect pathogens of Entomophthoromycotina, possess divergent sets of SSPs consistent with their phylogeny (Figure 4). Third, convergent ecologies do not share orthologous SSPs. Instead, they possess a unique repertoire of SSP OGs, a finding consistent with non-orthologous solutions to common ecological guilds.

Concluding remarks
These data and analyses revealed that the secretome composition of zygomycete fungi is primarily influenced by phylogeny with trophic mode and ecological guild accounting for a smaller impact. The early evolution of terrestrial fungi secretome content is strongly influenced by adaptation to plant versus non-plant hosts and substrates. This diversification is characterized by ancient losses of enzymatic diversity in the MRCAs and more derived lineage-specific expansions, that resulted in the ll OPEN ACCESS iScience 25, 104840, August 19, 2022 9 iScience Article diversification of PCWDE in Mucoromycotina, chitinases, and peptidases in Mortierellomycotina, and independent expansions of chitinases and proteases within Zoopagomycota. Regardless of ecology, there is a trade-off between SSPs and digestive enzymes, suggesting that either the ecologies of these organisms are incompletely known, or that SSPs are involved in important functions associated with saprotrophic ecologies, or both.

Limitations of the study
With the most extensive genome sampling on zygomycetes up to date, this work systematically investigated how genealogy, genome/secretome composition, and fungal lifestyle were interconnected with each other. Although we took similarity-based approaches to annotate the secretome, we noted there were potential disparities between sequence similarity and function similarity for some proteins. In addition, just like any other large-scale studies, the application of the same threshold for different proteins from different fungal lineages (e.g., the same e-value cut-offs for all peptidases) might also iScience Article introduce inaccuracies in the secretome annotation. Although the life-style annotation of the sampled fungi was based on the authors' best knowledge and cross-validated by mycologists with expertise in different fungal groups, what fungi do in the natural environment is challenging to study. Hence the life-style annotation included in this study is likely a simplified and partial reflection of their real ecological functions.

STAR+METHODS
Detailed methods are provided in the online version of this paper and include the following:

Material availability
This study did not generate new unique reagents.

Data and code availability
The accession numbers for the genome assemblies are listed in Table S1.
This paper does not report original code.
Any additional information of the data is available from the lead contact upon request.

EXPERIMENTAL MODELS AND SUBJECT DETAILS
No experimental model organism was used in this study.

METHOD DETAILS
Taxon sampling, trophic mode coding, and genome sequencing, assembly and annotation We sampled 182 fungal genomes, with representative species from all the subphyla in Kingdom Fungi, including 98 Mucoromycota and 34 Zoopagomycota genomes (Table S1). We included three non-fungal Opisthokont genomes (Monosiga brevicollis, Capsaspora owczarzki and Drosophila melanogaster) as the outgroups for phylogenetic reconstruction.
We coded the sampled zygomycete fungi according to their trophic mode and ecological guild (Tables S2  and S3). The trophic mode categories included symbiotroph, pathotroph, and saprotroph, while the ecological guilds included amoeba parasite, animal pathogen, animal symbiont, arbuscular mycorrhiza, dung saprotroph, ectomycorrhiza, endophyte, mycoparasite, plant pathogen, saprotroph, & 'unknown' (1). We adopted two coding schemes -1) Each genome coded by its primary trophic mode or ecological guild; 2) a binary coding method with '0' for absence and '1' for presence of a particular ecological guild. In this scheme, we allowed a fungus to have multiple ecological guilds when evidence suggests this possibility.
We newly sequenced 29 genomes in this work (Table S1). The genome of Podila verticillata (=Mortierella verticillata) NRRL 6337 was sequenced by the Broad Institute using Sanger and Illumina technology. Three libraries were constructed for Sanger sequencing, 4 and 10 kb plasmids and 40 kb Fosmids, and paired-end reads were generated. Illumina sequencing was performed on an Genome Analyzer II to generate 101 base paired reads. Reads were assembled using the Arachne Hybrid Assembler v. Nov. 2010 (Jaffe et al., 2003). The assembly was annotated using the Broad Institute fungal pipeline .
The remaining 28 genomes were sequenced by the Joint Genome Institution. Out of these 28 genomes, the genomes of Cunninghamella echinulata NRRL 1382 and Dissophora ornata CBS 347.77 were sequencing using Illumina platform, for which 100 ng of DNA was sheared to 300 bp using the Covaris LE220 and size selected using SPRI beads (Beckman Coulter). The fragments were treated with end-repair, A-tailing, and ligation of Illumina compatible adapters (IDT, Inc) using the KAPA-Illumina library creation kit (KAPA biosystems). qPCR was used to determine the concentration of the libraries and were sequenced on the Illumina Hiseq. The prepared libraries were quantified using KAPA Biosystems' next-generation sequencing library qPCR kit and run on a Roche LightCycler 480 real-time PCR instrument. The quantified libraries were then prepared for sequencing on the Illumina HiSeq sequencing platform utilizing a TruSeq paired-end cluster kit, v4. Sequencing of the flowcell was performed on the Illumina HiSeq2500 sequencer using HiSeq TruSeq SBS sequencing kits, v4, following a 2 3 150 indexed run recipe.
All other genomes were sequenced using the PacBio platform, for which, 1-6 ug of genomic DNA was sheared to 10kb or >10kb using Covaris g-Tubes. The sheared DNA was treated with exonuclease to remove single-stranded ends and DNA damage repair mix followed by end repair and ligation of blunt adapters using SMRTbell Template Prep Kit 1.0 (Pacific Biosciences). All libraries were then purified with AMPure PB beads. After libraries were prepared, PacBio Sequencing primer was then annealed to the SMRTbell template library and Version P6 sequencing polymerase was bound to them. Then, the prepared SMRTbell template libraries were then sequenced on a Pacific Biosciences RSII sequencer using Version C4 chemistry and 1 3 240 sequencing movie run times. However, for C. cucurbitarum, H. radiatus, H. pulchrum and C. mojavensis, sequencing polymerase was bound to them using Sequel Binding kit 2.0 and then the prepared SMRTbell template libraries were sequenced on a Pacific Biosystems' Sequel sequencer using v3 sequencing primer, 1M v2 SMRT cells, and Version 2.1 sequencing chemistry with 1 3 360 & 1 3 600 sequencing movie run times. For PacBio-sequenced genomes, a variety of methods were used for assembly as new tools became available. All genomes were assembled using Falcon (Chin et al., 2016). FinisherSC (Lam et al., 2015), was then used in most cases to improve the assembly, followed by polishing with either Quiver (Chin et al., 2013). or Arrow (https://github.com/PacificBiosciences/GenomicConsensus). See Table S1 for details on tools used for assembly, finishing and polishing.
All transcriptomes were sequenced with Illumina. For all lineages except D. ornata CBS 347.77 and H. pulchrum RSA 2064, stranded cDNA libraries were generated using the Illumina Truseq Stranded mRNA Library Prep kit. mRNA was purified from 1 ug of total RNA using magnetic beads containing poly-T oligos. mRNA was fragmented and reversed transcribed using random hexamers and SSII (Invitrogen) followed by second strand synthesis. The fragmented cDNA was treated with end-pair, A-tailing, adapter ligation, and 8 cycles of PCR. The prepared libraries were quantified using KAPA Biosystems' next-generation sequencing library qPCR kit and run on a Roche LightCycler 480 real-time PCR instrument. The quantified libraries were then prepared for sequencing on the Illumina HiSeq sequencing platform utilizing a TruSeq paired-end cluster kit, v4. Sequencing of the flowcell was performed on the Illumina HiSeq2000 sequencer using HiSeq TruSeq SBS sequencing kits, v4, following a 2 3 150 indexed run recipe.
For D. ornata CBS 347.77 and H. pulchrum RSA 2064, plate-based RNA sample prep was performed on the PerkinElmer Sciclone NGS robotic liquid handling system using Illumina's TruSeq Stranded mRNA HT sample prep kit utilizing poly-A selection of mRNA following the protocol outlined by Illumina in their user guide (https://support.illumina.com/sequencing/sequencing_kits/truseq-stranded-mrna.html). For D. ornata, the quantified libraries were sequenced as described earlier but using HiSeq2500. For H. pulchrum, sequencing of the flowcell was performed on the Illumina NovaSeq sequencer using NovaSeq XP V1 reagent kits, S4 flowcell, following a 2 3 150 indexed run recipe.
All reads were filtered and trimmed for quality. BBDuk (https://sourceforge.net/projects/bbmap/) was used to remove artifact sequences by kmer matching (kmer = 25), allowing 1 mismatch. Detected artifacts were trimmed from the 3 0 end of the reads. RNA spike-in reads, PhiX reads and reads containing any Ns were removed. Quality trimming was performed using the phred trimming method set at Q6. Finally, following trimming, reads under the length threshold were removed (minimum length 25 bases or 1/3 of the original read length -whichever is longer). Filtered reads were assembled de novo using Trinity (Grabherr et al., 2011) version 2.1.1.1 or 2.3.2. The assemblies and annotations of newly sequenced genomes are available at the fungal genome portal MycoCosm (https://mycocosm.jgi.doe.gov) (Grigoriev et al., 2014) and in the DDBJ/EMBL/GenBank repository (https://www.ncbi.nlm.nih.gov; Table S1).

Phylogenetic reconstruction
The phylogenetic reconstruction was based on a dataset comprising 434 protein markers that have previously proved useful for higher level phylogenetic analysis of fungi (Beaudet et al., 2018;Chang et al., 2019) (https://github.com/1KFG/Phylogenomics_HMMs) (12). The retrieval and subsequent alignment of the markers were performed using the pipeline PHYLING (http://github.com/stajichlab/PHYling_unified) (Stajich, 2020) with the default settings. Phylogenetic reconstruction was performed on the concatenated ll OPEN ACCESS iScience Article

Reconstruction of digestive enzyme evolution
The sequences for each digestive enzyme family were aligned using MAFFT 7.402 (Katoh and Standley, 2013) (L-INS-i, maxiterate = 1000). Individual enzyme trees were inferred using RAXML v.8.0.26 (Stamatakis, 2014) (-f a option, 100 bootstrap replicates and a PROTGAMMALG model). The evolution of enzyme families was reconstructed by reconciling the enzyme trees and the best RAxML species tree using NOTUNG 2.9 (Stolzer et al., 2012). To take into account the uncertainty in the enzyme trees, we allowed the branches with <70% ML support to be rearranged to obtain the most parsimonious solution during the reconciliation.

ll
OPEN ACCESS