Search

The complete genome of a novel mycovirus, Colletotrichum curcumae narnvirus 1 (CcNV1), derived from the phytopathogenic fungus Colletotrichum curcumae strain 780-2T was sequenced and analyzed. The full sequence of CcNV1 is 3,374 bp in length and contains a single large open reading frame (ORF) encoding an RNA-dependent RNA polymerase (RdRp) of 1,087 amino acids with a molecular mass of 124.2 kDa that shares the closest similarity with that of Monilinia narnavirus H (53.02% identity). RdRp phylogeny analysis showed that CcNV1 is a new member of the genus Betanarnavirus within family Narnaviridae. This is the first report of a novel narnavirus infecting the phytopathogenic fungus C. curcumae, the causal agent of leaf blight of Curcuma wenyujin.

All tissue development and replenishment relies upon the breaking of symmetries leading to the morphological and operational differentiation of progenitor cells into more specialized cells. One of the main engines driving this process is the Notch signal transduction pathway, a ubiquitous signalling system found in the vast majority of metazoan cell types characterized to date. Broadly speaking, Notch receptor activity is governed by a balance between two processes: 1) intercellular Notch transactivation triggered via interactions between receptors and ligands expressed in neighbouring cells; 2) intracellular cis inhibition caused by ligands binding to receptors within the same cell. Additionally, recent reports have also unveiled evidence of cis activation. Whilst context-dependent Notch receptor clustering has been hypothesized, to date, Notch signalling has been assumed to involve an interplay between receptor and ligand monomers. In this study, we demonstrate biochemically, through a mutational analysis of DLL4, both in vitro and in tissue culture cells, that Notch ligands can efficiently self-associate. We found that the membrane proximal EGF-like repeat of DLL4 was necessary and sufficient to promote oligomerization/dimerization. Mechanistically, our experimental evidence supports the view that DLL4 ligand dimerization is specifically required for cis-inhibition of Notch receptor activity. To further substantiate these findings, we have adapted and extended existing ordinary differential equation-based models of Notch signalling to take account of the ligand dimerization-dependent cis-inhibition reported here. Our new model faithfully recapitulates our experimental data and improves predictions based upon published data. Collectively, our work favours a model in which net output following Notch receptor/ligand binding results from ligand monomer-driven Notch receptor transactivation (and cis activation) counterposed by ligand dimer-mediated cis-inhibition.Author summaryThe growth and maintenance of tissues is a fundamental characteristic of metazoan life, controlled by a highly conserved core of cell signal transduction networks. One such pathway, the Notch signalling system, plays a unique role in these phenomena by orchestrating the generation of the phenotypic and genetic asymmetries which underlie tissue development and remodeling. At the molecular level, it achieves this via two specific types of receptor/ligand interaction: intercellular binding of receptors and ligands expressed in neighbouring cells, which triggers receptor activation (transactivation); intracellular receptor/ligand binding within the same cell which blocks receptor activation (cis inhibition). Together, these counterposed mechanisms determine the strength, the direction and the specificity of Notch signalling output. Whilst, the basic mechanisms of receptor transactivation have been delineated in some detail, the precise nature of cis inhibition has remained enigmatic. Through a combination of experimental approaches and computational modelling, in this study, we present a new model of Notch signalling in which ligand monomers promote Notch receptor transactivation, whereas cis inhibition is induced optimally via ligand dimers. This is the first model to include a concrete molecular distinction, in terms of ligand configuration, between the main branches of Notch signalling. Our model faithfully recapitulates both our presented experimental results as well as the recently published work of others, and provides a novel perspective for understanding Notch-regulated biological processes such as embryo development and angiogenesis.

Powdery mildew severely affects several important crops and cash plants. Disruption of

Most models of COVID-19 are implemented at a single micro or macro scale, ignoring the interplay between immune response, viral dynamics, individual infectiousness and epidemiological contact networks. Here we develop a data-driven model linking the within-host viral dynamics to the between-host transmission dynamics on a multilayer contact network to investigate the potential factors driving transmission dynamics and to inform how school closures and antiviral treatment can influence the epidemic. Using multi-source data, we initially determine the viral dynamics and estimate the relationship between viral load and infectiousness. Then, we embed the viral dynamics model into a four-layer contact network and formulate an agent-based model to simulate between-host transmission. The results illustrate that the heterogeneity of immune response between children and adults and between vaccinated and unvaccinated infections can produce different transmission patterns. We find that school closures play a significant effect on mitigating the pandemic as more adults get vaccinated and the virus mutates. If enough infected individuals are diagnosed by testing before symptom onset and then treated quickly, the transmission can be effectively curbed. Our multiscale model reveals the critical role played by younger individuals and antiviral treatment with testing in controlling the epidemic.

The

Alternative splicing (AS) and alternative polyadenylation (APA) contribute significantly to the regulation of gene expression in higher eukaryotes. Their biological impact in filamentous fungi, however, is largely unknown. Here we combine PacBio Isoform-Sequencing and strand-specific RNA-sequencing of multiple tissues and mutant characterization to reveal the landscape and regulation of AS and APA in Fusarium graminearum. We generated a transcript annotation comprising 51 617 isoforms from 17 189 genes. In total, 4997 and 11 133 genes are alternatively spliced and polyadenylated, respectively. Majority of the AS events alter coding sequences. Unexpectedly, the AS transcripts containing premature-termination codons are not sensitive to nonsense-mediated messenger RNA decay. Unlike in yeasts and animals, distal APA sites have strong signals, but proximal APA isoforms are highly expressed in F. graminearum. The 3'-end processing factors FgRNA15, FgHRP1, and FgFIP1 play roles in promoting proximal APA site usage and intron splicing. A genome-wide increase in intron inclusion and distal APA site usage and downregulation of the spliceosomal and 3'-end processing factors were observed in older and quiescent tissues, indicating intron inclusion and 3'-untranslated region lengthening as novel mechanisms in regulating aging and dormancy in fungi. This study provides new insights into the complexity and regulation of AS and APA in filamentous fungi.

Alternative splicing (AS) and alternative polyadenylation (APA) of pre-mRNAs contribute greatly to transcriptome complexity and gene expression regulation in higher eukaryotes. Their biological impact in filamentous fungi, however, has been poorly studied. Here we combine PacBio Isoform Sequencing and strand-specific RNA-Seq of multiple tissues together with mutant characterization to reveal the landscape, complexity and regulation of AS and APA in the filamentous plant pathogenic fungus Fusarium graminearum. We updated the reference genome and generated a comprehensive annotation comprising 51,617 transcript isoforms from 17,189 genes. Majority of the transcripts represent novel isoforms, including 2,998 undiscovered protein-coding genes. In total, 42.7% of multi-exonic genes and 64.8% of genes have AS and APA isoforms, respectively, suggesting AS and APA increase previously unrecognized transcriptome complexity in fungi. Nonsense-mediated mRNA decay factor FgUPF1 may not degrade AS transcripts with premature-stop codons but regulate ribosome biogenesis. Distal polyadenylation sites have a strong signal but proximal polyadenylation isoforms are high expressed. The core 3’-end processing factors FgRNA15, FgHRP1, and FgFIP1 play important roles in promoting proximal polyadenylation site usage and also intron splicing. Genome-wide increase in the abundance of transcripts with retained introns and long 3’-UTRs and downregulation of the spliceosomal and 3’-end processing factors are found in older tissues and quiescent conidia, indicating that intron retention and 3’-UTR lengthening may be a transcriptional signature of aging and dormancy in fungi. Overall, our study generates a comprehensive full-length transcript annotation for F. graminearum and provides new insights into the complexity and regulation of transcriptome in filamentous fungi.

Here, a novel umbra-like mycovirus, ‘Phoma matteucciicola RNA virus 2’ (PmRV2), isolated from Phoma matteucciicola strain HNQH1 in Hainan province of China, was sequenced and analyzed. The complete genomic sequence of PmRV2 is 3,460 nucleotides (nts) with a GC content of 56.71%. Sequence analysis of PmRV2 indicated that the presence of two noncontiguous open reading frames (ORFs) encoding a hypothetical protein and a RNA-dependent RNA polymerase (RdRp), respectively. PmRV2 contains a metal-binding ‘GDN’ triad in Motif C of RdRp while most + ssRNA mycoviruses contained a ‘GDD’ motif in the same region. Additionally, a BLASTp search showed that the RdRp amino acid sequence of PmRV2 was most closely related to the RdRp of Macrophomina phaseolina umbra-like virus 1 (50.72% identity) and Erysiphe necator umbra-like virus 2 (44.84% identity). Phylogenetic analysis indicated that PmRV2 grouped together with Erysiphe necator umbra-like virus 2 (EnUlV2) within the recently proposed family of ‘Mycotombusviridae’.

A novel hexa-segmented double-stranded RNA (dsRNA) mycovirus was isolated and characterized from the filamentous phytopathogenic fungus Diaporthe pseudophoenicicola and was named Diaporthe pseudophoenicicola chrysovirus 1 (DpCV1). The full-length cDNAs of dsRNA1-6 were 3335, 3030, 3039, 2980, 963, and 780 bp, respectively. Sequence analysis indicated the presence of nine open reading frames (ORFs) in the DpCV1 genome. ORF1 in dsRNA1 putatively encoded the RNA-dependent RNA polymerase (RdRp) and ORF3 in dsRNA2 encoded a capsid protein (CP). The seven remaining ORFs, ORF2 in dsRNA2, ORF4 in dsRNA3, ORF6, seven in dsRNA4, ORF8 in dsRNA5, and ORF9 in dsRNA6, encoded proteins with unknown functions. Phylogenetic analysis revealed that DpCV1 is closely related to members of the cluster I group within the family Chrysoviridae but formed a separate clade. Importantly, all the six segments of DpCV1 were cured successfully through single spore isolation to obtain the isogenic virus-free strains. DpCV1 can confer hypovirulence to the fungal host of Diaporthe pseudophoenicicola. Compared with the virus-free strain, WC02 harbouring the DpCV1 is more sensitive to fungicide prochloraz. Furthermore, the cell wall of DpCV1 infected strain was loose and enlarged. This is the first report of a hexa-segmented tentative chrysovirus in D. pseudophoenicicola.

The wheat head blight fungus Fusarium graminearum has two highly similar beta-tubulin genes with overlapping functions during vegetative growth but only TUB1 is important for sexual reproduction. To better understand their functional divergence during ascosporogenesis, in this study we characterized the sequence elements important for stage-specific functions of TUB1. Deletion of TUB1 blocked the late but not initial stages of perithecium formation. Perithecia formed by tub1 mutant had limited ascogenous hyphae and failed to develop asci. Silencing of TUB1 by MSUD also resulted in defects in ascospore formation. Interestingly, the 3'-UTR of TUB1 was dispensable for growth but essential for its function during sexual reproduction. RIP mutations that specifically affected Tub1 functions during sexual reproduction also were identified in two ascospore progeny. Furthermore, site-directed mutagenesis showed that whereas the non-editable mutations at three A-to-I RNA editing sites had no effects, the N347D (not T362D or I368V) edited mutation affected ascospore development. In addition, the F167Y, but not E198K or F200Y, mutation in TUB1 conferred tolerance to carbendazim and caused a minor defect in sexual reproduction. Taken together, our data indicate TUB1 plays an essential role in ascosporogenesis and sexual-specific functions of TUB1 require stage-specific RNA processing and Tub1 expression.

Infectious diseases attack humans from time to time and threaten the lives and survival of people all around the world. An important strategy to prevent the spatial spread of infectious diseases is to restrict population travel. With the reduction of the epidemic situation, when and where travel restrictions can be lifted, and how to organize orderly movement patterns become critical and fall within the scope of this study. We define a novel diffusion distance derived from the estimated mobility network, based on which we provide a general model to describe the spatiotemporal spread of infectious diseases with a random diffusion process and a deterministic drift process of the population. We consequently develop a multi-source data fusion method to determine the population flow in epidemic areas. In this method, we first select available subregions in epidemic areas, and then provide solutions to initiate new travel flux among these subregions. To verify our model and method, we analyze the multi-source data from mainland China and obtain a new travel flux triggering scheme in the selected 29 cities with the most active population movements in mainland China. The testable predictions in these selected cities show that reopening the borders in accordance with our proposed travel flux will not cause a second outbreak of COVID-19 in these cities. The finding provides a methodology of re-triggering travel flux during the weakening spread stage of the epidemic., Author summary Human infectious diseases spread from their origins to other places with population movements. In order to curb the spatial spread of infectious diseases, many countries and regions may introduce some travel restrictions when the epidemic is severe, and reopen the borders as the epidemic eases. This process involves some important issues such as the start and end time of travel restrictions, the geographical scope of the implementation of the exit strategy, and the allowable passenger flow on traffic lines. Here, we integrate multi-source data with a mathematical model, and consequently develop a new method to determine the travel flux in epidemic areas. As an application, we use this method to calculate when and where the travel restrictions targeting COVID-19 in China in early 2020 could be lifted, and how to optimize passenger flow along the traffic lines among the reopened cities. The testable predictions indicate that the population flow in accordance with our proposed movement pattern will not cause a resurgent outbreak of COVID-19 in the cities studied.

The filamentous ascomycete Fusarium graminearum contains two β-tubulin genes TUB1 and TUB2 that differ in functions during vegetative growth and sexual reproduction. To further characterize their functional relationship, in this study we determined the co-localization of Tub1 and Tub2 and assayed their expression levels in different mutants and roles in DON production. Tub1 co-localized with Tub2 to the same regions of microtubules in conidia, hyphae, and ascospores. Whereas deletion of TUB1 had no obvious effect on the transcription of TUB2 and two α-tubulin genes (TUB4 and TUB5), the tub2 mutant was up-regulated in TUB1 transcription. To assay their protein expression levels, polyclonal antibodies that could specifically detect four α- and β-tubulin proteins were generated. Western blot analyses showed that the abundance of Tub1 proteins was increased in tub2 but reduced in tub4 and tub5 mutants. Interestingly, protein expression of Tub4 and Tub5 was decreased in the tub1 mutant in comparison with the wild type, despite a lack of obvious changes in their transcription. In contrast, deletion of TUB2 had no effect on translation of TUB4 and TUB5. Ectopic expression of Tub2-mCherry partially recovered the growth defect of the tub1 mutant but did not rescue its defect in sexual reproduction. Expression of Tub1-GFP in the tub2 mutant also partially rescued its defects in vegetative growth, suggesting that disturbance in the balance of α- and β-tubulins contributes to mutant defects. The tub2 but not tub1 mutant was almost blocked in DON biosynthesis. Expression of TRI genes, toxisome formation, and DON-related cellular differentiation were significantly reduced in the tub2 mutant. Overall, our results showed that Tub1 and Tub2 share similar subcellular localization and have overlapping functions during vegetative growth but they differ in functions in DON production and ascosporogenesis in F. graminearum.

It has been observed that air pollution greatly affects respiratory infection and generates public health problem, but there are many challenges to quantifying the dynamics of air pollution and evaluating its impact on respiratory infections. A periodic Filippov system describing the state-dependent control strategy for air pollution, described by air quality index (AQI), is proposed. We theoretically analyze the non-autonomous Filippov subsystem for variation of AQI by converting into an autonomous Filippov system via increasing dimension of the system. We obtain that there is a unique periodic solution which is globally asymptotically stable. In particular, it shows that AQI stabilizes at either one of the periodic solutions of the free and control systems or a new periodic solution induced by the on-off control strategies, depending on the threshold levels. Then, a periodic system for respiratory infection with AQI-embedded transmission probability is formulated to examine the influence of air pollution on respiratory infection. We further obtain the control reproduction number (basic reproduction number with control strategies) of this periodic respiratory infection model. It is shown that the respiratory infection will go extinct if the control reproduction number is less than unity, while uniformly persists for larger than unity. The case study showed that the estimated threshold level coincides with the actual threshold which launches the traffic limitation measure in Xi'an, indicating the untimely density-dependent traffic limitation measure was ineffective in improving air quality. Numerically studies indicate that increasing the threshold level leads to an increase in the maximum value of the unique periodic solution for AQI and the control reproduction number for the AQI embedded SEIS model. These findings emphasize the importance of suitable threshold level to trigger interventions and suggest that untimely implementing control strategy may not effectively control air population.

Trichothecene mycotoxins, such as deoxynivalenol (DON) produced by the fungal pathogen, Fusarium graminearum, are not only important for plant infection but are also harmful to human and animal health. Trichothecene targets the ribosomal protein Rpl3 that is conserved in eukaryotes. Hence, a self-defense mechanism must exist in DON-producing fungi. It is reported that TRI (trichothecene biosynthesis) 101 and TRI12 are two genes responsible for self-defense against trichothecene toxins in Fusarium. In this study, however, we found that simultaneous disruption of TRI101 and TRI12 has no obvious influence on DON resistance upon exogenous DON treatment in F. graminearum, suggesting that other mechanisms may be involved in self-defense. By using RNA-seq, we identified 253 genes specifically induced in DON-treated cultures compared with samples from cultures treated or untreated with cycloheximide, a commonly used inhibitor of eukaryotic protein synthesis. We found that transporter genes are significantly enriched in this group of DON-induced genes. Of those genes, 15 encode major facilitator superfamily transporters likely involved in mycotoxin efflux. Significantly, we found that genes involved in the metabolism of gamma-aminobutyric acid (GABA), a known inducer of DON production in F. graminearum, are significantly enriched among the DON-induced genes. The GABA biosynthesis gene PROLINE UTILIZATION 2-2 (PUT2-2) is downregulated, while GABA degradation genes are upregulated at least twofold upon treatment with DON, resulting in decreased levels of GABA. Taken together, our results suggest that transporters influencing DON efflux are important for self-defense and that GABA mediates the balance of DON production and self-defense in F. graminearum.

This study systematically identified adenosine to inosine (A-to-I) editing sites in Neurospora crassa and showed the existence of stage-specific editing events at different sexual stages. Unlike in humans, fungal A-to-I editing mainly occurred in coding regions and caused nonsynonymous changes that significantly increased proteome complexity. In general, nonsynonymous editing sites in Neurospora are adaptive and favored by positive selection. RNA editing enables stage-specific functions or expression of proteins important for different sexual developmental processes. Some editing events are well conserved and may affect genes important for other genetic and epigenetic phenomena occurring during sexual reproduction. Overall, our results provide insights into the complex regulation of sexual development and reveal the role of A-to-I editing for adaptive evolution in Neurospora.

Although fungi lack adenosine deaminase acting on RNA (ADAR) enzymes, adenosine to inosine (A-to-I) RNA editing was reported recently in Fusarium graminearum during sexual reproduction. In this study, we profiled the A-to-I editing landscape and characterized its functional and adaptive properties in the model filamentous fungus Neurospora crassa. A total of 40,677 A-to-I editing sites were identified, and approximately half of them displayed stage-specific editing or editing levels at different sexual stages. RNA-sequencing analysis with the Δstc-1 and Δsad-1 mutants confirmed A-to-I editing occurred before ascus development but became more prevalent during ascosporogenesis. Besides fungal-specific sequence and secondary structure preference, 63.5% of A-to-I editing sites were in the coding regions and 81.3% of them resulted in nonsynonymous recoding, resulting in a significant increase in the proteome complexity. Many genes involved in RNA silencing, DNA methylation, and histone modifications had extensive recoding, including sad-1, sms-3, qde-1, and dim-2. Fifty pseudogenes harbor premature stop codons that require A-to-I editing to encode full-length proteins. Unlike in humans, nonsynonymous editing events in N. crassa are generally beneficial and favored by positive selection. Almost half of the nonsynonymous editing sites in N. crassa are conserved and edited in Neurospora tetrasperma. Furthermore, hundreds of them are conserved in F. graminearum and had higher editing levels. Two unknown genes with editing sites conserved between Neurospora and Fusarium were experimentally shown to be important for ascosporogenesis. This study comprehensively analyzed A-to-I editing in N. crassa and showed that RNA editing is stage-specific and generally adaptive, and may be functionally related to repeat induced point mutation and meiotic silencing by unpaired DNA.

Fusarium head blight caused by Fusarium graminearum is an important disease of wheat and barley worldwide. In a previous study on functional characterization of the F. graminearum kinome, one protein kinase gene important for virulence is orthologous to SCH9 that is functionally related to the cAMP-PKA and TOR pathways in the budding yeast. In this study, we further characterized the functions of FgSCH9 in F. graminearum and its ortholog in Magnaporthe oryzae. The ΔFgsch9 mutant was slightly reduced in growth rate but significantly reduced in conidiation, DON production, and virulence on wheat heads and corn silks. It had increased tolerance to elevated temperatures but became hypersensitive to oxidative, hyperosmotic, cell wall, and membrane stresses. The ΔFgsch9 deletion also had conidium morphology defects and produced smaller conidia. These results suggest that FgSCH9 is important for stress responses, DON production, conidiogenesis, and pathogenesis in F. graminearum. In the rice blast fungus Magnaporthe oryzae, the ΔMosch9 mutant also was defective in conidiogenesis and pathogenesis. Interestingly, it also produced smaller conidia and appressoria. Taken together, our data indicate that the SCH9 kinase gene may have a conserved role in regulating conidium size and plant infection in phytopathogenic ascomycetes.