amPEPpy 1.0: a portable and accurate antimicrobial peptide prediction tool.

SUMMARY: Antimicrobial peptides (AMPs) are promising alternative antimicrobial agents. Currently, however, portable, user-friendly and efficient methods for predicting AMP sequences from genome-scale data are not readily available. Here we present amPEPpy, an open-source, multi-threaded command-line application for predicting AMP sequences using a random forest classifier. AVAILABILITY AND IMPLEMENTATION: amPEPpy is implemented in Python 3 and is freely available through GitHub (https://github.com/tlawrence3/amPEPpy). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Ectomycorrhizal Fungus Laccaria bicolor Produces Lipochitooligosaccharides and Uses the Common Symbiosis Pathway to Colonize Populus Roots.

Mycorrhizal fungi form mutualistic associations with the roots of most land plants and provide them with mineral nutrients from the soil in exchange for fixed carbon derived from photosynthesis. The common symbiosis pathway (CSP) is a conserved molecular signaling pathway in all plants capable of associating with arbuscular mycorrhizal fungi. It is required not only for arbuscular mycorrhizal symbiosis but also for rhizobia-legume and actinorhizal symbioses. Given its role in such diverse symbiotic associations, we hypothesized that the CSP also plays a role in ectomycorrhizal associations. We showed that the ectomycorrhizal fungus Laccaria bicolor produces an array of lipochitooligosaccharides (LCOs) that can trigger both root hair branching in legumes and, most importantly, calcium spiking in the host plant Populus in a CASTOR/POLLUX-dependent manner. Nonsulfated LCOs enhanced lateral root development in Populus in a calcium/calmodulin-dependent protein kinase (CCaMK)-dependent manner, and sulfated LCOs enhanced the colonization of Populus by L. bicolor Compared with the wild-type Populus, the colonization of CASTOR/POLLUX and CCaMK RNA interference lines by L. bicolor was reduced. Our work demonstrates that similar to other root symbioses, L. bicolor uses the CSP for the full establishment of its mutualistic association with Populus.

Towards engineering ectomycorrhization into switchgrass bioenergy crops via a lectin receptor-like kinase.

Soil-borne microbes can establish compatible relationships with host plants, providing a large variety of nutritive and protective compounds in exchange for photosynthesized sugars. However, the molecular mechanisms mediating the establishment of these beneficial relationships remain unclear. Our previous genetic mapping and whole-genome resequencing studies identified a gene deletion event of a Populus trichocarpa lectin receptor-like kinase gene PtLecRLK1 in Populus deltoides that was associated with poor-root colonization by the ectomycorrhizal fungus Laccaria bicolor. By introducing PtLecRLK1 into a perennial grass known to be a non-host of L. bicolor, switchgrass (Panicum virgatum L.), we found that L. bicolor colonizes ZmUbipro-PtLecRLK1 transgenic switchgrass roots, which illustrates that the introduction of PtLecRLK1 has the potential to convert a non-host to a host of L. bicolor. Furthermore, transcriptomic and proteomic analyses on inoculated-transgenic switchgrass roots revealed genes/proteins overrepresented in the compatible interaction and underrepresented in the pathogenic defence pathway, consistent with the view that pathogenic defence response is down-regulated during compatible interaction. Metabolomic profiling revealed that root colonization in the transgenic switchgrass was associated with an increase in N-containing metabolites and a decrease in organic acids, sugars, and aromatic hydroxycinnamate conjugates, which are often seen in the early steps of establishing compatible interactions. These studies illustrate that PtLecRLK1 is able to render a plant susceptible to colonization by the ectomycorrhizal fungus L. bicolor and shed light on engineering mycorrhizal symbiosis into a non-host to enhance plant productivity and fitness on marginal lands.

Arabidopsis C-terminal binding protein ANGUSTIFOLIA modulates transcriptional co-regulation of MYB46 and WRKY33.

The apparent antagonism between salicylic acid (SA) and jasmonic acid (JA)/ethylene (ET) signalling resulting in trade-offs between defence against (hemi)biotrophic and necrotrophic pathogens has been widely described across multiple plant species. However, the underlying mechanism remains to be fully established. The molecular and cellular functions of ANGUSTIFOLIA (AN) were characterised, and its role in regulating the pathogenic response was studied in Arabidopsis. We demonstrated that AN, a plant homologue of mammalian C-TERMINAL BINDING PROTEIN (CtBP), antagonistically regulates plant resistance to the hemibiotrophic pathogen Pseudomonas syringae and the necrotrophic pathogen Botrytis cinerea. Consistent with phenotypic observations, transcription of genes involved in SA and JA/ET pathways was antagonistically regulated by AN. By interacting with another nuclear protein TYROSYL-DNA PHOSPHODIESTERASE1 (TDP1), AN imposes transcriptional repression on MYB46, encoding a transcriptional activator of PHENYLALANINE AMMONIA-LYASE (PAL) genes which are required for SA biosynthesis, while releasing TDP1-imposed transcriptional repression on WRKY33, a master regulator of the JA/ET signalling pathway. These findings demonstrate that transcriptional co-regulation of MYB46 and WRKY33 by AN mediates the coordination of SA and JA/ET pathways to optimise defences against (hemi)biotrophic and necrotrophic pathogens.

Fungal Endophytes of Populus trichocarpa Alter Host Phenotype, Gene Expression, and Rhizobiome Composition.

Mortierella and Ilyonectria genera include common species of soil fungi that are frequently detected as root endophytes in many plants, including Populus spp. However, the ecological roles of these and other endophytic fungi with respect to plant growth and function are still not well understood. The functional ecology of two key taxa from the P. trichocarpa rhizobiome, M. elongata PMI93 and I. europaea PMI82, was studied by coupling forest soil bioassays with environmental metatranscriptomics. Using soil bioassay experiments amended with fungal inoculants, M. elongata was observed to promote the growth of P. trichocarpa. This response was cultivar independent. In contrast, I. europaea had no visible effect on P. trichocarpa growth. Metatranscriptomic studies revealed that these fungi impacted rhizophytic and endophytic activities in P. trichocarpa and induced shifts in soil and root microbial communities. Differential expression of core genes in P. trichocarpa roots was observed in response to both fungal species. Expression of P. trichocarpa genes for lipid signaling and nutrient uptake were upregulated, and expression of genes associated with gibberellin signaling were altered in plants inoculated with M. elongata, but not I. europaea. Upregulation of genes for growth promotion, downregulation of genes for several leucine-rich repeat receptor kinases, and alteration of expression of genes associated with plant defense responses (e.g., jasmonic acid, salicylic acid, and ethylene signal pathways) also suggest that M. elongata manipulates plant defenses while promoting plant growth.

Russulaceae: a new genomic dataset to study ecosystem function and evolutionary diversification of ectomycorrhizal fungi with their tree associates.

The family Russulaceae is considered an iconic lineage of mostly mushroom-forming basidiomycetes due to their importance as edible mushrooms in many parts of the world, and their ubiquity as ectomycorrhizal symbionts in both temperate and tropical forested biomes. Although much research has been focused on this group, a comprehensive or cohesive synthesis by which to understand the functional diversity of the group has yet to develop. Interest in ectomycorrhizal fungi, of which Russulaceae is a key lineage, is prodigious due to the important roles they play as plant root mutualists in ecosystem functioning, global carbon sequestration, and a potential role in technology development toward environmental sustainability. As one of the most species-diverse ectomycorrhizal lineages, the Russulaceae has recently been the focus of a dense sampling and genome sequencing initiative with the Joint Genome Institute aimed at untangling their functional roles and testing whether functional niche specialization exists for independent lineages of ectomycorrhizal fungi. Here we present a review of important studies on this group to contextualize what we know about its members' evolutionary history and ecosystem functions, as well as to generate hypotheses establishing the Russulaceae as a valuable experimental system.

Integrated proteomics and metabolomics suggests symbiotic metabolism and multimodal regulation in a fungal-endobacterial system.

Many plant-associated fungi host endosymbiotic endobacteria with reduced genomes. While endobacteria play important roles in these tri-partite plant-fungal-endobacterial systems, the active physiology of fungal endobacteria has not been characterized extensively by systems biology approaches. Here, we use integrated proteomics and metabolomics to characterize the relationship between the endobacterium Mycoavidus sp. and the root-associated fungus Mortierella elongata. In nitrogen-poor media, M. elongata had decreased growth but hosted a large and growing endobacterial population. The active endobacterium likely extracted malate from the fungal host as the primary carbon substrate for energy production and biosynthesis of phospho-sugars, nucleobases, peptidoglycan and some amino acids. The endobacterium obtained nitrogen by importing a variety of nitrogen-containing compounds. Further, nitrogen limitation significantly perturbed the carbon and nitrogen flows in the fungal metabolic network. M. elongata regulated many pathways by concordant changes on enzyme abundances, post-translational modifications, reactant concentrations and allosteric effectors. Such multimodal regulations may be a general mechanism for metabolic modulation.

Genome-wide analysis of lectin receptor-like kinases in Populus.

BACKGROUND: Receptor-like kinases (RLKs) belong to a large protein family with over 600 members in Arabidopsis and over 1000 in rice. Among RLKs, the lectin receptor-like kinases (LecRLKs) possess a characteristic extracellular carbohydrate-binding lectin domain and play important roles in plant development and innate immunity. There are 75 and 173 LecRLKs in Arabidopsis and rice, respectively. However, little is known about LecRLKs in perennial woody plants. RESULTS: Here we report the genome-wide analysis of classification, domain architecture and expression of LecRLKs in the perennial woody model plant Populus. We found that the LecRLK family has expanded in Populus to a total of 231, including 180 G-type, 50 L-type and 1 C-type LecRLKs. Expansion of the Populus LecRLKs (PtLecRLKs) occurred partially through tandem duplication. Based on domain architecture and orientation features, we classified PtLecRLKs into eight different classes. RNA-seq-based transcriptomics analysis revealed diverse expression patterns of PtLecRLK genes among leaves, stems, roots, buds and reproductive tissues and organs. CONCLUSIONS: This study offers a comprehensive view of LecRLKs in the perennial woody model plant Populus and provides a foundation for functional characterization of this important family of receptor-like kinases.

Extensive gene flow over Europe and possible speciation over Eurasia in the ectomycorrhizal basidiomycete Laccaria amethystina complex.

Biogeographical patterns and large-scale genetic structure have been little studied in ectomycorrhizal (EM) fungi, despite the ecological and economic importance of EM symbioses. We coupled population genetics and phylogenetic approaches to understand spatial structure in fungal populations on a continental scale. Using nine microsatellite markers, we characterized gene flow among 16 populations of the widespread EM basidiomycete Laccaria amethystina over Europe (i.e. over 2900 km). We also widened our scope to two additional populations from Japan (10(4) km away) and compared them with European populations through microsatellite markers and multilocus phylogenies, using three nuclear genes (NAR, G6PD and ribosomal DNA) and two mitochondrial ribosomal genes. European L. amethystina populations displayed limited differentiation (average F(ST) = 0.041) and very weak isolation by distance (IBD). This panmictic European pattern may result from effective aerial dispersal of spores, high genetic diversity in populations and mutualistic interactions with multiple hosts that all facilitate migration. The multilocus phylogeny based on nuclear genes confirmed that Japanese and European specimens were closely related but clustered on a geographical basis. By using microsatellite markers, we found that Japanese populations were strongly differentiated from the European populations (F(ST) = 0.416), more than expected by extrapolating the European pattern of IBD. Population structure analyses clearly separated the populations into two clusters, i.e. European and Japanese clusters. We discuss the possibility of IBD in a continuous population (considering some evidence for a ring species over the Northern Hemisphere) vs. an allopatric speciation over Eurasia, making L. amethystina a promising model of intercontinental species for future studies.

Survey and analysis of simple sequence repeats in the Laccaria bicolor genome, with development of microsatellite markers.

It is becoming clear that simple sequence repeats (SSRs) play a significant role in fungal genome organization, and they are a large source of genetic markers for population genetics and meiotic maps. We identified SSRs in the Laccaria bicolor genome by in silico survey and analyzed their distribution in the different genomic regions. We also compared the abundance and distribution of SSRs in L. bicolor with those of the following fungal genomes: Phanerochaete chrysosporium, Coprinopsis cinerea, Ustilago maydis, Cryptococcus neoformans, Aspergillus nidulans, Magnaporthe grisea, Neurospora crassa and Saccharomyces cerevisiae. Using the MISA computer program, we detected 277,062 SSRs in the L. bicolor genome representing 8% of the assembled genomic sequence. Among the analyzed basidiomycetes, L. bicolor exhibited the highest SSR density although no correlation between relative abundance and the genome sizes was observed. In most genomes the short motifs (mono- to trinucleotides) were more abundant than the longer repeated SSRs. Generally, in each organism, the occurrence, relative abundance, and relative density of SSRs decreased as the repeat unit increased. Furthermore, each organism had its own common and longest SSRs. In the L. bicolor genome, most of the SSRs were located in intergenic regions (73.3%) and the highest SSR density was observed in transposable elements (TEs; 6,706 SSRs/Mb). However, 81% of the protein-coding genes contained SSRs in their exons, suggesting that SSR polymorphism may alter gene phenotypes. Within a L. bicolor offspring, sequence polymorphism of 78 SSRs was mainly detected in non-TE intergenic regions. Unlike previously developed microsatellite markers, these new ones are spread throughout the genome; these markers could have immediate applications in population genetics.

Phylogenetic diversity of 200+ isolates of the ectomycorrhizal fungus Cenococcum geophilum associated with Populus trichocarpa soils in the Pacific Northwest, USA and comparison to globally distributed representatives.

The ectomycorrhizal fungal symbiont Cenococcum geophilum is of high interest as it is globally distributed, associates with many plant species, and has resistance to multiple environmental stressors. C. geophilum is only known from asexual states but is often considered a cryptic species complex, since extreme phylogenetic divergence is often observed within nearly morphologically identical strains. Alternatively, C. geophilum may represent a highly diverse single species, which would suggest cryptic but frequent recombination. Here we describe a new isolate collection of 229 C. geophilum isolates from soils under Populus trichocarpa at 123 collection sites spanning a ~283 mile north-south transect in Western Washington and Oregon, USA (PNW). To further understanding of the phylogenetic relationships within C. geophilum, we performed maximum likelihood and Bayesian phylogenetic analyses to assess divergence within the PNW isolate collection, as well as a global phylogenetic analysis of 789 isolates with publicly available data from the United States, Japan, and European countries. Phylogenetic analyses of the PNW isolates revealed three distinct phylogenetic groups, with 15 clades that strongly resolved at >80% bootstrap support based on a GAPDH phylogeny and one clade segregating strongly in two principle component analyses. The abundance and representation of PNW isolate clades varied greatly across the North-South range, including a monophyletic group of isolates that spanned nearly the entire gradient at ~250 miles. A direct comparison between the GAPDH and ITS rRNA gene region phylogenies, combined with additional analyses revealed stark incongruence between the ITS and GAPDH gene regions, consistent with intra-species recombination between PNW isolates. In the global isolate collection phylogeny, 34 clades were strongly resolved using Maximum Likelihood and Bayesian approaches (at >80% MLBS and >0.90 BPP respectively), with some clades having intra- and intercontinental distributions. Together these data are highly suggestive of divergence within multiple cryptic species, however additional analyses such as higher resolution genotype-by-sequencing approaches are needed to distinguish potential species boundaries and the mode and tempo of recombination patterns.

Biosynthesis and characterization of deuterated chitosan in filamentous fungus and yeast.

Deuterated chitosan was produced from the filamentous fungus Rhizopus oryzae, cultivated with deuterated glucose in H2O medium, without the need for conventional chemical deacetylation. After extraction and purification, the chemical composition and structure were determined by Fourier-transform infrared spectroscopy (FTIR), nuclear magnetic resonance (NMR), and small-angle neutron scattering (SANS). 13C NMR experiments provided additional information about the position of the deuterons in the glucoseamine backbone. The NMR spectra indicated that the deuterium incorporation at the non-exchangeable hydrogen positions of the aminoglucopyranosyl ring in the C3 - C5 positions was at least 60-80 %. However, the C2 position was deuterated at a much lower level (6%). Also, SANS showed that the structure of deuterated chitosan was very similar compared to the non-deuterated counterpart. The most abundant radii of the protiated and deuterated chitosan fibers were 54 Å and 60 Å, respectively, but there is a broader distribution of fiber radii in the protiated chitosan sample. The highly deuterated, soluble fungal chitosan described here can be used as a model material for studying chitosan-enzyme complexes for future neutron scattering studies. Because the physical behavior of non-deuterated fungal chitosan mimicked that of shrimp shell chitosan, the methods presented here represent a new approach to producing a high quality deuterated non-animal-derived aminopolysaccharide for studying the structure-function association of biocomposite materials in drug delivery, tissue engineering and other bioactive chitosan-based composites.

Microbe to Microbiome: A Paradigm Shift in the Application of Microorganisms for Sustainable Agriculture.

Light, water and healthy soil are three essential natural resources required for agricultural productivity. Industrialization of agriculture has resulted in intensification of cropping practices using enormous amounts of chemical pesticides and fertilizers that damage these natural resources. Therefore, there is a need to embrace agriculture practices that do not depend on greater use of fertilizers and water to meet the growing demand of global food requirements. Plants and soil harbor millions of microorganisms, which collectively form a microbial community known as the microbiome. An effective microbiome can offer benefits to its host, including plant growth promotion, nutrient use efficiency, and control of pests and phytopathogens. Therefore, there is an immediate need to bring functional potential of plant-associated microbiome and its innovation into crop production. In addition to that, new scientific methodologies that can track the nutrient flux through the plant, its resident microbiome and surrounding soil, will offer new opportunities for the design of more efficient microbial consortia design. It is now increasingly acknowledged that the diversity of a microbial inoculum is as important as its plant growth promoting ability. Not surprisingly, outcomes from such plant and soil microbiome studies have resulted in a paradigm shift away from single, specific soil microbes to a more holistic microbiome approach for enhancing crop productivity and the restoration of soil health. Herein, we have reviewed this paradigm shift and discussed various aspects of benign microbiome-based approaches for sustainable agriculture.

Advancing How We Learn from Biodesign to Mitigate Risks with Next-Generation Genome Engineering.

In the last decade, the unprecedented simplicity and flexibility of the CRISPR-Cas system has made it the dominant transformative tool in gene and genome editing. However, this democratized technology is both a boon and a bane, for which we have yet to understand the full potential to investigate and rewrite genomes (also named "genome biodesign"). Rapid CRISPR advances in a range of applications in basic research, agriculture, and clinical applications pose new risks and raise several biosecurity concerns. In such a fast-moving field of research, we emphasize the importance of properly communicating the quality and accuracy of results and recommend new reporting requirements for results derived from next-generation genome engineering.

Insight into a highly polymorphic endophyte isolated from the roots of the halophytic seepweed Suaeda salsa: Laburnicola rhizohalophila sp. nov. (Didymosphaeriaceae, Pleosporales).

We surveyed root endophytic fungi of the coastal halophyte Suaeda salsa and detected a population of a novel species that we described here as Laburnicola rhizohalophila sp. nov. No sexual sporulating structure was observed. Instead, it produced a large amount of thalloconidia, 0-1 transverse septa, hyaline to darkly pigmented, often peanut-shaped and sometimes dumbbell-shaped, both ends enlarged with numerous oil droplets inside the hyphal cells. Surprisingly, a high degree of phenotypic and physiological intraspecific variation (e.g., salinity tolerance, growth under different carbon:nitrogen ratios, and carbon utilization pattern) was recorded. The inoculation test indicated that the isolates could successfully infect host roots and form microsclerotia-like structures in cortical cells, a typical trait of dark septate endophytes (DSEs). Furthermore, most isolates were shown to promote host seedling growth. To evaluate conspecificity and infer its phylogenetic affinity, multiloci data including nuclear rRNA loci (ITS1 and 2, partial 28S), partial RNA Polymerase II second-largest subunit (rpb2), and partial translation elongation factor-1α (tef1) were characterized. Genealogical concordance phylogenetic species recognition (GCPSR) detected a genetically isolated clade of L. rhizohalophila within the Pleosporales in the Didymosphaeriaceae. Maximum likelihood phylogenetic reconstruction revealed that the endophytic fungus was genetically close to Laburnicoladactylidis but separated by a relatively long genetic distance. Our work highlights that the pleosporalean taxa might represent an underexplored reservoir of root DSEs.

Lipo-chitooligosaccharides as regulatory signals of fungal growth and development.

Lipo-chitooligosaccharides (LCOs) are signaling molecules produced by rhizobial bacteria that trigger the nodulation process in legumes, and by some fungi that also establish symbiotic relationships with plants, notably the arbuscular and ecto mycorrhizal fungi. Here, we show that many other fungi also produce LCOs. We tested 59 species representing most fungal phyla, and found that 53 species produce LCOs that can be detected by functional assays and/or by mass spectroscopy. LCO treatment affects spore germination, branching of hyphae, pseudohyphal growth, and transcription in non-symbiotic fungi from the Ascomycete and Basidiomycete phyla. Our findings suggest that LCO production is common among fungi, and LCOs may function as signals regulating fungal growth and development.

Plant Biosystems Design Research Roadmap 1.0.

Human life intimately depends on plants for food, biomaterials, health, energy, and a sustainable environment. Various plants have been genetically improved mostly through breeding, along with limited modification via genetic engineering, yet they are still not able to meet the ever-increasing needs, in terms of both quantity and quality, resulting from the rapid increase in world population and expected standards of living. A step change that may address these challenges would be to expand the potential of plants using biosystems design approaches. This represents a shift in plant science research from relatively simple trial-and-error approaches to innovative strategies based on predictive models of biological systems. Plant biosystems design seeks to accelerate plant genetic improvement using genome editing and genetic circuit engineering or create novel plant systems through de novo synthesis of plant genomes. From this perspective, we present a comprehensive roadmap of plant biosystems design covering theories, principles, and technical methods, along with potential applications in basic and applied plant biology research. We highlight current challenges, future opportunities, and research priorities, along with a framework for international collaboration, towards rapid advancement of this emerging interdisciplinary area of research. Finally, we discuss the importance of social responsibility in utilizing plant biosystems design and suggest strategies for improving public perception, trust, and acceptance.

Increasing access to microfluidics for studying fungi and other branched biological structures.

BACKGROUND: Microfluidic systems are well-suited for studying mixed biological communities for improving industrial processes of fermentation, biofuel production, and pharmaceutical production. The results of which have the potential to resolve the underlying mechanisms of growth and transport in these complex branched living systems. Microfluidics provide controlled environments and improved optical access for real-time and high-resolution imaging studies that allow high-content and quantitative analyses. Studying growing branched structures and the dynamics of cellular interactions with both biotic and abiotic cues provides context for molecule production and genetic manipulations. To make progress in this arena, technical and logistical barriers must be overcome to more effectively deploy microfluidics in biological disciplines. A principle technical barrier is the process of assembling, sterilizing, and hydrating the microfluidic system; the lack of the necessary equipment for the preparatory process is a contributing factor to this barrier. To improve access to microfluidic systems, we present the development, characterization, and implementation of a microfluidics assembly and packaging process that builds on self-priming point-of-care principles to achieve "ready-to-use microfluidics." RESULTS: We present results from domestic and international collaborations using novel microfluidic architectures prepared with a unique packaging protocol. We implement this approach by focusing primarily on filamentous fungi; we also demonstrate the utility of this approach for collaborations on plants and neurons. In this work we (1) determine the shelf-life of ready-to-use microfluidics, (2) demonstrate biofilm-like colonization on fungi, (3) describe bacterial motility on fungal hyphae (fungal highway), (4) report material-dependent bacterial-fungal colonization, (5) demonstrate germination of vacuum-sealed Arabidopsis seeds in microfluidics stored for up to 2 weeks, and (6) observe bidirectional cytoplasmic streaming in fungi. CONCLUSIONS: This pre-packaging approach provides a simple, one step process to initiate microfluidics in any setting for fungal studies, bacteria-fungal interactions, and other biological inquiries. This process improves access to microfluidics for controlling biological microenvironments, and further enabling visual and quantitative analysis of fungal cultures.

Terpene Synthase Genes Originated from Bacteria through Horizontal Gene Transfer Contribute to Terpenoid Diversity in Fungi.

Fungi are successful eukaryotes of wide distribution. They are known as rich producers of secondary metabolites, especially terpenoids, which are important for fungi-environment interactions. Horizontal gene transfer (HGT) is an important mechanism contributing to genetic innovation of fungi. However, it remains unclear whether HGT has played a role in creating the enormous chemical diversity of fungal terpenoids. Here we report that fungi have acquired terpene synthase genes (TPSs), which encode pivotal enzymes for terpenoid biosynthesis, from bacteria through HGT. Phylogenetic analysis placed the majority of fungal and bacterial TPS genes from diverse taxa into two clades, indicating ancient divergence. Nested in the bacterial TPS clade is a number of fungal TPS genes that are inferred as the outcome of HGT. These include a monophyletic clade of nine fungal TPS genes, designated as BTPSL for bacterial TPS-like genes, from eight species of related entomopathogenic fungi, including seven TPSs from six species in the genus Metarhizium. In vitro enzyme assays demonstrate that all seven BTPSL genes from the genus Metarhizium encode active enzymes with sesquiterpene synthase activities of two general product profiles. By analyzing the catalytic activity of two resurrected ancestral BTPSLs and one closely related bacterial TPS, the trajectory of functional evolution of BTPSLs after HGT from bacteria to fungi and functional divergence within Metarhizium could be traced. Using M. brunneum as a model species, both BTPSLs and typical fungal TPSs were demonstrated to be involved in the in vivo production of terpenoids, illustrating the general importance of HGT of TPS genes from bacteria as a mechanism contributing to terpenoid diversity in fungi.