The pursuit of genetic gain in agricultural crops through the application of machine-learning to genomic prediction

Current practice in agriculture applies genomic prediction to assist crop breeding in the analysis of genetic marker data. Genomic selection methods typically use linear mixed models, but using machine-learning may provide further potential for improved selection accuracy, or may provide additional information. Here we describe SelectML, an automated pipeline for testing and comparing the performance of a range of linear mixed model and machine-learning-based genomic selection methods. We demonstrate the use of SelectML on an in silico-generated marker dataset which simulated a randomly-sampled (mixed) and an unevenly-sampled (unbalanced) population, comparing the relative performance of various methods included in SelectML on the two datasets. Although machine-learning based methods performed similarly overall to linear mixed models, they performed worse on the mixed dataset and marginally better on the unbalanced dataset, being more affected than linear mixed models by the imposed sampling bias. SelectML can assist in the training, comparison, and selection of genomic selection models, and is available from https://github.com/darcyabjones/selectml.

Association mapping combined with whole genome sequencing data reveals candidate causal variants for sclerotinia stem rot resistance in Brassica napus

Canola (Brassica napus) yield can be significantly reduced by the disease sclerotinia stem rot (SSR), which is caused by Sclerotinia sclerotiorum, a necrotrophic fungal pathogen with an unusually large host range. Breeding cultivars that are physiologically resistant to SSR is desirable to enhance crop productivity. However, the development of resistant varieties has proved challenging due to the highly polygenic nature of S. sclerotiorum resistance. Here, we identified regions of the B. napus genome associated with SSR resistance using data from a previous study by association mapping. We then validated their contribution to resistance in a follow-up screen. This follow-up screen also confirmed high levels of SSR resistance in several genotypes from the previous study. Using publicly available whole genome sequencing data for a panel of 83 B. napus genotypes, we identified non-synonymous polymorphisms linked to the SSR resistance loci. A qPCR analysis showed that two of the genes containing these polymorphisms were transcriptionally responsive to S. sclerotiorum infection. In addition, we provide evidence that homologues of three of the candidate genes contribute to resistance in the model Brassicaceae species Arabidopsis thaliana. The identification of resistant germplasm and candidate genomic loci associated with resistance are important findings that may be exploited by breeders to improve the genetic resistance of canola varieties.

Small RNAs from the plant pathogenic fungus Sclerotinia sclerotiorum highlight host candidate genes associated with quantitative disease resistance

Fungal plant pathogens secrete effector proteins and metabolites to cause disease. Additionally, some species transfer small RNAs (sRNAs) into plant cells to silence host mRNAs through complementary base pairing and suppress plant immunity. The fungus Sclerotinia sclerotiorum infects over 600 plant species, but little is known about the molecular processes that govern interactions with its many hosts. In particular, evidence for the production of sRNAs by S. sclerotiorum during infection is lacking. We sequenced sRNAs produced by S. sclerotiorum in vitro and during infection of two host species, Arabidopsis thaliana and Phaseolus vulgaris. We found that S. sclerotiorum produces at least 374 distinct highly abundant sRNAs during infection, mostly originating from repeat-rich plastic genomic regions. We predicted the targets of these sRNAs in A. thaliana and found that these genes were significantly more down-regulated during infection than the rest of the genome. Predicted targets of S. sclerotiorum sRNAs in A. thaliana were enriched for functional domains associated with plant immunity and were more strongly associated with quantitative disease resistance in a genome-wide association study (GWAS) than the rest of the genome. Mutants in A. thaliana predicted sRNA target genes SERK2 and SNAK2 were more susceptible to S. sclerotiorum than wild-type, suggesting that S. sclerotiorum sRNAs may contribute to the silencing of immune components in plants. The prediction of fungal sRNA targets in plant genomes can be combined with other global approaches, such as GWAS, to assist in the identification of plant genes involved in quantitative disease resistance.

Prediction of pathogenicity genes involved in adaptation to a lupin host in the fungal pathogens Botrytis cinerea and Sclerotinia sclerotiorum via comparative genomics

Background

Narrow-leafed lupin is an emerging crop of significance in agriculture, livestock feed and human health food. However, its susceptibility to various diseases is a major obstacle towards increased adoption. Sclerotinia sclerotiorum and Botrytis cinerea – both necrotrophs with broad host-ranges - are reported among the top 10 lupin pathogens. Whole-genome sequencing and comparative genomics are useful tools to discover genes responsible for interactions between pathogens and their hosts.

Results

Genomes were assembled for one isolate of B. cinerea and two isolates of S. sclerotiorum, which were isolated from either narrow-leafed or pearl lupin species. Comparative genomics analysis between lupin-derived isolates and others isolated from alternate hosts was used to predict between 94 to 98 effector gene candidates from among their respective non-conserved gene contents.

Conclusions

Detection of minor differences between relatively recently-diverged isolates, originating from distinct regions and with hosts, may highlight novel or recent gene mutations and losses resulting from host adaptation in broad host-range fungal pathogens.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5774-2) contains supplementary material, which is available to authorized users.

The complete genome sequence of the phytopathogenic fungus Sclerotinia sclerotiorum reveals insights into the genome architecture of broad host range pathogens

Sclerotinia sclerotiorum is a phytopathogenic fungus with over 400 hosts including numerous economically important cultivated species. This contrasts many economically destructive pathogens that only exhibit a single or very few hosts. Many plant pathogens exhibit a “two-speed” genome. So described because their genomes contain alternating gene rich, repeat sparse and gene poor, repeat-rich regions. In fungi, the repeat-rich regions may be subjected to a process termed repeat-induced point mutation (RIP). Both repeat activity and RIP are thought to play a significant role in evolution of secreted virulence proteins, termed effectors. We present a complete genome sequence of S. sclerotiorum generated using Single Molecule Real-Time Sequencing technology with highly accurate annotations produced using an extensive RNA sequencing data set. We identified 70 effector candidates and have highlighted their in planta expression profiles. Furthermore, we characterized the genome architecture of S. sclerotiorum in comparison to plant pathogens that exhibit “two-speed” genomes. We show that there is a significant association between positions of secreted proteins and regions with a high RIP index in S. sclerotiorum but we did not detect a correlation between secreted protein proportion and GC content. Neither did we detect a negative correlation between CDS content and secreted protein proportion across the S. sclerotiorum genome. We conclude that S. sclerotiorum exhibits subtle signatures of enhanced mutation of secreted proteins in specific genomic compartments as a result of transposition and RIP activity. However, these signatures are not observable at the whole-genome scale.

Codon optimization underpins generalist parasitism in fungi

The range of hosts that parasites can infect is a key determinant of the emergence and spread of disease. Yet, the impact of host range variation on the evolution of parasite genomes remains unknown. Here, we show that codon optimization underlies genome adaptation in broad host range parasites. We found that the longer proteins encoded by broad host range fungi likely increase natural selection on codon optimization in these species. Accordingly, codon optimization correlates with host range across the fungal kingdom. At the species level, biased patterns of synonymous substitutions underpin increased codon optimization in a generalist but not a specialist fungal pathogen. Virulence genes were consistently enriched in highly codon-optimized genes of generalist but not specialist species. We conclude that codon optimization is related to the capacity of parasites to colonize multiple hosts. Our results link genome evolution and translational regulation to the long-term persistence of generalist parasitism.

Unraveling incompatibility between wheat and the fungal pathogen Zymoseptoria tritici through apoplastic proteomics

BACKGROUND

Hemibiotrophic fungal pathogen Zymoseptoria tritici causes severe foliar disease in wheat. However, current knowledge of molecular mechanisms involved in plant resistance to Z. tritici and Z. tritici virulence factors is far from being complete. The present work investigated the proteome of leaf apoplastic fluid with emphasis on both host wheat and Z. tritici during the compatible and incompatible interactions.

RESULTS

The proteomics analysis revealed rapid host responses to the biotrophic growth, including enhanced carbohydrate metabolism, apoplastic defenses and stress, and cell wall reinforcement, might contribute to resistance. Compatibility between the host and the pathogen was associated with inactivated plant apoplastic responses as well as fungal defenses to oxidative stress and perturbation of plant cell wall during the initial biotrophic stage, followed by the strong induction of plant defenses during the necrotrophic stage. To study the role of anti-oxidative stress in Z. tritici pathogenicity in depth, a YAP1 transcription factor regulating antioxidant expression was deleted and showed the contribution to anti-oxidative stress in Z. tritici, but was not required for pathogenicity. This result suggests the functional redundancy of antioxidants in the fungus.

CONCLUSIONS

The data demonstrate that incompatibility is probably resulted from the proteome-level activation of host apoplastic defenses as well as fungal incapability to adapt to stress and interfere with host cell at the biotrophic stage of the interaction.

Transcriptome and metabolite profiling of the infection cycle of Zymoseptoria tritici on wheat reveals a biphasic interaction with plant immunity involving differential pathogen chromosomal contributions and a variation on the hemibiotrophic lifestyle definition

The hemibiotrophic fungus Zymoseptoria tritici causes Septoria tritici blotch disease of wheat (Triticum aestivum). Pathogen reproduction on wheat occurs without cell penetration, suggesting that dynamic and intimate intercellular communication occurs between fungus and plant throughout the disease cycle. We used deep RNA sequencing and metabolomics to investigate the physiology of plant and pathogen throughout an asexual reproductive cycle of Z. tritici on wheat leaves. Over 3,000 pathogen genes, more than 7,000 wheat genes, and more than 300 metabolites were differentially regulated. Intriguingly, individual fungal chromosomes contributed unequally to the overall gene expression changes. Early transcriptional down-regulation of putative host defense genes was detected in inoculated leaves. There was little evidence for fungal nutrient acquisition from the plant throughout symptomless colonization by Z. tritici, which may instead be utilizing lipid and fatty acid stores for growth. However, the fungus then subsequently manipulated specific plant carbohydrates, including fructan metabolites, during the switch to necrotrophic growth and reproduction. This switch coincided with increased expression of jasmonic acid biosynthesis genes and large-scale activation of other plant defense responses. Fungal genes encoding putative secondary metabolite clusters and secreted effector proteins were identified with distinct infection phase-specific expression patterns, although functional analysis suggested that many have overlapping/redundant functions in virulence. The pathogenic lifestyle of Z. tritici on wheat revealed through this study, involving initial defense suppression by a slow-growing extracellular and nutritionally limited pathogen followed by defense (hyper) activation during reproduction, reveals a subtle modification of the conceptual definition of hemibiotrophic plant infection.

Trichothecenes: reference materials and method validation

The frequent contamination of food and feed with trichothecene mycotoxins, the high consumption of these products, and the potential risk associated herewith, has led to an increasing public awareness and therefore to the establishment of measures to control trichothecene contamination. The analytical difficulty and the economic importance of controlling trichothecenes in food and feed support the need for certified reference materials (CRMs) and validated methods. They form invaluable tools to ensure comparability and traceability in analytical measurements and are very useful for the implementation of written standards, legislation/regulations and laboratory accreditation. The present paper provides an overview of previous work, current strategies and prospectives for the production of CRMs and validation of analytical methods in the field of trichothecene analysis. Additional information is given on methodological demands, normative frameworks and commonly accepted procedures.