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

Main menu

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

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Ecological and Evolutionary Science

Comparative Genomics Reveals Prophylactic and Catabolic Capabilities of Actinobacteria within the Fungus-Farming Termite Symbiosis

Robert Murphy, René Benndorf, Z. Wilhelm de Beer, John Vollmers, Anne-Kristin Kaster, Christine Beemelmanns, Michael Poulsen
Gary Sawers, Editor
Robert Murphy
aUniversity of Copenhagen, Department of Biology, Section for Ecology and Evolution, Copenhagen East, Denmark
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Robert Murphy
René Benndorf
bLeibniz Institute for Natural Product Research and Infection Biology, Hans Knöll Institute, Jena, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Z. Wilhelm de Beer
cDepartment of Microbiology and Plant Pathology, Forestry and Agriculture Biotechnology Institute, University of Pretoria, Pretoria, South Africa
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Z. Wilhelm de Beer
John Vollmers
dInstitute for Biological Interfaces (IBG 5), Karlsruhe Institute of Technology, Eggenstein-Leopoldshafen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for John Vollmers
Anne-Kristin Kaster
dInstitute for Biological Interfaces (IBG 5), Karlsruhe Institute of Technology, Eggenstein-Leopoldshafen, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Anne-Kristin Kaster
Christine Beemelmanns
bLeibniz Institute for Natural Product Research and Infection Biology, Hans Knöll Institute, Jena, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Christine Beemelmanns
Michael Poulsen
aUniversity of Copenhagen, Department of Biology, Section for Ecology and Evolution, Copenhagen East, Denmark
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael Poulsen
Gary Sawers
Martin Luther University of Halle-Wittenberg Institute of Biology/Microbiology
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphere.01233-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Actinobacteria, one of the largest bacterial phyla, are ubiquitous in many of Earth’s ecosystems and often act as defensive symbionts with animal hosts. Members of the phylum have repeatedly been isolated from basidiomycete-cultivating fungus-farming termites that maintain a monoculture fungus crop on macerated dead plant substrate. The proclivity for antimicrobial and enzyme production of Actinobacteria make them likely contributors to plant decomposition and defense in the symbiosis. To test this, we analyzed the prophylactic (biosynthetic gene cluster [BGC]) and metabolic (carbohydrate-active enzyme [CAZy]) potential in 16 (10 existing and six new genomes) termite-associated Actinobacteria and compared these to the soil-dwelling close relatives. Using antiSMASH, we identified 435 BGCs, of which 329 (65 unique) were similar to known compound gene clusters, while 106 were putatively novel, suggesting ample prospects for novel compound discovery. BGCs were identified among all major compound categories, including 26 encoding the production of known antimicrobial compounds, which ranged in activity (antibacterial being most prevalent) and modes of action that might suggest broad defensive potential. Peptide pattern recognition analysis revealed 823 (43 unique) CAZymes coding for enzymes that target key plant and fungal cell wall components (predominantly chitin, cellulose, and hemicellulose), confirming a substantial degradative potential of these bacteria. Comparison of termite-associated and soil-dwelling bacteria indicated no significant difference in either BGC or CAZy potential, suggesting that the farming termite hosts may have coopted these soil-dwelling bacteria due to their metabolic potential but that they have not been subject to genome change associated with symbiosis.

IMPORTANCE Actinobacteria have repeatedly been isolated in fungus-farming termites, and our genome analyses provide insights into the potential roles they may serve in defense and for plant biomass breakdown. These insights, combined with their relatively higher abundances in fungus combs than in termite gut, suggest that they are more likely to play roles in fungus combs than in termite guts. Up to 25% of the BGCs we identify have no similarity to known clusters, indicating a large potential for novel chemistry to be discovered. Similarities in metabolic potential of soil-dwelling and termite-associated bacteria suggest that they have environmental origins, but their consistent presence with the termite system suggests their importance for the symbiosis.

INTRODUCTION

Actinobacteria is a Gram-positive bacterial phylum that represents one of the largest bacterial clades. Members of the Actinobacteria are adapted to a range of environmental conditions and produce a variety of extracellular enzymes and natural products, many of the latter serving as drugs or drug leads (1). Actinobacteria associate with a diverse set of eukaryotic hosts, particularly insects, serving protective roles through the production of antimicrobials (2–7). Well-known examples include the European beewolves (genus Philanthus) that host antifungal-producing species in the antennae to protect their larvae from fungal infection (8, 9), and members of the genus Pseudonocardia in New World fungus-farming ants that help defend from specialized mycoparasites of the ants’ fungal mutualism (10–12). Actinobacteria have also repeatedly been isolated from the Old World fungus-farming termite symbiosis (13–15), but their potential symbiotic roles have remained elusive.

Termites in the subfamily Macrotermitinae (Termitidae: Blattodea) cultivate basidiomycete fungi in the genus Termitomyces (Agaricales: Lyophyllaceae) as their sole food source, using plant material foraged on by the termites (16). The termites, through the symbiosis, manage to fully utilize plant substrates (17) and maintain monoculture fungal farms in densely populated colonies of often millions of termite workers without apparent problems with infectious disease (18). Older workers collect plant biomass (Fig. 1A.1), which is brought back to the nest, where younger workers ingest the substrate along with asexual spores of Termitomyces produced in specialized nodules (Fig. 1A.2) and deposit the fecal matter as fresh comb. After this first gut passage, Termitomyces grows on the macerated plant substrate within the fungus comb (fungus garden [Fig. 1A.3/4a]). When the plant material is fully utilized, older workers ingest and digest the mature fungus comb (Fig. 1A.4c), after which all organic material is essentially utilized (19). A plethora of bacterial symbionts have been identified from fungus-farming termites, many of which facilitate metabolism of plant and fungal biomass (e.g., Alistipes, Bacteroidetes [20]) and defense against antagonists (e.g., Bacillus [21]). Recent studies have demonstrated the consistent presence of Actinobacteria in both termite guts and fungus combs, and the propensities for enzyme and antimicrobial production make them promising candidate plant decomposers and defensive symbionts.

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

Isolation sources and sites in fungus-farming termite colonies where Actinobacteria could contribute to symbiosis. (A) Biomass route through the fungus-growing termite symbiosis and schematic representation of nest structure. (B) Isolation location by bacterial genera; colors mirror the corresponding location in panel A. Ac, Actinomadura; Am, Amycolatopsis; Lu, Luteimicrobium; My, Mycolicibacterium; No, Nocardia; St, Streptomyces. (C) Relative abundances (box plot showing medians, 1st and 3rd quartiles, and all data points) of the six bacterial genera, to which our genomes belong, in previously published 16S rRNA amplicon sequencing community data sets from 10 M. natalensis colonies, totaling 10 comb and 86 gut samples (14). These were determined by BLASTn of the 16S rRNA in the termite-associated Actinobacteria genomes against a database of the amplicon sequencing community data. Successful hits were set at genus level (95%) and E value ≤ 0.01. Relative abundances were on average significantly higher in combs than guts (P = 0.0004; Mann-Whitney U test). Note that DNA extraction protocols tend to bias against Gram-positive bacteria; thus, actual relative abundances are conceivably higher than reported (42).

Actinobacterial contributions to defense and plant biomass decomposition in this symbiosis are most likely to occur within combs for three main reasons (15, 22, 23). First, the initial gut passage is rapid (24) and involves minimal plant decomposition (17; but see reference 25 regarding lignin cleavage), so enzyme contributions from Actinobacteria during gut passage are unlikely to be important. Second, although gut bacteria would be in a prime position to suppress competitors or antagonists of Termitomyces or the termites present in the plant substrate (26), recent findings in Macrotermes bellicosus refuted this “gut sanitation” hypothesis, as potential fungal antagonists appear to pass the gut unharmed and indeed enter fungus combs (27). Lastly, Actinobacteria are consistently present and more abundant in fungus combs (14) than guts (Fig. 1C). Thus, Actinobacteria contributions are conceivably in fungus combs, while their presence in termite guts is more likely due to ingestion from the plant substrate, soil, or the comb (14, 28, 29).

To elucidate putative functions of Actinobacteria in the symbiosis, we comparatively analyze whole-genome sequences of 16 isolates from the fungus-farming termite Macrotermes natalensis and characterize their antimicrobial and carbohydrate-active enzyme (CAZyme) potential. These isolates were obtained from either termite workers (cuticle, gut) or fungus combs (Fig. 1B). To place them phylogenetically, and to provide indications of genome adaptations to symbioses, we compare their genomes with those of the closely related soil-dwelling free-living counterparts and analyze their prophylactic and metabolic potential through whole-genome targeted mining.

RESULTS

Phylogenetic placement and genome quality.We first performed multilocus sequence typing (MLST) analyses based on 122 genes to obtain phylogenetic placement of the genomes, and most termite-associated isolates showed short phylogenetic distance to closely related soil species with reference genomes in RefSeq (determined by BLASTn of the 16S rRNA sequence) (Fig. 2). This suggests that the isolates are not specific to being in symbiosis with the termites but more likely originate from the mound or surrounding soil. The quality of the genomes, while varying, was generally very good, with high completeness, high L50 (length lower limit of contigs making up 50% of the assembly values), long largest contigs, and low N50 (number of contigs making up 50% of your assembly) values (Fig. 2; see also Table S1 in the supplemental material).

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

Phylogenetic placement, genome characteristics, CAZymes, and BGCs of the 16 termite-associated Actinobacteria. (A) Consensus maximum-likelihood distance tree based on 122 orthologous gene trees with nonparametric bootstrapping, placing the 16 termite-associated bacteria (seen in bold) alongside closely related soil-dwelling relatives. (B) N50, L50, and genome completeness (see Table S1 for the full results). N50 = number of contigs making 50% of the assembly; L50, length lower limit of contigs making up 50% of the assembly. (C) Log10-transformed abundance of CAZymes identified across the genomes. GH, glycoside hydrolase; GT, glycoside transferase; CBM, carbohydrate-binding module; CE, carbohydrate esterase; AA, auxiliary activities; PL, polysaccharide lyase. For full results, see Table S3. (D) Log10-transformed abundance of BGCs identified across the genomes. NRPS, nonribosomal peptide synthase; PKS, polyketide synthase; hybrid, NRPS-PKS hybrid; terpenes; and RiPPs, ribosomally synthesized and posttranslationally modified peptides. For the full results, see Table S2.

TABLE S1

Termite-associated Actinobacteria genome metadata and statistics. N50, number of contigs making 50% of the assembly; L50, length lower limit of contigs making up 50% of the assembly. Download Table S1, XLSX file, 0.02 MB.
Copyright © 2021 Murphy et al.

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

TABLE S2

(Sheet 1) Predicted BGCs identified by antiSMASH with a similarity to a reported BGC in the MIBiG database of 50% or higher. The class of BGC, the similar compound identified in the MIBiG database, and whether it has documented antimicrobial activity are detailed. The number relating a sample to a BGC indicates the occurrences of this BGC being found in that sample. (Sheet 2) All BGCs identified by antiSMASH, both novel and predicted (regardless of similarity). Download Table S2, XLSX file, 0.05 MB.
Copyright © 2021 Murphy et al.

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

TABLE S3

Overview results of HotPep-identified CAZymes. CAZy prediction summary indicating the EC number and corresponding enzyme name of CAZymes identified by HotPep, along with the enzyme substrate, category of enzyme (hydrolase, oxidoreductases, etc.), and whether the enzyme is catabolic or anabolic in nature. The number relating a sample to a specific enzyme is the frequency count of that enzyme in that sample. Download Table S3, XLSX file, 0.04 MB.
Copyright © 2021 Murphy et al.

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

Biosynthetic potential.A total of 435 biosynthetic gene clusters were identified via whole-genome mining by antiSMASH (30, 31) from termite-associated Actinobacteria, including 329 with a predicted similar gene cluster in the MIBiG database (32) and 106 putatively novel biosynthetic gene clusters (BGCs) (Table S2). Novel is defined here as no (i.e., 0%) homology to any BGCs in the MIBiG database. Of the 329 BGCs, 116 (65 unique) BGCs were ≥50% similar to a gene cluster in the MIBiG database. A similarity of ≥50% suggests likely similar function to the putative homologue in the MIBiG database. There were 24 unique BGCs encoding compounds with and 42 without known antimicrobial activity (Fig. 3; Table S2), based on the current literature. Unique is defined as a distinct cluster identified by antiSMASH. The total number of BGCs encoding the production of compounds with (25) and without (91) known antimicrobial activity differs considerably from the unique numbers (Table S2). This is likely due to the higher frequency of unique BGCs encoding enzymes that produce compounds without antimicrobial activity being present multiple times across the genomes (Fig. 3) or to the antimicrobial activities of many BGCs having yet to be discovered.

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

The 16 termite-associated Actinobacteria code for diverse BGCs, some of which are shared across isolates. Circular representation of the 112 identified BGCs, based on known biosynthetic gene clusters identified by antiSMASH with 50% or greater similarity to gene clusters in the MIBiG database. Each segment represents a unique BGC or strain with a link between them indicating that the BGC was identified in the strain. Names of BGCs with known antimicrobial activity are highlighted in bold. For the full results, see Table S4. (B) Structures of three recently identified secondary metabolites. The T1PKS-derived natural products natalamycin A and efomycin K were identified in the broad-range antifungal strain M56, although their role in the symbiosis remains unknown (59, 60). The T1PKS-derived macrotermycin A, identified from M39, has selective antifungal activity against Pseudoxylaria antagonists in the symbiosis (22).

TABLE S4

All essential genes identified with ARTS and their frequency of marker occurrence (duplication, proximity to a BGC, horizontal gene transfer, and whether they are a known target). Essential genes are ordered into categories. First, “Most likely resistance factor or target” indicating that in an individual genome for this gene, three or more markers were positive, one being proximity to a BGC. Second, “Likely resistance target or factor,” requiring the same as the previous category, just without proximity to a BGC being one of the markers. Lastly, “Unlikely resistance target or factor,” where fewer than three markers were identified in association with this essential gene in any singular genome. Download Table S4, XLSX file, 0.3 MB.
Copyright © 2021 Murphy et al.

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

One BGC that encodes the production of a compound with known activity and 15 BGCs encoding production of compounds without known activity were represented more than once across the isolates; for example, the BGC responsible for the production of ectoine was identified in 12 of the 16 genomes. Only melanin and alkylresorcinol were identified multiple times within a single isolate (twice in M41, RB17, and RB110-1 [henceforth known as RB110] and twice in RB33, respectively). With most BGCs being present in only one isolate, we see the genetic basis for the production of a diverse range of compounds with antimicrobial activities and potentially various modes and ranges of action. Ten of the BGCs encoding the production of natural products with antimicrobial activity have documented antifungal activities, 23 have bacteriostatic or bactericidal effects, two have antiviral effects, and several have antiprotozoal, anticoccidial, and general antiparasitic effects (Table S2). Many modes and ranges of action remained unknown, suggesting ample potential for discoveries of novel bioactivities.

Free-living Actinobacteria displayed similar BGC capacities and distributions, with 21 unique BGCs coding for enzymes that produce compounds with known activity and 41 unique BGCs for enzymes without known activity. The total number of BGCs coding for enzymes that produce compounds with and without known activity was comparable to those of termite-associated genomes (260 with versus 86 without). Three BGCs (melanin, alkylresorcinol, and fluostatins M to Q) were identified more than once in a single genome. However, origin of a genome (i.e., termite-associated versus free living) was not significant in predicting BGC class frequency in generalized linear models (glm) when comparing residual deviances (glm, P = 0.798) (residual deviance = 306.65 [df = 280] versus residual deviance = 306.59 [df = 279]).

A Dunn test for multiple comparisons, aiming to identify different frequencies of BGC class across termite-associated isolates, showed evidence for terpene enrichment compared to every other BGC class, apart from comparing terpenes to nonribosomal peptide synthases (NRPSs). There was no evidence for enrichment of NRPSs compared to other BGC classes. However, apparent terpene enrichment may be an artifact of contig assembly length, which could cause enrichment of terpene clusters due to their shorter lengths than other major BGC classes (Fig. S2) (33). Four of the five most abundant BGCs that passed our 50% threshold indeed encoded terpenes (Table S2). The BGC encoding the production of geosmin, a sesquiterpene alcohol that contributes to the distinctive smell of soil (34), alone accounts for 10 of the 112 BGCs. This is substantial given that multiple occurrences of BGCs are infrequent (only five were observed four or more times). The tricyclic sesquiterpene albaflavenone was the most frequently encoded compound with known antimicrobial activity. Albaflavenone has shown activity against Bacillus subtilis and is conserved in Streptomyces (35, 36).

FIG S2

Length (bp) of BGCs with a similarity score of 50% or greater to a known BGC in the MIBiG database (the same BGCs as in Fig. 3). Download FIG S2, TIF file, 1.8 MB.
Copyright © 2021 Murphy et al.

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

To explore the potential of unidentified antimicrobials present in the genomes, we utilized Antibiotic Resistance Target Seeker 2 (ARTS2), which identified 5,298 (422 unique) essential genes in the 16 termite-associated Actinobacteria (Table S3) through protein family (from the TIGRFAM database) homology to a reference set of core genes from related taxa (37, 38). For a gene to be accepted as a probable resistance factor or target, we determined it would need three of the four associated metadata markers (duplication, proximity to a BGC, known resistance factor/target, and horizontal gene transfer). Following this, 94 (52 unique) genes were categorized as probable resistance factors or targets (Table S3). Among these, 85 (50 unique) were located in close proximity to a BGC identified by antiSMASH. The frequency of metadata markers was not uniform across the genes that passed our threshold (Table S3). Duplication and horizontal gene transfer were the most prevalent markers for potential resistance genes or targets, with a mean relative frequency of 0.344 (±0.0535) and 0.426 (±0.0457), respectively. Comparative analysis showed that the addition of origin (termite-associated versus free-living) did not affect the residual deviance of a generalized linear model, indicating no difference between origins (glm, P = 0.061) (residual deviance = 322.37 [df = 459] versus residual deviance = 318.85 [df = 458]). Consistent with this, no apparent differences were observed between the frequency of markers by origin (Fig. S1; Table S3). This suggests that the antimicrobial potential is overall very similar between origins, consistent with the antiSMASH finding of comparable number and activities of BGCs across the genomes.

FIG S1

Frequency of essential gene markers produced by ARTS to indicate said gene being a resistance factor or resistant target. Each point is an identified essential gene, with three or more markers, observed in at least one termite-associated bacterial genome. The y axis is the frequency of how often that marker is positive in that specific essential gene across all genomes in which the essential gene was identified. Both axes have been jittered to improve visualization. “FGT” denotes that strains are fungus-growing termite-associated, and “free” indicates that strains are closely related soil-dwelling free-living bacteria. Download FIG S1, TIF file, 1.8 MB.
Copyright © 2021 Murphy et al.

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

Metabolic potential of Actinobacteria.HotPep identified 5,211 CAZymes across the 16 termite-associated genomes through homology to peptide patterns, of which 1,539 could be identified to EC number level, with 130 unique enzymes, comprised of 93 catabolic and 37 anabolic enzymes (Table S4). Ninety-seven of these were present in more than one genome (Fig. 4). Hydrolases accounted for 92% of all catabolic enzyme types, with 74 unique and 1,108 total enzymes (Fig. 4; Table S3). The predicted substrate targets of these were dominated by chitin, cellulose, and then hemicellulose, vital plant and fungal cell wall components (20) (Fig. 4; Table S3). Seven unique enzymes targeting chitin were observed 240 times, of which chitinase (EC 3.2.1.14) accounted for 65%. Notably, strains were often isolated on chitin medium, so this could bias this analysis. Enzymes cleaving cellulose and hemicellulose were observed on five and 10 unique occasions, with a total of 189 and 157 hits, respectively (Table S4). Cellulase and beta-glucosidase accounted for 40.7% and 41.7% of the cellulose-degrading enzymes, respectively, while acetylxylan esterase and endo-1,4-beta-xylanase accounted for 21% and 29% of hemicellulose-degrading enzymes, respectively. Termite-associated strains were not significantly different from free-living strains in their predicted degradative potential (glm, P = 0.253) (residual deviance = 1,020.3 [df = 1,409] versus 1,020.8 [df = 1,410]).

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

Diverse catabolic CAZymes are distributed across the 16 termite-associated Actinobacteria. Circular representation of the 43 unique catabolic plant and fungal cell wall-targeting CAZymes (total 823) sorted by target substrate. Each individual segment represents a unique CAZy gene or strain with a link between them indicating the identification of the CAZy gene in the strain. For the full results, see Table S4.

DISCUSSION

Actinobacteria are consistently present as low-abundance members of M. natalensis fungus combs, and to a lesser extent in gut microbiomes (Fig. 1C) (13). However, despite their low abundance, precedents from other work support that low-abundant taxa can be of significant importance (39–41). Furthermore, our current perception of their abundances in the termite symbiosis is likely to be underestimated, DNA extraction techniques tend to be biased against Gram-positive bacteria (42). Our analyses identified BGCs and CAZymes that may be of importance for the symbiosis, potentially with Actinobacteria as mutualistic symbionts, and our findings led us to three main conclusions. First, the Actinobacteria genomes contain numerous BGCs coding for natural products with potential antimicrobial activity (5, 43), of which a large proportion have yet to be elucidated chemically. This suggests ample potential for discovery of novel antimicrobials. Second, the genomes contain genes for enzymes targeting all plant and fungal substrates identified in fungus-growing termite gut community-level analyses (20). Lastly, although Actinobacteria represent a very small set of the complex microbiomes present within termite guts and fungus combs, their large enzymatic potential is consistent with their well-established roles in natural ecosystems as efficient biomass degraders (16, 20, 44–47) and prolific producers of antimicrobial compounds (5, 43). In line with this, the closely related free-living Actinobacteria show similar metabolic and enzymatic potential, thus making horizontal acquisition highly likely. In light of these findings, we discuss the potential for Actinobacteria serving roles in associations with farming termite hosts.

The persistent similarities in genomes irrespective of their origin being from termite colonies or free living suggest that the termite-obtained Actinobacteria do not display unique genomic signatures that point to specialized biomass decomposition or antimicrobial production. This could suggest that they may randomly colonize fungus combs or termites from the substrate harvested by the termites, from ingested soil (14, 28, 29, 48), or from contact with the surrounding soil. Although the net positive (mutualistic) benefit to the termites is currently not known, their presence implies potential effects in the comb environment, and considering that they are not removed by the termites, it makes it unlikely that they negatively impact the termite-Termitomyces association. The termites could then have coopted members of the phylum as symbionts because of their degradative or biosynthetic potential, without specific selection leading to apparent fungal genomic changes. This would mirror many other ectosymbiotic associations that have arisen because of symbiont (or host) metabolic capacities, but without genome erosion that is mostly apparent in strictly vertically transmitted endosymbionts (49–52). Furthermore, it appears to be only a subset of Actinobacteria in the surrounding soil that are capable of colonizing the symbiosis. Actinobacteria is one of the most abundant soil phyla (53–55), yet its members remain low in abundance within the symbiosis. To distinguish whether these are consistent colonizers of fungus-growing termites, we need to examine their presence in surrounding soils and within the symbiosis across nests.

If Actinobacteria symbiotically contribute to plant biomass decomposition by lignocellulolytic enzymes, this is most likely to occur within the comb environment (15, 22, 23) for two reasons. First, the predigestive gut passage of the substrate for fungiculture is comparably fast (Fig. 1A.1), and the relatively low abundance of the bacteria makes it unlikely that they play a major role at this stage of decomposition, despite their extensive cellulose- and hemicellulose-degrading capabilities. This is consistent with previous findings that cellulose, hemicellulose, and lignin-rich material enter the fungus comb in M. natalensis (17), within which Termitomyces dominates as the plant biomass degrader (16, 17, 56). Lignin may in part be cleaved during gut passage in the fungus-farming termite species Odontotermes formosanus (25), but although lignin-degrading enzymes such as laccase (EC 1.10.3.2) are extensively present in the Actinobacteria phylum (57), these occur infrequently in our isolates. Detailed gene expression analyses of combs to targeting both Termitomyces and bacterial communities remain to be performed but could provide insights into their relative contributions.

Maintaining monocultural fungus farming for years without other fungi present until colonies are compromised by external factors or succumb after the royal pair dies (58) requires highly effective defense functions. Actinobacteria have been proposed to play roles in facilitating these disease-free conditions through antimicrobial production (13, 22, 23, 59, 60). Our genome analyses of the termite-associated strains corroborate this promising defensive potential and detail a diverse set of NRPS-, PKS-, and terpene-related BGCs that are putative or known antimicrobials. The dominant compounds encoded by BGCs in both termite-associated and free-living genomes are common to Actinobacteria and not yet known to be antimicrobial (61–63). The potential broad spectrum of the diverse set of BGCs, supported by knowledge of compounds with known ranges and modes of action, may suggest an unselective antimicrobial landscape that could contribute to suppressing a series of targets. The most plausible site of Actinobacteria contributions is again within the fungus comb environment within which potential competitor or antagonist fungi might thrive (27).

One of our most striking findings was the high abundances of mycolytic enzymes targeting key fungal cell wall subunits encoded in the genomes. While termite gut microbes harbor large fungal biomass degradation potential (20, 25), the presence of these enzymes in predominantly comb-residing bacteria may appear counterintuitive. This could play a nutritional role; however, digestion of fungus material is not expected within fungus combs before ingestion by the termites. Alternatively, chitinase production could act as a layer of antifungal defense within combs if applied directly to sites of fungal infection by the termites, potentially targeting non-Termitomyces fungi. However, a perhaps more plausible explanation is that chitinases are ubiquitous in the comb and indiscriminately degrade fungal call wall within the comb, resulting in a potential cost of association for the termites and Termitomyces. This, however, would then potentially be mitigated by their low abundance and the persistent seeding and high abundance of Termitomyces or be outweighed by the putative antimicrobial benefits that Actinobacteria may provide.

The variety of antifungals identified within this study combined with the in vitro insights from previous bioassay reports (15) and tests of several specific isolates (M39, M56, RB108, and RB29 [13, 22, 23, 59, 60]) further supports unselective nest defense. Although a diverse range of fungal genera enter fungus combs after plant substrate passes the termite gut (27), their extremely low abundance within combs implies very effective growth suppression. Thus, further experimentation would be required to determine if seeding rate and general abundance of Termitomyces are what allows for monoculture with minimal presence of other fungi. Counter to this is identification of macrotermycins in Amycolatopsis M39, a putatively selective inhibitor of the stowaway fungus Pseudoxylaria (22), suggesting that a combination of targeted and untargeted approaches may be employed.

Our analyses of 16 termite-associated actinobacterial genomes also aimed to elucidate the genetic potential for the production of novel natural products and revealed a surprising number of BGCs with low/no homologies to previously reported BGCs. Our addition of ARTS2 resulted in the identification of an additional 50 essential genes that pass the three-marker threshold, with one of those markers being close proximity to an identified BGC. This remarkable finding speaks for the high potential to identify not only new bacterial diversity (64–66) but also novel natural products. Further work on both gene expression and compound production within termite colonies as well as biochemical studies to characterize the underlying gene clusters and their products is needed to further elaborate the roles of Actinobacteria within the fungus-farming termite system.

MATERIALS AND METHODS

Genomes, assembly, and annotation.Ten of the Actinobacteria genomes (RB29, RB68, RB99, M39, RB20, RB56, M41, M56, RB5, and RB17) were previously isolated from workers (cuticle, gut) or fungus combs (Fig. 1B) of the fungus-farming termite species Macrotermes natalensis, and genomes were downloaded from NCBI (see Table S1 in the supplemental material). We complement this set of genomes from the literature with five additional isolates that we sequence as part of this work; these have been deposited at NCBI (Table S1). The strains occupy the Actinomycetales order, spanning three suborders, six families, and six genera (Table S1).

Three of the six new strains (RB22, RB24, RB33, and RB108) were isolated from worker guts (RB22, RB24, and RB33) and two from fungus comb (RB108 and RB110-1 [referred to here as RB110]) of M. natalensis, and genomic DNA was extracted from cultures grown in ISP2 broth on a rotary shaker at 150 rpm at 30°C. Cells were harvested, and genomic DNA was extracted using the GenJet genomic DNA purification kit (Thermo Scientific, catalog no. K0721) following the manufacturer’s instructions, except with the following changes: lysozyme and proteinase K treatments were each extended to 40 min. For the Illumina library prep, genomic DNA was sheared using a Covaris S220 sonication device (Covaris Inc., MA, USA). Sequencing was performed using NEBNext Ultra II (New England Biolabs, Frankfurt, Germany) paired-end libraries on a MiSeq sequencer (RB22, RB24, RB33, and RB108) or NovaSeq sequencer (RB110). Quality trimming and adapter clipping were performed using Trimmomatic v.0.36 (67). Additional rounds of adapter clipping and filtering of low-complexity reads were performed using BBDuk of the BBTools package (68) v.36.84 (https://sourceforge.net/projects/bbmap/) and cutadapt (69) v.1.13 (https://github.com/marcelm/cutadapt). Overlapping read pairs were merged using FLASH (70) v.1.2.11 (https://ccb.jhu.edu/software/FLASH/). High-molecular-weight (HMW) DNA for PacBio-based whole-genome sequencing of RB110 was extracted using the NucleoBond HMW DNA kit (Macherey-Nagel).

Genomic DNA of Streptomyces sp. strain M32, isolated from an M. natalensis worker cuticle, was obtained from a 50-ml overnight culture at 30°C in ISP2. Extraction of DNA was performed using the GenElute bacterial genomic DNA kit (NA2100; Sigma-Aldrich), and the genome was sequenced at the Harvard Medical School Biopolymers Facility using the HiSeq2000 flow cell (Illumina CASAVA 1.8.2). Sequencing was performed using Illumina TruSeq 50-bp paired-end libraries on a HiSeq 2000 sequencer.

Because the new and downloaded genomes were compiled from different studies, assembly and annotation techniques and versions vary (see Table S1 for the full details). All Illumina MiSeq genomes (RB5 to -108) were assembled using SPAdes (71) v3.10.1 (https://github.com/ablab/spades), except RB29, for which v3.6.29 was used. M32 and M39 (22) were assembled with the A5 pipeline (72) v.20120518 (https://sourceforge.net/p/ngopt/wiki/A5PipelineREADME/). All assemblies were annotated with Prokka (73), with M56 and RB5 to -108 (except RB20, RB22, and RB56) utilizing v1.11 and the rest utilizing v1.12 beta. RB110 was hybrid assembled following a custom procedure. PacBio long reads were de novo assembled via Canu (74) v2.1.1 (https://github.com/marbl/canu) with default settings. This assembly was polished with the NextPolish pipeline v1.3.1 (75) using the quality-controlled (to a Phred score of 30) and merged Illumina NovaSeq and PacBio Sequel reads; see Text S1 in the supplemental material for details. The final assembly was annotated using Prokka (73) v1.14.5. RB22, RB24, RB33, and RB110 assemblies and annotated GenBank files were deposited to Zenodo (https://doi.org/10.5281/zenodo.4302144).

TEXT S1

Scripts used in biosynthetic gene cluster analysis and CAZyme and phylogenetic analysis, including statistical analysis. Consult the README file for specifics on how to run each analysis pipeline. All large-scale computation was performed on an HPC Conda environment manager with Slurm queue systems. Download Text S1, TXT file, 0.04 MB.
Copyright © 2021 Murphy et al.

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

To determine if termite-associated Actinobacteria differ genetically from free-living counterparts, we downloaded 14 high-quality genomes of the closely related free-living strains (based on BLASTn of the 16S rRNA sequence) from the NCBI RefSeq database (Table S1). All genomes of free-living isolates were annotated with Prokka (73) v1.14.5.

Phylogenetic placement of strains.Universal single-copy orthologous genes for all genomes were identified with BUSCO v4 (76) with default settings and autolineage flag engaged. Counts for all 130 genes in the BUSCO bacterial data set were determined, and multi-Fasta files were generated for ortholog genes present in three or more genomes and then aligned with Clustal Omega v1.2.4 (77). Phylogenies from the multilocus sequence typing (MLST) using 122 orthologs were generated using RAxML-NG v0.9.9 (78), employing the –all mode, GTR+G model, and a seed of 2. Branch support based on bootstrapping and transfer distance were obtained, before combining the gene trees into an unrooted species tree using ASTRAL-Pro v1.12 (79).

Biosynthetic potential and resistance mechanisms.Identification of the biosynthetic potential of Actinobacteria was carried out using antiSMASH v.5 (30, 31). The known clusters, subclusters, general, full HMMER, and pfam2go flags were utilized to generate a maximum output for a single run. Custom scripts were created to extract the known clusters information, that being biosynthetic gene clusters (BGCs) in which a similar compound was identified in the MIBiG database (https://mibig.secondarymetabolites.org/) along with the class of BGC (e.g., NRPSs, PKSs, terpenes, RiPPs [ribosomally synthesized and posttranslationally modified peptides], etc.). Only BGCs with ≥50% similarity score to the MIBiG database were retained to ensure likely similar functionality between homologous BGCs. A description of each compound was obtained via manual searching of the literature, including the mode of action and range of activity for compounds known to have antimicrobial properties. Comparative analysis of biosynthetic gene cluster class was undertaken by modeling the effects of isolate origin, either termite-associated or most closely related soil-dwelling free-living, on the occurrence of BGC class through Poisson distributed generalized linear models (glm) (from base R v4.0 [80]) with and without origin.

To explore the potential of unknown antimicrobial activity present in the genomes, we utilized Antibiotic Resistance Target Seeker (ARTS) version 2 (38). Prokka-annotated GenBank files with sequences for all genomes were run through the tool, which automates target-directed genome mining and builds on the concepts set out by Wright et al. (82) to exploit self-protection mechanisms inherent with the production of antimicrobially active compounds (37). This allowed searches for known resistance factors/targets and identification of potentially novel factor/targets through screening for duplications of essential genes, proximity to detected biosynthetic gene clusters (through antiSMASH), and whether a gene is likely to be horizontally transferred. These markers allow for inference of the essential gene being associated with self-protection and, thus, indicate antimicrobial activity of a gene-encoded natural product (37, 38). Essential genes are possible drug targets defined on a basis of ubiquity within a semibroad phylogenetic scope. Using Poisson-distributed generalized linear models, the results of ARTS2 were used to identify if origin (termite-associated versus soil-dwelling free-living) affected resistance factor/target marker counts associated with essential genes in a manner similar to the antiSMASH result modeling.

Metabolic potential.Identification of carbohydrate-active enzymes (CAZymes) was carried out via Homology to Peptide Pattern utilized by the HotPep tool (https://sourceforge.net/projects/hotpep/) using CAZyme PPRY patterns v.1. Following this, all identified CAZymes with an EC number to a specific enzyme were annotated with an enzyme name using the ExPASy enzyme.dat database (ftp.expasy.org/databases/enzyme) accessed on 17 February 2020. The overarching categories (lyases, transferases, hydrolases, etc.) of each enzyme were determined by taking the highest level of the corresponding EC number and searching the Brenda online database (https://brenda-enzymes.org/). The substrate of each enzyme was determined via manual searching of the Brenda and KEGG databases (http://www.genome.jp/kegg/). We comparatively analyzed summaries of counts of enzymes by target substrate by modeling the effect of origin (termite-associated versus free-living) using nested negative binomial distributed generalized linear models, and with substrate as a constant predictor variable in both models (glm.nb from MASS v7.3-51.6 [81]).

Data availability.All genomes are publicly available, either from previous publications, submission to the NCBI database under accession number JAEKDS000000000.1 for RB110, or Zenodo at https://doi.org/10.5281/zenodo.4302144 (RB22, RB24, RB33, M32, and RB108). All other accession numbers are provided in Table S1.

TABLE S5

Overview of all BGCs previously identified from Actinobacteria isolated from the fungus-growing termite symbiosis. Download Table S5, XLSX file, 0.02 MB.
Copyright © 2021 Murphy et al.

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

ACKNOWLEDGMENTS

We thank Suzanne G. Schmidt, Nick Bos, Veronica M. Sinotte, Kasun H. Bodawatta, Nils Peereboom, and Benjamin H. Conlon for comments on a previous draft of the manuscript.

This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germanýs Excellence Strategy – EXC 2051 – Project-ID 390713860 and BE 4799/3-1 to C.B. and a Villum Foundation Young Investigator grant (VKR10101) and a European Research Council Consolidator grant (771349) to M.P.

FOOTNOTES

    • Received December 16, 2020.
    • Accepted February 1, 2021.
  • Copyright © 2021 Murphy et al.

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

REFERENCES

  1. 1.↵
    1. Newman DJ,
    2. Cragg GM
    . 2016. Natural products as sources of new drugs from 1981 to 2014. J Nat Prod 79:629–661. doi:10.1021/acs.jnatprod.5b01055.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Mahmoud HM,
    2. Kalendar AA
    . 2016. Coral-associated Actinobacteria: diversity, abundance, and biotechnological potentials. Front Microbiol 7:204. doi:10.3389/fmicb.2016.00204.
    OpenUrlCrossRef
  3. 3.↵
    1. Egerton S,
    2. Culloty S,
    3. Whooley J,
    4. Stanton C,
    5. Ross RP
    . 2018. The gut microbiota of marine fish. Front Microbiol 9:873. doi:10.3389/fmicb.2018.00873.
    OpenUrlCrossRef
  4. 4.↵
    1. Jami M,
    2. Ghanbari M,
    3. Kneifel W,
    4. Domig KJ
    . 2015. Phylogenetic diversity and biological activity of culturable Actinobacteria isolated from freshwater fish gut microbiota. Microbiol Res 175:6–15. doi:10.1016/j.micres.2015.01.009.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Schrey SD,
    2. Erkenbrack E,
    3. Früh E,
    4. Fengler S,
    5. Hommel K,
    6. Horlacher N,
    7. Schulz D,
    8. Ecke M,
    9. Kulik A,
    10. Fiedler HP,
    11. Hampp R,
    12. Tarkka MT
    . 2012. Production of fungal and bacterial growth modulating secondary metabolites is widespread among mycorrhiza-associated streptomycetes. BMC Microbiol 12:164. doi:10.1186/1471-2180-12-164.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Hirsch AM,
    2. Valdés M
    . 2010. Micromonospora: an important microbe for biomedicine and potentially for biocontrol and biofuels. Soil Biol Biochem 42:536–542. doi:10.1016/j.soilbio.2009.11.023.
    OpenUrlCrossRef
  7. 7.↵
    1. Loqman S,
    2. Barka EA,
    3. Clément C,
    4. Ouhdouch Y
    . 2009. Antagonistic actinomycetes from Moroccan soil to control the grapevine gray mold. World J Microbiol Biotechnol 25:81–91. doi:10.1007/s11274-008-9864-6.
    OpenUrlCrossRef
  8. 8.↵
    1. Kroiss J,
    2. Kaltenpoth M,
    3. Schneider B,
    4. Schwinger MG,
    5. Hertweck C,
    6. Maddula RK,
    7. Strohm E,
    8. Svatos A
    . 2010. Symbiotic streptomycetes provide antibiotic combination prophylaxis for wasp offspring. Nat Chem Biol 6:261–263. doi:10.1038/nchembio.331.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    1. Kaltenpoth M,
    2. Göttler W,
    3. Herzner G,
    4. Strohm E
    . 2005. Symbiotic bacteria protect wasp larvae from fungal infestation. Curr Biol 15:475–479. doi:10.1016/j.cub.2004.12.084.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    1. Currie CR,
    2. Scott JA,
    3. Summerbell RC,
    4. Malloch D
    . 1999. Fungus-growing ants use antibiotic-producing bacteria to control garden parasites. Nature 398:701–704. doi:10.1038/19519.
    OpenUrlCrossRefWeb of Science
  11. 11.↵
    1. Currie CR,
    2. Bot ANM,
    3. Boomsma JJ
    . 2003. Experimental evidence of a tripartite mutualism: bacteria protect ant fungus gardens from specialized parasites. Oikos 101:91–102. doi:10.1034/j.1600-0706.2003.12036.x.
    OpenUrlCrossRefWeb of Science
  12. 12.↵
    1. Poulsen M,
    2. Cafaro MJ,
    3. Erhardt DP,
    4. Little AEF,
    5. Gerardo NM,
    6. Tebbets B,
    7. Klein BS,
    8. Currie CR
    . 2010. Variation in Pseudonocardia antibiotic defence helps govern parasite-induced morbidity in Acromyrmex leaf-cutting ants. Environ Microbiol Rep 2:534–540. doi:10.1111/j.1758-2229.2009.00098.x.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Benndorf R,
    2. Guo H,
    3. Sommerwerk E,
    4. Weigel C,
    5. Garcia-Altares M,
    6. Martin K,
    7. Hu H,
    8. Küfner M,
    9. de Beer ZW,
    10. Poulsen M,
    11. Beemelmanns C
    . 2018. Natural products from actinobacteria associated with fungus-growing termites. Antibiotics 7:83. doi:10.3390/antibiotics7030083.
    OpenUrlCrossRef
  14. 14.↵
    1. Otani S,
    2. Hansen LH,
    3. Sørensen SJ,
    4. Poulsen M
    . 2016. Bacterial communities in termite fungus combs are comprised of consistent gut deposits and contributions from the environment. Microb Ecol 71:207–220. doi:10.1007/s00248-015-0692-6.
    OpenUrlCrossRef
  15. 15.↵
    1. Visser AA,
    2. Nobre T,
    3. Currie CR,
    4. Aanen DK,
    5. Poulsen M
    . 2012. Exploring the potential for actinobacteria as defensive symbionts in fungus-growing termites. Microb Ecol 63:975–985. doi:10.1007/s00248-011-9987-4.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Poulsen M,
    2. Hu H,
    3. Li C,
    4. Chen Z,
    5. Xu L,
    6. Otani S,
    7. Nygaard S,
    8. Nobre T,
    9. Klaubauf S,
    10. Schindler PM,
    11. Hauser F,
    12. Pan H,
    13. Yang Z,
    14. Sonnenberg ASM,
    15. Wilhelm De Beer Z,
    16. Zhang Y,
    17. Wingfield MJ,
    18. De Grimmelikhuijzen CJP,
    19. Vries RP,
    20. Korb J,
    21. Aanen DK,
    22. Wang J,
    23. Boomsma JJ,
    24. Zhang G
    . 2014. Complementary symbiont contributions to plant decomposition in a fungus-farming termite. Proc Natl Acad Sci U S A 111:14500–14505. doi:10.1073/pnas.1319718111.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. da Costa RR,
    2. Hu H,
    3. Pilgaard B,
    4. Sabine SM,
    5. Schückel J,
    6. Pedersen KSK,
    7. Kračun SK,
    8. Busk PK,
    9. Harholt J,
    10. Sapountzis P,
    11. Lange L,
    12. Aanen DK,
    13. Poulsen M
    . 2017. Enzyme activities at different stages of plant biomass decomposition in three species of fungus-growing termites. Appl Environ Microbiol 84:e01815-17. doi:10.1128/AEM.01815-17.
    OpenUrlCrossRef
  18. 18.↵
    1. Otani S,
    2. Challinor VL,
    3. Kreuzenbeck NB,
    4. Kildgaard S,
    5. Krath Christensen S,
    6. Larsen LLM,
    7. Aanen DK,
    8. Rasmussen SA,
    9. Beemelmanns C,
    10. Poulsen M
    . 2019. Disease-free monoculture farming by fungus-growing termites. Sci Rep 9:8819. doi:10.1038/s41598-019-45364-z.
    OpenUrlCrossRef
  19. 19.↵
    1. Leuthold RH,
    2. Badertscher S,
    3. Imboden H
    . 1989. The inoculation of newly formed fungus comb with Termitomyces in Macrotermes colonies (Isoptera, Macrotermitinae). Insectes Soc 36:328–338. doi:10.1007/BF02224884.
    OpenUrlCrossRefWeb of Science
  20. 20.↵
    1. Hu H,
    2. da Costa RR,
    3. Pilgaard B,
    4. Schiøtt M,
    5. Lange L,
    6. Poulsen M
    . 2019. Fungiculture in termites is associated with a mycolytic gut bacterial community. mSphere 4:e00165-19. doi:10.1128/mSphere.00165-19.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Um S,
    2. Fraimout A,
    3. Sapountzis P,
    4. Oh DC,
    5. Poulsen M
    . 2013. The fungus-growing termite Macrotermes natalensis harbors bacillaene-producing Bacillus sp. that inhibit potentially antagonistic fungi. Sci Rep 3:3250. doi:10.1038/srep03250.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Beemelmanns C,
    2. Ramadhar TR,
    3. Kim KH,
    4. Klassen JL,
    5. Cao S,
    6. Wyche TP,
    7. Hou Y,
    8. Poulsen M,
    9. Bugni TS,
    10. Currie CR,
    11. Clardy J
    . 2017. Macrotermycins A-D, glycosylated macrolactams from a termite-associated Amycolatopsis sp. M39. Org Lett 19:1000–1003. doi:10.1021/acs.orglett.6b03831.
    OpenUrlCrossRef
  23. 23.↵
    1. Guo H,
    2. Benndorf R,
    3. Leichnitz D,
    4. Klassen JL,
    5. Vollmers J,
    6. Görls H,
    7. Steinacker M,
    8. Weigel C,
    9. Dahse HM,
    10. Kaster AK,
    11. de Beer ZW,
    12. Poulsen M,
    13. Beemelmanns C
    . 2017. Isolation, biosynthesis and chemical modifications of rubterolones A–F: rare tropolone alkaloids from Actinomadura sp. 5-2. Chemistry 23:9338–9345. doi:10.1002/chem.201701005.
    OpenUrlCrossRef
  24. 24.↵
    1. Sieber R,
    2. Leuthold RH
    . 1981. Behavioural elements and their meaning in incipient laboratory colonies of the fungus-growing termite Macrotermes michaelseni (Isoptera: Macrotermitinae). Insectes Soc 28:371–382. doi:10.1007/BF02224194.
    OpenUrlCrossRef
  25. 25.↵
    1. Li H,
    2. Yelle DJ,
    3. Li C,
    4. Yang M,
    5. Ke J,
    6. Zhang R,
    7. Liu Y,
    8. Zhu N,
    9. Liang S,
    10. Mo X,
    11. Ralph J,
    12. Currie CR,
    13. Mo J
    . 2017. Lignocellulose pretreatment in a fungus-cultivating termite. Proc Natl Acad Sci U S A 114:4709–4714. doi:10.1073/pnas.1618360114.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Nobre T,
    2. Aanen DK
    . 2012. Fungiculture or termite husbandry? The ruminant hypothesis. Insects 3:307–323. doi:10.3390/insects3010307.
    OpenUrlCrossRef
  27. 27.↵
    1. Bos N,
    2. Guimarães L,
    3. Palenzuela R,
    4. Renelies-Hamilton J,
    5. Maccario L,
    6. Silue SK,
    7. Koné NA,
    8. Poulsen M
    . 2020. You don’t have the guts: a diverse set of fungi survive passage through Macrotermes bellicosus termite guts. BMC Evol Biol 20:163. doi:10.1186/s12862-020-01727-z.
    OpenUrlCrossRef
  28. 28.↵
    1. Janssen PH
    . 2006. Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl Environ Microbiol 72:1719–1728. doi:10.1128/AEM.72.3.1719-1728.2006.
    OpenUrlFREE Full Text
  29. 29.↵
    1. Otani S,
    2. Mikaelyan A,
    3. Nobre T,
    4. Hansen LH,
    5. Koné NA,
    6. Sørensen SJ,
    7. Aanen DK,
    8. Boomsma JJ,
    9. Brune A,
    10. Poulsen M
    . 2014. Identifying the core microbial community in the gut of fungus-growing termites. Mol Ecol 23:4631–4644. doi:10.1111/mec.12874.
    OpenUrlCrossRef
  30. 30.↵
    1. Medema MH,
    2. Blin K,
    3. Cimermancic P,
    4. De Jager V,
    5. Zakrzewski P,
    6. Fischbach MA,
    7. Weber T,
    8. Takano E,
    9. Breitling R
    . 2011. antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res 39(Web Server issue):W339–W346. doi:10.1093/nar/gkr466.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Blin K,
    2. Shaw S,
    3. Steinke K,
    4. Villebro R,
    5. Ziemert N,
    6. Lee SY,
    7. Medema MH,
    8. Weber T
    . 2019. AntiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res 47:W81–W87. doi:10.1093/nar/gkz310.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Kautsar SA,
    2. Blin K,
    3. Shaw S,
    4. Navarro-Muñoz JC,
    5. Terlouw BR,
    6. Van Der Hooft JJJ,
    7. Van Santen JA,
    8. Tracanna V,
    9. Suarez Duran HG,
    10. Pascal Andreu V,
    11. Selem-Mojica N,
    12. Alanjary M,
    13. Robinson SL,
    14. Lund G,
    15. Epstein SC,
    16. Sisto AC,
    17. Charkoudian LK,
    18. Collemare J,
    19. Linington RG,
    20. Weber T,
    21. Medema MH
    . 2020. MIBiG 2.0: a repository for biosynthetic gene clusters of known function. Nucleic Acids Res 48(D1):D454–D458. doi:10.1093/nar/gkz882.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Dickschat JS
    . 2019. Bacterial diterpene biosynthesis. Angew Chem Int Ed Engl 58:15964–15976. doi:10.1002/anie.201905312.
    OpenUrlCrossRef
  34. 34.↵
    1. Yamada Y,
    2. Kuzuyama T,
    3. Komatsu M,
    4. Shin-Ya K,
    5. Omura S,
    6. Cane DE,
    7. Ikeda H
    . 2015. Terpene synthases are widely distributed in bacteria. Proc Natl Acad Sci U S A 112:857–862. doi:10.1073/pnas.1422108112.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Moody SC,
    2. Zhao B,
    3. Lei L,
    4. Nelson DR,
    5. Mullins JGL,
    6. Waterman MR,
    7. Kelly SL,
    8. Lamb DC
    . 2012. Investigating conservation of the albaflavenone biosynthetic pathway and CYP170 bifunctionality in streptomycetes. FEBS J 279:1640–1649. doi:10.1111/j.1742-4658.2011.08447.x.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Gürtler H,
    2. Pedersen R,
    3. Anthoni U,
    4. Christophersen C,
    5. Nielsen PH,
    6. Wellington EM,
    7. Pedersen C,
    8. Bock K
    . 1994. Albaflavenone, a sesquiterpene ketone with a zizaene skeleton produced by a streptomycete with a new rope morphology. J Antibiot (Tokyo) 47:434–439. doi:10.7164/antibiotics.47.434.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Alanjary M,
    2. Kronmiller B,
    3. Adamek M,
    4. Blin K,
    5. Weber T,
    6. Huson D,
    7. Philmus B,
    8. Ziemert N
    . 2017. The Antibiotic Resistant Target Seeker (ARTS), an exploration engine for antibiotic cluster prioritization and novel drug target discovery. Nucleic Acids Res 45:W42–W48. doi:10.1093/nar/gkx360.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Mungan MD,
    2. Alanjary M,
    3. Blin K,
    4. Weber T,
    5. Medema MH,
    6. Ziemert N
    . 2020. ARTS 2.0: feature updates and expansion of the Antibiotic Resistant Target Seeker for comparative genome mining. Nucleic Acids Res 48(W1):W546–W552. doi:10.1093/nar/gkaa374.
    OpenUrlCrossRef
  39. 39.↵
    1. Jousset A,
    2. Bienhold C,
    3. Chatzinotas A,
    4. Gallien L,
    5. Gobet A,
    6. Kurm V,
    7. Küsel K,
    8. Rillig MC,
    9. Rivett DW,
    10. Salles JF,
    11. Van Der Heijden MGA,
    12. Youssef NH,
    13. Zhang X,
    14. Wei Z,
    15. Hol GWH
    . 2017. Where less may be more: how the rare biosphere pulls ecosystems strings. ISME J 11:853–862. doi:10.1038/ismej.2016.174.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Benjamino J,
    2. Lincoln S,
    3. Srivastava R,
    4. Graf J
    . 2018. Low-abundant bacteria drive compositional changes in the gut microbiota after dietary alteration. Microbiome 6:86. doi:10.1186/s40168-018-0469-5.
    OpenUrlCrossRef
  41. 41.↵
    1. Dawson W,
    2. Hör J,
    3. Egert M,
    4. van Kleunen M,
    5. Pester M
    . 2017. A small number of low-abundance bacteria dominate plant species-specific responses during rhizosphere colonization. Front Microbiol 8:975. doi:10.3389/fmicb.2017.00975.
    OpenUrlCrossRef
  42. 42.↵
    1. Ketchum RN,
    2. Smith EG,
    3. Vaughan GO,
    4. Phippen BL,
    5. McParland D,
    6. Al-Mansoori N,
    7. Carrier TJ,
    8. Burt JA,
    9. Reitzel AM
    . 2018. DNA extraction method plays a significant role when defining bacterial community composition in the marine invertebrate Echinometra mathaei. Front Mar Sci 5:255. doi:10.3389/fmars.2018.00255.
    OpenUrlCrossRef
  43. 43.↵
    1. Sharma P,
    2. Thakur D
    . 2020. Antimicrobial biosynthetic potential and diversity of culturable soil actinobacteria from forest ecosystems of Northeast India. Sci Rep 10:4104. doi:10.1038/s41598-020-60968-6.
    OpenUrlCrossRef
  44. 44.↵
    1. Lacombe-Harvey MÈ,
    2. Brzezinski R,
    3. Beaulieu C
    . 2018. Chitinolytic functions in actinobacteria: ecology, enzymes, and evolution. Appl Microbiol Biotechnol 102:7219–7230. doi:10.1007/s00253-018-9149-4.
    OpenUrlCrossRef
  45. 45.↵
    1. Anderson I,
    2. Abt B,
    3. Lykidis A,
    4. Klenk H-P,
    5. Kyrpides N,
    6. Ivanova N
    . 2012. Genomics of aerobic cellulose utilization systems in actinobacteria. PLoS One 7:e39331. doi:10.1371/journal.pone.0039331.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Trojanowski J,
    2. Haider K,
    3. Sundman V
    . 1977. Decomposition of 14C-labelled lignin and phenols by a Nocardia sp. Arch Microbiol 114:149–153. doi:10.1007/BF00410776.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Petrosyan P,
    2. Luz-Madrigal A,
    3. Huitrón C,
    4. Flores ME
    . 2002. Characterization of a xylanolytic complex from Streptomyces sp. Biotechnol Lett 24:1473–1476. doi:10.1023/A:1019879200497.
    OpenUrlCrossRef
  48. 48.↵
    1. Makonde HM,
    2. Mwirichia R,
    3. Osiemo Z,
    4. Boga HI,
    5. Klenk HP
    . 2015. 454 Pyrosequencing-based assessment of bacterial diversity and community structure in termite guts, mounds and surrounding soils. Springerplus 4:471. doi:10.1186/s40064-015-1262-6.
    OpenUrlCrossRef
  49. 49.↵
    1. Kjeldsen KU,
    2. Bataillon T,
    3. Pinel N,
    4. De Mita S,
    5. Lund MB,
    6. Panitz F,
    7. Bendixen C,
    8. Stahl DA,
    9. Schramm A
    . 2012. Purifying selection and molecular adaptation in the genome of Verminephrobacter, the heritable symbiotic bacteria of earthworms. Genome Biol Evol 4:307–315. doi:10.1093/gbe/evs014.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Toft C,
    2. Andersson SGE
    . 2010. Evolutionary microbial genomics: insights into bacterial host adaptation. Nat Rev Genet 11:465–475. doi:10.1038/nrg2798.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Rio RVM,
    2. Lefevre C,
    3. Heddi A,
    4. Aksoy S
    . 2003. Comparative genomics of insect-symbiotic bacteria: influence of host environment on microbial genome composition. Appl Environ Microbiol 69:6825–6832. doi:10.1128/aem.69.11.6825-6832.2003.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Brown BP,
    2. Wernegreen JJ
    . 2019. Genomic erosion and extensive horizontal gene transfer in gut-associated Acetobacteraceae. BMC Genomics 20:472. doi:10.1186/s12864-019-5844-5.
    OpenUrlCrossRef
  53. 53.↵
    1. Cong J,
    2. Yang Y,
    3. Liu X,
    4. Lu H,
    5. Liu X,
    6. Zhou J,
    7. Li D,
    8. Yin H,
    9. Ding J,
    10. Zhang Y
    . 2015. Analyses of soil microbial community compositions and functional genes reveal potential consequences of natural forest succession. Sci Rep 5:10007. doi:10.1038/srep10007.
    OpenUrlCrossRef
  54. 54.↵
    1. Schrempf H
    . 2013. Actinobacteria within soils: capacities for mutualism, symbiosis and pathogenesis. FEMS Microbiol Lett 342:77–78. doi:10.1111/1574-6968.12147.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Hill P,
    2. Krištůfek V,
    3. Dijkhuizen L,
    4. Boddy C,
    5. Kroetsch D,
    6. van Elsas JD
    . 2011. Land use intensity controls actinobacterial community structure. Microb Ecol 61:286–302. doi:10.1007/s00248-010-9752-0.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. da Costa RR,
    2. Hu H,
    3. Li H,
    4. Poulsen M
    . 2019. Symbiotic plant biomass decomposition in fungus-growing termites. Insects 10:87. doi:10.3390/insects10040087.
    OpenUrlCrossRef
  57. 57.↵
    1. Janusz G,
    2. Pawlik A,
    3. Sulej J,
    4. Swiderska-Burek U,
    5. Jarosz-Wilkolazka A,
    6. Paszczynski A
    . 2017. Lignin degradation: microorganisms, enzymes involved, genomes analysis and evolution. FEMS Microbiol Rev 41:941–962. doi:10.1093/femsre/fux049.
    OpenUrlCrossRef
  58. 58.↵
    1. Batra LR,
    2. Batra SWT
    . 1966. Fungus-growing termites of tropical India and associated fungi. J Kans Entomol Soc 39:725–738.
    OpenUrl
  59. 59.↵
    1. Klassen JL,
    2. Lee SR,
    3. Poulsen M,
    4. Beemelmanns C,
    5. Kim KH
    . 2019. Efomycins K and L from a termite-associated Streptomyces sp. M56 and their putative biosynthetic origin. Front Microbiol 10:1739. doi:10.3389/fmicb.2019.01739.
    OpenUrlCrossRef
  60. 60.↵
    1. Kim KH,
    2. Ramadhar TR,
    3. Beemelmanns C,
    4. Cao S,
    5. Poulsen M,
    6. Currie CR,
    7. Clardy J
    . 2014. Natalamycin A, an ansamycin from a termite-associated Streptomyces sp. Chem Sci 5:4333–4338. doi:10.1039/c4sc01136h.
    OpenUrlCrossRefPubMed
  61. 61.↵
    1. Seipke RF,
    2. Loria R
    . 2009. Hopanoids are not essential for growth of Streptomyces scabies 87-22. J Bacteriol 191:5216–5223. doi:10.1128/JB.00390-09.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Czech L,
    2. Hermann L,
    3. Stöveken N,
    4. Richter AA,
    5. Höppner A,
    6. Smits SHJ,
    7. Heider J,
    8. Bremer E
    . 2018. Role of the extremolytes ectoine and hydroxyectoine as stress protectants and nutrients: genetics, phylogenomics, biochemistry, and structural analysis. Genes (Basel) 9:177. doi:10.3390/genes9040177.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Klausen C,
    2. Nicolaisen MH,
    3. Strobel BW,
    4. Warnecke F,
    5. Nielsen JL,
    6. Jørgensen NOG
    . 2005. Abundance of actinobacteria and production of geosmin and 2-methylisoborneol in Danish streams and fish ponds. FEMS Microbiol Ecol 52:265–278. doi:10.1016/j.femsec.2004.11.015.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    1. Benndorf R,
    2. Martin K,
    3. Küfner M,
    4. De Beer ZW,
    5. Vollmers J,
    6. Kaster AK,
    7. Beemelmanns C
    . 2020. Actinomadura rubteroloni sp. nov. and Actinomadura macrotermitis sp. nov., isolated from the gut of the fungus growing-termite Macrotermes natalensis. Int J Syst Evol Microbiol 70:5255–5262. doi:10.1099/ijsem.0.004403.
    OpenUrlCrossRef
  65. 65.↵
    1. Benndorf R,
    2. Schwitalla JW,
    3. Martin K,
    4. De Beer ZW,
    5. Vollmers J,
    6. Kaster AK,
    7. Poulsen M,
    8. Beemelmanns C
    . 2020. Nocardia macrotermitis sp. nov. and Nocardia aurantia sp. nov., isolated from the gut of the fungus-growing termite Macrotermes natalensis. Int J Syst Evol Microbiol 70:5226–5234. doi:10.1099/ijsem.0.004398.
    OpenUrlCrossRef
  66. 66.↵
    1. Schwitalla JW,
    2. Benndorf R,
    3. Martin K,
    4. Vollmers J,
    5. Kaster A-K,
    6. de Beer ZW,
    7. Poulsen M,
    8. Beemelmanns C
    . 2020. Streptomyces smaragdinus sp. nov., isolated from the gut of the fungus growing-termite Macrotermes natalensis. Int J Syst Evol Microbiol 70:5806–5811. doi:10.1099/ijsem.0.004478.
    OpenUrlCrossRef
  67. 67.↵
    1. Bolger AM,
    2. Lohse M,
    3. Usadel B
    . 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. doi:10.1093/bioinformatics/btu170.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Bushnell B
    . 2021. BBMap download. https://sourceforge.net/projects/bbmap/.
  69. 69.↵
    1. Martin M
    . 2011. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 17:10–12. doi:10.14806/ej.17.1.200.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Magoč T,
    2. Salzberg SL
    . 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27:2957–2963. doi:10.1093/bioinformatics/btr507.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Nurk S,
    2. Bankevich A,
    3. Antipov D,
    4. Gurevich A,
    5. Korobeynikov A,
    6. Lapidus A,
    7. Prjibelsky A,
    8. Pyshkin A,
    9. Sirotkin A,
    10. Sirotkin Y,
    11. Stepanauskas R,
    12. McLean J,
    13. Lasken R,
    14. Clingenpeel SR,
    15. Woyke T,
    16. Tesler G,
    17. Alekseyev MA,
    18. Pevzner PA
    . 2013. Assembling genomes and mini-metagenomes from highly chimeric reads. In Deng M, Jiang R, Sun F, Zhang X (ed), Research in computational molecular biology. RECOMB 2013. Lecture notes in computer science, vol 7821. Springer, Berlin, Germany.
  72. 72.↵
    1. Tritt A,
    2. Eisen JA,
    3. Facciotti MT,
    4. Darling AE
    . 2012. An integrated pipeline for de novo assembly of microbial genomes. PLoS One 7:e42304. doi:10.1371/journal.pone.0042304.
    OpenUrlCrossRefPubMed
  73. 73.↵
    1. Seemann T
    . 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069. doi:10.1093/bioinformatics/btu153.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Koren S,
    2. Walenz BP,
    3. Berlin K,
    4. Miller JR,
    5. Bergman NH,
    6. Phillippy AM
    . 2017. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res 27:722–736. doi:10.1101/gr.215087.116.
    OpenUrlAbstract/FREE Full Text
  75. 75.↵
    1. Hu J,
    2. Fan J,
    3. Sun Z,
    4. Liu S
    . 2020. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics 36:2253–2255. doi:10.1093/bioinformatics/btz891.
    OpenUrlCrossRef
  76. 76.↵
    1. Seppey M,
    2. Manni M,
    3. Zdobnov EM
    . 2019. BUSCO: assessing genome assembly and annotation completeness. Methods Mol Biol 1962:227–245. doi:10.1007/978-1-4939-9173-0_14.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Sievers F,
    2. Wilm A,
    3. Dineen D,
    4. Gibson TJ,
    5. Karplus K,
    6. Li W,
    7. Lopez R,
    8. McWilliam H,
    9. Remmert M,
    10. Soding J,
    11. Thompson JD,
    12. Higgins DG
    . 2011. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7:539–539. doi:10.1038/msb.2011.75.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    1. Kozlov AM,
    2. Darriba D,
    3. Flouri T,
    4. Morel B,
    5. Stamatakis A
    . 2019. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35:4453–4455. doi:10.1093/bioinformatics/btz305.
    OpenUrlCrossRefPubMed
  79. 79.↵
    1. Zhang C,
    2. Scornavacca C,
    3. Molloy E,
    4. Mirarab S
    . 2019. ASTRAL-Pro: quartet-based species tree inference despite paralogy. bioRxiv doi:10.1101/2019.12.12.874727.
    OpenUrlAbstract/FREE Full Text
  80. 80.↵
    R Development Core Team. 2018. A language and environment for statistical computing. 4.0. R Foundation for Statistical Computing, Vienna, Austria.
  81. 81.↵
    1. Venables WN,
    2. Ripley BD
    . 2002. Modern applied statistics with S (Google eBook). http://www.stats.ox.ac.uk/pub/MASS4/.
  82. 82.↵
    1. Wright GD
    . 2007. The antibiotic resistome: the nexus of chemical and genetic diversity. Nat Rev Microbiol 5:175–186. doi:10.1038/nrmicro1614.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Comparative Genomics Reveals Prophylactic and Catabolic Capabilities of Actinobacteria within the Fungus-Farming Termite Symbiosis
Robert Murphy, René Benndorf, Z. Wilhelm de Beer, John Vollmers, Anne-Kristin Kaster, Christine Beemelmanns, Michael Poulsen
mSphere Mar 2021, 6 (2) e01233-20; DOI: 10.1128/mSphere.01233-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

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

Enter multiple addresses on separate lines or separate them with commas.
Comparative Genomics Reveals Prophylactic and Catabolic Capabilities of Actinobacteria within the Fungus-Farming Termite Symbiosis
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Comparative Genomics Reveals Prophylactic and Catabolic Capabilities of Actinobacteria within the Fungus-Farming Termite Symbiosis
Robert Murphy, René Benndorf, Z. Wilhelm de Beer, John Vollmers, Anne-Kristin Kaster, Christine Beemelmanns, Michael Poulsen
mSphere Mar 2021, 6 (2) e01233-20; DOI: 10.1128/mSphere.01233-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

biosynthetic gene clusters
carbohydrate-active enzymes
Macrotermitinae
Streptomyces
Actinobacteria
Actinomadura
Amycolatopsis
Luteimicrobium
Mycolicibacterium
Nocardia
antimicrobial

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

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

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

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

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

Online ISSN: 2379-5042