Advanced backcross QTL analysis and comparative mapping with RIL QTL studies and GWAS provide an overview of QTL and marker haplotype diversity for resistance to Aphanomyces root rot in pea (Pisum sativum)

Aphanomyces euteiches is the most damaging soilborne pea pathogen in France. Breeding of pea resistant varieties combining a diversity of quantitative trait loci (QTL) is a promising strategy considering previous research achievements in dissecting polygenic resistance to A. euteiches. The objective of this study was to provide an overview of the diversity of QTL and marker haplotypes for resistance to A. euteiches, by integrating a novel QTL mapping study in advanced backcross (AB) populations with previous QTL analyses and genome-wide association study (GWAS) using common markers. QTL analysis was performed in two AB populations derived from the cross between the susceptible spring pea variety "Eden" and the two new sources of partial resistance "E11" and "LISA". The two AB populations were genotyped using 993 and 478 single nucleotide polymorphism (SNP) markers, respectively, and phenotyped for resistance to A. euteiches in controlled conditions and in infested fields at two locations. GWAS and QTL mapping previously reported in the pea-Aphanomyces collection and from four recombinant inbred line (RIL) populations, respectively, were updated using a total of 1,850 additional markers, including the markers used in the Eden x E11 and Eden x LISA populations analysis. A total of 29 resistance-associated SNPs and 171 resistance QTL were identified by GWAS and RIL or AB QTL analyses, respectively, which highlighted 10 consistent genetic regions confirming the previously reported QTL. No new consistent resistance QTL was detected from both Eden x E11 and Eden x LISA AB populations. However, a high diversity of resistance haplotypes was identified at 11 linkage disequilibrium (LD) blocks underlying consistent genetic regions, especially in 14 new sources of resistance from the pea-Aphanomyces collection. An accumulation of favorable haplotypes at these 11 blocks was confirmed in the most resistant pea lines of the collection. This study provides new SNP markers and rare haplotypes associated with the diversity of Aphanomyces root rot resistance QTL investigated, which will be useful for QTL pyramiding strategies to increase resistance levels in future pea varieties.

Effects of nitrogen fertilization and a commercial arbuscular mycorrhizal fungal inoculant on root rot and agronomic production of pea and lentil crops

In the Canadian prairies, pulse crops such as field pea (Pisum sativum L.) and lentil (Lens culinaris L.) are economically important and widely grown. However, in recent years, root rot, caused by a variety of fungal and oomycete pathogens, including Aphanomyces euteiches, has become a limiting factor on yield. In this study, we examined the impacts of nitrogen (N) fertilization and a commercial arbuscular mycorrhizal fungal (AMF) inoculant on pea and lentil plant health and agronomic production at three locations in Saskatchewan: Swift Current, Indian Head and Melfort. The AMF inoculation had no impact on root rot severity, and therefore is not considered a reliable method to manage root rot in pea and lentil. In contrast, N fertilization led to reductions in root rot in Swift Current, but not the other two sites. However, N fertilization did reduce nodulation. When both pea and lentil are considered, the abundance of A. euteiches in soil increased from pre-seeding to mid-bloom. A negative correlation between soil pH and disease severity was also observed. The high between-site variability highlights the importance of testing root rot mitigation strategies under multiple soil conditions to develop site-specific recommendations. Use of N fertilizer as a root rot management strategy merits further exploration, including investigation into its interactions with other management strategies, soil properties, and costs and benefits.

Evaluation of pea genotype PI180693 partial resistance towards aphanomyces root rot in commercial pea breeding

The cultivation of vining pea (Pisum sativum) faces a major constraint with root rot diseases, caused by a complex of soil-borne pathogens including the oomycetes Aphanomyces euteiches and Phytophtora pisi. Disease resistant commercial varieties are lacking but the landrace PI180693 is used as a source of partial resistance in ongoing pea breeding programs. In this study, the level of resistance and their interaction with A. euteiches virulence levels of six new back-crossed pea breeding lines, deriving from the cross between the susceptible commercial cultivar Linnea and PI180693, were evaluated for their resistance towards aphanomyces root rot in growth chamber and green house tests. Resistance towards mixed infections by A. euteiches and P. pisi and commercial production traits were evaluated in field trials. In growth chamber trials, pathogen virulence levels had a significant effect on plant resistance, as resistance was more consistent against A. euteiches strains exhibiting high or intermediate virulence compared with lowly virulent strains. In fact, line Z1701-1 showed to be significantly more resistant than both parents when inoculated with a lowly virulent strain. In two separate field trials in 2020, all six breeding lines performed equally well as the resistant parent PI180693 at sites only containing A. euteiches, as there were no differences in disease index. In mixed infections, PI180693 exhibited significantly lower disease index scores than Linnea. However, breeding lines displayed higher disease index scores compared with PI180693, indicating higher susceptibility towards P. pisi. Data on seedling emergence from the same field trials suggested that PI180693 was particularly sensitive towards seed decay/damping off disease caused by P. pisi. Furthermore, the breeding lines performed equally well as Linnea in traits important for green pea production, again emphasizing the commercial potential. In summary, we show that the resistance from PI180693 interacts with virulence levels of the pathogen A. euteiches and is less effective towards root rot caused by P. pisi. Our results show the potential use of combining PI180693 partial resistance against aphanomyces root rot with commercially favorable breeding traits in commercial breeding programs.

A fine-tuned defense at the pea root caps: Involvement of border cells and arabinogalactan proteins against soilborne diseases

Plants have to cope with a myriad of soilborne pathogens that affect crop production and food security. The complex interactions between the root system and microorganisms are determinant for the whole plant health. However, the knowledge regarding root defense responses is limited as compared to the aerial parts of the plant. Immune responses in roots appear to be tissue-specific suggesting a compartmentalization of defense mechanisms in these organs. The root cap releases cells termed root "associated cap-derived cells" (AC-DCs) or "border cells" embedded in a thick mucilage layer forming the root extracellular trap (RET) dedicated to root protection against soilborne pathogens. Pea (Pisum sativum) is the plant model used to characterize the composition of the RET and to unravel its function in root defense. The objective of this paper is to review modes of action of the RET from pea against diverse pathogens with a special focus on root rot disease caused by Aphanomyces euteiches, one of the most widely occurring and large-scale pea crop diseases. The RET, at the interface between the soil and the root, is enriched in antimicrobial compounds including defense-related proteins, secondary metabolites, and glycan-containing molecules. More especially arabinogalactan proteins (AGPs), a family of plant extracellular proteoglycans belonging to the hydroxyproline-rich glycoproteins were found to be particularly present in pea border cells and mucilage. Herein, we discuss the role of RET and AGPs in the interaction between roots and microorganisms and future potential developments for pea crop protection.

Identification of Novel Genes Associated with Partial Resistance to Aphanomyces Root Rot in Field Pea by BSR-Seq Analysis

Aphanomyces root rot, caused by Aphanomyces euteiches, causes severe yield loss in field pea (Pisum sativum). The identification of a pea germplasm resistant to this disease is an important breeding objective. Polygenetic resistance has been reported in the field pea cultivar '00-2067'. To facilitate marker-assisted selection (MAS), bulked segregant RNA-seq (BSR-seq) analysis was conducted using an F8 RIL population derived from the cross of 'Carman' × '00-2067'. Root rot development was assessed under controlled conditions in replicated experiments. Resistant (R) and susceptible (S) bulks were constructed based on the root rot severity in a greenhouse study. The BSR-seq analysis of the R bulks generated 44,595,510~51,658,688 reads, of which the aligned sequences were linked to 44,757 genes in a reference genome. In total, 2356 differentially expressed genes were identified, of which 44 were used for gene annotation, including defense-related pathways (jasmonate, ethylene and salicylate) and the GO biological process. A total of 344.1 K SNPs were identified between the R and S bulks, of which 395 variants were located in 31 candidate genes. The identification of novel genes associated with partial resistance to Aphanomyces root rot in field pea by BSR-seq may facilitate efforts to improve management of this important disease.

Faba bean root exudates alter pea root colonization by the oomycete Aphanomyces euteiches at early stages of infection

Aphanomyces euteiches is an oomycete pathogen that causes the pea root rot. We investigated the potential role of early belowground defense in pea (susceptible plant) and faba bean (tolerant plant) at three days after inoculation. Pea and faba bean were inoculated with A. euteiches zoospores. Root colonization was examined. Root exudates from pea and faba bean were harvested and their impact on A. euteiches development were assessed by using in vitro assays. A. euteiches root colonization and the influence of the oomycete inoculation on specialized metabolites patterns and arabinogalactan protein (AGP) concentration of root exudates were also determined. In faba bean root, A. euteiches colonization was very low as compared with that of pea. Whereas infected pea root exudates have a positive chemotaxis index (CI) on zoospores, faba bean exudate CI was negative suggesting a repellent effect. While furanoacetylenic compounds were only detected in faba bean exudates, AGP concentration was specifically increased in pea.This work showed that early in the course of infection, host susceptibility to A. euteiches is involved via a plant-species specific root exudation opening new perspectives in pea root rot disease management.

Untangling the Pea Root Rot Complex Reveals Microbial Markers for Plant Health

Plant health is recognised as a key element to ensure global food security. While plant breeding has substantially improved crop resistance against individual pathogens, it showed limited success for diseases caused by the interaction of multiple pathogens such as root rot in pea (Pisum sativum L.). To untangle the causal agents of the pea root rot complex and determine the role of the plant genotype in shaping its own detrimental or beneficial microbiome, fungal and oomycete root rot pathogens, as well as previously identified beneficials, i.e., arbuscular mycorrhizal fungi (AMF) and Clonostachys rosea, were qPCR quantified in diseased roots of eight differently resistant pea genotypes grown in four agricultural soils under controlled conditions. We found that soil and pea genotype significantly determined the microbial compositions in diseased pea roots. Despite significant genotype x soil interactions and distinct soil-dependent pathogen complexes, our data revealed key microbial taxa that were associated with plant fitness. Our study indicates the potential of fungal and oomycete markers for plant health and serves as a precedent for other complex plant pathosystems. Such microbial markers can be used to complement plant phenotype- and genotype-based selection strategies to improve disease resistance in one of the world's most important pulse crops of the world.

Effects of Lentil Genotype on the Colonization of Beneficial Trichoderma Species and Biocontrol of Aphanomyces Root Rot

Trichoderma species are opportunistic plant symbionts that are common in the root and rhizosphere ecosystems. Many Trichoderma species may enhance plant growth, nutrient acquisition, and disease resistance, and for these reasons, they are widely used in agriculture as biofertilizers or biocontrol agents. Host plant genotype and other microorganisms, such as root pathogens, may influence the efficacy of Trichoderma inoculants. Aphanomyces euteiches is an important soil-borne oomycete in western Canada that causes root rot in legume crops such as lentil and pea, and there is not yet any significantly resistant varieties or effective treatments available to control the disease. In this study, the composition of root-associated fungal communities and the abundance of Trichoderma species, T. harzianum strain T-22 and T. virens strain G41, was determined in the roots of eight Lens genotypes based on internal transcribed spacer (ITS) Illumina MiSeq paired-end sequencing, both in the presence and the absence of the root rot pathogen Aphanomyces euteiches. Biocontrol effects of T. harzianum on A. euteiches was also examined. Significant genotypic variations were observed in the composition of root-associated fungal communities and the abundance of the different Trichoderma species in the lentil roots. The presence of A. euteiches altered the composition of Trichoderma found associated to the lentil genotypes. Biocontrol of A. euteiches by T. harzianum T22 species was observed in vitro and positive correlations between the abundance of Trichoderma and plant root and shoot biomass were observed in vivo. These findings revealed that lentil genotype and infection by the phytopathogen A. euteiches greatly influenced the colonization of root-associated fungi and the abundance of the Trichoderma species, as well as the effect on plant growth promotion. The multipartite interactions observed among lentil genotypes, Trichoderma species and A. euteiches suggest possibilities to select compatible host-beneficial microbe combinations in lentil breeding programs and to develop application strategies to harness the beneficial effects of Trichoderma inoculants in sustainable crop production systems.

Dissecting the Genetic Architecture of Aphanomyces Root Rot Resistance in Lentil by QTL Mapping and Genome-Wide Association Study

Lentil (Lens culinaris Medikus) is an important source of protein for people in developing countries. Aphanomyces root rot (ARR) has emerged as one of the most devastating diseases affecting lentil production. In this study, we applied two complementary quantitative trait loci (QTL) analysis approaches to unravel the genetic architecture underlying this complex trait. A recombinant inbred line (RIL) population and an association mapping population were genotyped using genotyping by sequencing (GBS) to discover novel single nucleotide polymorphisms (SNPs). QTL mapping identified 19 QTL associated with ARR resistance, while association mapping detected 38 QTL and highlighted accumulation of favorable haplotypes in most of the resistant accessions. Seven QTL clusters were discovered on six chromosomes, and 15 putative genes were identified within the QTL clusters. To validate QTL mapping and genome-wide association study (GWAS) results, expression analysis of five selected genes was conducted on partially resistant and susceptible accessions. Three of the genes were differentially expressed at early stages of infection, two of which may be associated with ARR resistance. Our findings provide valuable insight into the genetic control of ARR, and genetic and genomic resources developed here can be used to accelerate development of lentil cultivars with high levels of partial resistance to ARR.

Polyphenolic Composition of Lentil Roots in Response to Infection by Aphanomyces euteiches

Polyphenols comprise the largest group of plant secondary metabolites and have critical roles in plant physiology and response to the biotic and abiotic environment. Changes in the content of polyphenols in the root extracts and root tissues of wild (Lens ervoides) and cultivated (Lens culinaris) lentil genotypes were examined in response to infection by Aphanomyces euteiches using liquid chromatography mass spectrometry (LC-MS). Genotype, infection and their interaction determined the composition of polyphenols in lentil roots. The levels of several polyphenols were lower in the root extract of the low-tannin genotype L. culinaris ZT-4 compared to L. ervoides L01-827A. Kaempferol derivatives including kaempferol dirutinoside and kaempferol 3-robinoside 7-rhamnoside were more concentrated in the healthy root tissues of L. ervoides L01-827A than in L. culinaris genotypes. Infection increased the concentration of kaempferol, apigenin, and naringenin in the root tissues of all genotypes, but had no effect on some polyphenols in the low-tannin genotype L. culinaris ZT-4. The concentrations of apigenin, naringenin, apigenin 4-glucoside, naringenin7-rutinoside, diosmetin, and hesperetin 7-rutinoside were higher in the infected root tissues of L. ervoides L01-827A compared with the L. culinaris genotypes. Organic acids including coumaric acid, vanillic acid, 4-aminosalicylic acid, 4-hydroxybenzoic acid, and 3,4-dihydroxybenzoic acid effectively suppressed the in-vitro hyphal growth of A. euteiches. Some of these bioactive polyphenols were more concentrated in roots of L. ervoides L01-827A but were low to undetectable in ZT-4. This study shows that genotypic differences exist in the composition of root polyphenols in lentil, and is related to the response to infection caused by A. euteiches. Polyphenols, particularly the organic acid content could be useful for selection and breeding of lentil genotypes that are resistant to Aphanomyces root rot (ARR) disease.

Is there genetic variation in mycorrhization of Medicago truncatula?

Differences in the plant's response among ecotypes or accessions are often used to identify molecular markers for the respective process. In order to analyze genetic diversity of Medicago truncatula in respect to interaction with the arbuscular mycorrhizal (AM) fungus Rhizophagus irregularis, mycorrhizal colonization was evaluated in 32 lines of the nested core collection representing the genetic diversity of the SARDI collection. All studied lines and the reference line Jemalong A17 were inoculated with R. irregularis and the mycorrhization rate was determined at three time points after inoculation. There were, however, no reliable and consistent differences in mycorrhization rates among all lines. To circumvent possible overlay of potential differences by use of the highly effective inoculum, native sandy soil was used in an independent experiment. Here, significant differences in mycorrhization rates among few of the lines were detectable, but the overall high variability in the mycorrhization rate hindered clear conclusions. To narrow down the number of lines to be tested in more detail, root system architecture (RSA) of in vitro-grown seedlings of all lines under two different phosphate (Pi) supply condition was determined in terms of primary root length and number of lateral roots. Under high Pi supply (100 µM), only minor differences were observed, whereas in response to Pi-limitation (3 µM) several lines exhibited a drastically changed number of lateral roots. Five lines showing the highest alterations or deviations in RSA were selected and inoculated with R. irregularis using two different Pi-fertilization regimes with either 13 mM or 3 mM Pi. Mycorrhization rate of these lines was checked in detail by molecular markers, such as transcript levels of RiTubulin and MtPT4. Under high phosphate supply, the ecotypes L000368 and L000555 exhibited slightly increased fungal colonization and more functional arbuscules, respectively. To address the question, whether capability for mycorrhizal colonization might be correlated to general invasion by microorganisms, selected lines were checked for infection by the root rot causing pathogen, Aphanoymces euteiches. The mycorrhizal colonization phenotype, however, did not correlate with the resistance phenotype upon infection with two strains of A. euteiches as L000368 showed partial resistance and L000555 exhibited high susceptibility as determined by quantification of A. euteiches rRNA within infected roots. Although there is genetic diversity in respect to pathogen infection, genetic diversity in mycorrhizal colonization of M. truncatula is rather low and it will be rather difficult to use it as a trait to access genetic markers.

Degenerate ITS7 primer enhances oomycete community coverage and PCR sensitivity to Aphanomyces species, economically important plant pathogens

Metagenomic analysis of oomycetes through deep amplicon sequencing has been conducted primarily using the ITS6-ITS7 primer set that targets the ITS1 region. While this primer set shows a perfect match to most oomycete taxa, ITS7 contains 3 mismatches to the corresponding binding site of plant pathogens within the genus Aphanomyces. Polymerase chain reaction (PCR) efficiency differs for taxa with uneven primer matching characteristics, which may explain why previous studies have detected this genus at low abundance. To overcome the impact of these mismatches on PCR sensitivity, the mismatched nucleotides were replaced with degenerate nucleotides. Oomycete communities from 35 soil samples collected from asymptomatic and root rot diseased sites in pea fields across Alberta were analyzed simultaneously using ITS6-ITS7 and ITS6-ITS7-a.e. (modified version of ITS7) primer sets on 1 Illumina MiSeq run. The number of high-quality reads obtained by ITS6-ITS7-a.e. was more than twice that of ITS6-ITS7. The relative abundance of Pythium spp. was reduced and Aphanomyces spp. increased. Aphanomyces cf. cladogamus and Aphanomyces euteiches were the second and third most abundant species, respectively, in the pea rhizosphere using the ITS7-a.e. primer, but were rare using the ITS7 primer. These results indicate that use of ITS7-a.e. provides a more accurate picture of oomycete communities than ITS7 by enhancing PCR sensitivity to Aphanomyces.

Systemic Induction of the Defensin and Phytoalexin Pisatin Pathways in Pea (Pisum sativum) against Aphanomyces euteiches by Acetylated and Nonacetylated Oligogalacturonides

Oligogalacturonides (OGs) are known for their powerful ability to stimulate the plant immune system but little is known about their mode of action in pea (Pisum sativum). In the present study, we investigated the elicitor activity of two fractions of OGs, with polymerization degrees (DPs) of 2-25, in pea against Aphanomyces euteiches. One fraction was nonacetylated (OGs - Ac) whereas the second one was 30% acetylated (OGs + Ac). OGs were applied by injecting the upper two rachises of the plants at three- and/or four-weeks-old. Five-week-old roots were inoculated with 10⁵ zoospores of A. euteiches. The root infection level was determined at 7, 10 and 14 days after inoculation using the quantitative real-time polymerase chain reaction (qPCR). Results showed significant root infection reductions namely 58, 45 and 48% in the plants treated with 80 µg OGs + Ac and 59, 56 and 65% with 200 µg of OGs - Ac. Gene expression results showed the upregulation of genes involved in the antifungal defensins, lignans and the phytoalexin pisatin pathways and a priming effect in the basal defense, SA and ROS gene markers as a response to OGs. The reduction of the efficient dose in OGs + Ac is suggesting that acetylation is necessary for some specific responses. Our work provides the first evidence for the potential of OGs in the defense induction in pea against Aphanomyces root rot.

Cross-Regulation between N Metabolism and Nitric Oxide (NO) Signaling during Plant Immunity

Plants are sessile organisms that have evolved a complex immune system which helps them cope with pathogen attacks. However, the capacity of a plant to mobilize different defense responses is strongly affected by its physiological status. Nitrogen (N) is a major nutrient that can play an important role in plant immunity by increasing or decreasing plant resistance to pathogens. Although no general rule can be drawn about the effect of N availability and quality on the fate of plant/pathogen interactions, plants' capacity to acquire, assimilate, allocate N, and maintain amino acid homeostasis appears to partly mediate the effects of N on plant defense. Nitric oxide (NO), one of the products of N metabolism, plays an important role in plant immunity signaling. NO is generated in part through Nitrate Reductase (NR), a key enzyme involved in nitrate assimilation, and its production depends on levels of nitrate/nitrite, NR substrate/product, as well as on L-arginine and polyamine levels. Cross-regulation between NO signaling and N supply/metabolism has been evidenced. NO production can be affected by N supply, and conversely NO appears to regulate nitrate transport and assimilation. Based on this knowledge, we hypothesized that N availability partly controls plant resistance to pathogens by controlling NO homeostasis. Using the Medicago truncatula/Aphanomyces euteiches pathosystem, we showed that NO homeostasis is important for resistance to this oomycete and that N availability impacts NO homeostasis by affecting S-nitrosothiol (SNO) levels and S-nitrosoglutathione reductase activity in roots. These results could therefore explain the increased resistance we noted in N-deprived as compared to N-replete M. truncatula seedlings. They open onto new perspectives for the studies of N/plant defense interactions.

MtNF-YA1, A Central Transcriptional Regulator of Symbiotic Nodule Development, Is Also a Determinant of Medicago truncatula Susceptibility toward a Root Pathogen

Plant NF-Y transcription factors control a wide array of biological functions enabling appropriate reproductive and developmental processes as well as adaptation to various abiotic and biotic environments. In Medicago truncatula, MtNF-YA1 was previously identified as a key determinant for nodule development and establishment of rhizobial symbiosis. Here, we highlight a new role for this protein in compatibility to Aphanomyces euteiches, a root pathogenic oomycete. The Mtnf-ya1-1 mutant plants showed better survival rate, reduced symptoms, and increased development of their root apparatus as compared to their wild-type (WT) background A17. MtNF-YA-1 was specifically up-regulated by A. euteiches in F83005.5, a highly susceptible natural accession of M. truncatula while transcript level remained stable in A17, which is partially resistant. The role of MtNF-YA1 in F83005.5 susceptibility was further documented by reducing MtNF-YA1 expression either by overexpression of the miR169q, a microRNA targeting MtNF-YA1, or by RNAi approaches leading to a strong enhancement in the resistance of this susceptible line. Comparative analysis of the transcriptome of WT and Mtnf-ya1-1 led to the identification of 1509 differentially expressed genes. Among those, almost 36 defense-related genes were constitutively expressed in Mtnf-ya1-1, while 20 genes linked to hormonal pathways were repressed. In summary, we revealed an unexpected dual role for this symbiotic transcription factor as a key player in the compatibility mechanisms to a pathogen.

High-density genome-wide association mapping implicates an F-box encoding gene in Medicago truncatula resistance to Aphanomyces euteiches

• The use of quantitative disease resistance (QDR) is a promising strategy for promoting durable resistance to plant pathogens, but genes involved in QDR are largely unknown. To identify genetic components and accelerate improvement of QDR in legumes to the root pathogen Aphanomyces euteiches, we took advantage of both the recently generated massive genomic data for Medicago truncatula and natural variation of this model legume. • A high-density (≈5.1 million single nucleotide polymorphisms (SNPs)) genome-wide association study (GWAS) was performed with both in vitro and glasshouse phenotyping data collected for 179 lines. • GWAS identified several candidate genes and pinpointed two independent major loci on the top of chromosome 3 that were detected in both phenotyping methods. Candidate SNPs in the most significant locus (σ(A)²= 23%) were in the promoter and coding regions of an F-box protein coding gene. Subsequent qRT-PCR and bioinformatic analyses performed on 20 lines demonstrated that resistance is associated with mutations directly affecting the interaction domain of the F-box protein rather than gene expression. • These results refine the position of previously identified QTL to specific candidate genes, suggest potential molecular mechanisms, and identify new loci explaining QDR against A. euteiches.