De novo assembly and annotation of three Leptosphaeria genomes using Oxford Nanopore MinION sequencing

Regulation of effector gene expression as concerted waves in Leptosphaeria maculans: a two-player game

Effector genes, encoding molecules involved in disease establishment, are concertedly expressed throughout the lifecycle of plant-pathogenic fungi. However, little is known about how effector gene expression is regulated. Since many effector genes are located in repeat-rich regions, the role of chromatin remodeling in their regulation was recently investigated, notably establishing that the repressive histone modification H3K9me3, deposited by KMT1, was involved in several fungal species including Leptosphaeria maculans. Nevertheless, previous data suggest that a second regulatory layer, probably involving a specific transcription factor (TF), might be required. In L. maculans, a Dothideomycete causing stem canker of oilseed rape, we identified the ortholog of Pf2, a TF belonging to the Zn2Cys6 fungal-specific family, and described as essential for pathogenicity and effector gene expression. We investigated its role together with KMT1, by inactivating and over-expressing LmPf2 in a wild-type strain and a ∆kmt1 mutant. Functional analyses of the corresponding transformants highlighted an essential role of LmPf2 in the establishment of pathogenesis and we found a major effect of LmPf2 on the induction of effector gene expression once KMT1 repression is lifted. Our results show, for the first time, a dual control of effector gene expression.

Location and timing govern tripartite interactions of fungal phytopathogens and host in the stem canker species complex

BACKGROUND

Leptosphaeria maculans "brassicae" (Lmb) and Leptosphaeria biglobosa "brassicae" (Lbb) make up a species complex involved in the stem canker (blackleg) disease of rapeseed (Brassica napus). They coinfect rapeseed together, from the early stage of infection on leaves to the final necrotic stage at the stem base, and both perform sexual crossings on plant residues. L. biglobosa is suggested to be a potential biocontrol agent against Lmb, but there has been no mechanistic investigation of the different types of interactions that may occur between the plant and the two fungal species.

RESULTS

We investigated the bi- or tripartite interaction mechanisms by (i) confronting Lmb and Lbb in culture conditions or during cotyledon infection, with different timing and/or spore concentration regimes, (ii) performing RNA-Seq experiments in vitro or on the kinetics of infection of cotyledons infected by Lmb and/or Lbb to evaluate the transcriptomic activity and the plant response when both fungal species are inoculated together. Lbb infection of B. napus cotyledons was typical of a necrotrophic behavior, with a very early setup of one pathogenicity program and very limited colonization of tissues. This contrasted with the complex succession of pathogenicity programs of the hemibiotroph Lmb. During simultaneous co-infection by both species, Lmb was strongly impacted in its growth and transcriptomic dynamics both in vitro and in planta, while Lbb was unaffected by the presence of Lmb. However, the drastic inhibition of Lmb growth by Lbb was ineffective in the case of delayed inoculation with Lbb or a lower amount of spores of Lbb compared to Lmb.

CONCLUSIONS

Our data suggest that Lmb growth inhibition by Lbb is the result of a combination of factors that may include competition for trophic resources, the generation by Lbb of an environment unsuitable for the lifecycle of Lmb or/and the effect on Lmb of plant defense responses induced by Lbb. It indicates that growth inhibition occurs in very specific conditions (i.e., co-inoculation at the same place of an equal amount of inoculum) that are unlikely to occur in the field where their coexistence does not prevent any species from completing their life cycle.

Deep amplicon sequencing reveals extensive allelic diversity in the erg11/CYP51 promoter and allows multi-population DMI fungicide resistance monitoring in the canola pathogen Leptosphaeria maculans

Continued use of fungicides provides a strong selection pressure towards strains with mutations to render these chemicals less effective. Previous research has shown that resistance to the demethylation inhibitor (DMI) fungicides, which target ergosterol synthesis, in the canola pathogen Leptosphaeria maculans has emerged in Australia and Europe. The change in fungicide sensitivity of individual isolates was found to be due to DNA insertions into the promoter of the erg11/CYP51 DMI target gene. Whether or not these were the only types of mutations and how prevalent they were in Australian populations was explored in the current study. New isolates with reduced DMI sensitivity were obtained from screens on DMI-treated plants, revealing eight independent insertions in the erg11 promoter. A novel deep amplicon sequencing approach applied to populations of ascospores fired from stubble identified an additional undetected insertion allele and quantified the frequencies of all known insertions, suggesting that, at least in the samples processed, the combined frequency of resistant alleles is between 0.0376% and 32.6%. Combined insertion allele frequencies positively correlated with population-level measures of in planta resistance to four different DMI treatments. Additionally, there was no evidence for erg11 coding mutations playing a role in conferring resistance in Australian populations. This research provides a key method for assessing fungicide resistance frequency in stubble-borne populations of plant pathogens and a baseline from which additional surveillance can be conducted in L. maculans. Whether or not the observed resistance allele frequencies are associated with loss of effective disease control in the field remains to be established.

The neighbouring genes AvrLm10A and AvrLm10B are part of a large multigene family of cooperating effector genes conserved in Dothideomycetes and Sordariomycetes

Fungal effectors (small-secreted proteins) have long been considered as species or even subpopulation-specific. The increasing availability of high-quality fungal genomes and annotations has allowed the identification of trans-species or trans-genera families of effectors. Two avirulence effectors, AvrLm10A and AvrLm10B, of Leptosphaeria maculans, the fungus causing stem canker of oilseed rape, are members of such a large family of effectors. AvrLm10A and AvrLm10B are neighbouring genes, organized in divergent transcriptional orientation. Sequence searches within the L. maculans genome showed that AvrLm10A/AvrLm10B belong to a multigene family comprising five pairs of genes with a similar tail-to-tail organization. The two genes, in a pair, always had the same expression pattern and two expression profiles were distinguished, associated with the biotrophic colonization of cotyledons and/or petioles and stems. Of the two protein pairs further investigated, AvrLm10A_like1/AvrLm10B_like1 and AvrLm10A_like2/AvrLm10B_like2, the second one had the ability to physically interact, similarly to what was previously described for the AvrLm10A/AvrLm10B pair, and cross-interactions were also detected for two pairs. AvrLm10A homologues were identified in more than 30 Dothideomycete and Sordariomycete plant-pathogenic fungi. One of them, SIX5, is an effector from Fusarium oxysporum f. sp. lycopersici physically interacting with the avirulence effector Avr2. We found that AvrLm10A/SIX5 homologues were associated with at least eight distinct putative effector families, suggesting that AvrLm10A/SIX5 is able to cooperate with different effectors. These results point to a general role of the AvrLm10A/SIX5 proteins as "cooperating proteins", able to interact with diverse families of effectors whose encoding gene is co-regulated with the neighbouring AvrLm10A homologue.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

Polymorphism of Avirulence Genes and Adaptation to Brassica Resistance Genes Is Gene-Dependent in the Phytopathogenic Fungus Leptosphaeria maculans

The fungal phytopathogen Leptosphaeria maculans, which causes stem canker (blackleg) of rapeseed (Brassica napus), is mainly controlled worldwide by genetic resistance, which includes major resistance genes (Rlm). This model is one of those for which the highest number of avirulence genes (AvrLm) has been cloned. In many systems, including the L. maculans-B. napus interaction, intense use of resistance genes exerts strong selection pressure on the corresponding avirulent isolates, and the fungi may rapidly escape resistance through various molecular events which modify the avirulence genes. In the literature, the study of polymorphism at avirulence loci is often focused on single genes under selection pressure. In this study, we investigate allelic polymorphism at 11 avirulence loci in a French population of 89 L. maculans isolates collected on a trap cultivar in four geographic locations in the 2017-2018 cropping season. The corresponding Rlm genes have been (i) used for a long time, (ii) recently used, or (iii) unused in agricultural practice. The sequence data generated indicate an extreme diversity of situations. For example, genes submitted to an ancient selection may have either been deleted in populations (AvrLm1) or replaced by a single-nucleotide mutated virulent version (AvrLm2, AvrLm5-9). Genes that have never been under selection may either be nearly invariant (AvrLm6, AvrLm10A, AvrLm10B), exhibit rare deletions (AvrLm11, AvrLm14), or display a high diversity of alleles and isoforms (AvrLmS-Lep2). These data suggest that the evolutionary trajectory of avirulence/virulence alleles is gene-dependent and independent of selection pressure in L. maculans. [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.

Unraveling the genomic complexity of sylvatic mosquitoes in changing Neotropical environments

Sylvatic New World mosquitoes (e.g. Old-growth Forest species) can transmit viruses among non-human primates. This could be a continuous source of viral cycling and spillover events from animals to humans, particularly in changing environments. However, most species of Neotropical sylvatic mosquitoes (genera Aedes, Haemagogus, and Sabethes), which include vector and non-vector species, currently lack genomic resources because there is no reliable and accurate approach for creating de novo reference genomes for these insects. This is a major knowledge gap in the biology of these mosquitoes, restricting our ability to predict and mitigate the emergence and spread of novel arboviruses in Neotropical regions. We discuss recent advances and potential solutions for generating hybrid de novo assemblies from vector and non-vector species using pools of consanguineous offspring. We also discussed research opportunities likely to emerge from these genomic resources.

Key Advances in the New Era of Genomics-Assisted Disease Resistance Improvement of Brassica Species

Disease resistance improvement remains a major focus in breeding programs as diseases continue to devastate Brassica production systems due to intensive cultivation and climate change. Genomics has paved the way to understand the complex genomes of Brassicas, which has been pivotal in the dissection of the genetic underpinnings of agronomic traits driving the development of superior cultivars. The new era of genomics-assisted disease resistance breeding has been marked by the development of high-quality genome references, accelerating the identification of disease resistance genes controlling both qualitative (major) gene and quantitative resistance. This facilitates the development of molecular markers for marker assisted selection and enables genome editing approaches for targeted gene manipulation to enhance the genetic value of disease resistance traits. This review summarizes the key advances in the development of genomic resources for Brassica species, focusing on improved genome references, based on long-read sequencing technologies and pangenome assemblies. This is further supported by the advances in pathogen genomics, which have resulted in the discovery of pathogenicity factors, complementing the mining of disease resistance genes in the host. Recognizing the co-evolutionary arms race between the host and pathogen, it is critical to identify novel resistance genes using crop wild relatives and synthetic cultivars or through genetic manipulation via genome-editing to sustain the development of superior cultivars. Integrating these key advances with new breeding techniques and improved phenotyping using advanced data analysis platforms will make disease resistance improvement in Brassica species more efficient and responsive to current and future demands.

CusProSe: a customizable protein annotation software with an application to the prediction of fungal secondary metabolism genes

We report here a new application, CustomProteinSearch (CusProSe), whose purpose is to help users to search for proteins of interest based on their domain composition. The application is customizable. It consists of two independent tools, IterHMMBuild and ProSeCDA. IterHMMBuild allows the iterative construction of Hidden Markov Model (HMM) profiles for conserved domains of selected protein sequences, while ProSeCDA scans a proteome of interest against an HMM profile database, and annotates identified proteins using user-defined rules. CusProSe was successfully used to identify, in fungal genomes, genes encoding key enzyme families involved in secondary metabolism, such as polyketide synthases (PKS), non-ribosomal peptide synthetases (NRPS), hybrid PKS-NRPS and dimethylallyl tryptophan synthases (DMATS), as well as to characterize distinct terpene synthases (TS) sub-families. The highly configurable characteristics of this application makes it a generic tool, which allows the user to refine the function of predicted proteins, to extend detection to new enzymes families, and may also be applied to biological systems other than fungi and to other proteins than those involved in secondary metabolism.

Indirect Evidence Based on Mating-Type Ratios for the Role of Sexual Reproduction in European and Chinese Populations of Plenodomus biglobosus (Blackleg of Oilseed Rape)

Blackleg (Phoma) disease, caused by the ascomycete fungi Plenodomus biglobosus and P. lingam, threatens oilseed rape (OSR; Brassica napus) crops internationally. In many parts of the world, both species co-occur, but in China only P. biglobosus has so far been reported. Plenodomus biglobosus reproduces asexually (pycnidiospores), but also sexually (pseudothecia-yielding ascospores), via a heterothallic mating system requiring MAT1-1 and MAT1-2 genotypes. However, the roles of airborne ascospore inoculum in driving blackleg disease outbreaks in China are less well understood compared to elsewhere in the world. This is despite the very different agronomic cropping practices in parts of China, in which paddy rice and OSR are often grown in rotation; OSR stubble is often submerged under water for long periods potentially affecting pseudothecial development. Here, we indirectly investigate the potential role of sexual reproduction by developing new polymerase chain reaction (PCR) -based mating-type diagnostics for P. biglobosus and subsequently screening an international collection of 59 European and 157 Chinese isolates. Overall, in both Europe and China, P. biglobosus mating types did not deviate from a 1:1 ratio, such as is generally thought to occur under frequency-dependent selection in sexually reproducing pathogen populations. Both mating types were balanced in all the individual European countries tested (Austria, France, Poland, UK). Conversely, in China, mating types were only balanced in the eastern region; in the northern and southwestern regions there were skewed ratios, more typical of predominantly asexual reproduction, towards MAT1-1 and MAT1-2, respectively. The implications of these findings and future research directions for improved understanding of P. biglobosus epidemiology on OSR, particularly in China, are considered.

"Late" effectors from Leptosphaeria maculans as tools for identifying novel sources of resistance in Brassica napus

The Dothideomycete Leptosphaeria maculans, causing stem canker (blackleg) of Brassica napus, secretes different cocktails of effectors at specific infection stages. Some effectors ("Late" effectors) are specifically produced during the long asymptomatic phase of stem colonization. By manipulating their expression so that they are overexpressed during cotyledon infection (OEC transformants of the fungus), we previously postulated that resistance genes operating in the stem may be involved in gene-for-gene relationship and thus contribute to quantitative disease resistance (QDR). Here, we selected 10 relevant new effector genes, and we generated OEC transformants to screen a collection of 130 B. napus genotypes, representative of the available diversity in the species. Five B. napus accessions showed a typical hypersensitive response when challenged with effectors LmSTEE98 or LmSTEE6826 at the cotyledon stage, and all belong to the semi-winter type of the diversity panel. In addition, five winter-type genotypes displayed an intermediate response to another late effector, LmSTEE7919. These new interactions now have to be genetically validated to check that they also correspond to gene-for-gene interactions. In all cases, they potentially provide novel resources, easy to breed for, and accounting for part of the quantitative resistance in a species for which we are currently facing limited resistance sources.

A new family of structurally conserved fungal effectors displays epistatic interactions with plant resistance proteins

Recognition of a pathogen avirulence (AVR) effector protein by a cognate plant resistance (R) protein triggers a set of immune responses that render the plant resistant. Pathogens can escape this so-called Effector-Triggered Immunity (ETI) by different mechanisms including the deletion or loss-of-function mutation of the AVR gene, the incorporation of point mutations that allow recognition to be evaded while maintaining virulence function, and the acquisition of new effectors that suppress AVR recognition. The Dothideomycete Leptosphaeria maculans, causal agent of oilseed rape stem canker, is one of the few fungal pathogens where suppression of ETI by an AVR effector has been demonstrated. Indeed, AvrLm4-7 suppresses Rlm3- and Rlm9-mediated resistance triggered by AvrLm3 and AvrLm5-9, respectively. The presence of AvrLm4-7 does not impede AvrLm3 and AvrLm5-9 expression, and the three AVR proteins do not appear to physically interact. To decipher the epistatic interaction between these L. maculans AVR effectors, we determined the crystal structure of AvrLm5-9 and obtained a 3D model of AvrLm3, based on the crystal structure of Ecp11-1, a homologous AVR effector candidate from Fulvia fulva. Despite a lack of sequence similarity, AvrLm5-9 and AvrLm3 are structural analogues of AvrLm4-7 (structure previously characterized). Structure-informed sequence database searches identified a larger number of putative structural analogues among L. maculans effector candidates, including the AVR effector AvrLmS-Lep2, all produced during the early stages of oilseed rape infection, as well as among effector candidates from other phytopathogenic fungi. These structural analogues are named LARS (for Leptosphaeria AviRulence and Suppressing) effectors. Remarkably, transformants of L. maculans expressing one of these structural analogues, Ecp11-1, triggered oilseed rape immunity in several genotypes carrying Rlm3. Furthermore, this resistance could be suppressed by AvrLm4-7. These results suggest that Ecp11-1 shares a common activity with AvrLm3 within the host plant which is detected by Rlm3, or that the Ecp11-1 structure is sufficiently close to that of AvrLm3 to be recognized by Rlm3.

An efficient strategy to control fungal diseases in the field is genetic control using resistant crop cultivars. Crop resistance mainly relies on gene-for-gene relationships between plant resistance (R) genes and pathogen avirulence (AVR) genes, as defined by Flor in the 1940s. However, such gene-for-gene relationships can increase in complexity over the course of plant-pathogen co-evolution. Resistance against the plant-pathogenic fungus Leptosphaeria maculans by Brassica napus and other Brassica species relies on the recognition of effector (AVR) proteins by R proteins; however, L. maculans produces an effector that suppresses a subset of these specific resistances. Using a protein structure approach, we revealed structural analogy between several of the resistance-triggering effectors, the resistance-suppressing effector, and effectors from other plant-pathogenic species in the Dothideomycetes and Sordariomycetes classes, defining a new family of effectors called LARS. Notably, cross-species expression of one LARS effector from Fulvia fulva, a pathogen of tomato, in L. maculans resulted in recognition by resistant cultivars of oilseed rape. These results highlight the need to integrate knowledge on effector structures to improve resistance management and to develop broad-spectrum resistances for multi-pathogen control of diseases.

A new family of structurally conserved fungal effectors displays epistatic interactions with plant resistance proteins

Recognition of a pathogen avirulence (AVR) effector protein by a cognate plant resistance (R) protein triggers a set of immune responses that render the plant resistant. Pathogens can escape this so-called Effector-Triggered Immunity (ETI) by different mechanisms including the deletion or loss-of-function mutation of the AVR gene, the incorporation of point mutations that allow recognition to be evaded while maintaining virulence function, and the acquisition of new effectors that suppress AVR recognition. The Dothideomycete Leptosphaeria maculans, causal agent of oilseed rape stem canker, is one of the few fungal pathogens where suppression of ETI by an AVR effector has been demonstrated. Indeed, AvrLm4-7 suppresses Rlm3- and Rlm9-mediated resistance triggered by AvrLm3 and AvrLm5-9, respectively. The presence of AvrLm4-7 does not impede AvrLm3 and AvrLm5-9 expression, and the three AVR proteins do not appear to physically interact. To decipher the epistatic interaction between these L. maculans AVR effectors, we determined the crystal structure of AvrLm5-9 and obtained a 3D model of AvrLm3, based on the crystal structure of Ecp11-1, a homologous AVR effector candidate from Fulvia fulva. Despite a lack of sequence similarity, AvrLm5-9 and AvrLm3 are structural analogues of AvrLm4-7 (structure previously characterized). Structure-informed sequence database searches identified a larger number of putative structural analogues among L. maculans effector candidates, including the AVR effector AvrLmS-Lep2, all produced during the early stages of oilseed rape infection, as well as among effector candidates from other phytopathogenic fungi. These structural analogues are named LARS (for Leptosphaeria AviRulence and Suppressing) effectors. Remarkably, transformants of L. maculans expressing one of these structural analogues, Ecp11-1, triggered oilseed rape immunity in several genotypes carrying Rlm3. Furthermore, this resistance could be suppressed by AvrLm4-7. These results suggest that Ecp11-1 shares a common activity with AvrLm3 within the host plant which is detected by Rlm3, or that the Ecp11-1 structure is sufficiently close to that of AvrLm3 to be recognized by Rlm3.

Two independent approaches converge to the cloning of a new Leptosphaeria maculans avirulence effector gene, AvrLmS-Lep2

Brassica napus (oilseed rape, canola) seedling resistance to Leptosphaeria maculans, the causal agent of blackleg (stem canker) disease, follows a gene-for-gene relationship. The avirulence genes AvrLmS and AvrLep2 were described to be perceived by the resistance genes RlmS and LepR2, respectively, present in B. napus 'Surpass 400'. Here we report cloning of AvrLmS and AvrLep2 using two independent methods. AvrLmS was cloned using combined in vitro crossing between avirulent and virulent isolates with sequencing of DNA bulks from avirulent or virulent progeny (bulked segregant sequencing). AvrLep2 was cloned using a biparental cross of avirulent and virulent L. maculans isolates and a classical map-based cloning approach. Taking these two approaches independently, we found that AvrLmS and AvrLep2 are the same gene. Complementation of virulent isolates with this gene confirmed its role in inducing resistance on Surpass 400, Topas-LepR2, and an RlmS-line. The gene, renamed AvrLmS-Lep2, encodes a small cysteine-rich protein of unknown function with an N-terminal secretory signal peptide, which is a common feature of the majority of effectors from extracellular fungal plant pathogens. The AvrLmS-Lep2/LepR2 interaction phenotype was found to vary from a typical hypersensitive response through intermediate resistance sometimes towards susceptibility, depending on the inoculation conditions. AvrLmS-Lep2 was nevertheless sufficient to significantly slow the systemic growth of the pathogen and reduce the stem lesion size on plant genotypes with LepR2, indicating the potential efficiency of this resistance to control the disease in the field.

High molecular weight DNA extraction methods lead to high quality filamentous ascomycete fungal genome assemblies using Oxford Nanopore sequencing

During the last two decades, whole-genome sequencing has revolutionized genetic research in all kingdoms, including fungi. More than 1000 fungal genomes have been submitted to sequence databases, mostly obtained through second generation short-read DNA sequencing. As a result, highly fragmented genome drafts have typically been obtained. However, with the emergence of third generation long-read DNA sequencing, the assembly challenge can be overcome and highly contiguous assemblies obtained. Such attractive results, however, are extremely dependent on the ability to extract highly purified high molecular weight (HMW) DNA. Extraction of such DNA is currently a significant challenge for all species with cell walls, not least fungi. In this study, four isolates of filamentous ascomycetes (Apiospora pterospermum, Aspergillus sp. (subgen. Cremei), Aspergillus westerdijkiae, and Penicillium aurantiogriseum) were used to develop extraction and purification methods that result in HMW DNA suitable for third generation sequencing. We have tested and propose two straightforward extraction methods based on treatment with either a commercial kit or traditional phenol-chloroform extraction both in combination with a single commercial purification method that result in high quality HMW DNA from filamentous ascomycetes. Our results demonstrated that using these DNA extraction methods and coverage, above 75 x of our haploid filamentous ascomycete fungal genomes result in complete and contiguous assemblies.

Hybrid Assembly Improves Genome Quality and Completeness of Trametes villosa CCMB561 and Reveals a Huge Potential for Lignocellulose Breakdown

Trametes villosa is a wood-decaying fungus with great potential to be used in the bioconversion of agro-industrial residues and to obtain high-value-added products, such as biofuels. Nonetheless, the lack of high-quality genomic data hampers studies investigating genetic mechanisms and metabolic pathways in T. villosa, hindering its application in industry. Herein, applying a hybrid assembly pipeline using short reads (Illumina HiSeq) and long reads (Oxford Nanopore MinION), we obtained a high-quality genome for the T. villosa CCMB561 and investigated its genetic potential for lignocellulose breakdown. The new genome possesses 143 contigs, N50 of 1,009,271 bp, a total length of 46,748,415 bp, 14,540 protein-coding genes, 22 secondary metabolite gene clusters, and 426 genes encoding Carbohydrate-Active enzymes. Our CAZome annotation and comparative genomic analyses of nine Trametes spp. genomes revealed T. villosa CCMB561 as the species with the highest number of genes encoding lignin-modifying enzymes and a wide array of genes encoding proteins for the breakdown of cellulose, hemicellulose, and pectin. These results bring to light the potential of this isolate to be applied in the bioconversion of lignocellulose and will support future studies on the expression, regulation, and evolution of genes, proteins, and metabolic pathways regarding the bioconversion of lignocellulosic residues.

A new avirulence gene of Leptosphaeria maculans, AvrLm14, identifies a resistance source in American broccoli (Brassica oleracea) genotypes

In many cultivated crops, sources of resistance to diseases are sparse and rely on introgression from wild relatives. Agricultural crops often are allopolyploids resulting from interspecific crosses between related species, which are sources of diversity for resistance genes. This is the case for Brassica napus (oilseed rape, canola), an interspecific hybrid between Brassica rapa (turnip) and Brassica oleracea (cabbage). B. napus has a narrow genetic basis and few effective resistance genes against stem canker (blackleg) disease, caused by the fungus Leptosphaeria maculans, are currently available. B. rapa diversity has proven to be a valuable source of resistance (Rlm, LepR) genes, while B. oleracea genotypes were mostly considered susceptible. Here we identified a new resistance source in B. oleracea genotypes from America, potentially effective against French L. maculans isolates under both controlled and field conditions. Genetic analysis of fungal avirulence and subsequent cloning and validation identified a new avirulence gene termed AvrLm14 and suggested a typical gene-for-gene interaction between AvrLm14 and the postulated Rlm14 gene. AvrLm14 shares all the usual characteristics of L. maculans avirulence genes: it is hosted in a genomic region enriched in transposable elements and heterochromatin marks H3K9me3, its expression is repressed during vegetative growth but shows a strong overexpression 5-9 days following cotyledon infection, and it encodes a small secreted protein enriched in cysteine residues with few matches in databases. Similar to the previously cloned AvrLm10-A, AvrLm14 contributes to reduce lesion size on susceptible cotyledons, pointing to a complex interplay between effectors promoting or reducing lesion development.

Genetic diversity and population structure of Leptosphaeria maculans isolates in Western Canada

Leptosphaeria maculans is a serious concern for canola production worldwide. For effective disease management, knowledge of the pathogen's genetic variability and population structure is a prerequisite. In this study, whole-genome sequencing was performed for 162 of 1590 L. maculans isolates collected in the years 2007-2008 and 2012-2014 in Western Canada. DNA variants in genome-wide and specific regions including avirulence (Avr) genes were characterized. A total of 31,870 high-quality polymorphic DNA variants were used to study L. maculans genetic diversity and population structure. Cluster analysis showed that 150 isolates were clustered into 2 main groups and 4 subgroups by DNA variants located in either Avr or small secreted protein-encoding genes and into 2 main groups and 6 subgroups by genome-wide variants. The analysis of nucleotide diversity and differentiation also confirmed genetic variation within a population and among populations. Principal component analysis with genome-wide variants showed that the isolates collected in 2012-2014 were more genetically diverse than those collected in 2007-2008. Population structure analysis discovered three distinct sub-populations. Although isolates from Saskatchewan and Alberta were of similar genetic composition, Manitoba isolates were highly diverse. Genome-wide association study detected DNA variants in genes AvrLm4-7, Lema_T86300, and Lema_T86310 associated with the years of collection.

An automated and combinative method for the predictive ranking of candidate effector proteins of fungal plant pathogens

Fungal plant-pathogens promote infection of their hosts through the release of 'effectors'-a broad class of cytotoxic or virulence-promoting molecules. Effectors may be recognised by resistance or sensitivity receptors in the host, which can determine disease outcomes. Accurate prediction of effectors remains a major challenge in plant pathology, but if achieved will facilitate rapid improvements to host disease resistance. This study presents a novel tool and pipeline for the ranking of predicted effector candidates-Predector-which interfaces with multiple software tools and methods, aggregates disparate features that are relevant to fungal effector proteins, and applies a pairwise learning to rank approach. Predector outperformed a typical combination of secretion and effector prediction methods in terms of ranking performance when applied to a curated set of confirmed effectors derived from multiple species. We present Predector ( https://github.com/ccdmb/predector ) as a useful tool for the ranking of predicted effector candidates, which also aggregates and reports additional supporting information relevant to effector and secretome prediction in a simple, efficient, and reproducible manner.

A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus

The control of stem canker disease of Brassica napus (rapeseed), caused by the fungus Leptosphaeria maculans is based largely on plant genetic resistance: single-gene specific resistance (Rlm genes) or quantitative, polygenic, adult-stage resistance. Our working hypothesis was that quantitative resistance partly obeys the gene-for-gene model, with resistance genes 'recognizing' fungal effectors expressed during late systemic colonization. Five LmSTEE (stem-expressed effector) genes were selected and placed under the control of the AvrLm4-7 promoter, an effector gene highly expressed at the cotyledon stage of infection, for miniaturized cotyledon inoculation test screening of a gene pool of 204 rapeseed genotypes. We identified a rapeseed genotype, 'Yudal', expressing hypersensitive response to LmSTEE98. The LmSTEE98-RlmSTEE98 interaction was further validated by inactivation of the LmSTEE98 gene with a CRISPR-Cas9 approach. Isolates with mutated versions of LmSTEE98 induced more severe stem symptoms than the wild-type isolate in 'Yudal'. This single-gene resistance was mapped in a 0.6 cM interval of the 'Darmor_bzh' × 'Yudal' genetic map. One typical gene-for-gene interaction contributes partly to quantitative resistance when L. maculans colonizes the stems of rapeseed. With numerous other effectors specific to stem colonization, our study provides a new route for resistance gene discovery, elucidation of quantitative resistance mechanisms and selection for durable resistance.

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

Large-scale transcriptomics to dissect 2 years of the life of a fungal phytopathogen interacting with its host plant

BACKGROUND

The fungus Leptosphaeria maculans has an exceptionally long and complex relationship with its host plant, Brassica napus, during which it switches between different lifestyles, including asymptomatic, biotrophic, necrotrophic, and saprotrophic stages. The fungus is also exemplary of "two-speed" genome organisms in the genome of which gene-rich and repeat-rich regions alternate. Except for a few stages of plant infection under controlled conditions, nothing is known about the genes mobilized by the fungus throughout its life cycle, which may last several years in the field.

RESULTS

We performed RNA-seq on samples corresponding to all stages of the interaction of L. maculans with its host plant, either alive or dead (stem residues after harvest) in controlled conditions or in field experiments under natural inoculum pressure, over periods of time ranging from a few days to months or years. A total of 102 biological samples corresponding to 37 sets of conditions were analyzed. We show here that about 9% of the genes of this fungus are highly expressed during its interactions with its host plant. These genes are distributed into eight well-defined expression clusters, corresponding to specific infection lifestyles or to tissue-specific genes. All expression clusters are enriched in effector genes, and one cluster is specific to the saprophytic lifestyle on plant residues. One cluster, including genes known to be involved in the first phase of asymptomatic fungal growth in leaves, is re-used at each asymptomatic growth stage, regardless of the type of organ infected. The expression of the genes of this cluster is repeatedly turned on and off during infection. Whatever their expression profile, the genes of these clusters are enriched in heterochromatin regions associated with H3K9me3 or H3K27me3 repressive marks. These findings provide support for the hypothesis that part of the fungal genes involved in niche adaptation is located in heterochromatic regions of the genome, conferring an extreme plasticity of expression.

CONCLUSION

This work opens up new avenues for plant disease control, by identifying stage-specific effectors that could be used as targets for the identification of novel durable disease resistance genes, or for the in-depth analysis of chromatin remodeling during plant infection, which could be manipulated to interfere with the global expression of effector genes at crucial stages of plant infection.

Current State and Future Directions of Genetics and Genomics of Endophytic Fungi for Bioprospecting Efforts

The bioprospecting of secondary metabolites from endophytic fungi received great attention in the 1990s and 2000s, when the controversy around taxol production from Taxus spp. endophytes was at its height. Since then, hundreds of reports have described the isolation and characterization of putative secondary metabolites from endophytic fungi. However, only very few studies also report the genetic basis for these phenotypic observations. With low sequencing cost and fast sample turnaround, genetics- and genomics-based approaches have risen to become comprehensive approaches to study natural products from a wide-range of organisms, especially to elucidate underlying biosynthetic pathways. However, in the field of fungal endophyte biology, elucidation of biosynthetic pathways is still a major challenge. As a relatively poorly investigated group of microorganisms, even in the light of recent efforts to sequence more fungal genomes, such as the 1000 Fungal Genomes Project at the Joint Genome Institute (JGI), the basis for bioprospecting of enzymes and pathways from endophytic fungi is still rather slim. In this review we want to discuss the current approaches and tools used to associate phenotype and genotype to elucidate biosynthetic pathways of secondary metabolites in endophytic fungi through the lens of bioprospecting. This review will point out the reported successes and shortcomings, and discuss future directions in sampling, and genetics and genomics of endophytic fungi. Identifying responsible biosynthetic genes for the numerous secondary metabolites isolated from endophytic fungi opens the opportunity to explore the genetic potential of producer strains to discover novel secondary metabolites and enhance secondary metabolite production by metabolic engineering resulting in novel and more affordable medicines and food additives.

A conserved Zn2Cys6 transcription factor, identified in a spontaneous mutant from in vitro passaging, is involved in pathogenicity of the blackleg fungus Leptosphaeria maculans

Continuous passaging in vitro can lead to the accumulation of changes in DNA sequence that potentially affect the properties of microbes, making them different from the original isolates. The identification of such genetic alterations is rare in fungi. A set of insertional mutants in the plant pathogenic fungus Leptosphaeria maculans, all derived from the same transformation experiment, had independent Agrobacterium T-DNA insertions and reduced pathogenicity on canola (Brassica napus). None of the insertions co-segregated in progeny from crosses with the reduction in pathogenicity. Genome sequences of three strains were analysed, and a mutation identified in a gene (ptf1, for pathogenicity-associated transcription factor 1) encoding a putative Zn₂(II)Cys₆ transcription factor. Homologs are found in other ascomycetes, and are required for pathogenicity by Fusarium graminearum, Fusarium oxysporum and Magnaporthe oryzae. The mutation in the L. maculans ptf1 gene co-segregates in progeny from crosses with the reduction in pathogenicity, a strain with an independent mutant allele isolated using CRISPR-Cas9 editing has reduced pathogenicity, and addition of wild type copies of the gene restores pathogenicity. Thus, this work defines a base pair substitution that occurred during in vitro passaging of a fungus that contributed to an attenuation of pathogenicity.

Population Genomics of Fungal Plant Pathogens and the Analyses of Rapidly Evolving Genome Compartments

Genome sequencing of fungal pathogens have documented extensive variation in genome structure and composition between species and in many cases between individuals of the same species. This type of genomic variation can be adaptive for pathogens to rapidly evolve new virulence phenotypes. Analyses of genome-wide variation in fungal pathogen genomes rely on high quality assemblies and methods to detect and quantify structural variation. Population genomic studies in fungi have addressed the underlying mechanisms whereby structural variation can be rapidly generated. Transposable elements, high mutation and recombination rates as well as incorrect chromosome segregation during mitosis and meiosis contribute to extensive variation observed in many species. We here summarize key findings in the field of fungal pathogen genomics and we discuss methods to detect and characterize structural variants including an alignment-based pipeline to study variation in population genomic data.

An Overview of Genomics, Phylogenomics and Proteomics Approaches in Ascomycota

Fungi are among the most successful eukaryotes on Earth: they have evolved strategies to survive in the most diverse environments and stressful conditions and have been selected and exploited for multiple aims by humans. The characteristic features intrinsic of Fungi have required evolutionary changes and adaptations at deep molecular levels. Omics approaches, nowadays including genomics, metagenomics, phylogenomics, transcriptomics, metabolomics, and proteomics have enormously advanced the way to understand fungal diversity at diverse taxonomic levels, under changeable conditions and in still under-investigated environments. These approaches can be applied both on environmental communities and on individual organisms, either in nature or in axenic culture and have led the traditional morphology-based fungal systematic to increasingly implement molecular-based approaches. The advent of next-generation sequencing technologies was key to boost advances in fungal genomics and proteomics research. Much effort has also been directed towards the development of methodologies for optimal genomic DNA and protein extraction and separation. To date, the amount of proteomics investigations in Ascomycetes exceeds those carried out in any other fungal group. This is primarily due to the preponderance of their involvement in plant and animal diseases and multiple industrial applications, and therefore the need to understand the biological basis of the infectious process to develop mechanisms for biologic control, as well as to detect key proteins with roles in stress survival. Here we chose to present an overview as much comprehensive as possible of the major advances, mainly of the past decade, in the fields of genomics (including phylogenomics) and proteomics of Ascomycota, focusing particularly on those reporting on opportunistic pathogenic, extremophilic, polyextremotolerant and lichenized fungi. We also present a review of the mostly used genome sequencing technologies and methods for DNA sequence and protein analyses applied so far for fungi.

Chromosome-level de novo assembly of Coprinopsis cinerea A43mut B43mut pab1-1 #326 and genetic variant identification of mutants using Nanopore MinION sequencing

The homokaryotic Coprinopsis cinerea strain A43mut B43mut pab1-1 #326 is a widely used experimental model for developmental studies in mushroom-forming fungi. It can grow on defined artificial media and complete the whole lifecycle within two weeks. The mutations in mating type factors A and B result in the special feature of clamp formation and fruiting without mating. This feature allows investigations and manipulations with a homokaryotic genetic background. Current genome assembly of strain #326 was based on short-read sequencing data and was highly fragmented, leading to the bias in gene annotation and downstream analyses. Here, we report a chromosome-level genome assembly of strain #326. Oxford Nanopore Technology (ONT) MinION sequencing was used to get long reads. Illumina short reads was used to polish the sequences. A combined assembly yield 13 chromosomes and a mitochondrial genome as individual scaffolds. The assembly has 15,250 annotated genes with a high synteny with the C. cinerea strain Okayama-7 #130. This assembly has great improvement on contiguity and annotations. It is a suitable reference for further genomic studies, especially for the genetic, genomic and transcriptomic analyses in ONT long reads. Single nucleotide variants and structural variants in six mutagenized and cisplatin-screened mutants could be identified and validated. A 66 bp deletion in Ras GTPase-activating protein (RasGAP) was found in all mutants. To make a better use of ONT sequencing platform, we modified a high-molecular-weight genomic DNA isolation protocol based on magnetic beads for filamentous fungi. This study showed the use of MinION to construct a fungal reference genome and to perform downstream studies in an individual laboratory. An experimental workflow was proposed, from DNA isolation and whole genome sequencing, to genome assembly and variant calling. Our results provided solutions and parameters for fungal genomic analysis on MinION sequencing platform.

Impact of a resistance gene against a fungal pathogen on the plant host residue microbiome: The case of the Leptosphaeria maculans-Brassica napus pathosystem

Oilseed rape residues are a crucial determinant of stem canker epidemiology as they support the sexual reproduction of the fungal pathogen Leptosphaeria maculans. The aim of this study was to characterize the impact of a resistance gene against L. maculans infection on residue microbial communities and to identify microorganisms interacting with this pathogen during residue degradation. We used near-isogenic lines to obtain healthy and infected host plants. The microbiome associated with the two types of plant residues was characterized by metabarcoding. A combination of linear discriminant analysis and ecological network analysis was used to compare the microbial communities and to identify microorganisms interacting with L. maculans. Fungal community structure differed between the two lines at harvest, but not subsequently, suggesting that the presence/absence of the resistance gene influences the microbiome at the base of the stem whilst the plant is alive, but that this does not necessarily lead to differential colonization of the residues by fungi. Direct interactions with other members of the community involved many fungal and bacterial amplicon sequence variants (ASVs). L. maculans appeared to play a minor role in networks, whereas one ASV affiliated to Plenodomus biglobosus (synonym Leptosphaeria biglobosa) from the Leptosphaeria species complex may be considered a keystone taxon in the networks at harvest. This approach could be used to identify and promote microorganisms with beneficial effects against residue-borne pathogens and, more broadly, to decipher the complex interactions between multispecies pathosystems and other microbial components in crop residues.

Comparative Genomic Analysis of Dactylonectria torresensis Strains from Grapevine, Soil and Weed Highlights Potential Mechanisms in Pathogenicity and Endophytic Lifestyle

The soil-borne fungus Dactylonectria torresensis is the most common causal agent of black-foot disease in Europe. However, there is a lack of understanding on how this fungus can provoke plant symptoms. In this study, we sequenced, annotated and analyzed the genomes of three isolates of D. torresensis collected from asymptomatic vine, weed and soil. Sequenced genomes were further compared to those of 27 fungal species including root and aerial pathogens, white rot degraders, indoor biodeterioration agents, saprotrophs, dark septate endophytes and mycorrhiza. Strains of D. torresensis present genomes with between 64 and 65 Mbp and with up to 18,548 predicted genes for each strain. Average Nucleotide Identity (ANI) shows that strains are different according to genome contents. Clusters of orthologous groups were compared, and clusters of genes related to necroses were particularly detected in all strains of D. torresensis (necrosis inducing peptides and proteins, and ethylene inducing peptides) as well as several genes involved in resistance against fungicides frequently used in viticulture such as copper. Interestingly, an expanded high number of genes related to carbohydrate-active enzymes were detected in each Dactylonectria strain, especially those related to glycoside hydrolases that could be involved in penetration of plant tissues or pathogenicity. An increased number of candidate genes for CAZyme classes AA9 and AA3-1 supports the ability of strains to efficiently degrade plant material. High numbers of genes of D. torresensis related to secretome and small secreted proteins were further characterized. Moreover, the presence of several gene clusters such as fujikurin-like genes was detected and were normally found in Fusariumfujikuroi, that have been linked to fungal pathogenicity. The phenotypes of the three strains investigated showed further difference in light response. We found that Dactylonectria strains have an increased number of photoreceptor encoding genes and we showed sequence alterations. Altogether, the results highlight several gene clusters present in D. torresensis strains that could be linked to endophytic lifestyle, pathogenicity, plant maceration and degradation of plant tissues as well as adaptation to soil contaminated with metals and metalloids and light response.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

Auxin biosynthesis in the phytopathogenic fungus Leptosphaeria maculans is associated with enhanced transcription of indole-3-pyruvate decarboxylase LmIPDC2 and tryptophan aminotransferase LmTAM1

Auxins are hormones that regulate growth and development in plants. Besides plants, various microorganisms also produce auxins. Here we investigate whether and how the phytopathogenic fungus Leptosphaeria maculans biosynthesizes auxins. We characterized the auxin profile of in vitro grown L. maculans. The culture was further supplied with the auxin biosynthetic-precursors tryptophan and tryptamine and gene expression and phytohormone content was analyzed. L. maculans in vitro produced IAA (indole-3-acetic acid) as the predominant auxin metabolite. IAA production could be further stimulated by supplying precursors. Expression of indole-3-pyruvate decarboxylase LmIPDC2, tryptophan aminotransferase LmTAM1 and nitrilase LmNIT1 genes was mainly upregulated after adding tryptophan and correlated with IAA production, suggesting that these genes are the key components of auxin biosynthesis in L. maculans. Tryptamine acted as a potent inducer of IAA production, though a pathway independent of LmIPDC2/LmTAM1 may be involved. Despite L. maculans being a rich source of bioactive IAA, the auxin metabolic profile of host plant Brassica napus was not altered upon infection. Exogenous IAA inhibited the growth of L. maculans in vitro when supplied in high concentration. Altogether, we showed that L. maculans is capable of IAA production and we have identified biosynthetic genes that were responsive to tryptophan treatment.

Complete, high-quality genomes from long-read metagenomic sequencing of two wolf lichen thalli reveals enigmatic genome architecture

Fungal genomes display incredible levels of complexity and diversity, and are exceptional study systems for genome evolution. Here we used the Oxford Nanopore MinION sequencing platform to generate high-quality fungal genomes from complex metagenomic samples of lichen thalli. We sequenced two wolf lichens using one flow cell per sample, generating 17.1 Gbps for Letharia lupina and 14.3 Gbps for Letharia columbiana. The resulting L. lupina genome is one of the most contiguous lichen genomes available to date, with 49.2 Mbp contained on 31 contigs. The L. columbiana genome, while less contiguous, is still relatively high quality, with 52.3 Mbp on a total of 161 contigs. Each thallus for both species contained multiple distinct haplotypes, a phenomenon that has rarely been empirically demonstrated. The Oxford Nanopore sequencing technologies are robust and effective when applied to complex symbioses, and have the potential to fundamentally transform our understanding of fungal genetics.

Mycovirus-Induced Hypervirulence of Leptosphaeria biglobosa Enhances Systemic Acquired Resistance to Leptosphaeria maculans in Brassica napus

Phoma stem canker (blackleg) is one of the most important diseases of winter oilseed rape (Brassica napus) worldwide and is caused by a complex that comprises at least two species: Leptosphaeria maculans and L. biglobosa. Screening a panel of field Leptosphaeria isolates from B. napus for the presence of mycoviruses revealed the presence of a novel double-stranded RNA quadrivirus in L. biglobosa and no viruses in L. maculans. Following elimination of the mycovirus, virus-infected and virus-free isogenic lines of L. biglobosa were created. A direct comparison of the growth and virulence of these isogenic lines illustrated that virus infection caused hypervirulence and resulted in induced systemic resistance toward L. maculans in B. napus following lower leaf preinoculation with the virus-infected isolate. Analysis of the plant transcriptome suggests that the presence of the virus leads to subtle alterations in metabolism and plant defenses. For instance, transcripts involved in carbohydrate and amino acid metabolism are enriched in plants treated with the virus-infected isolate, while pathogenesis-related proteins, chitinases and WRKY transcription factors are differentially expressed. These results illustrate the potential for deliberate inoculation of plants with hypervirulent L. biglobosa to decrease the severity of Phoma stem canker later in the growing season.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

A rapid approach to profiling diverse fungal communities using the MinION™ nanopore sequencer

The Oxford Nanopore Technologies MinION™ sequencer holds the capability to generate long amplicon reads; however, only a small amount of information is available regarding methodological approaches and the ability to identify a broad diversity of fungal taxa. To assess capabilities, three fungal mock communities were sequenced, each of which had varying ratios of 16 taxa. The data were processed through our selected pipeline. The MinION recovered all mock community members, when mixed at equal ratios. When a taxon was represented at a lower ratio, it was not recovered or decreased in relative abundance. Despite high error rates, highly accurate consensus sequences can be derived. This methodological approach identified all mock community taxa, demonstrating the MinION can be used as a practical alternative to profile fungal communities.

RaGOO: fast and accurate reference-guided scaffolding of draft genomes

We present RaGOO, a reference-guided contig ordering and orienting tool that leverages the speed and sensitivity of Minimap2 to accurately achieve chromosome-scale assemblies in minutes. After the pseudomolecules are constructed, RaGOO identifies structural variants, including those spanning sequencing gaps. We show that RaGOO accurately orders and orients 3 de novo tomato genome assemblies, including the widely used M82 reference cultivar. We then demonstrate the scalability and utility of RaGOO with a pan-genome analysis of 103 Arabidopsis thaliana accessions by examining the structural variants detected in the newly assembled pseudomolecules. RaGOO is available open source at https://github.com/malonge/RaGOO .