Search

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

The wheat pathogen Parastagonospora nodorum has emerged as a model necrotrophic fungal species with growing genomic resources. Recent population-level pan-genome studies were leveraged to provide novel insights into pathogen evolution and effector-like gene contents relevant to local crop disease outbreaks. In this study, we examined 156 isolates representing a regional population from the Western Australian (WA) wheat-belt region, and 17 internationally sourced isolates. We observed a highly diverse local population, within which were numerous small and highly similar clusters of isolates from hotter and drier regions. Pan-genome assembly and orthologous gene datasets resulted in 3579 predicted effector candidates, 2291 of which exhibited presence-absence variation (PAV) across the population, and 1362 were specific to WA isolates. There was an abundance of mutations (including repeat-induced point mutation (RIP)), distributed in ‘hot-spots’ within the pan-genomic landscape that were rich in effector candidates. Three characterised effector loci (ToxA, Tox1 and Tox3) were located within sub- telomeric regions of lower diversity, but were nestled within larger high-diversity regions. RIP was widespread across the genome, but non-synonymous RIP-like mutations were strongly selected against. These improved bioinformatic resources for P. nodorum, represent progressive advancements in fungal pan-genomics, with a view towards supporting region- specific surveillance of host-pathogen interactions.

Phytophthora cinnamomi is a pathogenic oomycete that causes plant dieback disease across a range of natural ecosystems and in many agriculturally important crops on a global scale. An annotated draft genome sequence and annotation is publicly available (JGI Mycocosm) and suggests 26,131 gene models. In this study, soluble mycelial, extracellular (secretome) and zoospore proteins of P. cinnamomi were exploited to refine the genome by correcting gene annotations and discovering novel genes. By implementing the diverse set of sub-proteomes into a generated proteogenomics pipeline, we were able to improve the P. cinnamomi genome. Liquid chromatography mass spectrometry was used to obtain high confidence peptides with spectral matching to both the annotated genome and a generated 6-frame translation. 2,764 annotations from the draft genome were confirmed by spectral matching. Using a proteogenomic pipeline, mass spectra were used to edit the P. cinnamomi genome and allowed identification of 23 new gene models and 60 edited gene features using high confidence peptides obtained by mass spectrometry, suggesting a rate of incorrect annotations of 3% of the detectable proteome. The novel features were further validated by total peptide support, alongside functional analysis including the use of Gene Ontology and functional domain identification. We demonstrated the use of spectral data in combination with our proteogenomics pipeline can be used to improve the genome of important plant diseases and identify biologically relevant missed genes. This study presents the first use of spectral data to edit and manually annotate an oomycete pathogen.

Background The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3. Results A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions. Conclusions We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

We report a fungal pan-genome study involving Parastagonospora spp., including 21 isolates of the wheat (Triticum aestivum) pathogen Parastagonospora nodorum, 10 of the grass-infecting Parastagonospora avenae, and 2 of a closely related undefined sister species. We observed substantial variation in the distribution of polymorphisms across the pan-genome, including repeat-induced point mutations, diversifying selection and gene gains and losses. We also discovered chromosome-scale inter and intraspecific presence/absence variation of some sequences, suggesting the occurrence of one or more accessory chromosomes or regions that may play a role in host-pathogen interactions. The presence of known pathogenicity effector loci SnToxA, SnTox1, and SnTox3 varied substantially among isolates. Three P. nodorum isolates lacked functional versions for all three loci, whereas three P. avenae isolates carried one or both of the SnTox1 and SnTox3 genes, indicating previously unrecognized potential for discovering additional effectors in the P. nodorum-wheat pathosystem. We utilized the pan-genomic comparative analysis to improve the prediction of pathogenicity effector candidates, recovering the three confirmed effectors among our top-ranked candidates. We propose applying this pan-genomic approach to identify the effector repertoire involved in other host-microbe interactions involving necrotrophic pathogens in the Pezizomycotina.

The traditional classification of fungal and oomycete phytopathogens into three classes – biotrophs, hemibiotrophs, or necrotrophs – is unsustainable. This study highlights multiple phytopathogen species for which these labels have been inappropriately applied. We propose a novel and reproducible classification based solely on genome-derived analysis of carbohydrate-active enzyme (CAZyme) gene content called CAZyme-Assisted Training And Sorting of -trophy (CATAStrophy). CATAStrophy defines four major divisions for species associated with living plants. These are monomertrophs (Mo) (corresponding to biotrophs), polymertrophs (P) (corresponding to necrotrophs), mesotrophs (Me) (corresponding to hemibiotrophs), and vasculartrophs (including species commonly described as wilts, rots, or anthracnoses). The Mo class encompasses symbiont, haustorial, and non-haustorial species. Me are divided into the subclasses intracellular and extracellular Me, and the P into broad and narrow host sub-classes. This gives a total of seven discrete plant-pathogenic classes. The classification provides insight into the properties of these species and offers a facile route to develop control measures for newly recognized diseases. Software for CATAStrophy is available online at https://github.com/ccdmb/catastrophy. We present the CATAStrophy method for the prediction of trophic phenotypes based on CAZyme gene content, as a complementary method to the traditional tripartite “biotroph–hemibiotroph–necrotroph” classifications that may encourage renewed investigation and revision within the fungal biology community.

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as “OcculterCut” (https://sourceforge.net/projects/occultercut, last accessed April 30, 2016), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighboring repeat-rich regions. Over time RIP perpetuates “two speeds” of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions.

Zymoseptoria tritici is a host-specific, necrotrophic pathogen of wheat. Infection by Z. tritici is characterized by its extended latent period, which typically lasts 2 wks, and is followed by extensive host cell death, and rapid proliferation of fungal biomass. This work characterizes the level of genomic variation in 13 isolates, for which we have measured virulence on 11 wheat cultivars with differential resistance genes. Between the reference isolate, IPO323, and the 13 Australian isolates we identified over 800,000 single nucleotide polymorphisms, of which ∼10% had an effect on the coding regions of the genome. Furthermore, we identified over 1700 probable presence/absence polymorphisms in genes across the Australian isolates using de novo assembly. Finally, we developed a gene tree sorting method that quickly identifies groups of isolates within a single gene alignment whose sequence haplotypes correspond with virulence scores on a single wheat cultivar. Using this method, we have identified

Background Pyrenophora tritici-repentis (Ptr) is a necrotrophic fungal pathogen that causes the major wheat disease, tan spot. We set out to provide essential genomics-based resources in order to better understand the pathogenicity mechanisms of this important pathogen. Results Here, we present eight new Ptr isolate genomes, assembled and annotated; representing races 1, 2 and 5, and a new race. We report a high quality Ptr reference genome, sequenced by PacBio technology with Illumina paired-end data support and optical mapping. An estimated 98% of the genome coverage was mapped to 10 chromosomal groups, using a two-enzyme hybrid approach. The final reference genome was 40.9 Mb and contained a total of 13,797 annotated genes, supported by transcriptomic and proteogenomics data sets. Conclusions Whole genome comparative analysis revealed major chromosomal segmental rearrangements and fusions, highlighting intraspecific genome plasticity in this species. Furthermore, the Ptr race classification was not supported at the whole genome level, as phylogenetic analysis did not cluster the ToxA producing isolates. This expansion of available Ptr genomics resources will directly facilitate research aimed at controlling tan spot disease. Electronic supplementary material The online version of this article (10.1186/s12864-018-4680-3) contains supplementary material, which is available to authorized users.

Fungal pathogen genomes can often be divided into core and accessory regions. Accessory regions ARs) may be comprised of either ARs (within core chromosomes (CCs) or wholly dispensable (accessory) chromosomes (ACs). Fungal ACs and ARs typically accumulate mutations and structural rearrangements more rapidly over time than CCs and many harbor genes relevant to host-pathogen interactions. These regions are of particular interest in plant pathology and include host-specific virulence factors and secondary metabolite synthesis gene clusters. This review outlines known ACs and ARs in fungal genomes, methods used for their detection, their common properties that differentiate them from the core genome, and what is currently known of their various roles in pathogenicity. Reports on the evolutionary processes generating and shaping AC and AR compartments are discussed, including repeat induced point mutation and breakage fusion bridge cycles. Previously ACs have been studied extensively within key genera, including Fusarium, Zymoseptoria, and Alternaria, but are growing in frequency of observation and perceived importance across a wider range of fungal species. Recent advances in sequencing technologies permit affordable genome assembly and resequencing of populations that will facilitate further discovery and routine screening of ACs.

Zymoseptoria tritici (syn. Mycosphaerella graminicola, Septoria tritici) is a haploid fungus belonging to the class Dothideomycetes. It is the causal agent of septoria leaf blotch – one of the world’s most significant diseases of wheat. Here we review the genomic and bioinformatic resources that have been generated for Z. tritici. These include the whole-genome reference assembly for isolate IPO323, genome resequencing of alternate isolates, mitochondrial genome sequences, transcriptome sequences and expression data, and annotations of gene structure and function. We also highlight important advances in our fundamental knowledge of genome evolution and its effects on adaptation and pathogenicity in Z. tritici that have been facilitated by these resources.

Stagonospora nodorum is an important wheat (Triticum aestivum) pathogen in many parts of the world, causing major yield losses. It was the first species in the large fungal Dothideomycete class to be genome sequenced. The reference genome sequence (SN15) has been instrumental in the discovery of genes encoding necrotrophic effectors that induce disease symptoms in specific host genotypes. Here we present the genome sequence of two further S. nodorum strains (Sn4 and Sn79) that differ in their effector repertoire from the reference. Sn79 is avirulent on wheat and produces no apparent effectors when infiltrated onto many cultivars and mapping population parents. Sn4 is pathogenic on wheat and has virulences not found in SN15. The new strains, sequenced with short-read Illumina chemistry, are compared with SN15 by a combination of mapping and de novo assembly approaches. Each of the genomes contains a large number of strain-specific genes, many of which have no meaningful similarity to any known gene. Large contiguous sections of the reference genome are absent in the two newly sequenced strains. We refer to these differences as “sectional gene absences.” The presence of genes in pathogenic strains and absence in Sn79 is added to computationally predicted properties of known proteins to produce a list of likely effector candidates. Transposon insertion was observed in the mitochondrial genomes of virulent strains where the avirulent strain retained the likely ancestral sequence. The study suggests that short-read enabled comparative genomics is an effective way to both identify new S. nodorum effector candidates and to illuminate evolutionary processes in this species.

The wheat pathogen Stagonospora nodorum, causal organism of the wheat disease Stagonospora nodorum blotch, has emerged as a model for the Dothideomycetes, a large fungal taxon that includes many important plant pathogens. The initial annotation of the genome assembly included 16 586 nuclear gene models. These gene models were used to design a microarray that has been interrogated with labelled transcripts from six cDNA samples: four from infected wheat plants at time points spanning early infection to sporulation, and two time points taken from growth in artificial media. Positive signals of expression were obtained for 12 281 genes. This represents strong corroborative evidence of the validity of these gene models. Significantly differential expression between the various time points was observed. When infected samples were compared with axenic cultures, 2882 genes were expressed at a higher level in planta and 3630 were expressed more highly in vitro. Similar numbers were differentially expressed between different developmental stages. The earliest time points in planta were particularly enriched in differentially expressed genes. A disproportionate number of the early expressed gene products were predicted to be secreted, but otherwise had no obvious sequence homology to functionally characterized genes. These genes are candidate necrotrophic effectors. We have focused attention on genes for carbohydrate metabolism and the specific biosynthetic pathways active during growth in planta. The analysis points to a very dynamic adjustment of metabolism during infection. Functional analysis of a gene in the coenzyme A biosynthetic pathway showed that the enzyme was dispensable for growth, indicating that a precursor is supplied by the plant.

Here, we report on the use of quantitative PCR (qPCR) to determine gene copy number in filamentous fungi. Using the sequenced dothideomycete Stagonospora nodorum, qPCR was used to unequivocally confirm the presence of single, two and three copy regions as predicted by in silico PCR. Further validation of the technique was demonstrated by verifying the copy numbers of introduced gene cassettes in previously characterised transformants of S. nodorum. Apart from increased sensitivity, this technique offers a high-throughput alternative to Southern blots for determining gene copy number, a

Rhizoctonia solani is a nectrotrophic fungal pathogen that causes billions of dollars of damage to agriculture worldwide and infects a broad host range including wheat, rice, potato and legumes. In this study we identify wheat genes that are differentially expressed in response to the R. solani isolate, AG8, using microarray technology. A significant number of wheat genes identified in this screen were involved in reactive oxygen species (ROS) production and redox regulation. Levels of ROS species were increased in wheat root tissue following R. solani infection as determined by Nitro Blue Tetrazolium (NBT), 3,3'-diaminobenzidine (DAB) and titanium sulphate measurements. Pathogen/ROS related genes from R. solani were also tested for expression patterns upon wheat infection. TmpL, a R. solani gene homologous to a gene associated with ROS regulation in Alternaria brassicicola, and OAH, a R. solani gene homologous to oxaloacetate acetylhydrolase which has been shown to produce oxalic acid in Sclerotinia sclerotiorum, were highly induced in R. solani when infecting wheat. We speculate that the interplay between the wheat and R. solani ROS generating proteins may be important for determining the outcome of the wheat/R. solani interaction.

We report the first gene-based linkage map of Lupinus angustifolius (narrow-leafed lupin) and its comparison to the partially sequenced genome of Medicago truncatula. The map comprises 382 loci in 20 major linkage groups, two triplets, three pairs and 11 unlinked loci and is 1,846 cM in length. The map was generated from the segregation of 163 RFLP markers, 135 gene-based PCR markers, 75 AFLP and 4 AFLP-derived SCAR markers in a mapping population of 93 recombinant inbred lines, derived from a cross between domesticated and wild-type parents. This enabled the mapping of five major genes controlling key domestication traits in L. angustifolius. Using marker sequence data, the L. angustifolius genetic map was compared to the partially completed M. truncatula genome sequence. We found evidence of conserved synteny in some regions of the genome despite the wide evolutionary distance between these legume species. We also found new evidence of widespread duplication within the L. angustifolius genome.

The plant-pathogenic fungus Mycosphaerella graminicola (asexual stage: Septoria tritici) causes septoria tritici blotch, a disease that greatly reduces the yield and quality of wheat. This disease is economically important in most wheat-growing areas worldwide and threatens global food production. Control of the disease has been hampered by a limited understanding of the genetic and biochemical bases of pathogenicity, including mechanisms of infection and of resistance in the host. Unlike most other plant pathogens, M. graminicola has a long latent period during which it evades host defenses. Although this type of stealth pathogenicity occurs commonly in Mycosphaerella and other Dothideomycetes, the largest class of plant-pathogenic fungi, its genetic basis is not known. To address this problem, the genome of M. graminicola was sequenced completely. The finished genome contains 21 chromosomes, eight of which could be lost with no visible effect on the fungus and thus are dispensable. This eight-chromosome dispensome is dynamic in field and progeny isolates, is different from the core genome in gene and repeat content, and appears to have originated by ancient horizontal transfer from an unknown donor. Synteny plots of the M. graminicola chromosomes versus those of the only other sequenced Dothideomycete, Stagonospora nodorum, revealed conservation of gene content but not order or orientation, suggesting a high rate of intra-chromosomal rearrangement in one or both species. This observed “mesosynteny” is very different from synteny seen between other organisms. A surprising feature of the M. graminicola genome compared to other sequenced plant pathogens was that it contained very few genes for enzymes that break down plant cell walls, which was more similar to endophytes than to pathogens. The stealth pathogenesis of M. graminicola probably involves degradation of proteins rather than carbohydrates to evade host defenses during the biotrophic stage of infection and may have evolved from endophytic ancestors., PLoS Genetics, 7 (6), ISSN:1553-7390, ISSN:1553-7404

Chickpea (Cicer arietinum) is the second most widely grown legume crop after soybean, accounting for a substantial proportion of human dietary nitrogen intake and playing a crucial role in food security in developing countries. We report the ∼738-Mb draft whole genome shotgun sequence of CDC Frontier, a kabuli chickpea variety, which contains an estimated 28,269 genes. Resequencing and analysis of 90 cultivated and wild genotypes from ten countries identifies targets of both breeding-associated genetic sweeps and breeding-associated balancing selection. Candidate genes for disease resistance and agronomic traits are highlighted, including traits that distinguish the two main market classes of cultivated chickpea - desi and kabuli. These data comprise a resource for chickpea improvement through molecular breeding and provide insights into both genome diversity and domestication. Copyright © 2013 Nature America, Inc.

The class Dothideomycetes is a recently defined taxon within the phylum Ascomycota, with approximately 20 000 member species. Many species are important phytopathogens employing diverse pathogenicity strategies. This review compares the currently available genome sequences, including their mitochondrial sequences. The genomes of Dothideomycetes exhibit remarkable plasticity, characterised by extensive rearrangement of gene order and orientation. Dothideomycetes also display the ability to take advantage of repeat induced point mutation (RIP), a fungal specific mechanism to protect their genomes against transposon invasion. RIP also affects areas of the genome adjacent to transposon insertion sites conferring the adaptive advantage of rapid generation of diversity in fungal pathogenicity effector genes. In addition, Dothideomycetes have the propensity for lateral gene transfer which can involve the transmission of pathogenicity effectors. The combination of these mechanisms has led to the ascendance of many Dothideomycetes as the dominant fungal pathogens of their cultivated host plant species.

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints., Leptosphaeria maculans is a plant pathogen that causes stem canker of oilseed rape. Rouxel et al. sequence and describe the key features of the L. maculans genome, including partitioning into AT-rich blocks that are enriched in effector genes and transposable elements affected by repeat-induced point mutation.

Background Repeat-induced point mutation (RIP) is a fungal genome defence mechanism guarding against transposon invasion. RIP mutates the sequence of repeated DNA and over time renders the affected regions unrecognisable by similarity search tools such as BLAST. Results DeRIP is a new software tool developed to predict the original sequence of a RIP-mutated region prior to the occurrence of RIP. In this study, we apply deRIP to the genome of the wheat pathogen Stagonospora nodorum SN15 and predict the origin of several previously uncharacterised classes of repetitive DNA. Conclusions Five new classes of transposon repeats and four classes of endogenous gene repeats were identified after deRIP. The deRIP process is a new tool for fungal genomics that facilitates the identification and understanding of the role and origin of fungal repetitive DNA. DeRIP is open-source and is available as part of the RIPCAL suite at http://www.sourceforge.net/projects/ripcal.

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a ‘gene for gene’ manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans., PLoS Pathogens, 6 (11), ISSN:1553-7374, ISSN:1553-7366

Stagonospora nodorum is a major necrotrophic fungal pathogen of wheat (Triticum aestivum) and a member of the Dothideomycetes, a large fungal taxon that includes many important plant pathogens affecting all major crop plant families. Here, we report the acquisition and initial analysis of a draft genome sequence for this fungus. The assembly comprises 37,164,227 bp of nuclear DNA contained in 107 scaffolds. The circular mitochondrial genome comprises 49,761 bp encoding 46 genes, including four that are intron encoded. The nuclear genome assembly contains 26 classes of repetitive DNA, comprising 4.5% of the genome. Some of the repeats show evidence of repeat-induced point mutations consistent with a frequent sexual cycle. ESTs and gene prediction models support a minimum of 10,762 nuclear genes. Extensive orthology was found between the polyketide synthase family in S. nodorum and Cochliobolus heterostrophus, suggesting an ancient origin and conserved functions for these genes. A striking feature of the gene catalog was the large number of genes predicted to encode secreted proteins; the majority has no meaningful similarity to any other known genes. It is likely that genes for host-specific toxins, in addition to ToxA, will be found among this group. ESTs obtained from axenic mycelium grown on oleate (chosen to mimic early infection) and late-stage lesions sporulating on wheat leaves were obtained. Statistical analysis shows that transcripts encoding proteins involved in protein synthesis and in the production of extracellular proteases, cellulases, and xylanases predominate in the infection library. This suggests that the fungus is dependant on the degradation of wheat macromolecular constituents to provide the carbon skeletons and energy for the synthesis of proteins and other components destined for the developing pycnidiospores.

Background The development of genetic markers is complex and costly in species with little pre-existing genomic information. Faba bean possesses one of the largest and least studied genomes among cultivated crop plants and no gene-based genetic maps exist. Gene-based orthologous markers allow chromosomal regions and levels of synteny to be characterised between species, reveal phylogenetic relationships and chromosomal evolution, and enable targeted identification of markers for crop breeding. In this study orthologous codominant cross-species markers have been deployed to produce the first exclusively gene-based genetic linkage map of faba bean (Vicia faba), using an F6 population developed from a cross between the lines Vf6 (equina type) and Vf27 (paucijuga type). Results Of 796 intron-targeted amplified polymorphic (ITAP) markers screened, 151 markers could be used to construct a comparative genetic map. Linkage analysis revealed seven major and five small linkage groups (LGs), one pair and 12 unlinked markers. Each LG was comprised of three to 30 markers and varied in length from 23.6 cM to 324.8 cM. The map spanned a total length of 1685.8 cM. A simple and direct macrosyntenic relationship between faba bean and Medicago truncatula was evident, while faba bean and lentil shared a common rearrangement relative to M. truncatula. One hundred and four of the 127 mapped markers in the 12 LGs, which were previously assigned to M. truncatula genetic and physical maps, were found in regions syntenic between the faba bean and M. truncatula genomes. However chromosomal rearrangements were observed that could explain the difference in chromosome numbers between these three legume species. These rearrangements suggested high conservation of M. truncatula chromosomes 1, 5 and 8; moderate conservation of chromosomes 2, 3, 4 and 7 and no conservation with M. truncatula chromosome 6. Multiple PCR amplicons and comparative mapping were suggestive of small-scale duplication events in faba bean. This study also provides a preliminary indication for finer scale macrosynteny between M. truncatula, lentil and faba bean. Markers originally designed from genes on the same M. truncatula BACs were found to be grouped together in corresponding syntenic areas in lentil and faba bean. Conclusion Despite the large size of the faba bean genome, comparative mapping did not reveal evidence for polyploidisation, segmental duplication, or significant rearrangements compared to M. truncatula, although a bias in the use of single locus markers may have limited the detection of duplications. Non-coding repetitive DNA or transposable element content provides a possible explanation for the difference in genome sizes. Similar patterns of rearrangements in faba bean and lentil compared to M. truncatula support phylogenetic studies dividing these species into the tribes Viceae and Trifoliae. However, substantial macrosynteny was apparent between faba bean and M. truncatula, with the exception of chromosome 6 where no orthologous markers were found, confirming previous investigations suggesting chromosome 6 is atypical. The composite map, anchored with orthologous markers mapped in M. truncatula, provides a central reference map for future use of genomic and genetic information in faba bean genetic analysis and breeding.

The wheat pathogen Zymoseptoria tritici possesses a large number of accessory chromosomes that may be present or absent in its genome. The genome of the reference isolate IPO323 has been assembled to a very high standard and contains 21 full length chromosome sequences, 8 of which represent accessory chromosomes. The IPO323 reference, when combined with low-cost next-generation sequencing and bioinformatics, can be used as a powerful tool to assess the presence or absence of accessory chromosomes. We present an outline of a range of bioinformatics techniques that can be applied to the analysis of presence–absence variation among accessory chromosomes across 13 novel isolates of Z. tritici.

Background The major proteins in lupin seeds are conglutins that have primary roles in supplying carbon, sulphur and nitrogen and energy for the germinating seedling. They fall into four families; α, β, γ and δ. Interest in these conglutins is growing as family members have been shown to have beneficial nutritional and pharmaceutical properties. Results An in-depth transcriptome and draft genome from the narrow-leafed lupin (NLL; Lupinus angustifolius) variety, Tanjil, were examined and 16 conglutin genes were identified. Using RNAseq data sets, the structure and expression of these 16 conglutin genes were analysed across eight lupin varieties from five lupin species. Phylogenic analysis suggest that the α and γ conglutins diverged prior to lupin speciation while β and δ members diverged both prior and after speciation. A comparison of the expression of the 16 conglutin genes was performed, and in general the conglutin genes showed similar levels of RNA expression among varieties within species, but quite distinct expression patterns between lupin species. Antibodies were generated against the specific conglutin families and immunoblot analyses were used to compare the levels of conglutin proteins in various tissues and during different stages of seed development in NLL, Tanjil, confirming the expression in the seed. This analysis showed that the conglutins were expressed highly at the mature seed stage, in all lupin species, and a range of polypeptide sizes were observed for each conglutin family. Conclusions This study has provided substantial information on the complexity of the four conglutin families in a range of lupin species in terms of their gene structure, phylogenetic relationships as well as their relative RNA and protein abundance during seed development. The results demonstrate that the majority of the heterogeneity of conglutin polypeptides is likely to arise from post-translational modification from a limited number of precursor polypeptides rather than a large number of different genes. Overall, the results demonstrate a high degree of plasticity for conglutin expression during seed development in different lupin species. Electronic supplementary material The online version of this article (doi:10.1186/s12870-015-0485-6) contains supplementary material, which is available to authorized users.