2,341
Views
0
CrossRef citations to date
0
Altmetric
Research Paper

Plant pathogens provide clues to the potential origin of bat white-nose syndrome Pseudogymnoascus destructans

ORCID Icon, ORCID Icon, , ORCID Icon & ORCID Icon
Pages 1020-1031 | Received 04 Feb 2022, Accepted 19 May 2022, Published online: 06 Jun 2022

ABSTRACT

White-nose syndrome has killed millions of bats, yet both the origins and infection strategy of the causative fungus, Pseudogymnoascus destructans, remain elusive. We provide evidence for a novel hypothesis that P. destructans emerged from plant-associated fungi and retained invasion strategies affiliated with fungal pathogens of plants. We demonstrate that P. destructans invades bat skin in successive biotrophic and necrotrophic stages (hemibiotrophic infection), a mechanism previously only described in plant fungal pathogens. Further, the convergence of hyphae at hair follicles suggests nutrient tropism. Tropism, biotrophy, and necrotrophy are often associated with structures termed appressoria in plant fungal pathogens; the penetrating hyphae produced by P. destructans resemble appressoria. Finally, we conducted a phylogenomic analysis of a taxonomically diverse collection of fungi. Despite gaps in genetic sampling of prehistoric and contemporary fungal species, we estimate an 88% probability the ancestral state of the clade containing P. destructans was a plant-associated fungus.

Introduction

Battling an emerging pathogen requires a thorough understanding of both its origin and the mechanisms by which it attacks the host. Often, these factors are overlooked in the early stages of research on an emerging pathogen as efforts are concentrated on disease surveillance, containment, and possible treatments. This pattern is especially apparent in the study of novel wildlife diseases where resources are limited, and where we lack an understanding of the basic biology of the pathogen and sometimes even the host. Perhaps no wildlife disease better exemplifies this pattern than white-nose syndrome (WNS), a disease with devastating effects on populations of cave-hibernating bats in North America [Citation1]. The causative agent of WNS, Pseudogymnoascus (Geomyces) destructans [Citation2,Citation3], was first isolated over a decade ago, but significant questions still remain about the ultimate origin of this fungus and the mechanism by which it invades the skin of hibernating bats.

The colloquial name “white-nose syndrome” () [Citation6] belies the primary site of infection for P. destructans: the skin of the wing membrane (). A visible infection of the muzzle with P. destructans is often used as an identifying field characteristic but infection of the muzzle can be absent in the presence of wing infection. An accurate diagnosis of WNS requires the identification of its novel pathognomonic skin lesions, which are best appreciated on the relatively hairless patagium of the wings. These unique and characteristic ‘cupping erosions’ progress to a potentially fatal, invasive, cutaneous ascomycosis as WNS progresses [Citation5].

Figure 1. Cave hibernating Myotis lucifugus with white-nose syndrome caused by Pseudogymnoascus destructions. (a) White patches seen grossly on the muzzle, ear, and wing membrane (white arrows) are aerial hyphae with conidia. Photograph used with permission from Gregory Turner, Pennsylvania game commission, USA. (b) Dark areas of contraction (white arrow) are due to the infection of wing membrane with P. destructans. The black arrow points to a region of more normal tissue of the wing membrane. (c) A bat wing illuminated by 385-nm ultraviolet light. The yellow fluorescence corresponds to intraepidermal aggregates of hyphae [Citation4]. (d) Periodic acid Schiff-stained section of wing membrane from a bat with white-nose syndrome. Packets of intraepidermal, magenta-stained hyphae of P. destructans (arrows), previously termed ‘cupping erosions’ [Citation5], cause the fluorescence in (c). Absence of inflammatory cell reaction is consistent with immune downregulation that occurs during hibernation.

Figure 1. Cave hibernating Myotis lucifugus with white-nose syndrome caused by Pseudogymnoascus destructions. (a) White patches seen grossly on the muzzle, ear, and wing membrane (white arrows) are aerial hyphae with conidia. Photograph used with permission from Gregory Turner, Pennsylvania game commission, USA. (b) Dark areas of contraction (white arrow) are due to the infection of wing membrane with P. destructans. The black arrow points to a region of more normal tissue of the wing membrane. (c) A bat wing illuminated by 385-nm ultraviolet light. The yellow fluorescence corresponds to intraepidermal aggregates of hyphae [Citation4]. (d) Periodic acid Schiff-stained section of wing membrane from a bat with white-nose syndrome. Packets of intraepidermal, magenta-stained hyphae of P. destructans (arrows), previously termed ‘cupping erosions’ [Citation5], cause the fluorescence in (c). Absence of inflammatory cell reaction is consistent with immune downregulation that occurs during hibernation.

The invasive nature of P. destructans contrasts with the behavior of typical dermatophytes of mammals (e.g. Trichophyton, Microsporum, Epidermophyton spp.), which are superficial parasites, more annoyance than pathogen (e.g. ringworm and athlete’s foot). Normal microflora of human skin, like the yeast Malassezia, can infect sebaceous glands and cause superficial infections under certain conditions. Likewise, if the skin surface is otherwise damaged, soil-associated fungi such as Madurella spp. and Pseudogymnoascus (Geomyces) pannorum can secondarily infect these penetrating wounds and “digest” animal tissue. Other fungal infections occur in hibernating bats that can cause patches of white fungal growth on skin that grossly resemble WNS, but histologically resemble typical dermatophytes with infection limited to the keratin layer of the skin without invasion [Citation7]. The important distinction is that, on their own, these dermatophytes cannot invade an intact skin surface. Indeed, most fungi that can cause pathology in animals exist primarily as saprobes and are opportunistic pathogens that do not require an animal host [Citation8]. P. destructans is unique among known mammal-associated cutaneous fungi in its ability to actively penetrate the normal, intact skin of its host.

Several other characteristics distinguish P. destructans from typical dermatophytes. P. destructans is an obligate psychrophile that inhabits the consistently cold hibernacula where bats overwinter [Citation9], while other dermatophytes tolerate the consistently high body temperatures of homeothermic mammals. Further, dermatophytes are nourished by non-living keratin on skin surface, which they access through secretion of keratinase [Citation10]. P. destructans has not been shown to hydrolyze keratin [Citation11,Citation12]. Instead, it produces lipases, proteinases, endopeptidases, collagenase, esterase, hemolysins, urease, and chitinase [Citation11,Citation13], pointing to a different nutrition strategy because these enzymes hydrolyze fat, protein, collagen, etc. Finally, it appears that P. destructans undergoes morphological shifts not seen in typical mammal dermatophytes during invasion. When found on the surface, hyphae of P. destructans are delicate with parallel walls and are non-fluorescent when exposed to longwave ultraviolet (UV) light. When found in the diagnostically characteristic cupping erosions, the hyphae change morphology and are wide and bulbous, with irregular non-parallel walls [Citation5]. Here, the distinct bundles of hyphae fluoresce under UV light () [Citation4]. The change in fluorescence is considered evidence of a physiological shift [Citation14], which is seen in fungal plant pathogens as invasion of their plant host progresses through various stages [Citation15]. Although P. destructans has been considered an unusual dermatophyte that was recently introduced to North America with devastating effect on hibernating bats, previous research has not made an explicit attempt to address the unusual pattern of skin invasion or the origin of P. destructans. Evidence would indicate it is unlikely P. destructans evolved from a typical mammalian dermatophyte.

Proteins associated with P. destructans sampled under different environmental conditions [Citation16] provide evidence that P. destructans may have originated from an environmental saprophyte. However, although it can survive in sediments when bats are not present, P. destructans is outcompeted by other fungi commonly found in hibernacula [Citation17,Citation18]. This is consistent with findings that P. destructans is less efficient than sympatric fungi at metabolizing plant derived cellulose substrates and remains of other microbes and invertebrates, including chitin, that are a primary nutrient sources for saprophytes [Citation13]. In addition, a comparative genome study reported a recent duplication of the high affinity nitrate transporter NRT2 in P. destructans that may reflect a more recent enhancement of P. destructans’ ability to use environmental nitrogen [Citation19]. These findings suggest that P. destructans is unlikely to be a primary saprophyte acting as an opportunistic pathogen, but rather a primary pathogen and an opportunistic saprophyte. However, none of this evidence rules out contributions from an ancestral environmental fungus, particularly an obligate cryophile, in the evolution of P. destructans.

Previous efforts to characterize P. destructans have left the question of its origins unanswered. Over the course of examining more than a hundred WNS-positive bats, we first recognized an unusual mechanism of skin invasion by P. destructans, previously undescribed in mammalian pathogens, and further recognized that invasion mechanisms might provide clues about the ultimate origin of P. destructans (C. Meteyer, U.S. Geological Survey, personal observation, [2008–2012]). Here, we detail these histopathological findings and add transmission electron microscopy and phylogenomic analyses to more explicitly address both the invasion mechanism and ultimate origins of P. destructans. These multiple lines of independent evidence point to plant-associated fungal pathogens as ancestors of P. destructans.

Results and Discussion

To investigate mechanisms of bat skin invasion by P. destructans, we analyzed microscopic details of the hyphae as they invaded wing membranes and muzzles in naturally infected bats. Our microscopic and electron microscopic studies of wing tissue elucidate two distinct stages of infection as P. destructans invades the skin of bats. The first stage of infection forms the epidermal bundles of irregular hyphae, previously termed “cupping erosions” (), that are so unique among fungal skin infections they continue to be the gold standard for diagnosing WNS, even with the development of molecular diagnostic tests [Citation5,Citation20]. Importantly, our studies reveal the cupping erosions are not actual erosions in the skin surface, but proliferation of hyphae within the intraepidermal space, confined above and below by epidermal cells without associated damage to the epidermis (). This initial stage of P. destructans infection is morphologically similar to the biotrophic strategy of infection used by fungal pathogens of plants [Citation22]. Biotrophic infection, the propagation of an invading organism within living host tissues without causing tissue damage, has been extensively studied in plant fungal pathogens; it has not been detailed as a strategy of fungal infections in animals.

Figure 2. Wing membrane. Histologic periodic acid Schiff stain (a, b, d) and transmission electron micrographic (TEM, c) sections of wing membrane from Myotis lucifugus infected with Pseudogymnoascus destructans causing white-nose syndrome. Fungal hyphae are stained magenta. Lack of inflammatory cell reaction is consistent with WNS in hibernating bats [Citation5]. (a) Pigmented cells of the epidermal layer (arrowheads) illustrate the intraepidermal boundary above and below dense packets of irregular P. destructans hyphae without epidermal damage [Citation21]; this is consistent with the biotrophic stage of invasion. (b) Intraepidermal proliferation of fungal hyphae demarcated above and below by pigmented epidermis (arrowheads). Arrows show hyphae penetrating the basement membrane of the epidermis initiating infection of the deeper dermis. This hyphal invasion resembles the penetration mechanism used by some types of appressoria. (c) TEM shows fungal hyphae (arrows) within the outer layer (stratum corneum) of the epidermis. Arrowheads indicate the layer of bat epidermal cells covering hyphae without morphologic change or evidence of damage to the host cells; this is consistent with biotrophic infection. (d) Progression of P. destructans infection into deeper dermis (arrows) with replacement of connective tissue of the wing membrane with fungal hyphae.

Figure 2. Wing membrane. Histologic periodic acid Schiff stain (a, b, d) and transmission electron micrographic (TEM, c) sections of wing membrane from Myotis lucifugus infected with Pseudogymnoascus destructans causing white-nose syndrome. Fungal hyphae are stained magenta. Lack of inflammatory cell reaction is consistent with WNS in hibernating bats [Citation5]. (a) Pigmented cells of the epidermal layer (arrowheads) illustrate the intraepidermal boundary above and below dense packets of irregular P. destructans hyphae without epidermal damage [Citation21]; this is consistent with the biotrophic stage of invasion. (b) Intraepidermal proliferation of fungal hyphae demarcated above and below by pigmented epidermis (arrowheads). Arrows show hyphae penetrating the basement membrane of the epidermis initiating infection of the deeper dermis. This hyphal invasion resembles the penetration mechanism used by some types of appressoria. (c) TEM shows fungal hyphae (arrows) within the outer layer (stratum corneum) of the epidermis. Arrowheads indicate the layer of bat epidermal cells covering hyphae without morphologic change or evidence of damage to the host cells; this is consistent with biotrophic infection. (d) Progression of P. destructans infection into deeper dermis (arrows) with replacement of connective tissue of the wing membrane with fungal hyphae.

As infection progresses, P. destructans invades beyond the epidermis into deeper layers of skin. This second stage of infection is initiated by hyphae that penetrate the basement membrane that separates the epidermis and dermis (). After penetration, hyphae proliferate in the dermis, hydrolyzing and replacing host tissue (). This stage of infection is similar to the necrotrophic strategy of plant fungal pathogens. Hemibiotrophic species, those that switch from biotrophic infection to necrotrophic infection, are well studied among plant pathogens [Citation23] but are unknown among fungal pathogens of vertebrates, including mammals. Invading hyphae seen in histologic sections of bats infected with P. destructans () resemble appressoria, which have also been well described in plant fungal pathogens. Appressoria represent a category of structures used for host invasion. They vary in morphology and mechanism, but have previously not been documented in fungal pathogens of mammals [Citation24].

Another well-known wildlife disease that has contributed to population declines, chytridiomycosis, serves as an informative contrast to demonstrate the uniqueness of P. destructans’ invasion strategy among vertebrate fungal pathogens. As P. destructans invades deeper layers of skin as it spreads along the wing membrane in hibernating bats, the disruption of the epidermis disturbs thermoregulation, water and electrolyte balance, and respiration, potentially leading to death [Citation25,Citation26]. The causative fungal pathogen of chytridiomycosis, Batrachochytrium dendrobatidis, also compromises the critical respiratory function of the skin in many species of frogs, but remains in the epidermis and does not invade the deeper dermis. This primitive, single-celled, non-filamentous fungus uses keratin as its nutritional source, encysts on the surface of the skin and develops a germ tube that penetrates a host cell to generate immature sporangia [Citation27]. As a surface infection, B. dendrobatidis survives on the superficial non-living keratin layer of the host which thickens in response to infection. As the keratin layer thickens (hyperkeratosis), the osmoregulation and respiratory functions of the amphibian skin is disrupted and can eventually lead to death [Citation27]. Thus, even though P. destructans and B. dendrobatidis lead to death through disruption of the physiology of the skin in their respective hosts, the infection strategies and the resulting pathology are very different.

Our microscopic and electron microscopic studies also provide evidence that the hyphae of P. destructans seem to detect and respond to nutritional sources similar to nutrient tropism in fungal pathogens of plants and insects. In WNS, P. destructans produces hyphae and conidia that cover the muzzle and create the titular “white nose”. The haired skin of the muzzle is endowed with abundant sebaceous glands [Citation28], which produce simple and complex lipids secreted through the sebaceous pore into the hair follicle (). In addition to intraepidermal invasion of the muzzle, similar to what occurs in tissues of the wing membrane, the hyphae of P. destructans also invade and replace sebaceous glands as the enzymes it produces (e.g. lipase, esterase, and proteinase) hydrolyze lipids, fatty acids, cholesterol esters, and triglycerides produced by these glands (). Many fungal pathogens of plants and insects have evolved nutrient sensing and signaling pathways [Citation29]. When these pathways are activated they induce growth of hyphae toward the nutrient source (tropism) and initiate the differentiation of specialized structures to access nutrients, including pseudohyphae, hyphal extensions, and appressoria [Citation29]. The term appressoria encompasses a wide range of morphologies and invasion mechanisms [Citation30–34]. Although appressoria or appressoria-like structures are widely described as crucial infection structures of some fungal pathogens of plants and insects [Citation35,Citation36], they have not been identified in vertebrate pathogens [Citation24].

Figure 3. Muzzle. Histologic sections of muzzle from Myotis lucifugus stained with hematoxylin and eosin (a), periodic acid Schiff stain (b, d), and transmission electron micrographic (TEM, c). (a) Normal muzzle tissue stained with hematoxylin and eosin. Arrowhead points to hair shaft. Black arrow points to sebaceous gland that secretes lipid material into the hair follicle. (b) Fungal hyphae of Pseudogymnoascus destructans (magenta) within and expanding from the follicle (arrowheads), invading and replacing lipid-laden sebaceous gland (arrow) and surrounding connective tissue. (c) TEM shows fungal hyphae (black arrows) entering the hair follicle identified by emerging hair shaft (white arrow). (d) Extensive invasion of P. destructans including hair follicle surrounding hair shaft (arrowheads) with replacement of associated sebaceous glands (arrows) and regional connective tissue.

Figure 3. Muzzle. Histologic sections of muzzle from Myotis lucifugus stained with hematoxylin and eosin (a), periodic acid Schiff stain (b, d), and transmission electron micrographic (TEM, c). (a) Normal muzzle tissue stained with hematoxylin and eosin. Arrowhead points to hair shaft. Black arrow points to sebaceous gland that secretes lipid material into the hair follicle. (b) Fungal hyphae of Pseudogymnoascus destructans (magenta) within and expanding from the follicle (arrowheads), invading and replacing lipid-laden sebaceous gland (arrow) and surrounding connective tissue. (c) TEM shows fungal hyphae (black arrows) entering the hair follicle identified by emerging hair shaft (white arrow). (d) Extensive invasion of P. destructans including hair follicle surrounding hair shaft (arrowheads) with replacement of associated sebaceous glands (arrows) and regional connective tissue.

Our transmission electron microscopic images of bat muzzles infected with P. destructans point to a similar tropism-induced growth of P. destructans hyphae toward a lipid nutrient source (sebaceous gland) as they are consistently seen at the openings of hair follicles () and within sebaceous glands in bats with “white noses”. In this nutrient rich environment, hyphae proliferate, hydrolyzing and obliterating regional tissue (). In addition to sebaceous glands, the epidermis covering the muzzle and wing membrane is also rich in lipids that bathe surface corneocytes in a complex extracellular lipid matrix [Citation37]. We propose that detection of lipids and other nutrients available in the intraepidermal spaces of bat skin and in sebaceous glands activates a similar nutrient signaling pathway in P. destructans that might lead to the formation of appressoria. Our transmission electron microscopic images were limited and did not allow us to conclusively characterize appressoria, but invading hyphae seen in histologic sections () are consistent with appressoria. Further in vitro and in vivo infection studies would be valuable to determine whether or not appressoria are induced in P. destructans to invade host tissues. A transcriptome study of P. destructans from infected bat wing membrane identified genes encoding enzymes involved in the remodeling of fungal cell walls [Citation38]. Cell wall remodeling is not specific to appressoria formation but it is necessary for their development [Citation39]. If confirmed, the development of appressoria or appressoria-like structures would provide additional evidence for the use of plant pathogen invasion strategies in a mammalian host.

To investigate the origin of P. destructans more thoroughly in light of its unique characteristics, we conducted a high resolution phylogenomic analyses of a broad inventory of ascomycete species. As previously shown in a more taxonomically focused phylogenomic reconstruction [Citation40], Pseudogymnoascus forms a well-supported sister group to the order Helotiales within the class Leotiomycetes, but with a long branch illustrating the relatively poor sampling of this lineage. To expand on this analysis, we broadened taxonomic coverage, making use of publicly available sequence data, and further collected information about the “lifestyle” associated with taxa for which information was available. The resulting phylogeny demonstrates the uniqueness of P. destructans in terms of host range among closely related species. The class Leotiomycetes is highly enriched with species known to be associated with plants in different ways: as pathogens, endophytes, or saprotrophs (e.g. Coleophoma spp., Cudoniella spp. and Lachnellula spp.) [Citation41,Citation42] (). For example, among the genera in the clade with P. destructans is Botrytis spp., which contain important plant pathogens [Citation43]. Using phylogenetic models of character evolution, we quantified the probability that the ancestral state of the clade containing P. destructans was associated with plants (as a pathogen or not) to be at least 88% (Table S-2).

Figure 4. Pseudogymnoascus species are phylogenetically related to plant-associated fungi. Maximum likelihood tree based on 1,146 families of orthologues. All displayed nodes are supported by a bootstrap value >95%. The number of sequences in each collapsed clade are indicated within parentheses. The tree is represented in its more complete form inFigure S-1. Within the Leotiomyceta clade, groups outside Pseudogymnoascus have been collapsed at the species level, with the corresponding number of strains indicated in parentheses. Color dots indicate known ecologies. The Leotiomyceta clade does not contain any representatives of the ecology groups for aquatic, lichen, mycoparasite of plant pathogen, or mammalian pathogens outside of the Pseudogymnoascus genus. Additional species have been discarded because they were closely related to species with more sequence data, a list of which can be found in Supplementary Table S1.

Figure 4. Pseudogymnoascus species are phylogenetically related to plant-associated fungi. Maximum likelihood tree based on 1,146 families of orthologues. All displayed nodes are supported by a bootstrap value >95%. The number of sequences in each collapsed clade are indicated within parentheses. The tree is represented in its more complete form inFigure S-1. Within the Leotiomyceta clade, groups outside Pseudogymnoascus have been collapsed at the species level, with the corresponding number of strains indicated in parentheses. Color dots indicate known ecologies. The Leotiomyceta clade does not contain any representatives of the ecology groups for aquatic, lichen, mycoparasite of plant pathogen, or mammalian pathogens outside of the Pseudogymnoascus genus. Additional species have been discarded because they were closely related to species with more sequence data, a list of which can be found in Supplementary Table S1.

Collectively, our morphologic and phylogenomic analyses provide evidence for our hypothesis that P. destructans evolved from an ancestor that was adapted to coexist not with mammals but with plants, with invasion strategies that resemble hemibiotrophic plant pathogens. In a different context, the similarity of P. destructans to plant pathogens was demonstrated by Pannkuk and colleagues [Citation44] who investigated an extracellular subtilisin-like serine protease produced by P. destructans. Although subtilisin serine proteases are produced by typical dermatophytes and promote skin adhesion during infection [Citation10], comparative amino-acid sequence analyses showed that the serine protease produced by P. destructans has a higher similarity to serine proteases of diverse fungal plant pathogens including Sclerotinia and Botrytis than dermatophyte pathogens like Arthroderma, Trichophyton, and Microsporum [Citation44].

Although P. pannorum, a widespread soil saprophyte, is considered a close contemporary to P. destructans, it was not included in our final phylogenomic analysis because of insufficient data (see methods for details). The comparison between these two fungi is interesting because of the close (presumed) relationship but pronounced physiological differences. P. pannorum is broadly thermotolerant allowing it to infect homeothermic mammals, including humans [Citation45]. It cannot, however, independently invade animal tissue, instead requiring preexisting damage to the skin for infection to occur. Conversely, P. destructans uses a novel approach that allows it to invade the intact skin of mammals while being restricted in its host range as an obligate psychrophile. Comparisons of these congeners may lead to a better understanding of the evolutionary history of the clade as a whole and may offer clues about the unique confluence of conditions that could facilitate a shift from plant-associated fungus to a cryophilic mammalian pathogen.

We can currently only speculate about what might have transpired millions of years ago that allowed P. destructans to adapt as a pathogen of hibernating bats using strategies that resemble those of plant pathogens, but there are hints in the literature as to where such an investigation might begin. A previous study comparing the evolutionary genomics of Eurasian bats and P. destructans suggests their ancestors could have temporally overlapped with the potential to co-evolve during the Cenozoic [Citation40]. During this era, Eurasia (and indeed most of the planet) would have been in the middle of the Cenozoic cooling, characterized by decreasing annual temperatures and increasing seasonality [Citation46,Citation47]. Such a climate would necessitate the evolution of strategies such as hibernation to allow mammals to survive cold periods when food was not available. A gradual shift to hibernation across geological time [Citation48] might have provided plant-associated ancestors of P. destructans an opportunity to intermittently and progressively adapt to a new host that was becoming thermally compatible for longer and longer periods as the body temperature of hibernating bats dropped within the ideal temperature range for P. destructans [Citation9]. Further, finding refuge in consistently cool caves to hibernate and avoid freezing could have brought bats in proximity to a cryophilic fungus like P. destructans. There are other possibilities, including the jump to hibernating bats through an intermediate host with labile body temperature, like seen in some extant tropical and subtropical mammals to this day [Citation49]. Still, the simplest explanation, and one warranting further exploration is the co-evolution of bats and P. destructans as bats adapted from mild to progressively more severe winters. Initially intermittent infections of short duration might have provided evolutionary pressure over a geological timescale for ancestors of P. destructans to progressively adapt virulence strategies to take advantage of hibernating bats as a new host niche. Initial infections of short duration during bouts of torpor might also have allowed ancestral bats to develop tolerance to P. destructans. Comparing immune effector proteins from populations exposed to P. destructans for different lengths of time [Citation50,Citation51] provides additional contemporary evidence for adaptation/tolerance and co-evolution.

This evolutionary scenario also provides a hypothetical explanation for the mild infections seen in Eurasian bats infected with P. destructans when compared to the severity of infections and death experienced by North American bats [Citation52]. Genetic diversity is higher among Eurasian isolates of P. destructans than North American isolates showing a longer evolutionary history in Eurasia and recent introduction with clonal expansion in North America [Citation3,Citation53,Citation54]. North American P. destructans isolates are more closely related to European isolates, whereas the genomes of isolates from these geographic regions are distinct from those isolated in China, supporting the introduction of P. destructans from Europe to North America [Citation53].

Geographical isolation of North American bats, and their recent ancestors, from P. destructans was breached with a recent transcontinental introduction of this unique pathogen [Citation53]. With patterns of host skin infection that resemble those of a plant pathogen, and recent introduction into North America, it would be expected that North American bats would have none of the tolerance found in Eurasian bats [Citation55–57]. Conversely, P. destructans would be perfectly adapted to infect and cause severe disease in North American bats. The result has been a predictable and rapid spread of P. destructans in hibernating bats across the continent, and unprecedented mortality associated with WNS [Citation1]. A similar scenario might have played out at least once before as well because mass accumulations of bat skeletal remains in Europe suggest a significant mortality event during the Pleistocene and Pliocene (5.4 million years ago − 11,700 years ago) that might have been associated with the introduction of a novel pathogen, proposed to be the initial introduction of P. destructans into Europe [Citation52].

Our study used multiple lines of evidence to provide a novel hypothesis about the origin of P. destructans. Specifically, we propose that P. destructans is adapted from plant-associated fungi, with invasion strategies compatible with ancestral contributions from plant pathogens. Importantly, the plant-pathogen hypothesis is testable. For example, transcriptome analysis can characterize significant physiological stages of invasion such as the gene transcriptome that functions as the “switch” from a biotrophic to necrotrophic infection [Citation23]. Similarly, in vitro systems exist for real-time imaging of invasion by fungal plant pathogens, including fluorescent contrast studies to track appressoria [Citation39]. Similar imaging and transcription studies could be applied to P. destructans.

We have proposed a novel, testable hypothesis that could act to bridge research between scientists studying the pathogenesis of plant fungi and scientists studying the pathogenesis of animal fungi. Research on WNS started from nothing about 15 years ago, but plant fungal pathogens have been the subject of intense research for decades. A trove of knowledge has been accumulated from plant-based research that might be applicable to WNS in bats. Taking advantage of this accumulated knowledge and recruiting scientists from the field of fungal pathogenesis in plants can accelerate the understanding of P. destructans and its novel host-pathogen relationship while potentially providing new paths to prevent WNS. Applying standard tools used in the study of plant fungal pathogens to P. destructans also has the potential to add a new model system to better understand mechanisms virulent fungal pathogens use to invade plant hosts. And finally, bridging research between scientists working with plant and animal pathogens will almost certainly lead to new approaches and insights into how new pathogens emerge from physical and biological environments and adapt to new hosts, including those in distant kingdoms.

Materials and methods

Microscopy

Histopathology that contributed to this paper included wing samples taken after full necropsy and examination of wings of 168 bats (47% positive for intraepidermal fungal proliferation without necrosis or inflammation – the hallmark of WNS) of 11 species (60% Myotis lucifugus, 16% Perimyotis subflavus) submitted to the U.S. Geological Survey (USGS) National Wildlife Health Center in Madison, Wisconsin, USA from 21 states between March 2009 and April 2012 [Citation4]. More details regarding the diagnostic findings for these bats can be found in [Citation13], including Table S1. We conducted histological and microscopic observations on bats submitted to the USGS National Wildlife Health Center for diagnosis and documentation of WNS. We followed standard protocols for histology [Citation58]. Briefly, we fixed sections of skin in 10% neutral buffered formalin, embedded them in paraffin, cut them to a thickness of 5 microns, placed them on glass slides, and stained with the periodic acid-Schiff (PAS). To conduct transmission electron microscopy, we fixed tissues in 10% neutral buffered formalin and then post-fixed samples of the patagium in glutaraldehyde. We embedded tissues in epon embedding resin and poststained thin sections with saturated aqueous uranyl acetate and Reynold’s lead citrate. Finally, we collected images with a 200kV JEOL 2100 LaB6, transmission electron microscope (Peabody, Massachusetts).

Phylogenomic analysis

We retrieved the protein gene catalog of P. destructans, version 20170919, from Mycocosm, the Joint Genome Institute (JGI) fungal genome database (https://mycocosm.jgi.doe.gov/Pseudest1/Pseudest1.home.html) [Citation59]. The complete proteome was then used as a query for a homology search in the National Center for Biotechnology Information (NCBI) non-redundant (nr) protein database. For computational efficiency, we downloaded the complete nr dataset from the NCBI server (http://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz) and used the DIAMOND program as a rapid alternate to BLAST [Citation60]. We structured the nr dataset into a DIAMOND database including taxonomic information, and then conducted the homology search using the “blastp” algorithm of DIAMOND, with the options -k0 (to report all significant alignments for each query sequence), -e 1e-6 (significance threshold of 1e-6 for reporting an alignment), -b4 (block size of 4, as automatically suggested by the program based on the computer configuration), and – taxonlist 4751 (to limit the search to fungal sequences). The search results were generated in tabular format, and we further split them per query sequence to facilitate downstream analyses. In cases where multiple high scoring pairs (HSP) were found for a reference sequence, we only retained the most significant one. All scripts and data are available at https://gitlab.gwdg.de/molsysevol/pseudogymnoascus-destructans-phylogeny/.

We then reconstructed the presence-absence pattern of each one of the 9,310 P. destructans genes using the taxonomic information produced in the DIAMOND database. Our database search revealed homologous sequences in 36,578 fungal species, so we generated a 9,310 x 36,578 table containing the number of homologous sequences found in each species for each gene. In this context, a table entry of zero indicates the gene was not found in a particular species, either because it does not have a homologous sequence (i.e. the corresponding homolog was not sequenced) or because the sequence was too divergent to lead to a similarity score higher than the threshold used for significance assessment.

As we are interested in speciation and not duplication events, we discarded paralogous sequences, in particular to minimize the problem of hidden paralogy. For each species, we set the occurrence number of the genes that were present in multiple copies to zero in the matrix. In order to minimize the final amount of missing data, we first only retain species that had at least 1,000 genes with an identified homolog in P. destructans, keeping a total of 1,257 species. A lack of quality sequence data for the congener, Pseudogymnoascus pannorum, precluded its inclusion beyond this step in this analysis (the cause of “snake fungal disease”, Ophidiomyces ophiodiicola, also lacked adequate sequence data to be included). Second, we kept only genes present in at least 80% of all species, reducing the number of genes to 1,146.

We next retrieved the amino-acid sequences of the selected species for each of the 1,146 selected genes. We aligned each family dataset using two distinct methods: Muscle [Citation61] and Clustal Omega [Citation62]. We compared the two resulting alignments using the bppAlnScore program from the Bio++ Program Suite [Citation63], and only selected positions recovered identically by the two methods for downstream analyses. We ran all analyses in parallel using the GNU parallel software [Citation64]. All 1,146 filtered alignments were concatenated, leading to a total alignment of 461,596 amino-acid positions.

We conducted a phylogenetic analysis in several steps because the large data matrix prohibited the inference of a maximum likelihood tree in a reasonable time. First, we used FastTree [Citation65] and the complete data set to build an initial maximum likelihood tree using the Le and Gascuel substitution model [Citation66]. In order to run a more thorough search using the RaxML-NG software [Citation67], we reduced the dataset by removing sequences that are highly similar. Using the PhySamp program [Citation68], we sampled sequences based on their similarity: when two sequences were more than 97% identical, the one with the lowest amount of missing data (unresolved character or gap) was kept and the other sequence discarded. The resulting similarity clusters and their representative members are listed in Supplementary Table S1, and the final sampled tree included 817 species. We used RaxML-NG software [Citation67] to estimate the phylogeny from the concatenated alignment by maximum likelihood (ML), following the protocol described in [Citation69]. We used an LG model of amino-acid substitutions [Citation66] with a four-classes gamma distribution of rates. We used the FastTree tree as starting tree and compared two models: one where all 1,146 genes followed the same model (1,632 parameters), and a partition model where each gene was allowed to have its own rate distribution and a distinctive branch scaling factor (3,922 parameters). This latter model led to a lower Akaike Information criterion (AIC) (397,373,191 vs. 398,450,754) and was selected for downstream analyses. The resulting tree is plotted in Supplementary Figure S1. We then focused on the Pezizomycotina clade. We pruned the subtree from the RaxML tree and selected the corresponding sequences. We then conducted a more thorough ML inference on this sub-dataset, using 20 distinct initial trees (10 random and 10 obtained by maximum parsimony) and retained the resulting estimated tree with the highest likelihood. To obtain support values for all branches in the tree, we further conducted a bootstrap analysis with 100 replicates, which was sufficient to get reliable bootstrap values based on the extended majority rule (MRE)-based bootstrapping test implemented in RaxML [Citation70]. We plotted the resulting phylogeny using the ggtree package for R [Citation71]. The data that support this study are openly available in pseudogymonoscus-destructans-phylogeny at https://gitlab.gwdg.de/molsysevol/

Life history trait annotations

We classified each species as one of the following: animal associated (including nematophagous, pathogen, coprophilous), plant associated (including saprotroph, endophyte, epiphyte, mycorrhizae), plant pathogen, mycoparasite, lichen, aquatic, soil. We first attempted to assign classifications based on descriptions from the JGI fungal genomics database. For species not well described in the database, we conducted a general literature search, and recorded species as missing data when no reliable information could be found. Many fungal species are associated with multiple environments or organisms, so we conservatively recorded a pathogen as animal associated if there is a record of any animal infection. All assigned life history traits are available in Figure S1.

Ancestral trait reconstruction

We inferred the characteristics of the P. destructans ancestor using models of discrete character evolution. We used the Pezizomycotina maximum likelihood tree and discarded species for which no life history trait data could be found in the literature. The resulting dataset contained 355 species and seven character states: “animal associated” (87 species), “aquatic” (7 species), “lichen” (3 species), “mycoparasite” (4 species), “plant associated” (114 species), “plant pathogen” (121 species) and “soil” (19 species). Because some character states were observed in low frequency, we also analyzed a reduced dataset restricted to the most frequent states: “animal associated”, “plant associated”, “plant pathogen” and “soil,” covering 341 species. We fit three models of character-state transition to each dataset: all changes equally likely (equal rate, ER), all rates different but symmetric changes equally likely (symmetric model, SYM), and all rates different (ARD). We used the “ace” function of the ‘ape’ package for R [Citation72] to fit models and calculate Akaike’s information criterion (AIC) for each model. To estimate the character state of the P. destructans ancestor, we used a stochastic mapping procedure [Citation73,Citation74] as implemented in the ‘phytools’ package for R [Citation75]. We sampled 1,000 maps of character changes from the posterior distribution, accounting for the uncertainty in the substitution matrix (argument Q = “msmc” in the ‘make.simmap’ function). We used the ‘describe.simmap’ function for ancestral state estimation of each model and the ‘aic.w’ function and AIC to average models (Supplementary Table S2). Scripts used for the analyses are openly available in https://gitlab.gwdg.de/molsysevol/pseudogymnoascus-destructans-phylogeny.

Author contributions

CUM conceived the study, performed analyses, wrote the original draft, performed histopathology, photography and visualization of images

JYD Formal analysis, investigation, resources, visualization and writing – review and editing

MKK Investigation, resources, visualization of the data and writing – review and editing

JGB Conceptualization and writing – review and editing

EHS Formal analysis, investigation, writing – review and editing

Supplemental material

Supplemental Material

Download Zip (690 KB)

Acknowledgements

We would like to thank the many field biologists from State, Regional, and Federal organizations that documented progression as P. destructans spread across the United States of America and affected regions in Canada. Special thanks to Kathryn Wesenberg for help with references and others at the National Wildlife Health Center for supporting the effort to define white-nose syndrome, and to Emmanuel Paradis for helping with the use of the R ‘ace’ function. EHS is grateful for support from the Canadian Institute for Advanced Research (CIFAR). Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Disclosure statement

No potential conflict of interest was reported by the authors.

Data availability statement

The data that support the findings of this study are openly available in Zenodo at DOI reference number 10.5281/zenodo. 5939906 at https://zenodo.org/record/5939906#.YjpiN-rMKUk.

Supplementary material

Supplemental data for this article can be accessed online at https://doi.org/10.1080/21505594.2022.2082139

Additional information

Funding

The authors reported there is no funding associated with the work featured in this article.

References

  • Hoyt JR, Kilpatrick AM, Langwig KE. Ecology and impacts of white-nose syndrome on bats. Nat Rev Microbiol. 2021;19(3):196–210.
  • Gargas A, Trest M, Christensen M, et al. Geomyces destructans sp. nov. associated with bat white-nose syndrome. Mycotaxon. 2009;108(1):147–154.
  • Minnis AM, Lindner DL. Phylogenetic evaluation of Geomyces and allies reveals no close relatives of Pseudogymnoascus destructans, comb. nov., in bat hibernacula of eastern North America. Fungal Biol. 2013;117(9):638–649.
  • Turner GG, Meteyer CU, Barton H, et al. Nonlethal screening of bat-wing skin with the use of ultraviolet fluorescence to detect lesions indicative of white-nose syndrome. J Wildl Dis. 2014;50(3):566–573. DOI:10.7589/2014-03-058
  • Meteyer CU, Buckles EL, Blehert DS, et al. Histopathologic criteria to confirm white-nose syndrome in bats. J Vet Diagn Investig. 2009;21(4):411–414. DOI:10.1177/104063870902100401
  • Blehert DS, Hicks AC, Behr M, et al. Bat white-nose syndrome: an emerging fungal pathogen? Science (Washington, DC). 2009;323:227.
  • Lorch JM, Minnis AM, Meteyer CU, et al. The fungus Trichophyton redellii sp. nov. causes skin infections that resemble white-nose syndrome of hibernating bats. J Wildl Dis. 2015;51(1):36–47. DOI:10.7589/2014-05-134
  • Berbee ML. The phylogeny of plant and animal pathogens in the Ascomycota. Physiol Mol Plant Pathol. 2001;59(4):165–187.
  • Verant ML, Boyles JG, Waldrep JW, et al. Temperature-Dependent growth of Geomyces destructans, the fungus that causes bat white-nose syndrome. PLoS One. 2012;7:e46280.
  • Chinnapun D. Virulence factors involved in pathogenicity of dermatophytes. Walailak J Sci Technol. 2015;12:573–580.
  • Raudabaugh DB, Miller AN. Nutritional capability of and substrate suitability for Pseudogymnoascus destructans, the causal agent of bat white-nose syndrome. PLoS One. 2013;8(10):e78300.
  • Chaturvedi V, Springer DJ, Behr MJ, et al. Morphological and molecular characterizations of psychrophilic fungus Geomyces destructans from New York bats with white nose syndrome (WNS). PLoS One. 2010;5(5):e10783. DOI:10.1371/journal.pone.0010783
  • Reynolds HT, Barton HA. Comparison of the white-nose syndrome agent Pseudogymnoascus destructans to cave-dwelling relatives suggests reduced saprotrophic enzyme activity. PLoS One. 2014;9(1):e86437.
  • Ye F, Li M, Zhu S, et al. Diagnosis of dermatophytosis using single fungus endogenous fluorescence spectrometry. Biomed Opt Express. 2018;9(6):2733–2742.
  • Cruz-Mireles N, Eseola AB, Osés-Ruiz M, et al. From appressorium to transpressorium—defining the morphogenetic basis of host cell invasion by the rice blast fungus. PLoS Path. 2021;17(7):e1009779.
  • Donaldson ME, Davy CM, Vanderwolf KJ, et al. Growth medium and incubation temperature alter the Pseudogymnoascus destructans transcriptome: implications in identifying virulence factors. Mycologia. 2018;110(2):300–315.
  • Reynolds HT, Ingersoll T, Barton HA. Modeling the environmental growth of Pseudogymnoascus destructans and its impact on the white-nose syndrome epidemic. J Wildl Dis. 2015;51(2):318–331.
  • Wilson MB, Held BW, Freiborg AH, et al. Resource capture and competitive ability of non-pathogenic Pseudogymnoascus spp. And P. destructans, the cause of white-nose syndrome in bats. PLoS One. 2017;12:e0178968.
  • Reynolds HT, Barton HA, Slot JC. Phylogenomic analysis supports a recent change in nitrate assimilation in the white-nose syndrome pathogen, Pseudogymnoascus destructans. Fungal Ecol. 2016;23:20–29.
  • McGuire LP, Turner JM, Warnecke L, et al. White-Nose syndrome disease severity and a comparison of diagnostic methods. Ecohealth. 2016;13(1):60–71. DOI:10.1007/s10393-016-1107-y
  • Meteyer CU, Valent M, Kashmer J, et al. Recovery of little brown bats (Myotis lucifugus) from natural infection with Geomyces destructans, white-nose syndrome. J Wildl Dis. 2011;47(3):618–626. DOI:10.7589/0090-3558-47.3.618
  • Koeck M, Hardham AR, Dodds PN. The role of effectors of biotrophic and hemibiotrophic fungi in infection. Cell Microbiol. 2011;13(12):1849–1857.
  • Chowdhury S, Basu A, Kundu S. Biotrophy-Necrotrophy switch in pathogen evoke differential response in resistant and susceptible sesame involving multiple signaling pathways at different phases. Sci Rep. 2017;7(1):17251.
  • Brun S, Silar P Convergent evolution of morphogenetic processes in fungi. In: Evolutionary biology–concepts, molecular and morphological evolution. Berlin, Heidelberg: Springer; 2010. p. 317–328.
  • Cryan PM, Meteyer CU, Boyles JG, et al. Wing pathology of white-nose syndrome in bats suggests life-threatening disruption of physiology. BMC Biol. 2010;8(1):1–8.
  • Verant ML, Meteyer CU, Speakman JR, et al. White-Nose syndrome initiates a cascade of physiologic disturbances in the hibernating bat host. BMC Physiol. 2014;14(1):1–11.
  • Van Rooij P, Martel A, Haesebrouck F, et al. Amphibian chytridiomycosis: a review with focus on fungus-host interactions. Vet Res. 2015;46(1):1–22.
  • Haffner M. The size of sebaceous glands in relation to the size of hair follicles on the heads of some small mammals (Insectivore, Chiroptera, Rodentia). Ann Anat. 1998;180(2):165–171.
  • Fuller KK, Rhodes JC. Protein kinase a and fungal virulence: a sinister side to a conserved nutrient sensing pathway. Virulence. 2012;3:109–121.
  • Tariq V, Jeffries P. Appressorium formation by Sclerotinia sclerotiorum: scanning electron microscopy. Trans Br Mycol Soc. 1984;82(4):645–651.
  • van Kan JA. Licensed to kill: the lifestyle of a necrotrophic plant pathogen. Trends Plant Sci. 2006;11(5):247–253.
  • Ikeda K, Park P, Nakayashiki H. Cell biology in phytopathogenic fungi during host infection: Commonalities and differences. J Gen Plant Pathol. 2019;85(3):163–173.
  • Demoor A, Silar P, Brun S. Appressorium: the breakthrough in Dikarya. J Fungi. 2019;5(3):72.
  • Schnee S, Rougeux E, Pezet R, et al. Evidence for constitutive cutinase in urgeminated conidia of Erysiphe necator Schwein. J Cytol Histol. 2013;4:197.
  • Duan Y, Wu H, Ma Z, et al. Scanning electron microscopy and histopathological observations of Beauveria bassiana infection of Colorado potato beetle larvae. Microb Pathog. 2017;111:435–439.
  • Hajek AE. Pathology and epizootiology of Entomophaga maimaiga infections in forest Lepidoptera. Microbiol Mol Biol Rev. 1999;63(4):814–835.
  • Brandner J, Zorn-Kruppa M, Yoshida T, et al. Epidermal tight junctions in health and disease. Tissue Barriers. 2015;3(1–2):e974451.
  • Reeder SM, Palmer JM, Prokkola JM, et al. Pseudogymnoascus destructans transcriptome changes during white-nose syndrome infections. Virulence. 2017;8(8):1695–1707.
  • Valent B, Khang CH. Recent advances in rice blast effector research. Curr Opin Plant Biol. 2010;13(4):434–441.
  • Palmer JM, Drees KP, Foster JT, et al. Extreme sensitivity to ultraviolet light in the fungal pathogen causing white-nose syndrome of bats. Nat Commun. 2018;9(1):1–10.
  • Crous PW, Groenewald JZ. They seldom occur alone. Fungal Biol. 2016;120(11):1392–1415.
  • Oguchi T. Studies on the species of Lachnellula in Hokkaido: their morphology, physiology, and pathogenicity. Bull Hokkaido for Exp Stn. 1981;19:187–246.
  • Dean R, Van Kan JA, Pretorius ZA, et al. The Top 10 fungal pathogens in molecular plant pathology. Mol Plant Pathol. 2012;13(4):414–430. DOI:10.1111/j.1364-3703.2011.00783.x
  • Pannkuk EL, Risch TS, Savary BJ. Isolation and identification of an extracellular subtilisin-like serine protease secreted by the bat pathogen Pseudogymnoascus destructans. PLoS One. 2015;10(3):e0120508.
  • Chaturvedi V, DeFiglio H, Chaturvedi S. Phenotype profiling of white-nose syndrome pathogen Pseudogymnoascus destructans and closely-related Pseudogymnoascus pannorum reveals metabolic differences underlying fungal lifestyles. F1000res. 2018;7:665.
  • Mosbrugger V, Utescher T, Dilcher DL. Cenozoic continental climatic evolution of Central Europe. Proc Natl Acad Sci USA. 2005;102(42):14964–14969.
  • Zachos JC, Dickens GR, Zeebe RE. An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. 2008;451(7176):279–283.
  • Lovegrove BG. The evolution of endothermy in Cenozoic mammals: a plesiomorphic‐apomorphic continuum. Biol Rev (Camb). 2012;87(1):128–162.
  • Canale CI, Levesque DL, Lovegrove BG. Tropical heterothermy: does the exception prove the rule or force a re-definition? Living in a Seasonal World. 2012:29–40.
  • Davy CM, Donaldson ME, Bandouchova H, et al. Transcriptional host–pathogen responses of Pseudogymnoascus destructans and three species of bats with white-nose syndrome. Virulence. 2020;11(1):781–794. DOI:10.1080/21505594.2020.1768018
  • Lilley TM, Prokkola JM, Blomberg AS, et al. Resistance is futile: RNA-sequencing reveals differing responses to bat fungal pathogen in Nearctic Myotis lucifugus and Palearctic Myotis myotis. Oecologia. 2019;191(2):295–309. DOI:10.1007/s00442-019-04499-6
  • Martínková N, Bačkor P, Bartonička T, et al. Increasing incidence of Geomyces destructans fungus in bats from the Czech Republic and Slovakia. PLoS One. 2010;5(11):e13853. DOI:10.1371/journal.pone.0013853
  • Drees KP, Lorch JM, Puechmaille SJ, et al. Phylogenetics of a fungal invasion: origins and widespread dispersal of white-nose syndrome. Mbio. 2017;8(6): e01941-17. doi:10.1128/mBio.01941-17.
  • Trivedi J, Lachapelle J, Vanderwolf KJ, et al. Fungus causing white-nose syndrome in bats accumulates genetic variability in North America with no sign of recombination. mSphere. 2017;2(4): e00271-17. doi:10.1128/mSphereDirect.00271-17.
  • Puechmaille SJ, Wibbelt G, Korn V, et al. Pan-European distribution of white-nose syndrome fungus (Geomyces destructans) not associated with mass mortality. PLoS One. 2011;6(4):e19167. DOI:10.1371/journal.pone.0019167
  • Warnecke L, Turner JM, Bollinger TK, et al. Inoculation of bats with European Geomyces destructans supports the novel pathogen hypothesis for the origin of white-nose syndrome. Proc Natl Acad Sci USA. 2012;109(18):6999–7003. DOI:10.1073/pnas.1200374109
  • Kovacova V, Zukal J, Bandouchova H, et al. White-Nose syndrome detected in bats over an extensive area of Russia. BMC Vet Res. 2018;14(1):1–9. DOI:10.1186/s12917-018-1521-1
  • Prophet EB. Laboratory methods in histotechnology. Amer Registry of Pathology; 1992.
  • Drees KP, Palmer JM, Sebra R, et al. Use of multiple sequencing technologies to produce a high-quality genome of the fungus Pseudogymnoascus destructans , the causative agent of bat white-nose syndrome. Genome Announc. 2016;4(3): e00445-16. doi:10.1128/genomeA.00445-16.
  • Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.
  • Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–1797.
  • Sievers F, Higgins DG. Clustal omega, accurate alignment of very large numbers of sequences. Multiple sequence alignment methods. Methods Mol Biol. 2014;1079:105–116.
  • Dutheil J, Boussau B. Non-Homogeneous models of sequence evolution in the bio++ suite of libraries and programs. BMC Evol Biol. 2008;8(1):1–12.
  • Tange O. GNU parallel: the command-line power tool. 2011.
  • Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490.
  • Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25(7):1307–1320.
  • Kozlov AM, Darriba D, Flouri T, et al. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics. 2019;35:4453–4455.
  • Guéguen L, Gaillard S, Boussau B, et al. Bio++: efficient extensible libraries and tools for computational molecular evolution. Mol Biol Evol. 2013;30(8):1745–1750. DOI:10.1093/molbev/mst097
  • Kozlov AM, Stamatakis A. Using RAxML-NG in practice. In: Scornavacca C; F Delsuc and N Galtier, editors. Phylogenomics in the Genomic Era: No commercial publisher | Authors open access book. 2020. pp. 1.3:11.3:1.3:25
  • Pattengale ND, Alipour M, Bininda-Emonds OR, et al. How many bootstrap replicates are necessary? J Comput Biol. 2010;17(3):337–354.
  • Yu G, Lam T-Y, Zhu H, et al. Two methods for mapping and visualizing associated data on phylogeny using ggtree. Mol Biol Evol. 2018;35(12):3041–3043.
  • Paradis E, Schliep K. Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2019;35(3):526–528.
  • Huelsenbeck JP, Nielsen R, Bollback JP. Stochastic mapping of morphological characters. Syst Biol. 2003;52(2):131–158.
  • Nielsen R, Huelsenbeck J. Mapping mutations on phylogenies. Syst Biol. 2002;51(5):729–739.
  • Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–223.