Genome-wide-association studies (GWAS) on the interaction between the cyst nematode <em>Heterodera schachtii </em>and<em> Arabidopsis thaliana</em>

Abstract

Plant-parasitic cyst nematodes are obligate sedentary parasites that infect the roots of a wide range of host plants. In research, the beet cyst nematode Heterodera schachtii infecting Arabidopsis thaliana is often used as a model. H. schachtii is sexually dimorphic and differentiation into male and female is strongly influenced by interactions with the environment. Under favorable conditions the unsexed infective juveniles (J2) tend to develop into females, whereas they rather form males under adverse conditions. After hatching from the egg, J2 enters the root tissue to find out a suitable infection site. Once a feeding site is selected, they become immobile and release specific proteinaceous effectors triggering the formation a syncytial nurse cell system that serves as the only nutrient supply throughout the nematode life cycle. To identify host loci underlying variation in susceptibility, we performed a genome-wide association study (GWAS) using different accessions of Arabidopsis that were exposed to H. schachtii infections. Three different susceptibility parameters were examined: average numbers of female per plant, average numbers of male per plant and their subsequent female to male sex ratio. Based on female- to -male sex ratio, we categorized all accessions into three groups: lowly susceptible, moderately susceptible and highly susceptible. GWAS revealed a novel QTL allele at chromosome four underlying potential candidates associated to sex ratio variations. Using the physical map of this locus, we selected five candidate genes from the potential QTL: GDSL-lipase (At4g18970), PPR-protein (pentatricopeptide repeat-protein; At4g18975), AtS40-3 (senescence-associated protein; At4g18980), XTH-29 (xyloglucan endotransglucosylase 29; At4g18990) and IWS-2 (interacts with SPT6; Atg19000). Finally, extreme bulks showing lowest and highest sex ratio were selected for detailed functional characterization of candidate genes. At first, we measured the expression of candidate genes among four extreme accessions showing lowest and highest sex ratio. This analysis revealed that out of the five, expression of two genes, PPR and AtS40-3 was strongly down-regulated in lowly susceptible accessions compared with highly susceptible accessions. In comparison, expression of GDSL-lipase, iWS2, and XTH29 was unchanged among all tested accessions. Together, these data suggest that changes in expression level of PPR and AtS40-3 might be correlated with variation in nematode sex ratio. AtS40-3 and PPR are divergent genes that are organized head-to-head (in opposite orientations) in a non-overlapping manner. Because of their specific regulation, we analyzed AtS40-3 and PPR for sequence variations in their common promoter region. Notably, we observed a consistent deletion in the promoter in lowly susceptible accessions, whereas no such deletions were present in highly susceptible accessions. Based on these data, we hypothesized that the deletions in the promoter region might be one of the reasons for their strong down-regulation in root tissues in lowly susceptible accessions. To further investigate this hypothesis, we performed an in silico analysis of the commonly deleted 28 bp promoter region from susceptible line Col-0 for the occurrence of important cis-acting element. Our analysis revealed the presence of some prominent transcription factor binding sequences including a typical TATA box. To understand the mechanism that regulates the differential promoter activity for AtS40-3 between lowly and highly susceptible lines, the 1500 bp promoter upstream of ATG for AtS40-3 from the low susceptibility accession Xan-1 (carrying a deletion of 28 nt) and the susceptible accession Col-0 (without deletion) was cloned and used to drive the expression of GFP. The expression of GFP was analyzed via a transient expression system in N. benthamiana through real-time qPCR. Indeed, we found a significantly reduced expression of GFP driven by the promoter cloned from the Xan-1 accession (low susceptibility to nematodes) compared with that from Col-0 (high susceptibility to nematodes). We obtained the T-DNA-inserted loss of function mutant lines of AtS40-3 and characterized these lines via nematode infection assays. Remarkably, a prominent increase in average number of males and decrease in average number of females was observed. Consequently, female-to-male sex ratio was also changed in Xan-1 compared with Col-0. Taken together, this study provides insights into genes and factors that influence the sexual variation in cyst nematode. Further exploitation of these factors may lead to development of novel resistance resources against nematode in crop plants.

Pflanzenpararsitäre Zystennematoden sind obligat biotrophe Parasiten, die Wurzeln zahlreicher Pflanzen befallen. Zu Forschungszwecken wird oft der Rübenzystennematode Heterodera schachtii als Parasit von Arabidopsis thaliana als Modell verwendet. Diese Art weist einen deutlichen Sexualdimorphismus auf. Die Differenzierung der Geschlechter ist stark von Umweltinteraktionen beeinflusst. Unter günstigen Bedingungen entwickeln sich die undifferenzierten Infektionslarven (J2) eher zu weiblichen Tieren, während sich unter ungünstigen Bedingungen aus ihnen eher Männchen bilden. Nach dem Schlupf aus dem Ei dringen die J2 auf der Suche nach einer geeigneten Infektionsstelle in die Wurzel ein. Nachdem sie eine geeignete Stelle gefunden haben, werden sie immobil und geben Effektorproteine ab, die die Bildung eines syncytialen Nährzellensystems auslösen, das als einzige Nahrungsquelle der Nematoden dient. Um Wirtsloci zu identifizieren, die im Zusammenhang mit der Anfälligkeit gegenüber H. schachtii stehen, wurden im Rahmen dieser Arbeit so genannte Genomweite Assoziationsstudien (GWAS) auf der Basis verschiedener Arabidopsisherkünfte durchgeführt. Dabei wurden drei verschiedene Parameter untersucht: Die durchschnittliche Anzahl von Weibchen pro Pflanze, die durchschnittliche Anzahl von Männchen pro Pflanze und das daraus folgende Geschlechterverhältnis. Anhand des Geschlechterverhältnisses wurden alle Herkünfte in drei Gruppen eingeteilt: Gering anfällig, mäßig anfällig, und hoch anfällig. Die GWAS Analyse erbrachte neuartiges QTL Allel auf Chromosom vier das mit Variabilität des Geschlechterverhältnisses assoziiert war. Auf der Basis der Genkartierung dieses Lokus wurden fünf Kandidatengene identifiziert: GDSL-Lipase (At4g18970), PPR-Protein (pentatricopeptide repeat-protein; At4g18975), AtS40-3 (senescence-associated protein; At4g18980), XTH-29 (xyloglucan endotransglucosylase 29; At4g18990) und IWS-2 (interacts with SPT6; Atg19000). Außerdem wurden zwei Extremgruppen von Herkünften mit dem geringsten und dem höchsten Geschlechterverhältnis gebildet, um eine funktionelle Charakterisierung der Kandidatengene durchzuführen. Hierzu wurde die Expression der Kandidatengene in vier Extremherkünften mit hohem und niedrigem Geschlechterverhältnis untersucht. Diese Untersuchung ergab, dass die Gene PPR (At4g18975) und AtS40-3 (At4g18980) in gering anfälligen Herkünften stark abreguliert waren. Im Vergleich dazu war die Expression der Gene GDSL-Lipase (At4g18970), iWS2 (At4g19000), and XTH29 (At4g18990) in diesen Herkünften unverändert. Diese Ergebnisse führten zu der Annahme, dass unterschiedliche Expressionslevel von PPR und AtS40-3 mit der Ausprägung des Geschlechterverhältnisses korrelieren könnten. Eine Sequenzanalyse der Promotorregion von PPR und AtS40-3 ergab, dass beide Gene Kopf and Kopf in gegensätzlicher Ausrichtung ohne Überlappung lokalisiert sind. Wegen ihrer spezifischen Regulation wurden die Promotersequenzen in den Extremherkünften untersucht. Es zeigte sich, dass die Promotorsequenzen von AtS40-3 und PPRder gering anfälligen Herkünfte Deletionen von 7-51 bp aufwiesen, während in den hoch anfälligen keine Deletionen zu finden waren. Infolge dessen verfolgten wir die Hypothese, dass die Deletionen im Zusammenhang mit der starken Abregulierung der Gene in gering anfälligen Herkünften stehen. Folglich führten wir eine in silico Analyse der zentralen 28/27 bp Promotordeletion durch, um wichtige Cis-Elemente zu identifizieren. Analyse ergab, dass die Deletion Transkriptionsfaktor-bindende Elemente, wie zum Beispiel eine TATA Box enthält. Um den Regulationsmechanismus für AtS40-3 in gering und hoch anfälligen Linien zu untersuchen, wurde die 1500 bp umfassende Promotorregion von AtS40-3 der gering anfälligen Linie Xan-1 (mit einer 28 bp Deletion) und der anfälligen Linie Col-0 kloniert und zur transienten Expression von GFP in N. benthamiana verwendet. In Real-time qPCR-Analysen konnte eine signifikant geringere Expression mit dem Promotor von N. benthamiana Xan-1 gemessen werden. Infektionsversuche an einer T-DNA Mutationslinie von AtS40-3 ergaben eine ähnliche Verschiebung der Geschlechterverhältnises bei der Entwicklung der Nematoden zu Gunsten von Männchen, wie sie bei Xan-1 beobachtet wurde. Die Ergebnisse zeigen, in welcher Weise die unterschiedliche Expression von Genen die Entwicklung von Nematoden beinflussen kann. Diese Erkenntnisse bilden eine wichtige Grundlage zu einer gezielten Züchtung von nematodenresistenten Pflanzen.