Gene make-up: rapid and massive intron gains after horizontal transfer of a bacterial α-amylase gene to Basidiomycetes

Horizontally transferred glycoside hydrolase 26 may aid hemipteran insects in plant tissue digestion

Glycoside hydrolases are enzymes that break down complex carbohydrates into simple sugars by catalyzing the hydrolysis of glycosidic bonds. There have been multiple instances of adaptive horizontal gene transfer of genes belonging to various glycoside hydrolase families from microbes to insects, as glycoside hydrolases can metabolize constituents of the carbohydrate-rich plant cell wall. In this study, we characterize the horizontal transfer of a gene from the glycoside hydrolase family 26 (GH26) from bacteria to insects of the order Hemiptera. Our phylogenies trace the horizontal gene transfer to the common ancestor of the superfamilies Pentatomoidea and Lygaeoidea, which include stink bugs and seed bugs. After horizontal transfer, the gene was assimilated into the insect genome as indicated by the gain of an intron, and a eukaryotic signal peptide. Subsequently, the gene has undergone independent losses and expansions in copy number in multiple lineages, suggesting an adaptive role of GH26s in some insects. Finally, we measured tissue-level gene expression of multiple stink bugs and the large milkweed bug using publicly available RNA-seq datasets. We found that the GH26 genes are highly expressed in tissues associated with plant digestion, especially in the principal salivary glands of the stink bugs. Our results are consistent with the hypothesis that this horizontally transferred GH26 was co-opted by the insect to aid in plant tissue digestion and that this HGT event was likely adaptive.

Two newly established and mutually related subfamilies GH13_48 and GH13_49 of the α-amylase family GH13

Currently, the main α-amylase family GH13 has been divided into 47 subfamilies in CAZy, with new subfamilies regularly emerging. The present in silico study was performed to highlight the groups, represented by the maltogenic amylase from Thermotoga neapolitana and the α-amylase from Haloarcula japonica, which are worth of creating their own new GH13 subfamilies. This enlarges functional annotation and thus allows more precise prediction of the function of putative proteins. Interestingly, those two share certain sequence features, e.g. the highly conserved cysteine in the second conserved sequence region (CSR-II) directly preceding the catalytic nucleophile, or the well-preserved GQ character of the end of CSR-VII. On the other hand, the two groups bear also specific and highly conserved positions that distinguish them not only from each other but also from representatives of remaining GH13 subfamilies established so far. For the T. neapolitana maltogenic amylase group, it is the stretch of residues at the end of CSR-V highly conserved as L-[DN]. The H. japonica α-amylase group can be characterized by a highly conserved [WY]-[GA] sequence at the end of CSR-II. Other specific sequence features include an almost fully conserved aspartic acid located directly preceding the general acid/base in CSR-III or well-preserved glutamic acid in CSR-IV. The assumption that these two groups represent two mutually related, but simultaneously independent GH13 subfamilies has been supported by phylogenetic analysis as well as by comparison of tertiary structures. The main α-amylase family GH13 has thus been expanded by two novel subfamilies GH13_48 and GH13_49. KEY POINTS: • In silico analysis of two groups of family GH13 members with characterized representatives • Identification of certain common, but also some specific sequence features in seven CSRs • Creation of two novel subfamilies-GH13_48 and GH13_49 within the CAZy database.

Functional Prokaryotic-Like Deoxycytidine Triphosphate Deaminases and Thymidylate Synthase in Eukaryotic Social Amoebae: Vertical, Endosymbiotic, or Horizontal Gene Transfer?

The de novo synthesis of deoxythymidine triphosphate uses several pathways: gram-negative bacteria use deoxycytidine triphosphate deaminase to convert deoxycytidine triphosphate into deoxyuridine triphosphate, whereas eukaryotes and gram-positive bacteria instead use deoxycytidine monophosphate deaminase to transform deoxycytidine monophosphate to deoxyuridine monophosphate. It is then unusual that in addition to deoxycytidine monophosphate deaminases, the eukaryote Dictyostelium discoideum has 2 deoxycytidine triphosphate deaminases (Dcd1Dicty and Dcd2Dicty). Expression of either DcdDicty can fully rescue the slow growth of an Escherichia coli dcd knockout. Both DcdDicty mitigate the hydroxyurea sensitivity of a Schizosaccharomyces pombe deoxycytidine monophosphate deaminase knockout. Phylogenies show that Dcd1Dicty homologs may have entered the common ancestor of the eukaryotic groups of Amoebozoa, Obazoa, Metamonada, and Discoba through an ancient horizontal gene transfer from a prokaryote or an ancient endosymbiotic gene transfer from a mitochondrion, followed by horizontal gene transfer from Amoebozoa to several other unrelated groups of eukaryotes. In contrast, the Dcd2Dicty homologs were a separate horizontal gene transfer from a prokaryote or a virus into either Amoebozoa or Rhizaria, followed by a horizontal gene transfer between them. ThyXDicty, the D. discoideum thymidylate synthase, another enzyme of the deoxythymidine triphosphate biosynthesis pathway, was suggested previously to be acquired from the ancestral mitochondria or by horizontal gene transfer from alpha-proteobacteria. ThyXDicty can fully rescue the E. coli thymidylate synthase knockout, and we establish that it was obtained by the common ancestor of social amoebae not from mitochondria but from a bacterium. We propose horizontal gene transfer and endosymbiotic gene transfer contributed to the enzyme diversity of the deoxythymidine triphosphate synthesis pathway in most social amoebae, many Amoebozoa, and other eukaryotes.

Functional Characterization of Recombinant Raw Starch Degrading α-Amylase from Roseateles terrae HL11 and Its Application on Cassava Pulp Saccharification

Exploring new raw starch-hydrolyzing α-amylases and understanding their biochemical characteristics are important for the utilization of starch-rich materials in bio-industry. In this work, the biochemical characteristics of a novel raw starch-degrading α-amylase (HL11 Amy) from Roseateles terrae HL11 was firstly reported. Evolutionary analysis revealed that HL11Amy was classified into glycoside hydrolase family 13 subfamily 32 (GH13_32). It contains four protein domains consisting of domain A, domain B, domain C and carbohydrate-binding module 20 (CMB20). The enzyme optimally worked at 50 °C, pH 4.0 with a specific activity of 6270 U/mg protein and 1030 raw starch-degrading (RSD) U/mg protein against soluble starch. Remarkably, HL11Amy exhibited activity toward both raw and gelatinized forms of various substrates, with the highest catalytic efficiency (kcat/Km) on starch from rice, followed by potato and cassava, respectively. HL11Amy effectively hydrolyzed cassava pulp (CP) hydrolysis, with a reducing sugar yield of 736 and 183 mg/g starch from gelatinized and raw CP, equivalent to 72% and 18% conversion based on starch content in the substrate, respectively. These demonstrated that HL11Amy represents a promising raw starch-degrading enzyme with potential applications in starch modification and cassava pulp saccharification.

Biotechnological potential and initial characterization of two novel sesquiterpene synthases from Basidiomycota Coniophora puteana for heterologous production of δ-cadinol

Background

Terpene synthases are versatile catalysts in all domains of life, catalyzing the formation of an enormous variety of different terpenoid secondary metabolites. Due to their diverse bioactive properties, terpenoids are of great interest as innovative ingredients in pharmaceutical and cosmetic applications. Recent advances in genome sequencing have led to the discovery of numerous terpene synthases, in particular in Basidiomycota like the wood rotting fungus Coniophora puteana, which further enhances the scope for the manufacture of terpenes for industrial purposes.

Results

In this study we describe the identification of two novel (+)-δ-cadinol synthases from C. puteana, Copu5 and Copu9. The sesquiterpene (+)-δ-cadinol was previously shown to exhibit cytotoxic activity therefore having an application as possible, new, and sustainably sourced anti-tumor agent. In an Escherichia coli strain, optimized for sesquiterpene production, titers of 225 mg l⁻¹ and 395 mg l⁻¹, respectively, could be achieved. Remarkably, both enzymes share the same product profile thereby representing the first two terpene synthases from Basidiomycota with identical product profiles. We solved the crystal structure of Copu9 in its closed conformation, for the first time providing molecular details of sesquiterpene synthase from Basidiomycota. Based on the Copu9 structure, we conducted structure-based mutagenesis of amino acid residues lining the active site, thereby altering the product profile. Interestingly, the mutagenesis study also revealed that despite the conserved product profiles of Copu5 and Copu9 different conformational changes may accompany the catalytic cycle of the two enzymes. This observation suggests that the involvement of tertiary structure elements in the reaction mechanism(s) employed by terpene synthases may be more complex than commonly expected.

Conclusion

The presented product selectivity and titers of Copu5 and Copu9 may pave the way towards a sustainable, biotechnological production of the potentially new bioactive (+)-δ-cadinol. Furthermore, Copu5 and Copu9 may serve as model systems for further mechanistic studies of terpenoid catalysis.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01791-8.

In Silico Analysis of Fungal and Chloride-Dependent α-Amylases within the Family GH13 with Identification of Possible Secondary Surface-Binding Sites

This study brings a detailed bioinformatics analysis of fungal and chloride-dependent α-amylases from the family GH13. Overall, 268 α-amylase sequences were retrieved from subfamilies GH13_1 (39 sequences), GH13_5 (35 sequences), GH13_15 (28 sequences), GH13_24 (23 sequences), GH13_32 (140 sequences) and GH13_42 (3 sequences). Eight conserved sequence regions (CSRs) characteristic for the family GH13 were identified in all sequences and respective sequence logos were analysed in an effort to identify unique sequence features of each subfamily. The main emphasis was given on the subfamily GH13_32 since it contains both fungal α-amylases and their bacterial chloride-activated counterparts. In addition to in silico analysis focused on eventual ability to bind the chloride anion, the property typical mainly for animal α-amylases from subfamilies GH13_15 and GH13_24, attention has been paid also to the potential presence of the so-called secondary surface-binding sites (SBSs) identified in complexed crystal structures of some particular α-amylases from the studied subfamilies. As template enzymes with already experimentally determined SBSs, the α-amylases from Aspergillus niger (GH13_1), Bacillus halmapalus, Bacillus paralicheniformis and Halothermothrix orenii (all from GH13_5) and Homo sapiens (saliva; GH13_24) were used. Evolutionary relationships between GH13 fungal and chloride-dependent α-amylases were demonstrated by two evolutionary trees—one based on the alignment of the segment of sequences spanning almost the entire catalytic TIM-barrel domain and the other one based on the alignment of eight extracted CSRs. Although both trees demonstrated similar results in terms of a closer evolutionary relatedness of subfamilies GH13_1 with GH13_42 including in a wider sense also the subfamily GH13_5 as well as for subfamilies GH13_32, GH13_15 and GH13_24, some subtle differences in clustering of particular α-amylases may nevertheless be observed.

Amyrel, a novel glucose-forming α-amylase from Drosophila with 4-α-glucanotransferase activity by disproportionation and hydrolysis of maltooligosaccharides

The α-amylase paralogue Amyrel present in true flies (Diptera Muscomorpha) has been classified as a glycoside hydrolase in CAZy family GH13 on the basis of its primary structure. Here we report that, in fact, Amyrel is currently unique amongst Animals as it possesses both the hydrolytic α-amylase activity (EC 3.2.1.1) and a 4-α-glucanotransferase (EC 2.4.1.25) transglycosylation activity. Amyrel reacts specifically on α-(1-4) glycosidic bonds of starch and related polymers but produces a complex mixture of maltooligosaccharides, in sharp contrast with canonical animal α-amylases. With model maltooligosaccharides G2 (maltose) to G7, the Amyrel reaction starts by a disproportionation leading to Gn-1 and Gn + 1 products, which become themselves substrates for new disproportionation cycles. As a result, all detectable odd- and even-numbered maltooligosaccharides at least up to G12 were observed. However, hydrolysis of these products proceeds simultaneously, as shown by p-nitrophenyl-tagged oligosaccharides and microcalorimetry, and upon prolonged reaction, glucose is the major end product followed by maltose. The main structural determinant of these atypical activities was found to be a Gly-His-Gly-Ala deletion in the so-called flexible loop bordering the active site. Indeed, engineering this deletion in pig pancreatic and D. melanogaster α-amylases results in reaction patterns similar to those of Amyrel. It is proposed that this deletion provides more freedom to the substrate for subsites occupancy and allows a less constrained action pattern resulting in versatile activities at the active site.

Comparative mitogenomics of Agaricomycetes: Diversity, abundance, impact and coding potential of putative open-reading frames

The mitochondrion is an organelle found in eukaryote organisms, and it is vital for different cellular pathways. The mitochondrion has its own DNA molecule and, because its genetic content is relatively conserved, despite the variation of size and structure, mitogenome sequences have been widely used as a promising molecular biomarker for taxonomy and evolution in fungi. In this study, the mitogenomes of two fungal species of Agaricomycetes class, Phellinotus piptadeniae and Trametes villosa, were assembled and annotated for the first time. We used these newly sequenced mitogenomes for comparative analyses with other 55 mitogenomes of Agaricomycetes available in public databases. Mitochondrial DNA (mtDNA) size and content are highly variable and non-coding and intronic regions, homing endonucleases (HEGs), and unidentified ORFs (uORFs) significantly contribute to the total size of the mitogenome. Furthermore, accessory genes (most of them as HEGs) are shared between distantly related species, most likely as a consequence of horizontal gene transfer events. Conversely, uORFs are only shared between taxonomically related species, most probably as a result of vertical evolutionary inheritance. Additionally, codon usage varies among mitogenomes and the GC content of mitochondrial features may be used to distinguish coding from non-coding sequences. Our results also indicated that transposition events of mitochondrial genes to the nuclear genome are not common. Despite the variation of size and content of the mitogenomes, mitochondrial genes seemed to be reliable molecular markers in our time-divergence analysis, even though the nucleotide substitution rates of mitochondrial and nuclear genomes of fungi are quite different. We also showed that many events of mitochondrial gene shuffling probably happened amongst the Agaricomycetes during evolution, which created differences in the gene order among species, even those of the same genus. Altogether, our study revealed new information regarding evolutionary dynamics in Agaricomycetes.

The emergence and evolution of intron-poor and intronless genes in intron-rich plant gene families

Eukaryotic genes can be classified into intronless (no introns), intron-poor (three or fewer introns per gene) or intron-rich. Early eukaryotic genes were mostly intron-rich, and their alternative splicing into multiple transcripts, giving rise to different proteins, might have played pivotal roles in adaptation and evolution. Interestingly, extant plant genomes contain many gene families with one or sometimes few sub-families with genes that are intron-poor or intronless, and it remains unknown when and how these intron-poor or intronless genes have originated and evolved, and what their possible functions are. In this study, we identified 33 such gene families that contained intronless and intron-poor sub-families. Intronless genes seemed to have first emerged in early land plant evolution, while intron-poor sub-families seemed first to have appeared in green algae. In contrast to intron-rich genes, intronless genes in intron-poor sub-families occurred later, and were subject to stronger functional constraints. Based on RNA-seq analyses in Arabidopsis and rice, intronless or intron-poor genes in AP2, EF-hand_7, bZIP, FAD_binding_4, STE_STE11, CAMK_CAMKL-CHK1 and C2 gene families were more likely to play a role in response to drought and salt stress, compared with intron-rich genes in the same gene families, whereas intronless genes in the B_lectin and S_locus_glycop gene family were more likely to participate in epigenetic processes and plant development. Understanding the origin and evolutionary trajectory, as well as the potential functions, of intronless and intron-poor sub-families provides further insight into plant genome evolution and the functional divergence of genes.

Bacterial-Like Nonribosomal Peptide Synthetases Produce Cyclopeptides in the Zygomycetous Fungus Mortierella alpina

Fungi are traditionally considered a reservoir of biologically active natural products. However, an active secondary metabolism has long not been attributed to early-diverging fungi such as Mortierella Here, we report on the biosynthesis of two series of cyclic pentapeptides, the malpicyclins and malpibaldins, as products of Mortierella alpina ATCC 32222. The molecular structures of malpicyclins were elucidated by high-resolution tandem mass spectrometry (HR-MS/MS), Marfey's method, and one-dimensional (1D) and 2D nuclear magnetic resonance (NMR) spectroscopy. In addition, malpibaldin biosynthesis was confirmed by HR-MS. Genome mining and comparative quantitative real-time PCR (qRT-PCR) expression analysis pointed at two pentamodular nonribosomal peptide synthetases (NRPSs), malpicyclin synthetase MpcA and malpibaldin synthetase MpbA, as candidate biosynthetic enzymes. Heterologous production of the respective adenylation domains and substrate specificity assays proved promiscuous substrate selection and confirmed their respective biosynthetic roles. In stark contrast to known fungal NRPSs, MpbA and MpcA contain bacterial-like dual epimerase/condensation domains allowing the racemization of enzyme-tethered l-amino acids and the subsequent incorporation of d-amino acids into the metabolites. Phylogenetic analyses of both NRPS genes indicated a bacterial origin and a horizontal gene transfer into the fungal genome. We report on the as-yet-unexplored nonribosomal peptide biosynthesis in basal fungi which highlights this paraphylum as a novel and underrated resource of natural products.IMPORTANCE Fungal natural compounds are industrially produced, with application in antibiotic treatment, cancer medications, and crop plant protection. Traditionally, higher fungi have been intensively investigated concerning their metabolic potential, but reidentification of already known compounds is frequently observed. Hence, alternative strategies to acquire novel bioactive molecules are required. We present the genus Mortierella as representative of the early-diverging fungi as an underestimated resource of natural products. Mortierella alpina produces two families of cyclopeptides, designated malpicyclins and malpibaldins, respectively, via two pentamodular nonribosomal peptide synthetases (NRPSs). These enzymes are much more closely related to bacterial than to other fungal NRPSs, suggesting a bacterial origin of these NRPS genes in Mortierella Both enzymes were biochemically characterized and are involved in as-yet-unknown biosynthetic pathways of natural products in basal fungi. Hence, this report establishes early-diverging fungi as prolific natural compound producers and sheds light on the origin of their biosynthetic capacity.

Extension of the taxonomic coverage of the family GH126 outside Firmicutes and in silico characterization of its non-catalytic terminal domains

The family GH126 is a family of glycoside hydrolases established in 2011. Officially, in the CAZy database, it counts ~ 1000 sequences originating solely from bacterial phylum Firmicutes. Two members, the proteins CPF_2247 from Clostridium perfringens and PssZ from Listeria monocytogenes have been characterized as a probable α-amylase and an exopolysaccharide-specific glycosidase, respectively; their three-dimensional structures being also solved as possessing catalytic (α/α)₆-barrel fold. Previously, based on a detailed in silico analysis, the seven conserved sequence regions (CSRs) were identified for the family along with elucidating basic evolutionary relationships within the family members. The present study represents a continuation study focusing on two particular aims: (1) to find out whether the taxonomic coverage of the family GH126 might be extended outside the Firmicutes and, if positive, to deliver those out-of-Firmicutes proteins with putting them into the context of the family; and (2) to identify the family members containing the N- and/or C-terminal extensions of their polypeptide chain, additional to the catalytic (α/α)₆-barrel domain, and perform the bioinformatics characterization of the extra domains. The main results could be summarized as follows: (1) 17 bacterial proteins caught by BLAST searches outside Firmicutes (especially from phyla Proteobacteria, Actinobacteria and Bacteroidetes) have been found and convincingly suggested as new family GH126 members; and (2) a thioredoxin-like fold and various leucine-rich repeat motifs identified by Phyre2 structure homology modelling have been recognized as extra domains occurring most frequently in the N-terminal extensions of family GH126 members possessing a modular organization.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus. Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Horizontal Transfer and Gene Loss Shaped the Evolution of Alpha-Amylases in Bilaterians

The subfamily GH13_1 of alpha-amylases is typical of Fungi, but it is also found in some unicellular eukaryotes (e.g., Amoebozoa, choanoflagellates) and non-bilaterian Metazoa. Since a previous study in 2007, GH13_1 amylases were considered ancestral to the Unikonts, including animals, except Bilateria, such that it was thought to have been lost in the ancestor of this clade. The only alpha-amylases known to be present in Bilateria so far belong to the GH13_15 and 24 subfamilies (commonly called bilaterian alpha-amylases) and were likely acquired by horizontal transfer from a proteobacterium. The taxonomic scope of Eukaryota genomes in databases has been greatly increased ever since 2007. We have surveyed GH13_1 sequences in recent data from ca. 1600 bilaterian species, 60 non-bilaterian animals and also in unicellular eukaryotes. As expected, we found a number of those sequences in non-bilaterians: Anthozoa (Cnidaria) and in sponges, confirming the previous observations, but none in jellyfishes and in Ctenophora. Our main and unexpected finding is that such fungal (also called Dictyo-type) amylases were also consistently retrieved in several bilaterian phyla: hemichordates (deuterostomes), brachiopods and related phyla, some molluscs and some annelids (protostomes). We discuss evolutionary hypotheses possibly explaining the scattered distribution of GH13_1 across bilaterians, namely, the retention of the ancestral gene in those phyla only and/or horizontal transfers from non-bilaterian donors.

A Robust Phylogenomic Time Tree for Biotechnologically and Medically Important Fungi in the Genera Aspergillus and Penicillium

Understanding the evolution of traits across technologically and medically significant fungi requires a robust phylogeny. Even though species in the Aspergillus and Penicillium genera (family Aspergillaceae, class Eurotiomycetes) are some of the most significant technologically and medically relevant fungi, we still lack a genome-scale phylogeny of the lineage or knowledge of the parts of the phylogeny that exhibit conflict among analyses. Here, we used a phylogenomic approach to infer evolutionary relationships among 81 genomes that span the diversity of Aspergillus and Penicillium species, to identify conflicts in the phylogeny, and to determine the likely underlying factors of the observed conflicts. Using a data matrix comprised of 1,668 genes, we found that while most branches of the phylogeny of the Aspergillaceae are robustly supported and recovered irrespective of method of analysis, a few exhibit various degrees of conflict among our analyses. Further examination of the observed conflict revealed that it largely stems from incomplete lineage sorting and hybridization or introgression. Our analyses provide a robust and comprehensive evolutionary genomic roadmap for this important lineage, which will facilitate the examination of the diverse technologically and medically relevant traits of these fungi in an evolutionary context.

The filamentous fungal family Aspergillaceae contains >1,000 known species, mostly in the genera Aspergillus and Penicillium. Several species are used in the food, biotechnology, and drug industries (e.g., Aspergillus oryzae and Penicillium camemberti), while others are dangerous human and plant pathogens (e.g., Aspergillus fumigatus and Penicillium digitatum). To infer a robust phylogeny and pinpoint poorly resolved branches and their likely underlying contributors, we used 81 genomes spanning the diversity of Aspergillus and Penicillium to construct a 1,668-gene data matrix. Phylogenies of the nucleotide and amino acid versions of this full data matrix as well as of several additional data matrices were generated using three different maximum likelihood schemes (i.e., gene-partitioned, unpartitioned, and coalescence) and using both site-homogenous and site-heterogeneous models (total of 64 species-level phylogenies). Examination of the topological agreement among these phylogenies and measures of internode certainty identified 11/78 (14.1%) bipartitions that were incongruent and pinpointed the likely underlying contributing factors, which included incomplete lineage sorting, hidden paralogy, hybridization or introgression, and reconstruction artifacts associated with poor taxon sampling. Relaxed molecular clock analyses suggest that Aspergillaceae likely originated in the lower Cretaceous and that the Aspergillus and Penicillium genera originated in the upper Cretaceous. Our results shed light on the ongoing debate on Aspergillus systematics and taxonomy and provide a robust evolutionary and temporal framework for comparative genomic analyses in Aspergillaceae. More broadly, our approach provides a general template for phylogenomic identification of resolved and contentious branches in densely genome-sequenced lineages across the tree of life.

Structure, evolution and diverse physiological roles of SWEET sugar transporters in plants

Present review describes the structure, evolution, transport mechanism and physiological functions of SWEETs. Their application using TALENs and CRISPR/CAS9 based genomic editing approach is discussed. Sugars Will Eventually be Exported Transporters (SWEET) proteins were first identified in plants as the novel family of sugar transporters which mediates the translocation of sugars across cell membranes. The SWEET family of sugar transporters is unique in terms of their structure which contains seven predicted transmembrane domains with two internal triple-helix bundles which possibly originate due to prokaryotic gene duplication. SWEETs perform diverse physiological functions such as pollen nutrition, nectar secretion, seed filling, phloem loading, and pathogen nutrition which we have discussed in the present review. We also discuss how transcriptional activator-like effector nucleases (TALENs) and CRISPR/CAS9 genome editing tools are used to engineer SWEET mutants which modulate pathogen resistance in plants and its applications in the field of agriculture. The expression of SWEETs promises to implement insights into many other cellular transport mechanisms. To conclude, the present review highlights the recent aspects which will further develop better understanding of molecular evolution, structure, and function of SWEET transporters in plants.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Eukaryotic Acquisition of a Bacterial Operon

Operons are a hallmark of bacterial genomes, where they allow concerted expression of functionally related genes as single polycistronic transcripts. They are rare in eukaryotes, where each gene usually drives expression of its own independent messenger RNAs. Here, we report the horizontal operon transfer of a siderophore biosynthesis pathway from relatives of Escherichia coli into a group of budding yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are expressed, exhibit eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer, several genetic changes occurred during subsequent evolution, including the gain of new transcription start sites that were sometimes within protein-coding sequences, acquisition of polyadenylation sites, structural rearrangements, and integration of eukaryotic genes into the cluster. We conclude that the genes were likely acquired as a unit, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly competitive, iron-limited environment.

Identification of sesquiterpene synthases from the Basidiomycota Coniophora puteana for the efficient and highly selective β-copaene and cubebol production in E. coli

Background

Terpenes are an important and extremely versatile class of secondary metabolites that are commercially used in the pharmaceutical, food and cosmetics sectors. Genome mining of different fungal collections has revealed the genetic basis for a steadily increasing number of putative terpene synthases without any detailed knowledge about their biochemical properties. The analysis and research of this rich genetic source provides a precious basis for the advancing biotechnological production of an almost endless number of valuable natural metabolites.

Results

Three annotated terpene synthases from the little investigated Basidiomycota Coniophora puteana were studied in this work. For biochemical characterization, the heterologous expression in E. coli was conducted leading to the identification of two sesquiterpene synthases capable of the highly selective generation of β-copaene and cubebol. These compounds are commercially used as food and flavor additives. The new enzymes show the highest reported product selectivity for their main compounds and therefore represent the first exclusive synthases for β-copaene (62% product selectivity) and cubebol (75% product selectivity) generation. In combination with an optimized heterologous microbial production system, we obtained product titers of 215 mg/L β-copaene and 497 mg/L cubebol.

Conclusion

The reported product selectivity and our generated terpene titers exceed all published biotechnological data regarding the production of β-copaene and cubebol. This represents a promising and economic alternative to extraction from natural plant sources and the associated complex product purification.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1010-z) contains supplementary material, which is available to authorized users.

Genome Mining of Non-Conventional Yeasts: Search and Analysis of MAL Clusters and Proteins

Genomic clustering of functionally related genes is rare in yeasts and other eukaryotes with only few examples available. Here, we summarize our data on a nontelomeric MAL cluster of a non-conventional methylotrophic yeast Ogataea (Hansenula) polymorpha containing genes for α-glucosidase MAL1, α-glucoside permease MAL2 and two hypothetical transcriptional activators. Using genome mining, we detected MAL clusters of varied number, position and composition in many other maltose-assimilating non-conventional yeasts from different phylogenetic groups. The highest number of MAL clusters was detected in Lipomyces starkeyi while no MAL clusters were found in Schizosaccharomyces pombe and Blastobotrys adeninivorans. Phylograms of α-glucosidases and α-glucoside transporters of yeasts agreed with phylogenesis of the respective yeast species. Substrate specificity of unstudied α-glucosidases was predicted from protein sequence analysis. Specific activities of Scheffersomycesstipitis α-glucosidases MAL7, MAL8, and MAL9 heterologously expressed in Escherichia coli confirmed the correctness of the prediction-these proteins were verified promiscuous maltase-isomaltases. α-Glucosidases of earlier diverged yeasts L. starkeyi, B. adeninivorans and S. pombe showed sequence relatedness with α-glucosidases of filamentous fungi and bacilli.

Sesquiterpene Synthase-3-Hydroxy-3-Methylglutaryl Coenzyme A Synthase Fusion Protein Responsible for Hirsutene Biosynthesis in Stereum hirsutum

The wood-rotting mushroom Stereum hirsutum is a known producer of a large number of namesake hirsutenoids, many with important bioactivities. Hirsutenoids form a structurally diverse and distinct class of sesquiterpenoids. No genes involved in hirsutenoid biosynthesis have yet been identified or their enzymes characterized. Here, we describe the cloning and functional characterization of a hirsutene synthase as an unexpected fusion protein of a sesquiterpene synthase (STS) with a C-terminal 3-hydroxy-3-methylglutaryl-coenzyme A (3-hydroxy-3-methylglutaryl-CoA) synthase (HMGS) domain. Both the full-length fusion protein and truncated STS domain are highly product-specific 1,11-cyclizing STS enzymes with kinetic properties typical of STSs. Complementation studies in Saccharomyces cerevisiae confirmed that the HMGS domain is also functional in vivo Phylogenetic analysis shows that the hirsutene synthase domain does not form a clade with other previously characterized sesquiterpene synthases from Basidiomycota. Comparative gene structure analysis of this hirsutene synthase with characterized fungal enzymes reveals a significantly higher intron density, suggesting that this enzyme may be acquired by horizontal gene transfer. In contrast, the HMGS domain is clearly related to other fungal homologs. This STS-HMGS fusion protein is part of a biosynthetic gene cluster that includes P450s and oxidases that are expressed and could be cloned from cDNA. Finally, this unusual fusion of a terpene synthase to an HMGS domain, which is not generally recognized as a key regulatory enzyme of the mevalonate isoprenoid precursor pathway, led to the identification of additional HMGS duplications in many fungal genomes, including the localization of HMGSs in other predicted sesquiterpenoid biosynthetic gene clusters.IMPORTANCE Hirsutenoids represent a structurally diverse class of bioactive sesquiterpenoids isolated from fungi. Identification of their biosynthetic pathways will provide access to this chemodiversity for the discovery and synthesis of molecules with new bioactivities. The identification and successful cloning of the previously elusive hirsutene synthase from the S. hirsutum provide important insights and strategies for biosynthetic gene discovery in Basidiomycota. The finding of a terpene synthase-HMGS fusion, the discovery of other sesquiterpenoid biosynthetic gene clusters with dedicated HMGS genes, and HMGS gene duplications in fungal genomes give new importance to the role of HMGS as a key regulatory enzyme in isoprenoid and sterol biosynthesis that should be exploited for metabolic engineering.

Isolation, identification and in silico analysis of alpha-amylase gene of Aspergillus niger strain CSA35 obtained from cassava undergoing spoilage

In this investigation, a gene (CDF_Amyl) encoding extracellular α-amylase in Aspergillus niger strain CSA35 associated with cassava spoilage was amplified using specific primers and characterized in silico. The gene had a partial nucleotide sequence of 968 bp and encoded a protein of 222 aa residues with a molecular weight and isoelectric point of 25.13 kDa and 4.17, respectively. Its catalytic site was located in the active site domain. BLASTp analysis showed that the protein primary sequence of the α-amylase gene had 98% and 99% homologies with the α-amylase of A. niger and A. oryzae RIB40, respectively. The gene is more closely related to α-amylase genes from fungi than to bacterial, plant, or animal α-amylase genes. Restriction mapping of the gene showed it can be digested with restriction enzymes like NcoI, PstI, SmaI, and BcLI among others but not with EcoRI and EcoRV. Its protein product had a hydrophobicity score of - 0.43 but no transmembrane helix. The CDF_Amyl protein was subcellularly localized in the secretory pathway, an indication of its release into extracellular space after secretion. Also, the 3D structure of the CDF-Amyl protein was barrel-shaped with domains characteristic of α-amylases. The encoded α-amylase Vmax is 6.90 U/mg protein and Km is 6.70 mg/ml. It was concluded that the unique characteristics of the CDF_Amyl gene and its deduced protein could find applications in biotechnological, food and pharmaceutical industries where cloning and further modification of this gene would be required for product development and improvement.

Purification of an alpha amylase from Aspergillus flavus NSH9 and molecular characterization of its nucleotide gene sequence

In this study, an alpha-amylase enzyme from a locally isolated Aspergillus flavus NSH9 was purified and characterized. The extracellular α-amylase was purified by ammonium sulfate precipitation and anion-exchange chromatography at a final yield of 2.55-fold and recovery of 11.73%. The molecular mass of the purified α-amylase was estimated to be 54 kDa using SDS-PAGE and the enzyme exhibited optimal catalytic activity at pH 5.0 and temperature of 50 °C. The enzyme was also thermally stable at 50 °C, with 87% residual activity after 60 min. As a metalloenzymes containing calcium, the purified α-amylase showed significantly increased enzyme activity in the presence of Ca²⁺ ions. Further gene isolation and characterization shows that the α-amylase gene of A. flavus NSH9 contained eight introns and an open reading frame that encodes for 499 amino acids with the first 21 amino acids presumed to be a signal peptide. Analysis of the deduced peptide sequence showed the presence of three conserved catalytic residues of α-amylase, two Ca²⁺-binding sites, seven conserved peptide sequences, and several other properties that indicates the protein belongs to glycosyl hydrolase family 13 capable of acting on α-1,4-bonds only. Based on sequence similarity, the deduced peptide sequence of A. flavus NSH9 α-amylase was also found to carry two potential surface/secondary-binding site (SBS) residues (Trp 237 and Tyr 409) that might be playing crucial roles in both the enzyme activity and also the binding of starch granules.

AS3MT-mediated tolerance to arsenic evolved by multiple independent horizontal gene transfers from bacteria to eukaryotes

Organisms have evolved the ability to tolerate toxic substances in their environments, often by producing metabolic enzymes that efficiently detoxify the toxicant. Inorganic arsenic is one of the most toxic and carcinogenic substances in the environment, but many organisms, including humans, metabolise inorganic arsenic to less toxic metabolites. This multistep process produces mono-, di-, and trimethylated arsenic metabolites, which the organism excretes. In humans, arsenite methyltransferase (AS3MT) appears to be the main metabolic enzyme that methylates arsenic. In this study, we examined the evolutionary origin of AS3MT and assessed the ability of different genotypes to produce methylated arsenic metabolites. Phylogenetic analysis suggests that multiple, independent horizontal gene transfers between different bacteria, and from bacteria to eukaryotes, increased tolerance to environmental arsenic during evolution. These findings are supported by the observation that genetic variation in AS3MT correlates with the capacity to methylate arsenic. Adaptation to arsenic thus serves as a model for how organisms evolve to survive under toxic conditions.

Evolutionary dynamics of triosephosphate isomerase gene intron location pattern in Metazoa: A new perspective on intron evolution in animals

Intron evolution, including its dynamics in the evolutionary transitions and diversification of eukaryotes, remains elusive. Inadequate taxon sampling due to data shortage, unclear phylogenetic framework, and inappropriate outgroup application might be among the causes. Besides, the integrity of all the introns within a gene was often neglected previously. Taking advantage of the ancient conserved triosephosphate isomerase gene (tim), the relatively robust phylogeny of Metazoa, and choanoflagellates as outgroup, the evolutionary dynamics of tim intron location pattern (ILP) in Metazoa was investigated. From 133 representative species of ten phyla, 30 types of ILPs were identified. A most common one, which harbors the maximum six intron positions, is deduced to be the common ancestral tim ILP of Metazoa, which almost had formed in their protozoan ancestor and was surprisingly retained and passed down till to each ancestors of metazoan phyla. In the subsequent animal diversification, it underwent different evolutionary trajectories: within Deuterostomia, it was almost completely retained only with changes in a few species with relatively recently fast-evolving histories, while within the rapidly radiating Protostomia, besides few but remarkable retention, it usually displayed extensive intron losses and a few gains. Therefore, a common ancestral exon-intron arrangement pattern of an animal gene is definitely discovered; besides the 'intron-rich view' of early animal genes being confirmed, the novel insight that high exon-intron re-arrangements of genes seem to be associated with the relatively recently rapid evolution of lineages/species/genomes but have no correlation with the ancient major evolutionary transitions in animal evolution, is revealed.

Transcriptomics analysis of hulless barley during grain development with a focus on starch biosynthesis

Hulless barley, with its unique nutritional value and potential health benefits, has increasingly attracted attentions in recent years. However, the transcription dynamics during hulless barley grain development is not well understood. In the present study, we investigated the transcriptome changes during barley grain development using Illumina paired-end RNA-sequencing. Two datasets of the developing grain transcriptomes from two barley landraces with the differential seed starch synthesis traits were generated, and comparative transcriptome approach in both genotypes was performed. The results showed that 38 differentially expressed genes (DEGs) were found co-modulated in both genotypes during the barley grain development. Of those, the proteins encoded by most of those DGEs were found, such as alpha-amylase-related proteins, lipid-transfer protein, homeodomain leucine zipper (HD-Zip), NUCLEAR FACTOR-Y, subunit B (NF-YBs), as well as MYB transcription factors. More interestingly, two genes Hvulgare_GLEAN_10012370 and Hvulgare_GLEAN_10021199 encoding SuSy, AGPase (Hvulgare_GLEAN_10033640 and Hvulgare_GLEAN_10056301), as well as SBE2b (Hvulgare_GLEAN_10018352) were found to significantly contribute to the regulatory mechanism during grain development in both genotypes. Moreover, six co-expression modules associated with specific biological processes or pathways (M1 to M6) were identified by consensus co-expression network. Significantly enriched pathways of those module genes showed difference in both genotypes. These results will expand our understanding of the complex molecular mechanism of starch synthesis during barley grain development.

A single amino-acid substitution toggles chloride dependence of the alpha-amylase paralog amyrel in Drosophila melanogaster and Drosophila virilis species

In animals, most α-amylases are chloride-dependent enzymes. A chloride ion is required for allosteric activation and is coordinated by one asparagine and two arginine side chains. Whereas the asparagine and one arginine are strictly conserved, the main chloride binding arginine is replaced by a glutamine in some rare instances, resulting in the loss of chloride binding and activation. Amyrel is a distant paralogue of α-amylase in Diptera, which was not characterized biochemically to date. Amyrel shows both substitutions depending on the species. In Drosophila melanogaster, an arginine is present in the sequence but in Drosophila virilis, a glutamine occurs at this position. We have investigated basic enzymological parameters and the dependence to chloride of Amyrel of both species, produced in yeast, and in mutants substituting arginine to glutamine or glutamine to arginine. We found that the amylolytic activity of Amyrel is about thirty times weaker than the classical Drosophila α-amylase, and that the substitution of the arginine by a glutamine in D. melanogaster suppressed the chloride-dependence but was detrimental to activity. In contrast, changing the glutamine into an arginine rendered D. virilis Amyrel chloride-dependent, and interestingly, significantly increased its catalytic efficiency. These results show that the chloride ion is not mandatory for Amyrel but stimulates the reaction rate. The possible phylogenetic origin of the arginine/glutamine substitution is also discussed.

Analysis on evolutionary relationship of amylases from archaea, bacteria and eukaryota

Amylase is one of the earliest characterized enzymes and has many applications in clinical and industrial settings. In biotechnological industries, the amylase activity is enhanced through modifying amylase structure and through cloning and expressing targeted amylases in different species. It is important to understand how engineered amylases can survive from generation to generation. This study used phylogenetic and statistical approaches to explore general patterns of amylases evolution, including 3118 α-amylases and 280 β-amylases from archaea, eukaryota and bacteria with fully documented taxonomic lineage. First, the phylogenetic tree was created to analyze the evolution of amylases with focus on individual amylases used in biofuel industry. Second, the average pairwise p-distance was computed for each kingdom, phylum, class, order, family and genus, and its diversity implies multi-time and multi-clan evolution. Finally, the variance was further partitioned into inter-clan variance and intra-clan variance for each taxonomic group, and they represent horizontal and vertical gene transfer. Theoretically, the results show a full picture on the evolution of amylases in manners of vertical and horizontal gene transfer, and multi-time and multi-clan evolution as well. Practically, this study provides the information on the surviving chance of desired amylase in a given taxonomic group, which may potentially enhance the successful rate of cloning and expression of amylase gene in different species.

Evaluation of the mechanisms of intron loss and gain in the social amoebae Dictyostelium

Background

Spliceosomal introns are a common feature of eukaryotic genomes. To approach a comprehensive understanding of intron evolution on Earth, studies should look beyond repeatedly studied groups such as animals, plants, and fungi. The slime mold Dictyostelium belongs to a supergroup of eukaryotes not covered in previous studies.

Results

We found 441 precise intron losses in Dictyostelium discoideum and 202 precise intron losses in Dictyostelium purpureum. Consistent with these observations, Dictyostelium discoideum was found to have significantly more copies of reverse transcriptase genes than Dictyostelium purpureum. We also found that the lost introns are significantly further from the 5′ end of genes than the conserved introns. Adjacent introns were prone to be lost simultaneously in Dictyostelium discoideum. In both Dictyostelium species, the exonic sequences flanking lost introns were found to have a significantly higher GC content than those flanking conserved introns. Together, these observations support a reverse-transcription model of intron loss in which intron losses were caused by gene conversion between genomic DNA and cDNA reverse transcribed from mature mRNA. We also identified two imprecise intron losses in Dictyostelium discoideum that may have resulted from genomic deletions. Ninety-eight putative intron gains were also observed. Consistent with previous studies of other lineages, the source sequences were found in only a small number of cases, with only two instances of intron gain identified in Dictyostelium discoideum.

Conclusions

Although they diverged very early from animals and fungi, Dictyostelium species have similar mechanisms of intron loss.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0567-y) contains supplementary material, which is available to authorized users.

Novel Introner-Like Elements in fungi Are Involved in Parallel Gains of Spliceosomal Introns

Spliceosomal introns are key components of the eukaryotic gene structure. Although they contributed to the emergence of eukaryotes, their origin remains elusive. In fungi, they might originate from the multiplication of invasive introns named Introner-Like Elements (ILEs). However, so far ILEs have been observed in six fungal species only, including Fulvia fulva and Dothistroma septosporum (Dothideomycetes), arguing against ILE insertion as a general mechanism for intron gain. Here, we identified novel ILEs in eight additional fungal species that are phylogenetically related to F. fulva and D. septosporum using PCR amplification with primers derived from previously identified ILEs. The ILE content appeared unique to each species, suggesting independent multiplication events. Interestingly, we identified four genes each containing two gained ILEs. By analysing intron positions in orthologues of these four genes in Ascomycota, we found that three ILEs had inserted within a 15 bp window that contains regular spliceosomal introns in other fungal species. These three positions are not the result of intron sliding because ILEs are newly gained introns. Furthermore, the alternative hypothesis of an inferred ancestral gain followed by independent losses contradicts the observed degeneration of ILEs. These observations clearly indicate three parallel intron gains in four genes that were randomly identified. Our findings suggest that parallel intron gain is a phenomenon that has been highly underestimated in ILE-containing fungi, and likely in the whole fungal kingdom.

Imprecise intron losses are less frequent than precise intron losses but are not rare in plants

Abstract

In this study, we identified 19 intron losses, including 11 precise intron losses (PILs), six imprecise intron losses (IILs), one de-exonization, and one exon deletion in tomato and potato, and 17 IILs in Arabidopsis thaliana. Comparative analysis of related genomes confirmed that all of the IILs have been fixed during evolution. Consistent with previous studies, our results indicate that PILs are a major type of intron loss. However, at least in plants, IILs are unlikely to be as rare as previously reported.

Reviewers

This article was reviewed by Jun Yu and Zhang Zhang. For complete reviews, see the Reviewers’ Reports section.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0056-7) contains supplementary material, which is available to authorized users.

Horizontal transfer of carbohydrate metabolism genes into ectomycorrhizal Amanita

The genus Amanita encompasses both symbiotic, ectomycorrhizal fungi and asymbiotic litter decomposers; all species are derived from asymbiotic ancestors. Symbiotic species are no longer able to degrade plant cell walls. The carbohydrate esterases family 1 (CE1s) is a diverse group of enzymes involved in carbon metabolism, including decomposition and carbon storage. CE1 genes of the ectomycorrhizal A. muscaria appear diverged from all other fungal homologues, and more similar to CE1s of bacteria, suggesting a horizontal gene transfer (HGT) event. In order to test whether AmanitaCE1s were acquired horizontally, we built a phylogeny of CE1s collected from across the tree of life, and describe the evolution of CE1 genes among Amanita and relevant lineages of bacteria. CE1s of symbiotic Amanita were very different from CE1s of asymbiotic Amanita, and are more similar to bacterial CE1s. The protein structure of one CE1 gene of A. muscaria matched a depolymerase that degrades the carbon storage molecule poly((R)-3-hydroxybutyrate) (PHB). Asymbiotic Amanita do not carry sequence or structural homologues of these genes. The CE1s acquired through HGT may enable novel metabolisms, or play roles in signaling or defense. This is the first evidence for the horizontal transfer of carbohydrate metabolism genes into ectomycorrhizal fungi.

Traversing the fungal terpenome

Fungi (Ascomycota and Basidiomycota) are prolific producers of structurally diverse terpenoid compounds. Classes of terpenoids identified in fungi include the sesqui-, di- and triterpenoids. Biosynthetic pathways and enzymes to terpenoids from each of these classes have been described. These typically involve the scaffold generating terpene synthases and cyclases, and scaffold tailoring enzymes such as e.g. cytochrome P450 monoxygenases, NAD(P)+ and flavin dependent oxidoreductases, and various group transferases that generate the final bioactive structures. The biosynthesis of several sesquiterpenoid mycotoxins and bioactive diterpenoids has been well-studied in Ascomycota (e.g. filamentous fungi). Little is known about the terpenoid biosynthetic pathways in Basidiomycota (e.g. mushroom forming fungi), although they produce a huge diversity of terpenoid natural products. Specifically, many trans-humulyl cation derived sesquiterpenoid natural products with potent bioactivities have been isolated. Biosynthetic gene clusters responsible for the production of trans-humulyl cation derived protoilludanes, and other sesquiterpenoids, can be rapidly identified by genome sequencing and bioinformatic methods. Genome mining combined with heterologous biosynthetic pathway refactoring has the potential to facilitate discovery and production of pharmaceutically relevant fungal terpenoids.

The complex intron landscape and massive intron invasion in a picoeukaryote provides insights into intron evolution

Genes in pieces and spliceosomal introns are a landmark of eukaryotes, with intron invasion usually assumed to have happened early on in evolution. Here, we analyze the intron landscape of Micromonas, a unicellular green alga in the Mamiellophyceae lineage, demonstrating the coexistence of several classes of introns and the occurrence of recent massive intron invasion. This study focuses on two strains, CCMP1545 and RCC299, and their related individuals from ocean samplings, showing that they not only harbor different classes of introns depending on their location in the genome, as for other Mamiellophyceae, but also uniquely carry several classes of repeat introns. These introns, dubbed introner elements (IEs), are found at novel positions in genes and have conserved sequences, contrary to canonical introns. This IE invasion has a huge impact on the genome, doubling the number of introns in the CCMP1545 strain. We hypothesize that each IE class originated from a single ancestral IE that has been colonizing the genome after strain divergence by inserting copies of itself into genes by intron transposition, likely involving reverse splicing. Along with similar cases recently observed in other organisms, our observations in Micromonas strains shed a new light on the evolution of introns, suggesting that intron gain is more widespread than previously thought.

Horizontal gene transfer of a bacterial insect toxin gene into the Epichloë fungal symbionts of grasses

Horizontal gene transfer is recognized as an important factor in genome evolution, particularly when the newly acquired gene confers a new capability to the recipient species. We identified a gene similar to the makes caterpillars floppy (mcf1 and mcf2) insect toxin genes in Photorhabdus, bacterial symbionts of nematodes, in the genomes of the Epichloë fungi, which are intercellular symbionts of grasses. Infection by Epichloë spp. often confers insect resistance to the grass hosts, largely due to the production of fungal alkaloids. A mcf-like gene is present in all of the Epichloë genome sequences currently available but in no other fungal genomes. This suggests the Epichloë genes were derived from a single lineage-specific HGT event. Molecular dating was used to estimate the time of the HGT event at between 7.2 and 58.8 million years ago. The mcf-like coding sequence from Epichloë typhina subsp. poae was cloned and expressed in Escherichia coli. E. coli cells expressing the Mcf protein were toxic to black cutworms (Agrotis ipsilon), whereas E. coli cells containing the vector only were non-toxic. These results suggest that the Epichloëmcf-like genes may be a component, in addition to the fungal alkaloids, of the insect resistance observed in Epichloë-infected grasses.

α-Amylase: an enzyme specificity found in various families of glycoside hydrolases

α-Amylase (EC 3.2.1.1) represents the best known amylolytic enzyme. It catalyzes the hydrolysis of α-1,4-glucosidic bonds in starch and related α-glucans. In general, the α-amylase is an enzyme with a broad substrate preference and product specificity. In the sequence-based classification system of all carbohydrate-active enzymes, it is one of the most frequently occurring glycoside hydrolases (GH). α-Amylase is the main representative of family GH13, but it is probably also present in the families GH57 and GH119, and possibly even in GH126. Family GH13, known generally as the main α-amylase family, forms clan GH-H together with families GH70 and GH77 that, however, contain no α-amylase. Within the family GH13, the α-amylase specificity is currently present in several subfamilies, such as GH13_1, 5, 6, 7, 15, 24, 27, 28, 36, 37, and, possibly in a few more that are not yet defined. The α-amylases classified in family GH13 employ a reaction mechanism giving retention of configuration, share 4-7 conserved sequence regions (CSRs) and catalytic machinery, and adopt the (β/α)8-barrel catalytic domain. Although the family GH57 α-amylases also employ the retaining reaction mechanism, they possess their own five CSRs and catalytic machinery, and adopt a (β/α)7-barrel fold. These family GH57 attributes are likely to be characteristic of α-amylases from the family GH119, too. With regard to family GH126, confirmation of the unambiguous presence of the α-amylase specificity may need more biochemical investigation because of an obvious, but unexpected, homology with inverting β-glucan-active hydrolases.

Biosynthesis of terpenoid natural products in fungi

Tens of thousands of terpenoid natural products have been isolated from plants and microbial sources. Higher fungi (Ascomycota and Basidiomycota) are known to produce an array of well-known terpenoid natural products, including mycotoxins, antibiotics, antitumor compounds, and phytohormones. Except for a few well-studied fungal biosynthetic pathways, the majority of genes and biosynthetic pathways responsible for the biosynthesis of a small number of these secondary metabolites have only been discovered and characterized in the past 5-10 years. This chapter provides a comprehensive overview of the current knowledge on fungal terpenoid biosynthesis from biochemical, genetic, and genomic viewpoints. Enzymes involved in synthesizing, transferring, and cyclizing the prenyl chains that form the hydrocarbon scaffolds of fungal terpenoid natural products are systematically discussed. Genomic information and functional evidence suggest differences between the terpenome of the two major fungal phyla--the Ascomycota and Basidiomycota--which will be illustrated for each group of terpenoid natural products.