Understanding the Effectiveness of Genomic Prediction in Tetraploid Potato.

Use of genomic prediction (GP) in tetraploid is becoming more common. Therefore, we think it is the right time for a comparison of GP models for tetraploid potato. GP models were compared that contrasted shrinkage with variable selection, parametric vs. non-parametric models and different ways of accounting for non-additive genetic effects. As a complement to GP, association studies were carried out in an attempt to understand the differences in prediction accuracy. We compared our GP models on a data set consisting of 147 cultivars, representing worldwide diversity, with over 39 k GBS markers and measurements on four tuber traits collected in six trials at three locations during 2 years. GP accuracies ranged from 0.32 for tuber count to 0.77 for dry matter content. For all traits, differences between GP models that utilised shrinkage penalties and those that performed variable selection were negligible. This was surprising for dry matter, as only a few additive markers explained over 50% of phenotypic variation. Accuracy for tuber count increased from 0.35 to 0.41, when dominance was included in the model. This result is supported by Genome Wide Association Study (GWAS) that found additive and dominance effects accounted for 37% of phenotypic variation, while significant additive effects alone accounted for 14%. For tuber weight, the Reproducing Kernel Hilbert Space (RKHS) model gave a larger improvement in prediction accuracy than explicitly modelling epistatic effects. This is an indication that capturing the between locus epistatic effects of tuber weight can be done more effectively using the semi-parametric RKHS model. Our results show good opportunities for GP in 4x potato.

Neofunctionalisation of the Sli gene leads to self-compatibility and facilitates precision breeding in potato.

Genetic gain in potato is hampered by the heterozygous tetraploid genome of cultivated potato. Converting potato into a diploid inbred-line based F1-hybrid crop provides a promising route towards increased genetic gain. The introduction of a dominant S-locus inhibitor (Sli) gene into diploid potato germplasm allows efficient generation of self-fertilized seeds and thus the development of potato inbred lines. Little is known about the structure and function of the Sli locus. Here we describe the mapping of Sli to a 12.6 kb interval on chromosome 12 using a recombinant screen approach. One of two candidate genes present in this interval shows a unique sequence that is exclusively present in self-compatible lines. We describe an expression vector that converts self-incompatible genotypes into self-compatible and a CRISPR-Cas9 vector that converts SC genotypes into SI. The Sli gene encodes an F-box protein that is specifically expressed in pollen from self-compatible plants. A 533 bp insertion in the promotor of that gene leads to a gain of function mutation, which overcomes self-pollen rejection.

Solyntus, the New Highly Contiguous Reference Genome for Potato (Solanum tuberosum).

With the rapid expansion of the application of genomics and sequencing in plant breeding, there is a constant drive for better reference genomes. In potato (Solanum tuberosum), the third largest food crop in the world, the related species S. phureja, designated "DM", has been used as the most popular reference genome for the last 10 years. Here, we introduce the de novo sequenced genome of Solyntus as the next standard reference in potato genome studies. A true Solanum tuberosum made up of 116 contigs that is also highly homozygous, diploid, vigorous and self-compatible, Solyntus provides a more direct and contiguous reference then ever before available. It was constructed by sequencing with state-of-the-art long and short read technology and assembled with Canu. The 116 contigs were assembled into scaffolds to form each pseudochromosome, with three contigs to 17 contigs per chromosome. This assembly contains 93.7% of the single-copy gene orthologs from the Solanaceae set and has an N50 of 63.7 Mbp. The genome and related files can be found at https://www.plantbreeding.wur.nl/Solyntus/ With the release of this research line and its draft genome we anticipate many exciting developments in (diploid) potato research.

Specific Hypersensitive Response-Associated Recognition of New Apoplastic Effectors from Cladosporium fulvum in Wild Tomato.

Tomato leaf mold disease is caused by the biotrophic fungus Cladosporium fulvum. During infection, C. fulvum produces extracellular small secreted protein (SSP) effectors that function to promote colonization of the leaf apoplast. Resistance to the disease is governed by Cf immune receptor genes that encode receptor-like proteins (RLPs). These RLPs recognize specific SSP effectors to initiate a hypersensitive response (HR) that renders the pathogen avirulent. C. fulvum strains capable of overcoming one or more of all cloned Cf genes have now emerged. To combat these strains, new Cf genes are required. An effectoromics approach was employed to identify wild tomato accessions carrying new Cf genes. Proteomics and transcriptome sequencing were first used to identify 70 apoplastic in planta-induced C. fulvum SSPs. Based on sequence homology, 61 of these SSPs were novel or lacked known functional domains. Seven, however, had predicted structural homology to antimicrobial proteins, suggesting a possible role in mediating antagonistic microbe-microbe interactions in planta. Wild tomato accessions were then screened for HR-associated recognition of 41 SSPs, using the Potato virus X-based transient expression system. Nine SSPs were recognized by one or more accessions, suggesting that these plants carry new Cf genes available for incorporation into cultivated tomato.

Novel Introner-Like Elements in fungi Are Involved in Parallel Gains of Spliceosomal Introns.

Spliceosomal introns are key components of the eukaryotic gene structure. Although they contributed to the emergence of eukaryotes, their origin remains elusive. In fungi, they might originate from the multiplication of invasive introns named Introner-Like Elements (ILEs). However, so far ILEs have been observed in six fungal species only, including Fulvia fulva and Dothistroma septosporum (Dothideomycetes), arguing against ILE insertion as a general mechanism for intron gain. Here, we identified novel ILEs in eight additional fungal species that are phylogenetically related to F. fulva and D. septosporum using PCR amplification with primers derived from previously identified ILEs. The ILE content appeared unique to each species, suggesting independent multiplication events. Interestingly, we identified four genes each containing two gained ILEs. By analysing intron positions in orthologues of these four genes in Ascomycota, we found that three ILEs had inserted within a 15 bp window that contains regular spliceosomal introns in other fungal species. These three positions are not the result of intron sliding because ILEs are newly gained introns. Furthermore, the alternative hypothesis of an inferred ancestral gain followed by independent losses contradicts the observed degeneration of ILEs. These observations clearly indicate three parallel intron gains in four genes that were randomly identified. Our findings suggest that parallel intron gain is a phenomenon that has been highly underestimated in ILE-containing fungi, and likely in the whole fungal kingdom.

Transcriptome sequencing uncovers the Avr5 avirulence gene of the tomato leaf mold pathogen Cladosporium fulvum.

The Cf-5 gene of tomato confers resistance to strains of the fungal pathogen Cladosporium fulvum carrying the avirulence gene Avr5. Although Cf-5 has been cloned, Avr5 has remained elusive. We report the cloning of Avr5 using a combined bioinformatic and transcriptome sequencing approach. RNA-Seq was performed on the sequenced race 0 strain (0WU; carrying Avr5), as well as a race 5 strain (IPO 1979; lacking a functional Avr5 gene) during infection of susceptible tomato. Forty-four in planta-induced C. fulvum candidate effector (CfCE) genes of 0WU were identified that putatively encode a secreted, small cysteine-rich protein. An expressed transcript sequence comparison between strains revealed two polymorphic CfCE genes in IPO 1979. One of these conferred avirulence to IPO 1979 on Cf-5 tomato following complementation with the corresponding 0WU allele, confirming identification of Avr5. Complementation also led to increased fungal biomass during infection of susceptible tomato, signifying a role for Avr5 in virulence. Seven of eight race 5 strains investigated escape Cf-5-mediated resistance through deletion of the Avr5 gene. Avr5 is heavily flanked by repetitive elements, suggesting that repeat instability, in combination with Cf-5-mediated selection pressure, has led to the emergence of race 5 strains deleted for the Avr5 gene.

Functional analysis of the conserved transcriptional regulator CfWor1 in Cladosporium fulvum reveals diverse roles in the virulence of plant pathogenic fungi.

Fungal Wor1-like proteins are conserved transcriptional regulators that are reported to regulate the virulence of several plant pathogenic fungi by affecting the expression of virulence genes. Here, we report the functional analysis of CfWor1, the homologue of Wor1 in Cladosporium fulvum. Δcfwor1 mutants produce sclerotium-like structures and rough hyphae, which are covered with a black extracellular matrix. These mutants do not sporulate and are no longer virulent on tomato. A CE.CfWor1 transformant that constitutively expresses CfWor1 produces fewer spores with altered morphology and is also reduced in virulence. RNA-seq and RT-qrtPCR analyses suggest that reduced virulence of Δcfwor1 mutants is due to global downregulation of transcription, translation and mitochondrial respiratory chain. The reduced virulence of the CE.CfWor1 transformant is likely due to downregulation of effector genes. Complementation of a non-virulent Δfosge1 (Wor1-homologue) mutant of Fusarium oxysporum f. sp. lycopersici with CfWor1 restored expression of the SIX effector genes in this fungus, but not its virulence. Chimeric proteins of CfWor1/FoSge1 also only partially restored defects of the Δfosge1 mutant, suggesting that these transcriptional regulators have functionally diverged. Altogether, our results suggest that CfWor1 primarily regulates development of C. fulvum, which indirectly affects the expression of a subset of virulence genes.

Pseudogenization in pathogenic fungi with different host plants and lifestyles might reflect their evolutionary past.

Pseudogenes are genes with significant homology to functional genes, but contain disruptive mutations (DMs) leading to the production of non- or partially functional proteins. Little is known about pseudogenization in pathogenic fungi with different lifestyles. Here, we report the identification of DMs causing pseudogenes in the genomes of the fungal plant pathogens Botrytis cinerea, Cladosporium fulvum, Dothistroma septosporum, Mycosphaerella fijiensis, Verticillium dahliae and Zymoseptoria tritici. In these fungi, we identified 1740 gene models containing 2795 DMs obtained by an alignment-based gene prediction method. The contribution of sequencing errors to DMs was minimized by analyses of resequenced genomes to obtain a refined dataset of 924 gene models containing 1666 true DMs. The frequency of pseudogenes varied from 1% to 5% in the gene catalogues of these fungi, being the highest in the asexually reproducing fungus C. fulvum (4.9%), followed by D. septosporum (2.4%) and V. dahliae (2.1%). The majority of pseudogenes do not represent recent gene duplications, but members of multi-gene families and unitary genes. In general, there was no bias for pseudogenization of specific genes in the six fungi. Single exceptions were those encoding secreted proteins, including proteases, which appeared more frequently pseudogenized in C. fulvum than in D. septosporum. Most pseudogenes present in these two phylogenetically closely related fungi are not shared, suggesting that they are related to adaptation to a different host (tomato versus pine) and lifestyle (biotroph versus hemibiotroph).

Automated alignment-based curation of gene models in filamentous fungi.

BACKGROUND: Automated gene-calling is still an error-prone process, particularly for the highly plastic genomes of fungal species. Improvement through quality control and manual curation of gene models is a time-consuming process that requires skilled biologists and is only marginally performed. The wealth of available fungal genomes has not yet been exploited by an automated method that applies quality control of gene models in order to obtain more accurate genome annotations. RESULTS: We provide a novel method named alignment-based fungal gene prediction (ABFGP) that is particularly suitable for plastic genomes like those of fungi. It can assess gene models on a gene-by-gene basis making use of informant gene loci. Its performance was benchmarked on 6,965 gene models confirmed by full-length unigenes from ten different fungi. 79.4% of all gene models were correctly predicted by ABFGP. It improves the output of ab initio gene prediction software due to a higher sensitivity and precision for all gene model components. Applicability of the method was shown by revisiting the annotations of six different fungi, using gene loci from up to 29 fungal genomes as informants. Between 7,231 and 8,337 genes were assessed by ABFGP and for each genome between 1,724 and 3,505 gene model revisions were proposed. The reliability of the proposed gene models is assessed by an a posteriori introspection procedure of each intron and exon in the multiple gene model alignment. The total number and type of proposed gene model revisions in the six fungal genomes is correlated to the quality of the genome assembly, and to sequencing strategies used in the sequencing centre, highlighting different types of errors in different annotation pipelines. The ABFGP method is particularly successful in discovering sequence errors and/or disruptive mutations causing truncated and erroneous gene models. CONCLUSIONS: The ABFGP method is an accurate and fully automated quality control method for fungal gene catalogues that can be easily implemented into existing annotation pipelines. With the exponential release of new genomes, the ABFGP method will help decreasing the number of gene models that require additional manual curation.

At the origin of spliceosomal introns: Is multiplication of introner-like elements the main mechanism of intron gain in fungi?

The recent discovery of introner-like elements (ILEs) in six fungal species shed new light on the origin of regular spliceosomal introns (RSIs) and the mechanism of intron gains. These novel spliceosomal introns are found in hundreds of copies, are longer than RSIs and harbor stable predicted secondary structures. Yet, they are prone to degeneration in sequence and length to become undistinguishable from RSIs, suggesting that ILEs are predecessors of most RSIs. In most fungi, other near-identical introns were found duplicated in lower numbers in the same gene or in unrelated genes, indicating that intron duplication is a widespread phenomenon. However, ILEs are associated with the majority of intron gains, suggesting that the other types of duplication are of minor importance to the overall gains of introns. Our data support the hypothesis that ILEs' multiplication corresponds to the main mechanism of intron gain in fungi.

The genomes of the fungal plant pathogens Cladosporium fulvum and Dothistroma septosporum reveal adaptation to different hosts and lifestyles but also signatures of common ancestry.

We sequenced and compared the genomes of the Dothideomycete fungal plant pathogens Cladosporium fulvum (Cfu) (syn. Passalora fulva) and Dothistroma septosporum (Dse) that are closely related phylogenetically, but have different lifestyles and hosts. Although both fungi grow extracellularly in close contact with host mesophyll cells, Cfu is a biotroph infecting tomato, while Dse is a hemibiotroph infecting pine. The genomes of these fungi have a similar set of genes (70% of gene content in both genomes are homologs), but differ significantly in size (Cfu >61.1-Mb; Dse 31.2-Mb), which is mainly due to the difference in repeat content (47.2% in Cfu versus 3.2% in Dse). Recent adaptation to different lifestyles and hosts is suggested by diverged sets of genes. Cfu contains an α-tomatinase gene that we predict might be required for detoxification of tomatine, while this gene is absent in Dse. Many genes encoding secreted proteins are unique to each species and the repeat-rich areas in Cfu are enriched for these species-specific genes. In contrast, conserved genes suggest common host ancestry. Homologs of Cfu effector genes, including Ecp2 and Avr4, are present in Dse and induce a Cf-Ecp2- and Cf-4-mediated hypersensitive response, respectively. Strikingly, genes involved in production of the toxin dothistromin, a likely virulence factor for Dse, are conserved in Cfu, but their expression differs markedly with essentially no expression by Cfu in planta. Likewise, Cfu has a carbohydrate-degrading enzyme catalog that is more similar to that of necrotrophs or hemibiotrophs and a larger pectinolytic gene arsenal than Dse, but many of these genes are not expressed in planta or are pseudogenized. Overall, comparison of their genomes suggests that these closely related plant pathogens had a common ancestral host but since adapted to different hosts and lifestyles by a combination of differentiated gene content, pseudogenization, and gene regulation.

Birth of new spliceosomal introns in fungi by multiplication of introner-like elements.

Spliceosomal introns are noncoding sequences that separate exons in eukaryotic genes and are removed from pre-messenger RNAs by the splicing machinery. Their origin has remained a mystery in biology since their discovery because intron gains seem to be infrequent in many eukaryotic lineages. Although a few recent intron gains have been reported, none of the proposed gain mechanisms can convincingly explain the high number of introns in present-day eukaryotic genomes. Here we report on particular spliceosomal introns that share high sequence similarity and are reminiscent of introner elements. These elements multiplied in unrelated genes of six fungal genomes and account for the vast majority of intron gains in these fungal species. Such introner-like elements (ILEs) contain all typical characteristics of regular spliceosomal introns (RSIs) but are longer and predicted to harbor more stable secondary structures. However, dating of multiplication events showed that they degenerate in sequence and length within 100,000 years to eventually become indistinguishable from RSIs. We suggest that ILEs not only account for intron gains in six fungi but also in ancestral eukaryotes to give rise to most RSIs by a yet unknown multiplication mechanism.

Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis.

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.

In silico miRNA prediction in metazoan genomes: balancing between sensitivity and specificity.

BACKGROUND: MicroRNAs (miRNAs), short approximately 21-nucleotide RNA molecules, play an important role in post-transcriptional regulation of gene expression. The number of known miRNA hairpins registered in the miRBase database is rapidly increasing, but recent reports suggest that many miRNAs with restricted temporal or tissue-specific expression remain undiscovered. Various strategies for in silico miRNA identification have been proposed to facilitate miRNA discovery. Notably support vector machine (SVM) methods have recently gained popularity. However, a drawback of these methods is that they do not provide insight into the biological properties of miRNA sequences. RESULTS: We here propose a new strategy for miRNA hairpin prediction in which the likelihood that a genomic hairpin is a true miRNA hairpin is evaluated based on statistical distributions of observed biological variation of properties (descriptors) of known miRNA hairpins. These distributions are transformed into a single and continuous outcome classifier called the L score. Using a dataset of known miRNA hairpins from the miRBase database and an exhaustive set of genomic hairpins identified in the genome of Caenorhabditis elegans, a subset of 18 most informative descriptors was selected after detailed analysis of correlation among and discriminative power of individual descriptors. We show that the majority of previously identified miRNA hairpins have high L scores, that the method outperforms miRNA prediction by threshold filtering and that it is more transparent than SVM classifiers. CONCLUSION: The L score is applicable as a prediction classifier with high sensitivity for novel miRNA hairpins. The L-score approach can be used to rank and select interesting miRNA hairpin candidates for downstream experimental analysis when coupled to a genome-wide set of in silico-identified hairpins or to facilitate the analysis of large sets of putative miRNA hairpin loci obtained in deep-sequencing efforts of small RNAs. Moreover, the in-depth analyses of miRNA hairpins descriptors preceding and determining the L score outcome could be used as an extension to miRBase entries to help increase the reliability and biological relevance of the miRNA registry.

Large-scale gene discovery in the septoria tritici blotch fungus Mycosphaerella graminicola with a focus on in planta expression.

The foliar disease septoria tritici blotch, caused by the fungus Mycosphaerella graminicola, is currently the most important wheat disease in Europe. Gene expression was examined under highly different conditions, using 10 expressed sequence tag libraries generated from M. graminicola isolate IPO323 using seven in vitro and three in planta growth conditions. To identify fungal clones in the interaction libraries, we developed a selection method based on hybridization with the entire genomic DNA of M. graminicola, to selectively enrich these libraries for fungal genes. Assembly of the 27,007 expressed sequence tags resulted in 9,190 unigenes, representing 5.2 Mb of the estimated 39-Mb genome size of M. graminicola. All libraries contributed significantly to the number of unigenes, especially the in planta libraries representing different stages of pathogenesis, which covered 15% of the library-specific unigenes. Even under presymptomatic conditions (5 days postinoculation), when fungal biomass is less than 5%, this method enabled us to efficiently capture fungal genes expressed during pathogenesis. Many of these genes were uniquely expressed in planta, indicating that in planta gene expression significantly differed from in vitro expression. Examples of gene discovery included a number of cell wall-degrading enzymes, a broad set of genes involved in signal transduction (n=11) and a range of ATP-binding cassette (n=20) and major facilitator superfamily transporter genes (n=12) potentially involved in protection against antifungal compounds or the secretion of pathogenicity factors. In addition, evidence is provided for a mycovirus in M. graminicola that is highly expressed under various stress conditions, in particular, under nitrogen starvation. Our analyses provide a unique window on in vitro and in planta gene expression of M. graminicola.

High-throughput bioinformatics with the Cyrille2 pipeline system.

BACKGROUND: Modern omics research involves the application of high-throughput technologies that generate vast volumes of data. These data need to be pre-processed, analyzed and integrated with existing knowledge through the use of diverse sets of software tools, models and databases. The analyses are often interdependent and chained together to form complex workflows or pipelines. Given the volume of the data used and the multitude of computational resources available, specialized pipeline software is required to make high-throughput analysis of large-scale omics datasets feasible. RESULTS: We have developed a generic pipeline system called Cyrille2. The system is modular in design and consists of three functionally distinct parts: 1) a web based, graphical user interface (GUI) that enables a pipeline operator to manage the system; 2) the Scheduler, which forms the functional core of the system and which tracks what data enters the system and determines what jobs must be scheduled for execution, and; 3) the Executor, which searches for scheduled jobs and executes these on a compute cluster. CONCLUSION: The Cyrille2 system is an extensible, modular system, implementing the stated requirements. Cyrille2 enables easy creation and execution of high throughput, flexible bioinformatics pipelines.

The use of multiple hierarchically independent gene ontology terms in gene function prediction and genome annotation.

The Gene Ontology (GO) is a widely used controlled vocabulary for the description of gene function. In this study we quantify the usage of multiple and hierarchically independent GO terms in the curated genome annotations of seven well-studied species. In most genomes, significant proportions (6-60%) of genes have been annotated with multiple and hierarchically independent terms. This may be necessary to attain adequate specificity of description. One noticeable exception is Arabidopsis thaliana, in which genes are much less frequently annotated with multiple terms (6-14%). In contrast, an analysis of the occurrence of InterPro hits in the proteomes of the seven species, followed by a mapping of the hits to GO terms, did not reveal an aberrant pattern for the A. thaliana genome. This study shows the widespread usage of multiple hierarchically independent GO terms in the functional annotation of genes. By consequence, probabilistic methods that aim to predict gene function automatically through integration of diverse genomic datasets, and that employ the GO, must be able to predict such multiple terms. We attribute the low frequency with which multiple GO terms are used in Arabidopsis to deviating practices in the genome annotation and curation process between communities of annotators. This may bias genome-scale comparisons of gene function between different species. GO term assignment should therefore be performed according to strictly similar rules and standards.