Mycena genomes resolve the evolution of fungal bioluminescence.

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

Whole-Genome Sequence Resource of Calonectria ilicicola, the Casual Pathogen of Soybean Red Crown Rot.

Calonectria ilicicola (anamorph: Cylindrocladium parasiticum) is a soilborne plant-pathogenic fungus with a broad host range, and it can cause red crown rot of soybean and Cylindrocladium black rot of peanut, which has become an emerging threat to crop production worldwide. Limited molecular studies have focused on Calonectria ilicicola and one of the possible difficulties is the lack of genomic resources. This study presents the first high quality and near-completed genome of C. ilicicola, using the Oxford Nanopore GridION sequencing platform. A total of 16 contigs were assembled and the genome of C. ilicicola isolate F018 was estimated to have 11 chromosomes. Currently, the C. ilicicola F018 genome represents the most contiguous assembly, which has the lowest contig number and the highest contig N50 among all Calonectria genome resources. Putative protein-coding sequences and secretory proteins were estimated to be 17,308 and 1,930 in the C. ilicicola F018 genome, respectively; and the prediction was close to other plant-pathogenic fungi, such as Fusarium species, within the Nectriaceae family. The availability of this high-quality genome resource is expected to facilitate research on fungal biology and genetics of C. ilicicola and to support advanced understanding of pathogen virulence and disease management.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Blast Fungal Genomes Show Frequent Chromosomal Changes, Gene Gains and Losses, and Effector Gene Turnover.

Pyricularia is a fungal genus comprising several pathogenic species causing the blast disease in monocots. Pyricularia oryzae, the best-known species, infects rice, wheat, finger millet, and other crops. As past comparative and population genomics studies mainly focused on isolates of P. oryzae, the genomes of the other Pyricularia species have not been well explored. In this study, we obtained a chromosomal-level genome assembly of the finger millet isolate P. oryzae MZ5-1-6 and also highly contiguous assemblies of Pyricularia sp. LS, P. grisea, and P. pennisetigena. The differences in the genomic content of repetitive DNA sequences could largely explain the variation in genome size among these new genomes. Moreover, we found extensive gene gains and losses and structural changes among Pyricularia genomes, including a large interchromosomal translocation. We searched for homologs of known blast effectors across fungal taxa and found that most avirulence effectors are specific to Pyricularia, whereas many other effectors share homologs with distant fungal taxa. In particular, we discovered a novel effector family with metalloprotease activity, distinct from the well-known AVR-Pita family. We predicted 751 gene families containing putative effectors in 7 Pyricularia genomes and found that 60 of them showed differential expression in the P. oryzae MZ5-1-6 transcriptomes obtained under experimental conditions mimicking the pathogen infection process. In summary, this study increased our understanding of the structural, functional, and evolutionary genomics of the blast pathogen and identified new potential effector genes, providing useful data for developing crops with durable resistance.

Co-opting sulphur-carrier proteins from primary metabolic pathways for 2-thiosugar biosynthesis.

Sulphur is an essential element for life and is ubiquitous in living systems. Yet how the sulphur atom is incorporated into many sulphur-containing secondary metabolites is poorly understood. For bond formation between carbon and sulphur in primary metabolites, the major ionic sulphur sources are the persulphide and thiocarboxylate groups on sulphur-carrier (donor) proteins. Each group is post-translationally generated through the action of a specific activating enzyme. In all reported bacterial cases, the gene encoding the enzyme that catalyses the carbon-sulphur bond formation reaction and that encoding the cognate sulphur-carrier protein exist in the same gene cluster. To study the production of the 2-thiosugar moiety in BE-7585A, an antibiotic from Amycolatopsis orientalis, we identified a putative 2-thioglucose synthase, BexX, whose protein sequence and mode of action seem similar to those of ThiG, the enzyme that catalyses thiazole formation in thiamine biosynthesis. However, no gene encoding a sulphur-carrier protein could be located in the BE-7585A cluster. Subsequent genome sequencing uncovered a few genes encoding sulphur-carrier proteins that are probably involved in the biosynthesis of primary metabolites but only one activating enzyme gene in the A. orientalis genome. Further experiments showed that this activating enzyme can adenylate each of these sulphur-carrier proteins and probably also catalyses the subsequent thiolation, through its rhodanese domain. A proper combination of these sulphur-delivery systems is effective for BexX-catalysed 2-thioglucose production. The ability of BexX to selectively distinguish sulphur-carrier proteins is given a structural basis using X-ray crystallography. This study is, to our knowledge, the first complete characterization of thiosugar formation in nature and also demonstrates the receptor promiscuity of the A. orientalis sulphur-delivery system. Our results also show that co-opting the sulphur-delivery machinery of primary metabolism for the biosynthesis of sulphur-containing natural products is probably a general strategy found in nature.

Functional characterization of the agtABCD and agtSR operons for 4-aminobutyrate and 5-aminovalerate uptake and regulation in Pseudomonas aeruginosa PAO1.

Growth of Pseudomonas aeruginosa on diamines cadaverine, putrescine, and diaminopropane requires the γ-glutamylation pathway to convert these diamines into δ-aminovalerate (AMV), γ-aminobutyrate (GABA), and ß-alanine. From DNA microarrays experiments the agtABCD operon (PA0603-0606) encoding components for an ABC transporter system was found inducible by exogenous AMV, GABA, and ß-alanine, but not by diamines. Induction of the agtABCD operon was abolished in the mutants of upstream agtS (PA0600) or agtR (PA0601) genes encoding the membrane-anchored sensor and the response regulator of a two-component regulatory system, respectively. Growth phenotype analysis supports the physiological functions of these agt genes on utilization of AMV and GABA. Through measurements of ß-galactosidase activities from an agtA::lacZ fusion, the requirement of a functional AgtS in control of the induction effect by exogenous AMV and GABA was further substantiated. The recombinant hexa-hisidine tagged agtR was constructed and purified to demonstrate its specific interactions with the agtA promoter region by electrophoretic mobility shift assays. In summary, this study establishes the functions of agtSR and agtABCD operons in AMV and GABA uptake, and provides a potential linkage between AMV/GABA metabolism and polymicrobial infection through the recently reported function of agtR in sensing of peptidoglycan shed by gram-positive bacteria (Korgaonkar et al., Proc Natl Acad Sci USA 110:1059-1064, 2013).

Decoding the genome of bloodsucking midge Forcipomyia taiwana (Diptera: Ceratopogonidae): Insights into odorant receptor expansion.

Biting midges, notably those within the Ceratopogonidae family, have long been recognized for their epidemiological significance, both as nuisances and vectors for disease transmission in vertebrates. Despite their impact, genomic insights into these insects, particularly beyond the Culicoides genus, remain limited. In this study, we assembled the Forcipomyia taiwana (Shiraki) genome, comprising 113 scaffolds covering 130.4 Mbps-with the longest scaffold reaching 7.6 Mbps and an N50 value of 2.6 Mbps-marking a pivotal advancement in understanding the genetic architecture of ceratopogonid biting midges. Phylogenomic analyses reveal a shared ancestry between F. taiwana and Culicoides sonorensis Wirth & Jones, dating back approximately 124 million years, and highlight a dynamic history of gene family expansions and contractions within the Ceratopogonidae family. Notably, a substantial expansion of the odorant receptor (OR) gene family was observed, which is crucial for the chemosensory capabilities that govern biting midges' interactions with their environment, including host seeking and oviposition behaviors. The distribution of OR genes across the F. taiwana genome displays notable clusters on scaffolds, indicating localized tandem gene duplication events. Additionally, several collinear regions were identified, hinting at segmental duplications, inversions, and translocations, contributing to the olfactory system's evolutionary complexity. Among the 156 ORs identified in F. taiwana, 134 are biting midge-specific ORs, distributed across three distinct clades, each exhibiting unique motif features that distinguish them from the others. Through weighted gene co-expression network analysis, we correlated distinct gene modules with sex and reproductive status, laying the groundwork for future investigations into the interplay between gene expression and adaptive behaviors in F. taiwana. In conclusion, our study not only highlights the unique olfactory repertoire of ceratopogonid biting midges but also sets the stage for future studies into the genetic underpinnings of their unique biological traits and ecological strategies.

Molecular characterization of PauR and its role in control of putrescine and cadaverine catabolism through the γ-glutamylation pathway in Pseudomonas aeruginosa PAO1.

Pseudomonas aeruginosa PAO1 grows on a variety of polyamines as the sole source of carbon and nitrogen. Catabolism of polyamines is mediated by the γ-glutamylation pathway, which is complicated by the existence of multiple homologous enzymes with redundant specificities toward different polyamines for a more diverse metabolic capacity in this organism. Through a series of markerless gene knockout mutants and complementation tests, specific combinations of pauABCD (polyamine utilization) genes were deciphered for catabolism of different polyamines. Among six pauA genes, expression of pauA1, pauA2, pauA4, and pauA5 was found to be inducible by diamines putrescine (PUT) and cadaverine (CAD) but not by diaminopropane. Activation of these promoters was regulated by the PauR repressor, as evidenced by constitutively active promoters in the pauR mutant. The activities of these promoters were further enhanced by exogenous PUT or CAD in the mutant devoid of all six pauA genes. The recombinant PauR protein with a hexahistidine tag at its N terminus was purified, and specific bindings of PauR to the promoter regions of most pau operons were demonstrated by electromobility shift assays. Potential interactions of PUT and CAD with PauR were also suggested by chemical cross-linkage analysis with glutaraldehyde. In comparison, growth on PUT was more proficient than that on CAD, and this observed growth phenotype was reflected in a strong catabolite repression of pauA promoter activation by CAD but was completely absent as reflected by activation by PUT. In summary, this study clearly establishes the function of PauR in control of pau promoters in response to PUT and CAD for their catabolism through the γ-glutamylation pathway.

Severe acute respiratory syndrome coronavirus 3C-like protease-induced apoptosis.

The pathogenesis of severe acute respiratory syndrome-associated coronavirus (SARS-CoV) is an important issue for the treatment and prevention of severe acute respiratory syndrome. Recently, SARS-CoV has been demonstrated to induce cell apoptosis in Vero-E6 cells. The possible role of SARS-CoV 3C-like protease (3CLpro) in virus-induced apoptosis is characterized in this study. Growth arrest and apoptosis via caspase-3 and caspase-9 activities were demonstrated in SARS-CoV 3CLpro -expressing human promonocyte cells. The fluorescence intensity of dihydrorhodamine 123 staining indicated that cellular reactive oxygen species were markedly increased in SARS-CoV 3CLpro -expressing cells. Moreover, in vivo signalling pathway assay indicated that 3CLpro increased the activation of the nuclear factor-kappa B-dependent reporter, but inhibited activator protein-1-dependent transcription. This finding is likely to be responsible for virus-induced apoptotic signalling.

Binding interaction of SARS coronavirus 3CL(pro) protease with vacuolar-H+ ATPase G1 subunit.

The pathogenesis of severe acute respiratory syndrome coronavirus (SARS-CoV) is an important issue for treatment and prevention of SARS. Recently, SARS-CoV 3CL(pro) protease has been implied to be possible relevance to SARS-CoV pathogenesis. In this study, we intended to identify potential 3CL(pro)-interacting cellular protein(s) using the phage-displayed human lung cDNA library. The vacuolar-H+ ATPase (V-ATPase) G1 subunit that contained a 3CL(pro) cleavage site-like motif was identified as a 3CL(pro)-interacting protein, as confirmed using the co-immunoprecipitation assay and the relative affinity assay. In addition, our result also demonstrated the cleavage of the V-ATPase G1 fusion protein and the immunoprecipitation of cellular V-ATPase G1 by the 3CL(pro). Moreover, loading cells with SNARF-1 pH-sensitive dye showed that the intracellular pH in 3CL(pro)-expressing cells was significantly lower as compared to mock cells.

Functional characterization of cellulases identified from the cow rumen fungus Neocallimastix patriciarum W5 by transcriptomic and secretomic analyses.

BACKGROUND: Neocallimastix patriciarum is one of the common anaerobic fungi in the digestive tracts of ruminants that can actively digest cellulosic materials, and its cellulases have great potential for hydrolyzing cellulosic feedstocks. Due to the difficulty in culture and lack of a genome database, it is not easy to gain a global understanding of the glycosyl hydrolases (GHs) produced by this anaerobic fungus. RESULTS: We have developed an efficient platform that uses a combination of transcriptomic and proteomic approaches to N. patriciarum to accelerate gene identification, enzyme classification and application in rice straw degradation. By conducting complementary studies of transcriptome (Roche 454 GS and Illumina GA IIx) and secretome (ESI-Trap LC-MS/MS), we identified 219 putative GH contigs and classified them into 25 GH families. The secretome analysis identified four major enzymes involved in rice straw degradation: ß-glucosidase, endo-1,4-ß-xylanase, xylanase B and Cel48A exoglucanase. From the sequences of assembled contigs, we cloned 19 putative cellulase genes, including the GH1, GH3, GH5, GH6, GH9, GH18, GH43 and GH48 gene families, which were highly expressed in N. patriciarum cultures grown on different feedstocks. CONCLUSIONS: These GH genes were expressed in Pichia pastoris and/or Saccharomyces cerevisiae for functional characterization. At least five novel cellulases displayed cellulytic activity for glucose production. One ß-glucosidase (W5-16143) and one exocellulase (W5-CAT26) showed strong activities and could potentially be developed into commercial enzymes.