Development of a Decision Support System for the Management of Mummy Berry Disease in Northwestern Washington

Mummy berry, caused by Monilinia vaccinii-corymbosi, is the most important disease of the northern highbush blueberry (Vaccinium corymbosum L.) in North America and can cause up to 70% yield losses in affected fields. A key event in the mummy berry disease cycle is the primary infection phase where ascospores are released by apothecia that infect emerging floral and vegetative tissues. Current management of mummy berry disease in northwestern Washington is predominantly reliant on the prevention of primary infections through prophylactic, calendar-based fungicide spray applications early in the growing season. To improve the understanding of risk during these periods and to help tailor management strategies, we developed a decision support system (DSS) based on field records spanning over five seasons and four locations in northwestern Washington. Environmental conditions across the region were highly uniform but different dynamics of apothecial development were observed under high- and low-management regimes. Based on our analysis, we suggest basing the initial iteration of the DSS on two sub-models. The first sub-model predicts the onset of apothecia based on chill-unit accumulation under high- and low-management regimes, and the second predicts primary infection risk, which provides opportunities to improve the timing of fungicide applications. The synoptic DSS proposed here is based on the current biological knowledge of the pathosystem and available data for the northwestern Washington region. We provide the analysis and the DSS implementation and evaluation as an open-source repository, providing opportunities for further improvements. Finally, we provide suggestions for future research and the operational efforts needed for improving the utility and accuracy of the mummy berry DSS.publishedVersio

Genetic Variation Among Natural Populations of Tilletia controversa and T. bromi

Isolates of Tilletia controversa and T. bromi were sampled from wheat and two Bromus species hosts, respectively, in the Pacific Northwest, and genetic variation within and among populations was determined. Fifty-one random amplified polymorphic DNA markers from eleven primers were treated as phenotypic 1 and 0 character state data to estimate similarities and analyze molecular variance (AMOVA) among populations and as putative genetic loci to carry out analyses of gene diversity. Phenotypic analysis of T. controversa and T. bromi isolates revealed two distinct clusters that were 37% similar. The T. bromi cluster was subdivided further into two groups, corresponding to host, with 40% similarity. Cluster analysis based on allele frequencies produced similar results and also supported two T. bromi groups based on host. No evidence of natural hybridization and introgression was detected between the T. controversa and T. bromi populations. Both AMOVA and gene diversity analyses detected moderate levels of differentiation among T. controversa populations, whereas T. bromi populations were highly differentiated. The level of genetic differentiation observed between the T. bromi populations on different Bromus species hosts supports the hypothesis that a high degree of host specificity exists in the wild grass-infecting smuts. We speculate that the higher level of genetic differentiation among the T. bromi populations compared with the T. controversa populations on wheat may be due to selection by a more genetically diverse host population

Diversity of Hypoviruses and Other Double-Stranded RNAs in Cryphonectria parasitica in North America

Double-stranded (ds) RNAs in Cryphonectria parasitica were randomly sampled from nine subpopulations in North America using an antibody-based detection system for dsRNA. dsRNA was detected in 166 (28%) of a total of 595 C. parasitica isolates sampled by immunoblotting. Incidence of dsRNA infection within subpopulations ranged from 0% in samples from New Hampshire and Ontario to 100% in County Line, MI. Most of the dsRNAs sampled were approximately 9 to 13 kb in size. dsRNAs from 72 isolates analyzed by probing Northern blots with 32P-labeled dsRNAs were in one of three hybridization groups. One hybridization group was widespread throughout eastern North America, being found in New York, New Jersey, Maryland, West Virginia, Kentucky, and Michigan. These dsRNAs hybridized to dsRNA from the previously described C. parasitica isolate SR2 from Maryland and are referred to as SR2-type dsRNAs. The second hybridization group was found almost exclusively in Michigan. The Michigan dsRNAs cross-hybridized to Cryphonectria hypovirus 3-GH2 (CHV3-GH2) and are referred to as CHV3-type dsRNAs.One dsRNA sampled from Kentucky hybridized to CHV3-type dsRNAs from Michigan. This dsRNA was probably derived from a fungal isolate that had been intentionally released for biological control at this same site 10 years previously and had become established in Kentucky. The third hybridization group was found only in New Jersey. These dsRNAs were much smaller than all other dsRNAs (3 and 5 kb) and were found in all 11 isolates that were probed; two of these isolates also had SR2-type dsRNA in mixed infections. None of the North American dsRNAs hybridized to CHV1 from Europe, even though CHV1 has been released in numerous locations in eastern North America for biological control of chestnut blight. Similarly, no dsRNAs hybridized to CHV2-NB58, a hypovirus found previously in New Jersey. Mixed infections of SR2-type and CHV3-type dsRNAs were found in 13 of 15 isolates from Frankfort, MI, while another nearby subpopulation (County Line) was infected with only CHV3-type dsRNAs. The distribution of dsRNA hybridization groups in C. parasitica thus presents a mixed picture, since one hybridization group is widespread, whereas two others are primarily restricted to smaller geographic areas

Variation in Tolerance and Virulence in the Chestnut Blight Fungus-Hypovirus Interaction

Chestnut blight, caused by the fungus Cryphonectria parasitica , has been effectively controlled with double-stranded RNA hypoviruses in Europe for over 40 years. The marked reduction in the virulence of C. parasitica by hypoviruses is a phenomenon known as hypovirulence. This virus-fungus pathosystem has become a model system for the study of biological control of fungi with viruses. We studied variation in tolerance to hypoviruses in fungal hosts and variation in virulence among virus isolates from a local population in Italy. Tolerance is defined as the relative fitness of a fungal individual when infected with hypoviruses (compared to being uninfected); virulence is defined for each hypovirus as the reduction in fitness of fungal hosts relative to virus-free hosts. Six hypovirus-infected isolates of C. parasitica were sampled from the population, and each hypovirus was transferred into six hypovirus-free recipient isolates. The resulting 36 hypovirus-fungus combinations were used to estimate genetic variation in tolerance to hypoviruses, in hypovirus virulence, and in virus-fungus interactions. Four phenotypes were evaluated for each virus-fungus combination to estimate relative fitness: (i) sporulation, i.e., the number of asexual spores (conidia) produced; (ii) canker area on field-inoculated chestnut trees, (iii) vertical transmission of hypoviruses into conidia, and (iv) conidial germination. Two-way analysis of variance (ANOVA) revealed significant interactions ( P < 0.001) between viruses and fungal isolates for sporulation and canker area but not for conidial germination or transmission. One-way ANOVA among hypoviruses (within each fungal isolate) and among fungal isolates (within each hypovirus) revealed significant genetic variation ( P < 0.01) in hypovirus virulence and fungal tolerance within several fungal isolates, and hypoviruses, respectively. These interactions and the significant genetic variation in several fitness characters indicate the potential for future evolution of these characters. However, biological control is unlikely to break down due to evolution of tolerance to hypoviruses in the fungus because the magnitudes of tolerance and interactions were relatively small

Contrasting Codon Usage Patterns and Purifying Selection at the Mating Locus in Putatively Asexual Alternaria Fungal Species

Sexual reproduction in heterothallic ascomycete fungi is controlled by a single mating-type locus called MAT1 with two alternate alleles or idiomorphs, MAT1-1 and MAT1-2 . These alleles lack sequence similarity and encode different transcriptional regulators. A large number of phytopathogenic fungi including Alternaria spp. are considered asexual, yet still carry expressed MAT1 genes. The molecular evolution of Alternaria MAT1 was explored using nucleotide diversity, nonsynonymous vs. synonymous substitution ( dn/ds ) ratios and codon usage statistics. Likelihood ratio tests of site-branch models failed to detect positive selection on MAT1-1-1 or MAT1-2-1 . Codon-site models demonstrated that both MAT1-1-1 and MAT1-2-1 are under purifying selection and significant differences in codon usage were observed between MAT1-1-1 and MAT1-2-1 . Mean GC content at the third position (GC3) and effective codon usage (ENC) were significantly different between MAT1-1-1 and MAT1-2-1 with values of 0.57 and 48 for MAT1-1-1 and 0.62 and 46 for MAT1-2-1 , respectively. In contrast, codon usage of Pleospora spp. (anamorph Stemphylium ), a closely related Dothideomycete genus, was not significantly different between MAT1-1-1 and MAT1-2-1 . The purifying selection and biased codon usage detected at the MAT1 locus in Alternaria spp. suggest a recent sexual past, cryptic sexual present and/or that MAT1 plays important cellular role(s) in addition to mating

A Novel Type Pathway-Specific Regulator and Dynamic Genome Environments of a Solanapyrone Biosynthesis Gene Cluster in the Fungus Ascochyta rabiei

Secondary metabolite genes are often clustered together and situated in particular genomic regions, like the subtelomere, that can facilitate niche adaptation in fungi. Solanapyrones are toxic secondary metabolites produced by fungi occupying different ecological niches. Full-genome sequencing of the ascomycete Ascochyta rabiei revealed a solanapyrone biosynthesis gene cluster embedded in an AT-rich region proximal to a telomere end and surrounded by Tc1/Mariner-type transposable elements. The highly AT-rich environment of the solanapyrone cluster is likely the product of repeat-induced point mutations. Several secondary metabolism-related genes were found in the flanking regions of the solanapyrone cluster. Although the solanapyrone cluster appears to be resistant to repeat-induced point mutations, a P450 monooxygenase gene adjacent to the cluster has been degraded by such mutations. Among the six solanapyrone cluster genes (sol1 to sol6), sol4 encodes a novel type of Zn(II)2Cys6 zinc cluster transcription factor. Deletion of sol4 resulted in the complete loss of solanapyrone production but did not compromise growth, sporulation, or virulence. Gene expression studies with the sol4 deletion and sol4-overexpressing mutants delimited the boundaries of the solanapyrone gene cluster and revealed that sol4 is likely a specific regulator of solanapyrone biosynthesis and appears to be necessary and sufficient for induction of the solanapyrone cluster genes. Despite the dynamic surrounding genomic regions, the solanapyrone gene cluster has maintained its integrity, suggesting important roles of solanapyrones in fungal biology

Predicting Ascospore Release of Monilinia vaccinii-corymbosi of Blueberry with Machine Learning

Mummy berry, caused by Monilinia vaccinii-corymbosi, causes economic losses of highbush blueberry in the U.S. Pacific Northwest (PNW). Apothecia develop from mummified berries overwintering on soil surfaces and produce ascospores that infect tissue emerging from floral and vegetative buds. Disease control currently relies on fungicides applied on a calendar basis rather than inoculum availability. To establish a prediction model for ascospore release, apothecial development was tracked in three fields, one in western Oregon and two in northwestern Washington in 2015 and 2016. Air and soil temperature, precipitation, soil moisture, leaf wetness, relative humidity and solar radiation were monitored using in-field weather stations and Washington State University's AgWeatherNet stations. Four modeling approaches were compared: logistic regression, multivariate adaptive regression splines, artificial neural networks, and random forest. A supervised learning approach was used to train the models on two data sets: training (70%) and testing (30%). The importance of environmental factors was calculated for each model separately. Soil temperature, soil moisture, and solar radiation were identified as the most important factors influencing ascospore release. Random forest models, with 78% accuracy, showed the best performance compared with the other models. Results of this research helps PNW blueberry growers to optimize fungicide use and reduce production costs

Development and Characterization of 37 Novel EST-SSR Markers in Pisum sativum

Premise of the study: Simple sequence repeat markers were developed based on expressed sequence tags (EST-SSR) and screened for polymorphism among 23 Pisum sativum individuals to assist development and refinement of pea linkage maps. In particular, the SSR markers were developed to assist in mapping of white mold disease resistance quantitative trait loci. Methods and Results: Primer pairs were designed for 46 SSRs identified in EST contiguous sequences assembled from a 454 pyrosequenced transcriptome of the pea cultivar, ‘LIFTER’. Thirty-seven SSR markers amplified PCR products, of which 11 (30%) SSR markers produced polymorphism in 23 individuals, including parents of recombinant inbred lines, with two to four alleles. The observed and expected heterozygosities ranged from 0 to 0.43 and from 0.31 to 0.83, respectively. Conclusions: These EST-SSR markers for pea will be useful for refinement of pea linkage maps, and will likely be useful for comparative mapping of pea and as tools for marker-based pea breeding

Rapid transcriptome characterization and parsing of sequences in a non-model host-pathogen interaction; pea-<it>Sclerotinia sclerotiorum</it>

Abstract Background White mold, caused by Sclerotinia sclerotiorum, is one of the most important diseases of pea (Pisum sativum L.), however, little is known about the genetics and biochemistry of this interaction. Identification of genes underlying resistance in the host or pathogenicity and virulence factors in the pathogen will increase our knowledge of the pea-S. sclerotiorum interaction and facilitate the introgression of new resistance genes into commercial pea varieties. Although the S. sclerotiorum genome sequence is available, no pea genome is available, due in part to its large genome size (~3500 Mb) and extensive repeated motifs. Here we present an EST data set specific to the interaction between S. sclerotiorum and pea, and a method to distinguish pathogen and host sequences without a species-specific reference genome. Results 10,158 contigs were obtained by de novo assembly of 128,720 high-quality reads generated by 454 pyrosequencing of the pea-S. sclerotiorum interactome. A method based on the tBLASTx program was modified to distinguish pea and S. sclerotiorum ESTs. To test this strategy, a mixture of known ESTs (18,490 pea and 17,198 S. sclerotiorum ESTs) from public databases were pooled and parsed; the tBLASTx method successfully separated 90.1% of the artificial EST mix with 99.9% accuracy. The tBLASTx method successfully parsed 89.4% of the 454-derived EST contigs, as validated by PCR, into pea (6,299 contigs) and S. sclerotiorum (2,780 contigs) categories. Two thousand eight hundred and forty pea ESTs and 996 S. sclerotiorum ESTs were predicted to be expressed specifically during the pea-S. sclerotiorum interaction as determined by homology search against 81,449 pea ESTs (from flowers, leaves, cotyledons, epi- and hypocotyl, and etiolated and light treated etiolated seedlings) and 57,751 S. sclerotiorum ESTs (from mycelia at neutral pH, developing apothecia and developing sclerotia). Among those ESTs specifically expressed, 277 (9.8%) pea ESTs were predicted to be involved in plant defense and response to biotic or abiotic stress, and 93 (9.3%) S. sclerotiorum ESTs were predicted to be involved in pathogenicity/virulence. Additionally, 142 S. sclerotiorum ESTs were identified as secretory/signal peptides of which only 21 were previously reported. Conclusions We present and characterize an EST resource specific to the pea-S. sclerotiorum interaction. Additionally, the tBLASTx method used to parse S. sclerotiorum and pea ESTs was demonstrated to be a reliable and accurate method to distinguish ESTs without a reference genome.</p