Pan-cancer analyses reveal cancer-type-specific fungal ecologies and bacteriome interactions.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

Key processes required for the different stages of fungal carnivory by a nematode-trapping fungus.

Nutritional deprivation triggers a switch from a saprotrophic to predatory lifestyle in soil-dwelling nematode-trapping fungi (NTF). In particular, the NTF Arthrobotrys oligospora secretes food and sex cues to lure nematodes to its mycelium and is triggered to develop specialized trapping devices. Captured nematodes are then invaded and digested by the fungus, thus serving as a food source. In this study, we examined the transcriptomic response of A. oligospora across the stages of sensing, trap development, and digestion upon exposure to the model nematode Caenorhabditis elegans. A. oligospora enacts a dynamic transcriptomic response, especially of protein secretion-related genes, in the presence of prey. Two-thirds of the predicted secretome of A. oligospora was up-regulated in the presence of C. elegans at all time points examined, and among these secreted proteins, 38.5% are predicted to be effector proteins. Furthermore, functional studies disrupting the t-SNARE protein Sso2 resulted in impaired ability to capture nematodes. Additionally, genes of the DUF3129 family, which are expanded in the genomes of several NTF, were highly up-regulated upon nematode exposure. We observed the accumulation of highly expressed DUF3129 proteins in trap cells, leading us to name members of this gene family as Trap Enriched Proteins (TEPs). Gene deletion of the most highly expressed TEP gene, TEP1, impairs the function of traps and prevents the fungus from capturing prey efficiently. In late stages of predation, we observed up-regulation of a variety of proteases, including metalloproteases. Following penetration of nematodes, these metalloproteases facilitate hyphal growth required for colonization of prey. These findings provide insights into the biology of the predatory lifestyle switch in a carnivorous fungus and provide frameworks for other fungal-nematode predator-prey systems.

Toward a Fully Resolved Fungal Tree of Life.

In this review, we discuss the current status and future challenges for fully elucidating the fungal tree of life. In the last 15 years, advances in genomic technologies have revolutionized fungal systematics, ushering the field into the phylogenomic era. This has made the unthinkable possible, namely access to the entire genetic record of all known extant taxa. We first review the current status of the fungal tree and highlight areas where additional effort will be required. We then review the analytical challenges imposed by the volume of data and discuss methods to recover the most accurate species tree given the sea of gene trees. Highly resolved and deeply sampled trees are being leveraged in novel ways to study fungal radiations, species delimitation, and metabolic evolution. Finally, we discuss the critical issue of incorporating the unnamed and uncultured dark matter taxa that represent the vast majority of fungal diversity.

The pan-genome of Aspergillus fumigatus provides a high-resolution view of its population structure revealing high levels of lineage-specific diversity driven by recombination.

Aspergillus fumigatus is a deadly agent of human fungal disease where virulence heterogeneity is thought to be at least partially structured by genetic variation between strains. While population genomic analyses based on reference genome alignments offer valuable insights into how gene variants are distributed across populations, these approaches fail to capture intraspecific variation in genes absent from the reference genome. Pan-genomic analyses based on de novo assemblies offer a promising alternative to reference-based genomics with the potential to address the full genetic repertoire of a species. Here, we evaluate 260 genome sequences of A. fumigatus including 62 newly sequenced strains, using a combination of population genomics, phylogenomics, and pan-genomics. Our results offer a high-resolution assessment of population structure and recombination frequency, phylogenetically structured gene presence-absence variation, evidence for metabolic specificity, and the distribution of putative antifungal resistance genes. Although A. fumigatus disperses primarily via asexual conidia, we identified extraordinarily high levels of recombination with the lowest linkage disequilibrium decay value reported for any fungal species to date. We provide evidence for 3 primary populations of A. fumigatus, with recombination occurring only rarely between populations and often within them. These 3 populations are structured by both gene variation and distinct patterns of gene presence-absence with unique suites of accessory genes present exclusively in each clade. Accessory genes displayed functional enrichment for nitrogen and carbohydrate metabolism suggesting that populations may be stratified by environmental niche specialization. Similarly, the distribution of antifungal resistance genes and resistance alleles were often structured by phylogeny. Altogether, the pan-genome of A. fumigatus represents one of the largest fungal pan-genomes reported to date including many genes unrepresented in the Af293 reference genome. These results highlight the inadequacy of relying on a single-reference genome-based approach for evaluating intraspecific variation and the power of combined genomic approaches to elucidate population structure, genetic diversity, and putative ecological drivers of clinically relevant fungi.

The population genetics of the causative agent of snake fungal disease indicate recent introductions to the USA.

Snake fungal disease (SFD; ophidiomycosis), caused by the pathogen Ophidiomyces ophiodiicola (Oo), has been documented in wild snakes in North America and Eurasia, and is considered an emerging disease in the eastern United States of America. However, a lack of historical disease data has made it challenging to determine whether Oo is a recent arrival to the USA or whether SFD emergence is due to other factors. Here, we examined the genomes of 82 Oo strains to determine the pathogen's history in the eastern USA. Oo strains from the USA formed a clade (Clade II) distinct from European strains (Clade I), and molecular dating indicated that these clades diverged too recently (approximately 2,000 years ago) for transcontinental dispersal of Oo to have occurred via natural snake movements across Beringia. A lack of nonrecombinant intermediates between clonal lineages in Clade II indicates that Oo has actually been introduced multiple times to North America from an unsampled source population, and molecular dating indicates that several of these introductions occurred within the last few hundred years. Molecular dating also indicated that the most common Clade II clonal lineages have expanded recently in the USA, with time of most recent common ancestor mean estimates ranging from 1985 to 2007 CE. The presence of Clade II in captive snakes worldwide demonstrates a potential mechanism of introduction and highlights that additional incursions are likely unless action is taken to reduce the risk of pathogen translocation and spillover into wild snake populations.

Diploid-dominant life cycles characterize the early evolution of Fungi.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Phylogenomic analysis of the Neocallimastigomycota: proposal of Caecomycetaceae fam. nov., Piromycetaceae fam. nov., and emended description of the families Neocallimastigaceae and Anaeromycetaceae.

The anaerobic gut fungi (AGF) represent a coherent phylogenetic clade within the Mycota. Twenty genera have been described so far. Currently, the phylogenetic and evolutionary relationships between AGF genera remain poorly understood. Here, we utilized 52 transcriptomic datasets from 14 genera to resolve AGF inter-genus relationships using phylogenomics, and to provide a quantitative estimate (amino acid identity, AAI) for intermediate rank assignments. We identify four distinct supra-genus clades, encompassing all genera producing polyflagellated zoospores, bulbous rhizoids, the broadly circumscribed genus Piromyces, and the Anaeromyces and affiliated genera. We also identify the genus Khoyollomyces as the earliest evolving AGF genus. Concordance between phylogenomic outputs and RPB1 and D1/D2 LSU, but not RPB2, MCM7, EF1α or ITS1, phylogenies was observed. We combine phylogenomic analysis and AAI outputs with informative phenotypic traits to propose accommodating 14/20 AGF genera into four families: Caecomycetaceae fam. nov. (encompassing the genera Caecomyces and Cyllamyces), Piromycetaceae fam. nov. (encompassing the genus Piromyces), emend the description of the family Neocallimastigaceae to encompass the genera Neocallimastix, Orpinomyces, Pecoramyces, Feramyces, Ghazallomyces, Aestipascuomyces and Paucimyces, as well as the family Anaeromycetaceae to include the genera Oontomyces, Liebetanzomyces and Capellomyces in addition to Anaeromyces. We refrain from proposing families for the deeply branching genus Khoyollomyces and for genera with uncertain position (Buwchfawromyces, Joblinomyces, Tahromyces, Agriosomyces and Aklioshbomyces) pending availability of additional isolates and sequence data; and these genera are designated as 'genera incertae sedis' in the order Neocallimastigales. Our results establish an evolutionary-grounded Linnaean taxonomic framework for the AGF, provide quantitative estimates for rank assignments, and demonstrate the utility of RPB1 as an additional informative marker in Neocallimastigomycota taxonomy.

The reference genome assembly of the bright cobblestone lichen, Acarospora socialis.

Acarospora socialis, the bright cobblestone lichen, is commonly found in southwestern North America. This charismatic yellow lichen is a species of key ecological significance as it is often a pioneer species in new environments. Despite their ecological importance virtually no research has been conducted on the genomics of A. socialis. To address this, we used long-read sequencing to generate the first high-quality draft genome of A. socialis. Lichen thallus tissue was collected from Pinkham Canyon in Joshua Tree National Park, California and deposited in the UC Riverside herbarium under accession #295874. The de novo assembly of the mycobiont partner of the lichen was generated from Pacific Biosciences HiFi long reads and Dovetail Omni-C chromatin capture data. After removing algal and bacterial contigs, the fungal genome was approximately 31.2 Mb consisting of 38 scaffolds with contig and scaffold N50 of 2.4 Mb. The BUSCO completeness score of the assembled genome was 97.5% using the Ascomycota gene set. Information on the genome of A. socialis is important for California conservation purposes given that this lichen is threatened in some places locally by wildfires due to climate change. This reference genome will be used for understanding the genetic diversity, population genomics, and comparative genomics of A. socialis species. Genomic resources for this species will support population and landscape genomics investigations, exploring the use of A. socialis as a bioindicator species for climate change, and in studies of adaptation by comparing populations that occur across aridity gradients in California.

Genomic diversity generated by a transposable element burst in a rice recombinant inbred population.

Genomes of all characterized higher eukaryotes harbor examples of transposable element (TE) bursts-the rapid amplification of TE copies throughout a genome. Despite their prevalence, understanding how bursts diversify genomes requires the characterization of actively transposing TEs before insertion sites and structural rearrangements have been obscured by selection acting over evolutionary time. In this study, rice recombinant inbred lines (RILs), generated by crossing a bursting accession and the reference Nipponbare accession, were exploited to characterize the spread of the very active Ping/mPing family through a small population and the resulting impact on genome diversity. Comparative sequence analysis of 272 individuals led to the identification of over 14,000 new insertions of the mPing miniature inverted-repeat transposable element (MITE), with no evidence for silencing of the transposase-encoding Ping element. In addition to new insertions, Ping-encoded transposase was found to preferentially catalyze the excision of mPing loci tightly linked to a second mPing insertion. Similarly, structural variations, including deletion of rice exons or regulatory regions, were enriched for those with break points at one or both ends of linked mPing elements. Taken together, these results indicate that structural variations are generated during a TE burst as transposase catalyzes both the high copy numbers needed to distribute linked elements throughout the genome and the DNA cuts at the TE ends known to dramatically increase the frequency of recombination.

Natural diversity in the predatory behavior facilitates the establishment of a robust model strain for nematode-trapping fungi.

Nematode-trapping fungi (NTF) are a group of specialized microbial predators that consume nematodes when food sources are limited. Predation is initiated when conserved nematode ascaroside pheromones are sensed, followed by the development of complex trapping devices. To gain insights into the coevolution of this interkingdom predator-prey relationship, we investigated natural populations of nematodes and NTF that we found to be ubiquitous in soils. Arthrobotrys species were sympatric with various nematode species and behaved as generalist predators. The ability to sense prey among wild isolates of Arthrobotrys oligospora varied greatly, as determined by the number of traps after exposure to Caenorhabditis elegans While some strains were highly sensitive to C. elegans and the nematode pheromone ascarosides, others responded only weakly. Furthermore, strains that were highly sensitive to the nematode prey also developed traps faster. The polymorphic nature of trap formation correlated with competency in prey killing, as well as with the phylogeny of A. oligospora natural strains, calculated after assembly and annotation of the genomes of 20 isolates. A chromosome-level genome assembly and annotation were established for one of the most sensitive wild isolates, and deletion of the only G-protein ß-subunit-encoding gene of A. oligospora nearly abolished trap formation. In summary, our study establishes a highly responsive A. oligospora wild isolate as a model strain for the study of fungus-nematode interactions and demonstrates that trap formation is a fitness character in generalist predators of the nematode-trapping fungus family.

Transcriptome and population structure of glassy-winged sharpshooters (Homalodisca vitripennis) with varying insecticide resistance in southern California.

BACKGROUND: Homalodisca vitripennis Germar, the glassy-winged sharpshooter, is an invasive insect in California and a critical threat to agriculture through its transmission of the plant pathogen, Xylella fastidiosa. Quarantine, broad-spectrum insecticides, and biological control have been used for population management of H. vitripennis since its invasion and subsequent proliferation throughout California. Recently wide-spread neonicotinoid resistance has been detected in populations of H. vitripennis in the southern portions of California's Central Valley. In order to better understand potential mechanisms of H. vitripennis neonicotinoid resistance, we performed RNA sequencing on wild-caught insecticide-resistant and relatively susceptible sharpshooters to profile their transcriptome and population structure. RESULTS: We identified 81 differentially expressed genes with higher expression in resistant individuals. The significant largest differentially expressed candidate gene linked to resistance status was a cytochrome P450 gene with similarity to CYP6A9. Furthermore, we observed an over-enrichment of GO terms representing functions supportive of roles in resistance mechanisms (cytochrome P450s, M13 peptidases, and cuticle structural proteins). Finally, we saw no evidence of broad-scale population structure, perhaps due to H. vitripennis' relatively recent introduction to California or due to the relatively small geographic scale investigated here. CONCLUSIONS: In this work, we characterized the transcriptome of insecticide-resistant and susceptible H. vitripennis and identified candidate genes that may be involved in resistance mechanisms for this species. Future work should seek to build on the transcriptome profiling performed here to confirm the role of the identified genes, particularly the cytochrome P450, in resistance in H. vitripennis. We hope this work helps aid future population management strategies for this and other species with growing insecticide resistance.

Single nucleotide polymorphisms and chromosomal copy number variation may impact the Sporothrix brasiliensis antifungal susceptibility and sporotrichosis clinical outcomes.

Feline-transmitted sporotrichosis has garnered attention due to the recent high incidence and the lack of efficient control in the epicenter of the epidemic, Rio de Janeiro, Brazil. Sporothrix brasiliensis is the major pathogen involved in feline-to-human sporotrichosis in Brazil and displays more virulent genotypes than the closely related species S. schenckii. Over the last two decades, several reports of antifungal-resistant strains have emerged. Sequencing and comparison analysis of the outbreak strains allowed us to observe that the azole non-wild-type S. brasiliensis strain CFP 1054 had significant chromosomal variations compared to wild-type strains. One of these variants includes a region of 231 Kb containing 75 duplicated genes, which were overrepresented for lipid and isoprenoid metabolism. We also identified an additional strain (CFP 1055) that was resistant to itraconazole and amphotericin B, which had a single nucleotide polymorphism in the tac1 gene. The patients infected with these two strains showed protracted clinical course and sequelae. Even though our sample size is modest, these results suggest the possibility of identifying specific point mutations and large chromosomal duplications potentially associated with antifungal resistance and clinical outcomes of sporotrichosis.

Comparative genomics reveals dynamic genome evolution in host specialist ectomycorrhizal fungi.

While there has been significant progress characterizing the 'symbiotic toolkit' of ectomycorrhizal (ECM) fungi, how host specificity may be encoded into ECM fungal genomes remains poorly understood. We conducted a comparative genomic analysis of ECM fungal host specialists and generalists, focusing on the specialist genus Suillus. Global analyses of genome dynamics across 46 species were assessed, along with targeted analyses of three classes of molecules previously identified as important determinants of host specificity: small secreted proteins (SSPs), secondary metabolites (SMs) and G-protein coupled receptors (GPCRs). Relative to other ECM fungi, including other host specialists, Suillus had highly dynamic genomes including numerous rapidly evolving gene families and many domain expansions and contractions. Targeted analyses supported a role for SMs but not SSPs or GPCRs in Suillus host specificity. Phylogenomic-based ancestral state reconstruction identified Larix as the ancestral host of Suillus, with multiple independent switches between white and red pine hosts. These results suggest that like other defining characteristics of the ECM lifestyle, host specificity is a dynamic process at the genome level. In the case of Suillus, both SMs and pathways involved in the deactivation of reactive oxygen species appear to be strongly associated with enhanced host specificity.

Parasitoid Jewel Wasp Mounts Multipronged Neurochemical Attack to Hijack a Host Brain.

The parasitoid emerald jewel wasp Ampulex compressa induces a compliant state of hypokinesia in its host, the American cockroach Periplaneta americana through direct envenomation of the central nervous system (CNS). To elucidate the biochemical strategy underlying venom-induced hypokinesia, we subjected the venom apparatus and milked venom to RNAseq and proteomics analyses to construct a comprehensive "venome," consisting of 264 proteins. Abundant in the venome are enzymes endogenous to the host brain, including M13 family metalloproteases, phospholipases, adenosine deaminase, hyaluronidase, and neuropeptide precursors. The amphipathic, alpha-helical ampulexins are among the most abundant venom components. Also prominent are members of the Toll/NF-κB signaling pathway, including proteases Persephone, Snake, Easter, and the Toll receptor ligand Spätzle. We find evidence that venom components are processed following envenomation. The acidic (pH∼4) venom contains unprocessed neuropeptide tachykinin and corazonin precursors and is conspicuously devoid of the corresponding processed, biologically active peptides. Neutralization of venom leads to appearance of mature tachykinin and corazonin, suggesting that the wasp employs precursors as a prolonged time-release strategy within the host brain post-envenomation. Injection of fully processed tachykinin into host cephalic ganglia elicits short-term hypokinesia. Ion channel modifiers and cytolytic toxins are absent in A. compressa venom, which appears to hijack control of the host brain by introducing a "storm" of its own neurochemicals. Our findings deepen understanding of the chemical warfare underlying host-parasitoid interactions and in particular neuromodulatory mechanisms that enable manipulation of host behavior to suit the nutritional needs of opportunistic parasitoid progeny.

Evolution of drug resistance in an antifungal-naive chronic Candida lusitaniae infection.

Management of the limited number of antimicrobials currently available requires the identification of infections that contain drug-resistant isolates and the discovery of factors that promote the evolution of drug resistance. Here, we report a single fungal infection in which we have identified numerous subpopulations that differ in their alleles of a single gene that impacts drug resistance. The diversity at this locus was markedly greater than the reported heterogeneity of alleles conferring antibiotic resistance in bacterial infections. Analysis of genomes from hundreds of Clavispora (Candida) lusitaniae isolates, through individual and pooled isolate sequencing, from a single individual with cystic fibrosis revealed at least 25 nonsynonymous mutations in MRR1, which encodes a transcription factor capable of inducing fluconazole (FLZ) resistance in Candida species. Isolates with high-activity Mrr1 variants were resistant to FLZ due to elevated expression of the MDR1-encoded efflux pump. We found that high Mrr1-regulated Mdr1 activity protected against host and bacterial factors, suggesting drug resistance can be selected for indirectly and perhaps explaining the Mrr1 heterogeneity in this individual who had no prior azole exposure. Regional analysis of C. lusitaniae populations from the upper and lower lobes of the right lung suggested intermingling of subpopulations throughout. Our retrospective characterization of sputum and lung populations by pooled sequencing found that alleles that confer FLZ resistance were a minority in each pool, possibly explaining why they were undetected before unsuccessful FLZ therapy. New susceptibility testing regimes may detect problematical drug-resistant subpopulations in heterogeneous single-species infections.

Clustering analysis of large-scale phenotypic data in the model filamentous fungus Neurospora crassa.

BACKGROUND: With 9730 protein-coding genes and a nearly complete gene knockout strain collection, Neurospora crassa is a major model organism for filamentous fungi. Despite this abundance of information, the phenotypes of these gene knockout mutants have not been categorized to determine whether there are broad correlations between phenotype and any genetic features. RESULTS: Here, we analyze data for 10 different growth or developmental phenotypes that have been obtained for 1168 N. crassa knockout mutants. Of these mutants, 265 (23%) are in the normal range, while 903 (77%) possess at least one mutant phenotype. With the exception of unclassified functions, the distribution of functional categories for genes in the mutant dataset mirrors that of the N. crassa genome. In contrast, most genes do not possess a yeast ortholog, suggesting that our analysis will reveal functions that are not conserved in Saccharomyces cerevisiae. To leverage the phenotypic data to identify pathways, we used weighted Partitioning Around Medoids (PAM) approach with 40 clusters. We found that genes encoding metabolic, transmembrane and protein phosphorylation-related genes are concentrated in subsets of clusters. Results from K-Means clustering of transcriptomic datasets showed that most phenotypic clusters contain multiple expression profiles, suggesting that co-expression is not generally observed for genes with shared phenotypes. Analysis of yeast orthologs of genes that co-clustered in MAPK signaling cascades revealed potential networks of interacting proteins in N. crassa. CONCLUSIONS: Our results demonstrate that clustering analysis of phenotypes is a promising tool for generating new hypotheses regarding involvement of genes in cellular pathways in N. crassa. Furthermore, information about gene clusters identified in N. crassa should be applicable to other filamentous fungi, including saprobes and pathogens.

Belowground impacts of alpine woody encroachment are determined by plant traits, local climate, and soil conditions.

Global climate and land use change are causing woody plant encroachment in arctic, alpine, and arid/semi-arid ecosystems around the world, yet our understanding of the belowground impacts of this phenomenon is limited. We conducted a globally distributed field study of 13 alpine sites across four continents undergoing woody plant encroachment and sampled soils from both woody encroached and nearby herbaceous plant community types. We found that woody plant encroachment influenced soil microbial richness and community composition across sites based on multiple factors including woody plant traits, site level climate, and abiotic soil conditions. In particular, root symbiont type was a key determinant of belowground effects, as Nitrogen-fixing woody plants had higher soil fungal richness, while Ecto/Ericoid mycorrhizal species had higher soil bacterial richness and symbiont types had distinct soil microbial community composition. Woody plant leaf traits indirectly influenced soil microbes through their impact on soil abiotic conditions, primarily soil pH and C:N ratios. Finally, site-level climate affected the overall magnitude and direction of woody plant influence, as soil fungal and bacterial richness were either higher or lower in woody encroached versus herbaceous soils depending on mean annual temperature and precipitation. All together, these results document global impacts of woody plant encroachment on soil microbial communities, but highlight that multiple biotic and abiotic pathways must be considered to scale up globally from site- and species-level patterns. Considering both the aboveground and belowground effects of woody encroachment will be critical to predict future changes in alpine ecosystem structure and function and subsequent feedbacks to the global climate system.

Tracking the genome-wide outcomes of a transposable element burst over decades of amplification.

To understand the success strategies of transposable elements (TEs) that attain high copy numbers, we analyzed two pairs of rice (Oryza sativa) strains, EG4/HEG4 and A119/A123, undergoing decades of rapid amplification (bursts) of the class 2 autonomous Ping element and the nonautonomous miniature inverted repeat transposable element (MITE) mPing Comparative analyses of whole-genome sequences of the two strain pairs validated that each pair has been maintained for decades as inbreds since divergence from their respective last common ancestor. Strains EG4 and HEG4 differ by fewer than 160 SNPs and a total of 264 new mPing insertions. Similarly, strains A119 and A123 exhibited about half as many SNPs (277) as new mPing insertions (518). Examination of all other potentially active TEs in these genomes revealed only a single new insertion out of ∼40,000 loci surveyed. The virtual absence of any new TE insertions in these strains outside the mPing bursts demonstrates that the Ping/mPing family gradually attains high copy numbers by maintaining activity and evading host detection for dozens of generations. Evasion is possible because host recognition of mPing sequences appears to have no impact on initiation or maintenance of the burst. Ping is actively transcribed, and both Ping and mPing can transpose despite methylation of terminal sequences. This finding suggests that an important feature of MITE success is that host recognition does not lead to the silencing of the source of transposase.

Probe-Based Multiplex Real-Time PCR as a Diagnostic Tool to Distinguish Distinct Fungal Symbionts Associated With Euwallacea kuroshio and Euwallacea whitfordiodendrus in California.

California has been invaded by two distinct Euwallacea spp. that vector unique plant pathogenic symbiotic fungi on multiple hosts and cause Fusarium dieback. The objective of this study was to develop multiplex real-time quantitative PCR assays using hydrolysis probes targeting the ß-tubulin gene to detect, distinguish, and quantify fungi associated with the polyphagous shot hole borer (PSHB; Euwallacea whitfordiodendrus, Fusarium euwallaceae, Graphium euwallaceae, and Paracremonium pembeum) as well as the Kuroshio shot hole borer (KSHB; Euwallacea kuroshio, Fusarium kuroshium, and Graphium kuroshium) from various sample types. Absolute quantification reaction efficiencies ranged from 88.2 to 104.3%, with a coefficient of determination >0.992 and a limit of detection of 100 copies µl-1 for all targets across both assays. Qualitative detection using the real-time assays on artificially inoculated avocado shoot extracts showed more sensitivity compared with conventional fungal isolation from wood. All symbiotic fungi, except P. pembeum, from PSHB and KSHB female heads were detectable and quantified. Field samples from symptomatic Platanus racemosa, Populus spp., and Salix spp. across 17 of 26 city parks were positively identified as PSHB and KSHB through detection of their symbiotic fungi, and both were found occurring together on five trees from three different park locations. The molecular assays presented here can be utilized to accurately identify fungi associated with these invasive pests in California.

Resolving the Mortierellaceae phylogeny through synthesis of multi-gene phylogenetics and phylogenomics.

Early efforts to classify Mortierellaceae were based on macro- and micromorphology, but sequencing and phylogenetic studies with ribosomal DNA (rDNA) markers have demonstrated conflicting taxonomic groupings and polyphyletic genera. Although some taxonomic confusion in the family has been clarified, rDNA data alone is unable to resolve higher level phylogenetic relationships within Mortierellaceae. In this study, we applied two parallel approaches to resolve the Mortierellaceae phylogeny: low coverage genome (LCG) sequencing and high-throughput, multiplexed targeted amplicon sequencing to generate sequence data for multi-gene phylogenetics. We then combined our datasets to provide a well-supported genome-based phylogeny having broad sampling depth from the amplicon dataset. Resolving the Mortierellaceae phylogeny into monophyletic groups led to the definition of 14 genera, 7 of which are newly proposed. Low-coverage genome sequencing proved to be a relatively cost-effective means of generating a well-resolved phylogeny. The multi-gene phylogenetics approach enabled much greater sampling depth and breadth than the LCG approach, but was unable to resolve higher-level organization of groups. We present this work to resolve some of the taxonomic confusion and provide a genus-level framework to empower future studies on Mortierellaceae diversity, biology, and evolution.