Exploring the genetics of lesion and nodal resistance in pea (Pisum sativum L.) to Sclerotinia sclerotiorum using genome-wide association studies and RNA-Seq

Application of bulk segregant RNA-Seq (BSR-Seq) and allele-specific primers to study soybean powdery mildew resistance

BACKGROUND

Powdery mildew (PM) is one of the important soybean diseases, and host resistance could practically contribute to soybean PM management. To date, only the Rmd locus on chromosome (Chr) 16 was identified through traditional QTL mapping and GWAS, and it remains unclear if the bulk segregant RNA-Seq (BSR-Seq) methodology is feasible to explore additional PM resistance that might exist in other varieties.

RESULTS

BSR-Seq was applied to contrast genotypes and gene expressions between the resistant bulk (R bulk) and the susceptible bulk (S bulk), as well as the parents. The ∆(SNP-index) and G' value identified several QTL and significant SNPs/Indels on Chr06, Chr15, and Chr16. Differentially expressed genes (DEGs) located within these QTL were identified using HISAT2 and Kallisto, and allele-specific primers (AS-primers) were designed to validate the accuracy of phenotypic prediction. While the AS-primers on Chr06 or Chr15 cannot distinguish the resistant and susceptible phenotypes, AS-primers on Chr16 exhibited 82% accuracy prediction with an additive effect, similar to the SSR marker Satt431.

CONCLUSIONS

Evaluation of additional AS-primers in the linkage disequilibrium (LD) block on Chr16 further confirmed the resistant locus, derived from the resistant parental variety 'Kaohsiung 11' ('KS11'), not only overlaps with the Rmd locus with unique up-regulated LRR genes (Glyma.16G213700 and Glyma.16G215100), but also harbors a down-regulated MLO gene (Glyma.16G145600). Accordingly, this study exemplified the feasibility of BSR-Seq in studying biotrophic disease resistance in soybean, and showed the genetic makeup of soybean variety 'KS11' comprising the Rmd locus and one MLO gene.

Characterization of SsHog1 and Shk1 Using Efficient Gene Knockout Systems through Repeated Protoplasting and CRISPR/Cas9 Ribonucleoprotein Approaches in Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is the causal agent of sclerotinia stem rot in over 400 plant species. In a previous study, the group III histidine kinase gene of S. sclerotiorum (Shk1) revealed its involvement in iprodione and fludioxonil sensitivity and osmotic stress. To further investigate the fungicide sensitivity associated with the high-osmolarity glycerol (HOG) pathway, we functionally characterized SsHog1, which is the downstream kinase of Shk1. To generate knockout mutants, split marker transformation combined with a newly developed repeated protoplasting method and CRISPR/Cas9 ribonucleoprotein (RNP) delivery approach were used. The pure SsHog1 and Shk1 knockout mutants showed reduced sensitivity to fungicides and increased sensitivity to osmotic stress. In addition, the SsHog1 knockout mutants demonstrated reduced virulence compared to Shk1 knockout mutants and wild-type. Our results indicate that the repeated protoplasting method and RNP approach can generate genetically pure homokaryotic mutants and SsHog1 is involved in osmotic adaptation, fungicide sensitivity, and virulence in S. sclerotiorum.

Colonization compatibility with Bacillus altitudinis confers soybean seed rot resistance

The plant microbiome and plant-associated bacteria are known to support plant health, but there are limited studies on seed and seedling microbiome to reveal how seed-associated bacteria may confer disease resistance. In this study, the application of antibiotics on soybean seedlings indicated that seed-associated bacteria were involved in the seed rot resistance against a soil-borne pathogen Calonectria ilicicola, but this resistance cannot be carried to withstand root rot. Using PacBio 16S rRNA gene full-length sequencing and microbiome analyses, 14 amplicon sequence variants (ASVs) including 2 ASVs matching to Bacillus altitudinis were found to be more abundant in the four most resistant varieties versus the four most susceptible varieties. Culture-dependent isolation obtained two B. altitudinis isolates that both exhibit antagonistic capability against six fungal pathogens. Application of B. altitudinis on the most resistant and susceptible soybean varieties revealed different colonization compatibility, and the seed rot resistance was restored in the five varieties showing higher bacterial colonization. Moreover, quantitative PCR confirmed the persistence of B. altitudinis on apical shoots till 21 days post-inoculation (dpi), but 9 dpi on roots of the resistant variety TN5. As for the susceptible variety HC, the persistence of B. altitudinis was only detected before 6 dpi on both shoots and roots. The short-term colonization of B. altitudinis on roots may explain the absence of root rot resistance. Collectively, this study advances the insight of B. altitudinis conferring soybean seed rot resistance and highlights the importance of considering bacterial compatibility with plant varieties and colonization persistence on plant tissues.

A High-Throughput Phenotyping Tool to Identify Field-Relevant Anthracnose Resistance in White Lupin

The seed- and air-borne pathogen Colletotrichum lupini, the causal agent of lupin anthracnose, is the most important disease in white lupin (Lupinus albus) worldwide and can cause total yield loss. The aims of this study were to establish a reliable high-throughput phenotyping tool to identify anthracnose resistance in white lupin germplasm and to evaluate a genomic prediction model, accounting for previously reported resistance quantitative trait loci, on a set of independent lupin genotypes. Phenotyping under controlled conditions, performing stem inoculation on seedlings, showed to be applicable for high throughput, and its disease score strongly correlated with field plot disease assessments (r = 0.95, P < 0.0001) and yield (r = -0.64, P = 0.035). Traditional one-row field disease phenotyping showed no significant correlation with field plot disease assessments (r = 0.31, P = 0.34) and yield (r = -0.45, P = 0.17). Genomically predicted resistance values showed no correlation with values observed under controlled or field conditions, and the parental lines of the recombinant inbred line population used for constructing the prediction model exhibited a resistance pattern opposite to that displayed in the original (Australian) environment used for model construction. Differing environmental conditions, inoculation procedures, or population structure may account for this result. Phenotyping a diverse set of 40 white lupin accessions under controlled conditions revealed eight accessions with improved resistance to anthracnose. The standardized area under the disease progress curves (sAUDPC) ranged from 2.1 to 2.8, compared with the susceptible reference accession with a sAUDPC of 3.85. These accessions can be incorporated into white lupin breeding programs. In conclusion, our data support stem inoculation-based disease phenotyping under controlled conditions as a time-effective approach to identify field-relevant resistance, which can now be applied to further identify sources of resistance and their underlying genetics.

Analysis and Identification of QTL for Resistance to Sclerotinia sclerotiorum in Pea (Pisum sativum L.)

White mold caused by Sclerotinia sclerotiorum is an important constraint to field pea (Pisum sativum L.) production worldwide. To transfer white mold resistance into an adapted background, and study the genetics of the disease, two recombinant inbred line (RIL) populations (PRIL17 and PRIL19) were developed by crossing two partially resistant plant introductions with two susceptible pea cultivars. PRIL17 (Lifter × PI240515), and PRIL19 (PI169603 × Medora) were evaluated for resistance to white mold by measuring lesion expansion inhibition (LEI) and nodal transmission inhibition (NTI) at 3, 7, and 14 days post inoculation (dpi) under controlled environmental conditions. Lesion expansion inhibition percentage (LEIP), survival rate (SR), and area under disease progress curves (AUDPC) were also calculated accordingly. Because of a positive correlation between LEI and NTI with height, short and long internode individuals of each population were analyzed separately to avoid any confounding effect of height to pathogen response. A total of 22 short genotypes demonstrated partial resistance based on at least two Porter's resistance criteria. Only two pea genotypes with partial resistance to white mold (PRIL19-18 and PRIL19-124) had both semi-leafless (afila) and short internode traits. Both the RIL populations were genotyped using genotyping by sequencing (GBS). For PRIL17 and PRIL19, genetic maps were constructed from a total of 1,967 and 1,196 single nucleotide polymorphism (SNP) and spanned over 1,494 cM and 1,415 cM representing seven and nine linkage groups, respectively. A consensus map constructed using data from both populations, had 1,486 unique SNPs over 2,461 cM belonging to seven linkage groups. Inclusive composite interval mapping (ICIM) identified thirteen quantitative trait loci (QTL) associated with white mold resistance traits in both populations. Three of them were co-located with height genes (a morphological trait that reduces infection risk and acts as disease avoidance) and the other ten QTL were associated with two forms of physiological resistance (seven for LEI and three for NTI) with LOD and r2 ranging from 3.0 to 28.5 and 5.1 to 64.3, respectively. The development of resistance lines, genetic dissection and identification of markers associated will help accelerate breeding efforts for white mold resistance using molecular breeding approaches.

Key Genes and Genetic Interactions of Plant-Pathogen Functional Modules in Poplar Infected by Marssonina brunnea

Marssonina brunnea, the causative pathogen of Marssonina leaf spot of poplars (MLSP), devastates poplar plantations by forming black spots on leaves and defoliating trees. Although MLSP has been studied for over 30 years, the key genes that function during M. brunnea infection and their effects on plant growth are poorly understood. Here, we used multigene association studies to investigate the effects of key genes in the plant-pathogen interaction pathway, as revealed by transcriptome analysis, on photosynthesis and growth in a natural population of 435 Populus tomentosa individuals. By analyzing transcriptomic changes during three stages of infection, we detected 628 transcription factor genes among the 7,611 differentially expressed genes that might be associated with basal defense responses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that transcriptomic changes across different stages of infection lead to the reprogramming of metabolic processes possibly related to defense activation. We identified 29,399 common single-nucleotide polymorphisms (SNPs) within 221 full-length genes in plant-pathogen interaction pathways that were significantly associated with photosynthetic and growth traits. We also detected 4,460 significant epistatic pairs associated with stomatal conductance, tree diameter, and tree height. Epistasis analysis uncovered significant interactions between 2,561 SNP-SNP pairs from different functional modules in the plant-pathogen interaction pathway, revealing possible genetic interactions. This analysis revealed many key genes that function during M. brunnea infection and their potential roles in mediating photosynthesis and plant growth, shedding light on genetic interactions between functional modules in the plant-pathogen interaction pathway.

A Sclerotinia sclerotiorum Transcription Factor Involved in Sclerotial Development and Virulence on Pea

Sclerotinia sclerotiorum is a plant-pathogenic ascomycete fungus and infects over 400 host plants, including pea (Pisum sativum L.). The fungus causes white mold on pea, and substantial yield loss is attributed to the disease. To improve white mold management, further understanding of S. sclerotiorum pathogenicity is crucial. In this study, 389 transcription factors (TFs) were mined from the complete genome sequence of S. sclerotiorum and their in planta expression patterns were determined in susceptible and partially resistant pea lines and compared to in vitro expression patterns on culture medium. One of the transcription factors was significantly induced in planta at 24 and 48 h postinfection compared to the expression in vitro This putative C6 transcription factor of S. sclerotiorum (SsC6TF1) was knocked down using a gene-silencing approach to investigate its functions in vegetative growth and sclerotial development as well as its virulence and pathogenicity in pea. While the SsC6TF1 knockdown mutants had hyphal growth rates identical to those of the wild-type strain and were capable of infection, the knockdown mutants produced no sclerotia or significantly fewer and smaller sclerotia on the culture medium and exhibited reduced virulence on both pea lines. This study profiled genome-wide expression for S. sclerotiorum transcription factors in planta and in vitro and functionally characterized a novel transcription factor, SsC6TF1, which positively regulates sclerotial development and virulence on pea. The finding provides molecular insights into S. sclerotiorum biology and interaction with pea and other economically important crops.IMPORTANCE White mold, caused by Sclerotinia sclerotiorum, is a destructive disease on important legume species such as soybean, dry bean, and pea. This study investigated expression levels of transcription factors in S. sclerotiorumin planta (pea lines) and in vitro (culture medium). One transcription factor displaying high expression in planta was found to be involved in sclerotial development and virulence on pea. This report provides a new understanding regarding transcription factors of S. sclerotiorum in development and virulence.