Summary Hybrid offspring can look very different from their parents, including having greatly increased or decreased fitness. In many plant species, conflicts between divergent elements of the immune system can cause hybrids to express autoimmunity, a generally deleterious syndrome known as hybrid necrosis. We are investigating multiple hybrid necrosis cases in Arabidopsis thaliana that are caused by allele-specific interactions between different variants at two unlinked resistance (R) gene clusters. One is the RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7) cluster, which encodes an intracellular nucleotide binding site-leucine rich repeat (NLR) immune receptors that confer strain-specific resistance to oomycetes. The other is the RESISTANCE TO POWDERY MILDEW 8 (RPW8)/HOMOLOG OF RPW8 (HR) locus, which encodes atypical resistance proteins that can confer broad-spectrum resistance to filamentous pathogens. There is extensive structural variation in the RPW8/HR cluster, both at the level of gene copy number and at the level of C-terminal protein repeats of unknown function. We demonstrate that the number of RPW8/HR repeats correlate, albeit in a complex manner, with the severity of hybrid necrosis when these alleles are combined with specific RPP7 variants. This observation suggests that gross structural differences, rather than individual amino acid polymorphisms, guide the genetic interaction between RPW8/HR and RPP7 alleles. We discuss these findings in light of the similarity of RPW8/HR proteins with pore-forming toxins, MLKL and HET-S, from mammals and fungi.

Disease is both among the most important selection pressures in nature and among the main causes of yield loss in agriculture. In plants, resistance to disease is often conferred by Nucleotide-binding Leucine-rich Repeat (NLR) proteins. These proteins function as intracellular immune receptors that recognize pathogen proteins and their effects on the plant. Consistent with evolutionarily dynamic interactions between plants and pathogens, NLRs are known to be encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define the majority of the Arabidopsis thaliana species-wide 'NLRome'. From NLR sequence enrichment and long-read sequencing of 65 diverse A. thaliana accessions, we infer that the pan-NLRome saturates with approximately 40 accessions. Despite the high diversity of NLRs, half of the pan-NLRome is present in most accessions. We chart the architectural diversity of NLR proteins, identify novel architectures, and quantify the selective forces that act on specific NLRs, domains, and positions. Our study provides a blueprint for defining the pan-NLRome of plant species.

As regulators of gene expression in multicellular organisms, microRNAs (miRNAs) are crucial for growth and development. While a plethora of factors involved in their biogenesis and action in Arabidopsis thaliana have been described, these processes and their fine-tuning are not fully understood in plants. Here, we used plants expressing an artificial miRNA target mimic (MIM) to screen for negative regulators of miR156 activity. We identified a new mutant allele of the F-box protein HAWAIIAN SKIRT (HWS; At3G61590), hws-5, as a suppressor of the MIM156-induced developmental and molecular phenotypes. In hws plants, levels of endogenous miRNAs are increased and their mRNA targets decreased. Plants constitutively expressing full-length HWS — but not a truncated version lacking the F-box domain — display morphological and molecular phenotypes resembling those of mutants defective in miRNA biogenesis and activity. In combination with such mutants, hws loses its delayed floral organ abscission (‘skirt’) phenotype, suggesting epistasis. Also, the overall hws transcriptome profile partially resembles well-known miRNA mutants hyl1-2 and se-3, indicating action in a common pathway. We thus propose HWS as a novel, F-box dependent regulator of miRNA biogenesis.

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves as a nexus for the trade-off between growth and pathogen resistance in wild populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, first identified in the natural accession Est-1, is found in over 10% of wild strains, even though it causes a clear fitness penalty under optimal growth conditions. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that the population genetic basis of ACD6 modulation is complex, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

The gold standard for studying natural selection and adaptation in the wild is to quantify lifetime fitness of individuals from natural populations that have been grown together in a common garden, or that have been reciprocally transplanted. By combining fitness values with species traits and genome sequences, one can infer selection coefficients at the genetic level. Here we present a rainfall-manipulation experiment with 517 whole-genome sequenced natural accessions of the plant Arabidopsis thaliana spanning the global distribution of the species. The experiments were conducted in two field stations in contrasting climates, in the Mediterranean and in Central Europe, where we built rainout shelters and simulated high and low rainfall. Using custom image analysis we quantified fitness- and phenology-related traits for 23,154 pots, which contained about 14,500 plants growing independently, and over 310,000 plants growing in small populations (max. 30 plants). This large field experiment dataset, which associates fitness and ecologically-relevant traits with genomes, will provide an important resource to test eco-evolutionary genetic theories and to understand the potential evolutionary impacts of future climates on an important plant model species.