Identification of Diverse Mycoviruses through Metatranscriptomics Characterization of the Viromes of Five Major Fungal Plant Pathogens.

UNLABELLED: Mycoviruses can have a marked effect on natural fungal communities and influence plant health and productivity. However, a comprehensive picture of mycoviral diversity is still lacking. To characterize the viromes of five widely dispersed plant-pathogenic fungi, Colletotrichum truncatum, Macrophomina phaseolina, Diaporthe longicolla, Rhizoctonia solani, and Sclerotinia sclerotiorum, a high-throughput sequencing-based metatranscriptomic approach was used to detect viral sequences. Total RNA and double-stranded RNA (dsRNA) from mycelia and RNA from samples enriched for virus particles were sequenced. Sequence data were assembled de novo, and contigs with predicted amino acid sequence similarities to viruses in the nonredundant protein database were selected. The analysis identified 72 partial or complete genome segments representing 66 previously undescribed mycoviruses. Using primers specific for each viral contig, at least one fungal isolate was identified that contained each virus. The novel mycoviruses showed affinity with 15 distinct lineages: Barnaviridae, Benyviridae, Chrysoviridae, Endornaviridae, Fusariviridae, Hypoviridae, Mononegavirales, Narnaviridae, Ophioviridae, Ourmiavirus, Partitiviridae, Tombusviridae, Totiviridae, Tymoviridae, and Virgaviridae More than half of the viral sequences were predicted to be members of the Mitovirus genus in the family Narnaviridae, which replicate within mitochondria. Five viral sequences showed strong affinity with three families (Benyviridae, Ophioviridae, and Virgaviridae) that previously contained no mycovirus species. The genomic information provides insight into the diversity and taxonomy of mycoviruses and coevolution of mycoviruses and their fungal hosts. IMPORTANCE: Plant-pathogenic fungi reduce crop yields, which affects food security worldwide. Plant host resistance is considered a sustainable disease management option but may often be incomplete or lacking for some crops to certain fungal pathogens or strains. In addition, the rising issues of fungicide resistance demand alternative strategies to reduce the negative impacts of fungal pathogens. Those fungus-infecting viruses (mycoviruses) that attenuate fungal virulence may be welcome additions for mitigation of plant diseases. By high-throughput sequencing of the RNAs from 275 isolates of five fungal plant pathogens, 66 previously undescribed mycoviruses were identified. In addition to identifying new potential biological control agents, these results expand the grand view of the diversity of mycoviruses and provide possible insights into the importance of intracellular and extracellular transmission in fungus-virus coevolution.

Molecular response to the pathogen Phytophthora sojae among ten soybean near isogenic lines revealed by comparative transcriptomics.

BACKGROUND: Phytophthora root and stem rot (PRR) of soybean, caused by Phytophthora sojae, is controlled by Rps genes. However, little is known regarding the Rps-induced molecular responses to P. sojae and how they actually overlap. We thus sequenced, analyzed, and compared the transcriptomes of 10 near isogenic lines (NILs), each with a unique Rps gene/allele, and the susceptible parent Williams, pre- and post-inoculation with the pathogen. RESULTS: A total of 4,330 differentially expressed genes (DEGs) were identified in Williams versus 2,014 to 5,499 DEGs in individual NILs upon inoculation with the pathogen. Comparisons of the DEGs between the NILs and Williams identified incompatible interaction genes (IIGs) and compatible interaction genes (CIGs). Hierarchical cluster and heatmap analyses consistently grouped the NILs into three clusters: Cluster I (Rps1-a), Cluster II (Rps1-b, 1-c and 1-k) and Cluster III (Rps3-a, 3-b, 3-c, 4, 5, and 6), suggesting an overlap in Rps-induced defense signaling among certain NILs. Gene ontology (GO) analysis revealed associations between members of the WRKY family and incompatible reactions and between a number of phytohormone signaling pathways and incompatible/compatible interactions. These associations appear to be distinguished according to the NIL clusters. CONCLUSIONS: This study characterized genes and multiple branches of putative regulatory networks associated with resistance to P. sojae in ten soybean NILs, and depicted functional "fingerprints" of individual Rps-mediated resistance responses through comparative transcriptomic analysis. Of particular interest are dramatic variations of detected DEGs, putatively involved in ethylene (ET)-, jasmonic acid (JA)-, (reactive oxygen species) ROS-, and (MAP-kinase) MAPK- signaling, among these soybean NILs, implicating their important roles of these signaling in differentiating molecular defense responses. We hypothesize that different timing and robustness in defense signaling to the same pathogen may be largely responsible for such variations.

Genetic architecture and evolution of the mating type locus in fusaria that cause soybean sudden death syndrome and bean root rot.

Fusarium tucumaniae is the only known sexually reproducing species among the seven closely related fusaria that cause soybean sudden death syndrome (SDS) or bean root rot (BRR). In a previous study, laboratory mating of F. tucumaniae yielded recombinant ascospore progeny but required two mating-compatible strains, indicating that it is heterothallic. To assess the reproductive mode of the other SDS and BRR fusaria, and their potential for mating, whole-genome sequences of two SDS and one BRR pathogen were analyzed to characterize their mating type (MAT) loci. This bioinformatic approach identified a MAT1-1 idiomorph in F. virguliforme NRRL 22292 and MAT1-2 idiomorphs in F. tucumaniae NRRL 34546 and F. azukicola NRRL 54364. Alignments of the MAT loci were used to design PCR primers within the conserved regions of the flanking genes APN1 and SLA2, which enabled primer walking to obtain nearly complete sequences of the MAT region for six MAT1-1 and five MAT1-2 SDS/BRR fusaria. As expected, sequences of the highly divergent 4.7 kb MAT1-1 and 3.7 kb MAT1-2 idiomorphs were unalignable. However, sequences of the respective idiomorphs and those that flank MAT1-1 and MAT1-2 were highly conserved. In addition to three genes at MAT1-1 (MAT1-1-1, MAT1-1-2, MAT1-1-3) and two at MAT1-2 (MAT1-2-1, MAT1-2-3), the MAT loci of the SDS/BRR fusaria also include a putative gene predicted to encode for a 252 amino acid protein of unknown function. Alignments of the MAT1-1-3 and MAT1-2-1 sequences were used to design a multiplex PCR assay for the MAT loci. This assay was used to screen DNA from 439 SDS/BRR isolates, which revealed that each isolate possessed MAT1-1 or MAT1-2, consistent with heterothallism. Both idiomorphs were represented among isolates of F. azukicola, F. brasiliense, F. phaseoli and F. tucumaniae, whereas isolates of F. virguliforme and F. cuneirostrum were only MAT1-1 and F. crassistipitatum were only MAT1-2. Finally, nucleotide sequence data from the RPB1 and RPB2 genes were used to date the origin of the SDS/BRR group, which was estimated to have occurred about 0.75 Mya (95% HPD interval: 0.27, 1.68) in the mid-Pleistocene, long before the domestication of the common bean or soybean.

Molecular mapping of two genes conferring resistance to Phytophthora sojae in a soybean landrace PI 567139B.

Phytophthora root and stem rot (PRR), caused by the soil-borne oomycete pathogen Phytophthora sojae, is one of the most destructive diseases of soybean. PRR can be effectively controlled by race-specific genes conferring resistance to P. sojae (Rps). However, the Rps genes are usually non-durable, as populations of P. sojae are highly diverse and quick to adapt, and can be overcome 8-15 years after deployment. Thus, it is important to identify novel Rps genes for development of resistant soybean cultivars. PI 567139B is a soybean landrace carrying excellent resistance to nearly all predominant P. sojae races in Indiana. A mapping population consisting of 245 F2 individuals and 403 F2:3 families was developed from a cross between PI 567139B and the susceptible cultivar 'Williams', and used to dissect the resistance carried by PI 567139B. We found that the resistance in PI 567139B was conferred by two independent Rps genes, designated RpsUN1 and RpsUN2. The former was mapped to a 6.5 cM region between SSR markers Satt159 and BARCSOYSSR_03_0250 that spans the Rps1 locus on chromosome 3, while the latter was mapped to a 3.0 cM region between BARCSOYSSR_16_1275 and Sat_144, approximately 3.0-3.4 cM upstream of Rps2 on chromosome 16. According to the 'Williams 82' reference genome sequence, both regions are highly enriched with NBS-LRR genes. Marker assisted resistance spectrum analyses of these genes with 16 isolates of P. sojae, in combination with the mapping results, suggested that RpsUN1 was likely to be a novel allele at the Rps1 locus, while RpsUN2 was more likely to be a novel Rps gene.

Heritability of Resistance to Verticillium Wilt in Alfalfa.

Verticillium wilt of alfalfa, caused by Verticillium albo-atrum, may reduce forage yields by up to 50% in alfalfa-producing areas of the northern United States and Canada. It has been suggested that cultivars require at least 60% resistant plants to afford maximum protection against disease. Our objective was to calculate heritability estimates of resistance to Verticillium wilt in alfalfa. Estimates were generated for two alfalfa populations developed from the cvs. Affinity + Z and Depend + EV. Heritability on a half-sib progeny means basis was calculated based on data from greenhouse pathogenicity tests. Estimates based on repeated experiments conducted for single years (2004 and 2005) were high for both populations, ranging from 0.86 to 0.92. The heritability estimate based on data collected over 2 years was 0.26 for Affinity + Z and 0.66 for Depend + EV. Disease was more severe in 2005 than in 2004. However, the Spearman rank correlation between mean disease severity index values for half-sib families over 2 years was positive and significant for both populations. Results of pathogenicity tests suggested that neither cultivar had resistance levels approaching 60%. The heritability estimates suggest that resistance levels in both Affinity + Z and Depend + EV could be improved further through selection.

Coordination of MicroRNAs, PhasiRNAs, and NB-LRR Genes in Response to a Plant Pathogen: Insights from Analyses of a Set of Soybean Rps Gene Near-Isogenic Lines.

Disease-related genes, particularly the nucleotide binding site (NB)-leucine-rich repeat (LRR) class of R plant genes can be triggered by microRNAs (miRNAs) to generate phased small interfering RNAs (phasiRNAs), which could reduce the transcript levels of their targets. However, how global changes in NB-LRR transcript levels coordinate with changes in miRNA and phasiRNA levels in defense responses remains largely unknown. Here, we investigated changes in the relative abundance of small RNAs (sRNAs), with a focus on miRNAs and phasiRNAs and their potential targets in response to the pathogen Phytophthora sojae in the susceptible soybean [Glycine max (L.) Merr.] 'Williams' and nine resistant near-isogenic lines (NILs), each carrying a unique resistance to P. sojae (Rps) gene. In total, 369 distinct miRNAs, including 78 new ones, were identified in the 10 soybean lines. The majority of miRNAs were downregulated by the pathogen. Of the 525 NB-LRR genes found in the soybean reference genome, 257 were predicted to be the targets of eight abundant miRNA families and 126 (dubbed phasi-NB-LRRs or pNLs) were predicted to have produced phasiRNAs. Upregulation of 15 phasi-NB-LRRs was associated with downregulation of their corresponding phasiRNAs in the NILs; these phasiRNAs were predicted to regulate 75 additional NB-LRRs in trans. In addition, we identified putative 24-nucleotide (nt) phasiRNAs from transposons, possibly representing a novel general epigenetic mechanism for regulation of transposon activity under biotic stresses. Together, these observations suggest that miRNAs and phasiRNAs play an important role in response to plant pathogens through complex, multiple layers of post-transcriptional regulation.

Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean.

The rhg1-b allele of soybean is widely used for resistance against soybean cyst nematode (SCN), the most economically damaging pathogen of soybeans in the United States. Gene silencing showed that genes in a 31-kilobase segment at rhg1-b, encoding an amino acid transporter, an α-SNAP protein, and a WI12 (wound-inducible domain) protein, each contribute to resistance. There is one copy of the 31-kilobase segment per haploid genome in susceptible varieties, but 10 tandem copies are present in an rhg1-b haplotype. Overexpression of the individual genes in roots was ineffective, but overexpression of the genes together conferred enhanced SCN resistance. Hence, SCN resistance mediated by the soybean quantitative trait locus Rhg1 is conferred by copy number variation that increases the expression of a set of dissimilar genes in a repeated multigene segment.