In this study, a thermophilic cellulolytic filamentous fungus Thermothelomyces heterothallica was engineered for direct production of cellobionic acid (CBA) from alkali-pretreated wheat straw without the addition of external cellulases. The CBA bioproduction pathway was established through sequential disruption of genes involved in cellobiose and CBA catabolism, including eight β-glucosidases, one cellobiose/cellobionate phosphorylase, and two putative cellobionate transporters using a ribonucleoprotein (RNP)-based CRISPR-Cas9 system. The resulting T. heterothallica strain, TH11, was able to produce cellobionate from cellobiose and pretreated wheat straw at a high yield. Further optimization of fermentation pH using 50 mM citrate buffer (pH 6.0) increased the final cellobionate concentration. Fermentation temperature optimization showed that 43 °C accelerated substrate conversion without significantly affecting final titers. The engineered strain, TH11 produces 64 mM of cellobionate in 5 days. The cellulose conversion was about 88
We report a draft genome assembly of Polystigma amygdalinum, the causal agent of almond red leaf blotch. DNA extracted from infected leaves was sequenced using PacBio HiFi, and host-derived reads were removed bioinformatically. The 238.7-Mb assembly (90.8% BUSCO completeness) is highly repetitive (82.3%) and unusually large for an ascomycete.
Phytophthora is a long-established, well-known, and globally important genus of plant pathogens. Phylogenetic evidence has shown that the biologically distinct, obligate biotrophic downy mildews evolved from Phytophthora at least twice. Because, cladistically, this renders Phytophthora "paraphyletic," it has been proposed that Phytophthora evolutionary clades be split into multiple genera (Crous et al. 2021; Runge et al. 2011; Thines 2023, 2024). In this letter, we review arguments for the retention of the generic name Phytophthora with a broad circumscription made by Brasier et al. (2022) and by many delegates at an open workshop organized by The American Phytopathological Society. We present our well-considered responses to the genus splitting proposals, both in general terms and in terms of the specific proposals for new genera, alongside new information regarding the biological properties and mode of origin of the Phytophthora clades. We consider that the proposals are mostly non-rigorous and not supported by the scientific evidence. Further, given (i) the apparent lack of any distinguishing biological characteristics (synapomorphies) between the Phytophthora clades; (ii) the fundamental monophyly of Phytophthora in the original Haeckelian sense (Haeckel 1877); (iii) the fact that paraphyly is not a justification for taxonomic splitting; and (iv) the considerable likely damage to effective scientific communication and disease management from an unnecessary breakup of the genus, we report that workshop delegates voted unanimously in favor of preserving the current generic concept and for seeking endorsement of this view by a working group of the International Commission on the Taxonomy of Fungi. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.
Tanoak ( Notholithocarpus densiflorus ) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks ( Lithocarpus spp .), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks ( Quercus ) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to ‘sudden oak death’ (SOD), a plant pathogen ( Phytophthora ramorum ) that has caused widespread mortality of tanoaks. Here, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using co-barcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore (ONT) long reads to 30X coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published until now, with a contig N50 of ∼1.2 Mb and a scaffold N50 of ∼2.1 Mb. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using co-barcoded data we were able to generate complete haplotype coverage of all 11 genomes.
Gluconic acid’s potential as a wheat straw pretreatment agent was studied at different concentrations (0.125–1 M) and temperatures (160–190 °C) for 30 min, followed by enzymatic hydrolysis. 0.125 M gluconic acid, 170 °C, yielded the highest xylose output, while 0.5 M gluconic acid at 190 °C yielded the best glucose yield. A fraction of gluconic acid decomposed during pretreatment. Detoxified hemicellulose hydrolysate from 0.125 M gluconate at 170 °C for 60 min showed promise for ethanol production. The gluconate contained in the detoxified hemicellulose hydrolysate can be fermented to ethanol along with other hemicellulose sugars present by Escherichia coli SL100. The ethanol yield from gluconate and sugars was about 90.4 ± 1.8
This study investigated cellobionate production from a lignocellulosic substrate using Neurospora crassa HL10. Utilizing NaOH-pretreated wheat straw as the substrate obviated the need for an exogenous redox mediator addition, as lignin contained in the pretreated wheat served as a natural mediator. The low laccase production by N. crassa HL10 on pretreated wheat straw caused slow cellobionate production, and exogenous laccase addition accelerated the process. Cycloheximide induced substantial laccase production in N. crassa HL10, enabling the strain to yield approximately 57 mM cellobionate from pretreated wheat straw (equivalent to 20 g/L cellulose), shortening the conversion time from 8 to 6 days. About 92% of the cellulose contained in the pretreated wheat straw is converted to cellobionate. In contrast to existing methods requiring pure cellobiose or cellulase enzymes, this process efficiently converts a low-cost feedstock into cellobionate at a high yield without enzyme or redox mediator supplementation.
Phytophthora species are oomycete plant pathogens that cause great economic and ecological impacts. The Phytophthora genus includes over 180 known species, infecting a wide range of plant hosts, including crops, trees, and ornamentals. We sequenced the genomes of 31 individual Phytophthora species and 24 individual transcriptomes to study genetic relationships across the genus. De novo genome assemblies revealed variation in genome sizes, numbers of predicted genes, and in repetitive element content across the Phytophthora genus. A genus-wide comparison evaluated orthologous groups of genes. Predicted effector gene counts varied across Phytophthora species by effector family, genome size, and plant host range. Predicted numbers of apoplastic effectors increased as the host range of Phytophthora species increased. Predicted numbers of cytoplasmic effectors also increased with host range but leveled off or decreased in Phytophthora species that have enormous host ranges. With extensive sequencing across the Phytophthora genus, we now have the genomic resources to evaluate horizontal gene transfer events across the oomycetes. Using a machine-learning approach to identify horizontally transferred genes with bacterial or fungal origin, we identified 44 candidates over 36 Phytophthora species genomes. Phylogenetic reconstruction indicates that the transfers of most of these 44 candidates happened in parallel to major advances in the evolution of the oomycetes and Phytophthora spp. We conclude that the 31 genomes presented here are essential for investigating genus-wide genomic associations in genus Phytophthora. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.
The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.
We report a novel production process for lactobionic acid (LBA) production using an engineered Neurospora crassa strain F5. The wild-type N. crassa strain produces cellobiose dehydrogenase (CDH) and uses lactose as a carbon source. N. crassa strain F5, which was constructed by deleting six out of the seven β-glucosidases in the wild type, showed a much slower lactose utilization rate and produced a much higher level of cellobiose dehydrogenase (CDH) than the wild type. Strain N. crassa F5 produced CDH and laccase simultaneously on the pretreated wheat straw with 3 µM of cycloheximide added as the laccase inducer. The deproteinized cheese whey was added directly to the shake flasks with the fungus present to achieve LBA production. Strain F5 produced about 37 g/L of LBA from 45 g/L of lactose in 27 h since deproteinized cheese whey addition. The yield of LBA from consumed lactose was about 85
Background Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations. Results Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce’s disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes. Conclusions This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.
Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. In the genome of Neurospora crassa, most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, some of which regulate self-recognition and define vegetative incompatibility groups. In contrast, the LSG-encoding proteins possess few to no domains that would help to identify potential functional roles. Possible functional roles of LSGs were further assessed by performing transcriptomic profiling in genetic mutants and in response to environmental alterations, as well as examining gene knockouts for phenotypes. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural, a wildfire-produced chemical that is a strong inducer of sexual development, and the structurally-related furan 5-hydroxymethyl furfural (HMF). Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that also regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. These observations encouraged further investigation of the roles of clustered lineage-specific and HET-domain genes in ecology and reproduction regulation in Neurospora, especially the regulation of the switch from the asexual growth to sexual reproduction, in response to dramatic environmental conditions changes.
As de novo elements of the genome, orphan genes have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. However, there is a lack of working hypotheses as to what role they play. Systematic investigation of the evolutionary history, genome structures, expression dynamics and regulation of orphan genes in well-studied models can shed light on their functions in reproductive isolation. The 670 orphan genes that we identified in the genome of the fungal model Neurospora crassa were aggregated adjacent to the telomeres, and over 63% of them formed clusters with 61% of the het -like genes, genes that regulate self-recognition and define vegetative compatibility groups. Nearly all orphan-het clusters are syntenic within closely related species, suggesting relatively recent colocalization. Analysis of transcriptomic data from N. crassa conditions associated with growth and reproduction reveals 342 orphan genes that are dynamically expressed during both asexual and sexual phases. Among these, 37% were detectably expressed in culture on common carbon resources, but 64% were detectably expressed in culture on unusual carbon sources such as furfural and HMF-wildfire-produced chemicals that are a strong inducer of sexual development in N. crassa . Expression of a significant portion of the orphan genes was sensitive to light and temperature, factors that regulate asexual and sexual reproduction, among other fungal activities. Coordinate expression in orphan-het gene clusters was detected during early hyphal branching. Furthermore, orphan genes and clustered het-like genes respond similarly to mutations in transcription factors adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the orphan genes was affected by a mating locus mutant. Functional interactions between orphan and het -like genes likely consolidate vegetative incompatibility during asexual reproduction, but possibly promote crossings between vegetative compatibility groups during sexual reproduction. Their involvement in the balance between genome homogeneity and heterogeneity in vegetative compatibility likely contribute to Neurospora speciation. Orphan genes are potential targets for the manipulation of fungal growth, a key aspect of both the control of fungal pathogenicity and the iterative improvement of fungal biotechnology.### Competing Interest StatementThe authors have declared no competing interest.
Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. However, there is a lack of working hypotheses as to how they might play that role. In the genome of Neurospora crassa , most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, which regulate self-recognition and define vegetative incompatibility groups. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural and HMF—wildfire-produced chemicals that are a strong inducer of sexual development. Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. Accordingly, we propose a gene-by-environment interaction model encouraging further investigation of the roles of LSGs and HET-domain genes in speciation in Neurospora . This gene-by-environment interaction model emphasizes the roles of the LSGs in response to genetic and environmental factors, leading to the regulation of the switch from the asexual growth and fusion, such that vegetative incompatibility governed by allorecognition promotes allelic homogeneity, sexual reproduction, and outbreeding, whereas VI repression and meiotic recombination promotes allelic polymorphism.
An engineered Neurospora crassa strain (HL10) naturally produces cellulases, cellobiose dehydrogenase (CDH), and heterologously expresses laccase. It can convert Avicel (cellulose) to cellobionic acid (cellobionate) without any enzyme addition at high yields. However, the addition of a catalytic amount of artificial redox mediators was required to achieve a high yield conversion. When a lignocellulosic substrate like wheat straw was used as the carbon source, it was found that adding exogenous artificial redox mediators did not improve cellobionate yield. We hypothesized that lignin and lignin degradation products were able to serve as redox mediators for the CDH-laccase conversion system. Experiments were designed to prove the hypothesis. When N. crassa HL10 was grown on Avicel, the addition of enzyme hydrolyzed lignin (EHL) can achieve a similar level of cellobionate production as with ABTS addition. The addition of EHL can increase the cellobionate production rate in the CDH-laccase bi-enzyme system using cellobiose as the substrate. The formation of lignin radicals, and the quenching of lignin radicals by reduced CDH, were verified by an electron paramagnetic resonance (EPR) experiment, providing further evidence that lignin radicals can serve as the electron acceptor of reduced CDH. Hence, lignocellulosic biomass is a self-sufficient substrate for the production of cellobionate. The conversion of the lignocellulosic substrate to cellobionate using N. crassa HL10 can be achieved without any exogenous enzyme addition or redox mediator addition.
Phytophthora species cause crown, root, and aerial canker diseases in diverse horticultural crops, resulting in a loss of yield, loss of quality, and plant death. Breeding programs are interested in developing genotypes resistant to these and other pathogens, but current screening methods for tree crops are time-consuming, potentially spanning many years. Here we present a rapid in vitro assay in which walnut shoots obtained from tissue culture can be challenged with Phytophthora zoospores and monitored for symptom development within a week. As a proof of concept, two cultivars, scion variety ‘Chandler’ and clonal rootstock ‘RX1’, with known different degrees of crown rot susceptibility, were inoculated in vitro with P. pini (formerly included in P. citricola complex), and results of the assay reproducibly matched the differences in disease susceptibility generally observed for these two walnut selections under greenhouse and orchard conditions. In addition, when the in vitro assay was used for inoculations with P. capsici, which is not pathogenic on walnut, no disease development was observable, indicating that the assay can discriminate between species of Phytophthora that are aggressive on walnut and those not known as walnut pathogens. Our results suggest that this in vitro shoot inoculation procedure may provide rapid assessments useful for prescreening resistance to Phytophthora in walnuts. Further testing of the assay is justified to determine whether it can resolve more subtle differences in resistance and whether it can be useful with other perennial hosts of Phytophthora.
Phytophthora ramorum is an invasive, broad host range pathogen that causes ramorum blight and sudden oak death in forest landscapes of western North America. In commercial nurseries, asymptomatic infections of nursery stock by P. ramorum and other Phytophthora species create unacceptable risk and complicate inspection and certification programs designed to prevent introduction and spread of these pathogens. In this study, we continue development of a volatile organic compound (VOC)-based test for detecting asymptomatic infections of P. ramorum in Rhododendron sp. We confirmed detection of P. ramorum from volatiles collected from asymptomatic root-inoculated Rhododendron plants in a nursery setting, finding that the VOC profile of infected plants is detectably different from that of healthy plants, when measured from both ambient VOC emissions and VOCs extracted from leaf material. Predicting infection status was successful from ambient volatiles, which had a mean area under the curve (AUC) value of 0.71 ± 0.17, derived from corresponding receiver operating characteristic curves from an extreme gradient boosting discriminant analysis. This finding compares with that of extracted leaf volatiles, which resulted in a lower AUC value of 0.51 ± 0.21. In a growth chamber, we contrasted volatile profiles of asymptomatic Rhododendron plants having roots infected with one of three pathogens: P. ramorum, P. cactorum, and Rhizoctonia solani. Each pathogen induced unique and measurable changes, but generally the infections reduced volatile emissions until 17 weeks after inoculation, when emissions trended upward relative to those of mock-inoculated controls. Forty-five compounds had significant differences compared with mock-inoculated controls in at least one host-pathogen combination.
Phosphites have been used to control Sudden Oak Death; however, their precise mode of action is not fully understood. To study the mechanism of action of phosphites, we conducted an inoculation experiment on two open-pollinated tanoak families, previously found to be partially resistant. Stems of treatment group individuals were sprayed with phosphite, and seven days later, distal leaves were inoculated with the Sudden Oak Death pathogen Phytophthora ramorum. Leaves from treated and untreated control plants were harvested before and seven days after inoculation, and transcriptomes of both host and pathogen were analyzed. We found that tanoak families differed in the presence of innate resistance (resistance displayed by untreated tanoak) and in the response to phosphite treatment. A set of expressed genes associated with innate resistance was found to overlap with an expressed gene set for phosphite-induced resistance. This observation may indicate that phosphite treatment increases the resistance of susceptible host plants. In addition, genes of the pathogen involved in detoxification were upregulated in phosphite-treated plants compared to phosphite-untreated plants. In summary, our RNA-Seq analysis supports a two-fold mode of action of phosphites, including a direct toxic effect on P. ramorum and an indirect enhancement of resistance in the tanoak host.
Phytophthora megakarya and P. palmivora are oomycete pathogens that cause black pod rot of cacao (Theobroma cacao), the most economically important disease on cacao globally. While P. palmivora is a cosmopolitan pathogen, P. megakarya, which is more aggressive on cacao than P. palmivora, has been reported only in West and Central Africa where it has been spreading and devastating cacao farms since the 1950s. In this study, we reconstructed the complete diploid genomes of multiple isolates of both species using single-molecule real-time sequencing. Thirty-one additional genotypes were sequenced to analyze inter- and intra-species genomic diversity. The P. megakarya genome is exceptionally large (222 Mbp) and nearly twice the size of P. palmivora (135 Mbp) and most known Phytophthora species (∼100 Mbp on average). Previous reports pointed toward a whole-genome duplication (WGD) in P. palmivora. In this study, we demonstrate that both species underwent independent and relatively recent WGD events. In P. megakarya we identified a unique combination of WGD and large-scale transposable element driven genome expansion, which places this genome in the upper range of Phytophthora genome sizes, as well as effector pools with 1,382 predicted RxLR effectors. Finally, this study provides evidence of adaptive evolution of effectors like RxLRs and Crinklers, and discusses the implications of effector expansion and diversification.
It is commonly assumed that asexual lineages are short-lived evolutionarily, yet many asexual organisms can generate genetic and phenotypic variation, providing an avenue for further evolution. Previous work on the asexual plant pathogen Phytophthora ramorum NA1 revealed considerable genetic variation in the form of Structural Variants (SVs). To better understand how SVs arise and their significance to the California NA1 population, we studied the evolutionary histories of SVs and the forest conditions associated with their emergence. Ancestral state reconstruction suggests that SVs arose by somatic mutations among multiple independent lineages, rather than by recombination. We asked if this unusual phenomenon of parallel evolution between isolated populations is transmitted to extant lineages and found that SVs persist longer in a population if their genetic background had a lower mutation load. Genetic parallelism was also found in geographically distant demes where forest conditions such as host density, solar radiation, and temperature, were similar. Parallel SVs overlap with genes involved in pathogenicity such as RXLRs and have the potential to change the course of an epidemic. By combining genomics and environmental data, we identified an unexpected pattern of repeated evolution in an asexual population and identified environmental factors potentially driving this phenomenon.
The NA1 clonal lineage of Phytophthora ramorum is responsible for sudden oak death, an epidemic that has devastated California coastal forest ecosystems. An NA1 isolate, Pr102, derived from coast live oak in California, was previously sequenced and reported with a 65Mb assembly containing 12 Mb of gaps in 2,576 scaffolds. Here, we report an improved 70-Mb genome in 1,512 scaffolds with 6,752 bp of gaps after incorporating PacBio P5-C3 long reads. This assembly contains 19,494 gene models (average gene length of 2,515 bp) compared with 16,134 genes (average gene length of 1,673 bp) in the previous version. We predicted 29 new RXLR genes and 76 new paralogs of a total 392 RXLR genes from this assembly. We predicted 35 CRN genes compared with 19 in an earlier version with six paralogs. Our long non-coding RNA prediction identified 255 candidates. This new resource will be invaluable for future evolution studies on the invasive plant pathogen.