Showing papers by "Detlef Weigel published in 2019"

Arabidopsis FLL2 promotes liquid–liquid phase separation of polyadenylation complexes

Abstract: An important component of cellular biochemistry is the concentration of proteins and nucleic acids in non-membranous compartments1,2. These biomolecular condensates are formed from processes that include liquid-liquid phase separation. The multivalent interactions necessary for liquid-liquid phase separation have been extensively studied in vitro1,3. However, the regulation of this process in vivo is poorly understood. Here we identify an in vivo regulator of liquid-liquid phase separation through a genetic screen targeting factors required for Arabidopsis RNA-binding protein FCA function. FCA contains prion-like domains that phase-separate in vitro, and exhibits behaviour in vivo that is consistent with phase separation. The mutant screen identified a functional requirement for FLL2, a coiled-coil protein, in the formation of FCA nuclear bodies. FCA reduces transcriptional read-through by promoting proximal polyadenylation at many sites in the Arabidopsis genome3,4. FLL2 was required to promote this proximal polyadenylation, but not the binding of FCA to target RNA. Ectopic expression of FLL2 increased the size and number of FCA nuclear bodies. Crosslinking with formaldehyde captured in vivo interactions between FLL2, FCA and the polymerase and nuclease modules of the RNA 3'-end processing machinery. These 3' RNA-processing components colocalized with FCA in the nuclear bodies in vivo, which indicates that FCA nuclear bodies compartmentalize 3'-end processing factors to enhance polyadenylation at specific sites. Our findings show that coiled-coil proteins can promote liquid-liquid phase separation, which expands our understanding of the principles that govern the in vivo dynamics of liquid-like bodies.

Natural selection on the Arabidopsis thaliana genome in present and future climates

Abstract: Through the lens of evolution, climate change is an agent of natural selection that forces populations to change and adapt, or face extinction. However, current assessments of the risk of biodiversity associated with climate change1 do not typically take into account how natural selection influences populations differently depending on their genetic makeup2. Here we make use of the extensive genome information that is available for Arabidopsis thaliana and measure how manipulation of the amount of rainfall affected the fitness of 517 natural Arabidopsis lines that were grown in Spain and Germany. This allowed us to directly infer selection along the genome3. Natural selection was particularly strong in the hot-dry location in Spain, where 63% of lines were killed and where natural selection substantially changed the frequency of approximately 5% of all genome-wide variants. A significant portion of this climate-driven natural selection of variants was predictable from signatures of local adaptation (R2 = 29–52%), as genetic variants that were found in geographical areas with climates more similar to the experimental sites were positively selected. Field-validated predictions across the species range indicated that Mediterranean and western Siberian populations—at the edges of the environmental limits of this species—currently experience the strongest climate-driven selection. With more frequent droughts and rising temperatures in Europe4, we forecast an increase in directional natural selection moving northwards from the southern end of Europe, putting many native A. thaliana populations at evolutionary risk. The predicted increase in frequency of droughts and rising temperatures in Europe will lead core populations of a temperate plant to an evolutionary dead-end unless they acquire genetic alleles that are present only in extreme edge Mediterranean, Scandinavian, or Siberian populations.

Multiple modes of convergent adaptation in the spread of glyphosate-resistant Amaranthus tuberculatus.

Abstract: The selection pressure exerted by herbicides has led to the repeated evolution of herbicide resistance in weeds. The evolution of herbicide resistance on contemporary timescales in turn provides an outstanding opportunity to investigate key questions about the genetics of adaptation, in particular the relative importance of adaptation from new mutations, standing genetic variation, or geographic spread of adaptive alleles through gene flow. Glyphosate-resistant Amaranthus tuberculatus poses one of the most significant threats to crop yields in the Midwestern United States, with both agricultural populations and herbicide resistance only recently emerging in Canada. To understand the evolutionary mechanisms driving the spread of resistance, we sequenced and assembled the A. tuberculatus genome and investigated the origins and population genomics of 163 resequenced glyphosate-resistant and susceptible individuals from Canada and the United States. In Canada, we discovered multiple modes of convergent evolution: in one locality, resistance appears to have evolved through introductions of preadapted US genotypes, while in another, there is evidence for the independent evolution of resistance on genomic backgrounds that are historically nonagricultural. Moreover, resistance on these local, nonagricultural backgrounds appears to have occurred predominantly through the partial sweep of a single haplotype. In contrast, resistant haplotypes arising from the Midwestern United States show multiple amplification haplotypes segregating both between and within populations. Therefore, while the remarkable species-wide diversity of A. tuberculatus has facilitated geographic parallel adaptation of glyphosate resistance, more recently established agricultural populations are limited to adaptation in a more mutation-limited framework.

Fine-Grained Analysis of Spontaneous Mutation Spectrum and Frequency in Arabidopsis thaliana.

Abstract: Mutations are the ultimate source of all genetic variation. However, few direct estimates of the contribution of mutation to molecular genetic variation are available. To address this issue, we first analyzed the rate and spectrum of mutations in the Arabidopsis thaliana reference accession after 25 generations of single-seed descent. We then compared the mutation profile in these mutation accumulation (MA) lines against genetic variation observed in the 1001 Genomes Project. The estimated haploid single nucleotide mutation (SNM) rate for A. thaliana is 6.95 × 10−9 (SE ± 2.68 × 10−10) per site per generation, with SNMs having higher frequency in transposable elements (TEs) and centromeric regions. The estimated indel mutation rate is 1.30 × 10−9 (±1.07 × 10−10) per site per generation, with deletions being more frequent and larger than insertions. Among the 1694 unique SNMs identified in the MA lines, the positions of 389 SNMs (23%) coincide with biallelic SNPs from the 1001 Genomes population, and in 289 (17%) cases the changes are identical. Of the 329 unique indels identified in the MA lines, 96 (29%) overlap with indels from the 1001 Genomes dataset, and 16 indels (5% of the total) are identical. These overlap frequencies are significantly higher than expected, suggesting that de novo mutations are not uniformly distributed and arise at polymorphic sites more frequently than assumed. These results suggest that high mutation rate potentially contributes to high polymorphism and low mutation rate to reduced polymorphism in natural populations providing insights of mutational inputs in generating natural genetic diversity.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Abstract: Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

Comparing Arabidopsis receptor kinase and receptor protein-mediated immune signaling reveals BIK1-dependent differences.

Abstract: Pattern recognition receptors (PRRs) sense microbial patterns and activate innate immunity against attempted microbial invasions. The leucine-rich repeat receptor kinases (LRR-RK) FLS2 and EFR, and the LRR receptor protein (LRR-RP) receptors RLP23 and RLP42, respectively, represent prototypical members of these two prominent and closely related PRR families. We conducted a survey of Arabidopsis thaliana immune signaling mediated by these receptors to address the question of commonalities and differences between LRR-RK and LRR-RP signaling. Quantitative differences in timing and amplitude were observed for several early immune responses, with RP-mediated responses typically being slower and more prolonged than those mediated by RKs. Activation of RLP23, but not FLS2, induced the production of camalexin. Transcriptomic analysis revealed that RLP23-regulated genes represent only a fraction of those genes differentially expressed upon FLS2 activation. Several positive and negative regulators of FLS2-signaling play similar roles in RLP23 signaling. Intriguingly, the cytoplasmic receptor kinase BIK1, a positive regulator of RK signaling, acts as a negative regulator of RP-type immune receptors in a manner dependent on BIK1 kinase activity. Our study unveiled unexpected differences in two closely related receptor systems and reports a new negative role of BIK1 in plant immunity.

Long-term balancing selection drives evolution of immunity genes in Capsella.

Abstract: Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

RPW8/HR repeats control NLR activation in Arabidopsis thaliana

Abstract: In many plant species, conflicts between divergent elements of the immune system, especially nucleotide-binding oligomerization domain-like receptors (NLR), can lead to hybrid necrosis. Here, we report deleterious allele-specific interactions between an NLR and a non-NLR gene cluster, resulting in not one, but multiple hybrid necrosis cases in Arabidopsis thaliana. The NLR cluster is RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7), which can confer strain-specific resistance to oomycetes. The non-NLR cluster is RESISTANCE TO POWDERY MILDEW 8 (RPW8) / HOMOLOG OF RPW8 (HR), which can confer broad-spectrum resistance to both fungi and oomycetes. RPW8/HR proteins contain at the N-terminus a potential transmembrane domain, followed by a specific coiled-coil (CC) domain that is similar to a domain found in pore-forming toxins MLKL and HET-S from mammals and fungi. C-terminal to the CC domain is a variable number of 21- or 14-amino acid repeats, reminiscent of regulatory 21-amino acid repeats in fungal HET-S. The number of repeats in different RPW8/HR proteins along with the sequence of a short C-terminal tail predicts their ability to activate immunity in combination with specific RPP7 partners. Whether a larger or smaller number of repeats is more dangerous depends on the specific RPW8/HR autoimmune risk variant.

AraPheno and the AraGWAS Catalog 2020: a major database update including RNA-Seq and knockout mutation data for Arabidopsis thaliana.

Abstract: Genome-wide association studies (GWAS) are integral for studying genotype-phenotype relationships and gaining a deeper understanding of the genetic architecture underlying trait variation. A plethora of genetic associations between distinct loci and various traits have been successfully discovered and published for the model plant Arabidopsis thaliana. This success and the free availability of full genomes and phenotypic data for more than 1,000 different natural inbred lines led to the development of several data repositories. AraPheno (https://arapheno.1001genomes.org) serves as a central repository of population-scale phenotypes in A. thaliana, while the AraGWAS Catalog (https://aragwas.1001genomes.org) provides a publicly available, manually curated and standardized collection of marker-trait associations for all available phenotypes from AraPheno. In this major update, we introduce the next generation of both platforms, including new data, features and tools. We included novel results on associations between knockout-mutations and all AraPheno traits. Furthermore, AraPheno has been extended to display RNA-Seq data for hundreds of accessions, providing expression information for over 28 000 genes for these accessions. All data, including the imputed genotype matrix used for GWAS, are easily downloadable via the respective databases.

Phytochrome activates the plastid-encoded RNA polymerase for chloroplast biogenesis via nucleus-to-plastid signaling.

Abstract: Light initiates chloroplast biogenesis by activating photosynthesis-associated genes encoded by not only the nuclear but also the plastidial genome, but how photoreceptors control plastidial gene expression remains enigmatic. Here we show that the photoactivation of phytochromes triggers the expression of photosynthesis-associated plastid-encoded genes (PhAPGs) by stimulating the assembly of the bacterial-type plastidial RNA polymerase (PEP) into a 1000-kDa complex. Using forward genetic approaches, we identified REGULATOR OF CHLOROPLAST BIOGENESIS (RCB) as a dual-targeted nuclear/plastidial phytochrome signaling component required for PEP assembly. Surprisingly, RCB controls PhAPG expression primarily from the nucleus by interacting with phytochromes and promoting their localization to photobodies for the degradation of the transcriptional regulators PIF1 and PIF3. RCB-dependent PIF degradation in the nucleus signals the plastids for PEP assembly and PhAPG expression. Thus, our findings reveal the framework of a nucleus-to-plastid anterograde signaling pathway by which phytochrome signaling in the nucleus controls plastidial transcription.

RST1 and RIPR connect the cytosolic RNA exosome to the Ski complex in Arabidopsis.

Abstract: The RNA exosome is a key 3'-5' exoribonuclease with an evolutionarily conserved structure and function. Its cytosolic functions require the co-factors SKI7 and the Ski complex. Here we demonstrate by co-purification experiments that the ARM-repeat protein RESURRECTION1 (RST1) and RST1 INTERACTING PROTEIN (RIPR) connect the cytosolic Arabidopsis RNA exosome to the Ski complex. rst1 and ripr mutants accumulate RNA quality control siRNAs (rqc-siRNAs) produced by the post-transcriptional gene silencing (PTGS) machinery when mRNA degradation is compromised. The small RNA populations observed in rst1 and ripr mutants are also detected in mutants lacking the RRP45B/CER7 core exosome subunit. Thus, molecular and genetic evidence supports a physical and functional link between RST1, RIPR and the RNA exosome. Our data reveal the existence of additional cytosolic exosome co-factors besides the known Ski subunits. RST1 is not restricted to plants, as homologues with a similar domain architecture but unknown function exist in animals, including humans.

Nonlinear phenotypic variation uncovers the emergence of heterosis in Arabidopsis thaliana

Abstract: Heterosis describes the phenotypic superiority of hybrids over their parents in traits related to agronomic performance and fitness. Understanding and predicting nonadditive inheritance such as heterosis is crucial for evolutionary biology as well as for plant and animal breeding. However, the physiological bases of heterosis remain debated. Moreover, empirical data in various species have shown that diverse genetic and molecular mechanisms are likely to explain heterosis, making it difficult to predict its emergence and amplitude from parental genotypes alone. In this study, we examined a model of physiological dominance initially proposed by Sewall Wright to explain the nonadditive inheritance of traits like metabolic fluxes at the cellular level. We evaluated Wright's model for two fitness-related traits at the whole-plant level, growth rate and fruit number, using 450 hybrids derived from crosses among natural accessions of A. thaliana. We found that allometric relationships between traits constrain phenotypic variation in a nonlinear and similar manner in hybrids and accessions. These allometric relationships behave predictably, explaining up to 75% of heterosis amplitude, while genetic distance among parents at best explains 7%. Thus, our findings are consistent with Wright's model of physiological dominance and suggest that the emergence of heterosis on plant performance is an intrinsic property of nonlinear relationships between traits. Furthermore, our study highlights the potential of a geometric approach of phenotypic relationships for predicting heterosis of major components of crop productivity and yield.

NCP activates chloroplast transcription by controlling phytochrome-dependent dual nuclear and plastidial switches.

Abstract: Phytochromes initiate chloroplast biogenesis by activating genes encoding the photosynthetic apparatus, including photosynthesis-associated plastid-encoded genes (PhAPGs). PhAPGs are transcribed by a bacterial-type RNA polymerase (PEP), but how phytochromes in the nucleus activate chloroplast gene expression remains enigmatic. We report here a forward genetic screen in Arabidopsis that identified NUCLEAR CONTROL OF PEP ACTIVITY (NCP) as a necessary component of phytochrome signaling for PhAPG activation. NCP is dual-targeted to plastids and the nucleus. While nuclear NCP mediates the degradation of two repressors of chloroplast biogenesis, PIF1 and PIF3, NCP in plastids promotes the assembly of the PEP complex for PhAPG transcription. NCP and its paralog RCB are non-catalytic thioredoxin-like proteins that diverged in seed plants to adopt nonredundant functions in phytochrome signaling. These results support a model in which phytochromes control PhAPG expression through light-dependent double nuclear and plastidial switches that are linked by evolutionarily conserved and dual-localized regulatory proteins.

The relationship between microbial biomass and disease in the Arabidopsis thaliana phyllosphere

Abstract: A central goal in microbiome research is to learn what distinguishes a healthy from a dysbiotic microbial community. Shifts in diversity and taxonomic composition are important indicators of dysbiosis, but a full understanding also requires knowledge of absolute microbial biomass. Simultaneous information on both microbiome composition and the quantity of its components can provide insight into microbiome function and disease state. Here we use shotgun metagenomics to simultaneously assess microbiome composition and microbial load in the phyllosphere of wild populations of the plant Arabidopsis thaliana. We find that wild plants vary substantially in the load of colonizing microbes, and that high loads are typically associated with the proliferation of single taxa, with only a few putatively pathogenic taxa achieving high abundances in the field. Our results suggest (i) that the inside of a plant leaf is on average sparsely colonized with an estimated two bacterial genomes per plant genome and an order of magnitude fewer eukaryotic microbial genomes, and (ii) that higher levels of microbial biomass often indicate successful colonization by pathogens. Lastly, our results show that load is a significant explanatory variable for loss of estimated Shannon diversity in phyllosphere microbiomes, implying that reduced diversity may be a significant predictor of microbial dysbiosis in a plant leaf.

The relationship between microbial population size and disease in the Arabidopsis thaliana phyllosphere

Abstract: A central goal in microbiome research is to learn what distinguishes a healthy from a dysbiotic microbial community. Shifts in diversity and taxonomic composition are important indicators of dysbiosis, but a full understanding also requires knowledge of absolute microbial population sizes. In addition to the number of microbial cells, information on taxonomic composition can provide important insight into microbiome function and disease state. Here we use shotgun metagenomics to simultaneously assess microbiome composition and microbial load in the phyllosphere of wild populations of the plant Arabidopsis thaliana. We find that wild plants vary substantially in the load of colonizing microbes, and that high loads are typically associated with the proliferation of single taxa, with only a few putatively pathogenic taxa achieving high abundances in the field. Our results suggest (i) that the inside of a plant leaf is on average sparsely colonized with an estimated two bacterial genomes per plant genome and an order of magnitude fewer eukaryotic microbial genomes, and (ii) that higher levels of microbial cells often indicate successful colonization by pathogens. Lastly, our results show that load is a significant explanatory variable for loss of estimated Shannon diversity in phyllosphere microbiomes, implying that reduced diversity may be a significant predictor of microbial dysbiosis in a plant leaf.

Combining whole genome shotgun sequencing and rDNA amplicon analyses to improve detection of microbe-microbe interaction networks in plant leaves

Abstract: Microorganisms from all domains of life establish associations with plants. Although some harm the plant, others antagonize pathogens or prime the plant immune system, acquire nutrients, tune plant hormone levels, or perform additional services. Most culture-independent plant microbiome research has focused on amplicon sequencing of 16S rDNA and/or the internal transcribed spacer (ITS) of rDNA loci, but the decreasing cost of high-throughput sequencing has made shotgun metagenome sequencing increasingly accessible. Here, we describe shotgun sequencing of 275 wild Arabidopsis thaliana leaf microbiomes from southwest Germany, with additional bacterial 16S rDNA and eukaryotic ITS1 amplicon data from 176 of these samples. The shotgun data were dominated by bacterial sequences, with eukaryotes contributing only a minority of reads. For shotgun and amplicon data, microbial membership showed weak associations with both site of origin and plant genotype, both of which were highly confounded in this dataset. There was large variation among microbiomes, with one extreme comprising samples of low complexity and a high load of microorganisms typical of infected plants, and the other extreme being samples of high complexity and a low microbial load. We use the metagenome data, which captures the ratio of bacterial to plant DNA in leaves of wild plants, to scale the 16S rDNA amplicon data such that they reflect absolute bacterial abundance. We show that this cost-effective hybrid strategy overcomes compositionality problems in amplicon data and leads to fundamentally different conclusions about microbiome community assembly.

The Arabidopsis thaliana pan-NLRome

Abstract: 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.

Multiple modes of convergent adaptation in the spread of glyphosate-resistant Amaranthus tuberculatus

Abstract: The selection pressure exerted by herbicides has led to the repeated evolution of resistance in weeds. The evolution of herbicide resistance on contemporary timescales provides an outstanding opportunity to investigate key open questions about the genetics of adaptation, in particular the relative importance of adaptation from new mutations, standing genetic variation, and geographic spread of adaptive alleles through gene flow. Glyphosate-resistant Amaranthus tuberculatus poses one of the most significant threats to crop yields in the midwestern United States, with both agricultural populations and resistance only recently emerging in Canada. To understand the evolutionary mechanisms driving the spread of resistance, we sequenced and assembled the A. tuberculatus genome and investigated the origins and population genomics of 163 resequenced glyphosate-resistant and susceptible individuals in Canada and the USA. In Canada, we discovered multiple modes of convergent evolution: in one locality, resistance appears to have evolved through introductions of preadapted US genotypes, while in another, there is evidence for the independent evolution of resistance on genomic backgrounds that are historically non-agricultural. Moreover, resistance on these local, non-agricultural backgrounds appears to have occurred predominantly through the partial sweep of a single amplification haplotype. In contrast, US genotypes and those in Canada introduced from the US show multiple amplification haplotypes segregating both between and within populations. Therefore, while the remarkable diversity of A. tuberculatus has facilitated geographic parallel adaptation of glyphosate resistance, different timescales of selection have favored either adaptation from standing variation or de novo mutation in certain parts of the range.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses

Abstract: The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

An ultra high-density Arabidopsis thaliana crossover map that refines the influences of structural variation and epigenetic features

Abstract: Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of over 17,000 COs between the Col-0 and Ler accessions of Arabidopsis thaliana. COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders, and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large and small-scale SVs, representing influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

EFFECTOR OF TRANSCRIPTION factors are novel plant-specific regulators associated with genomic DNA methylation in Arabidopsis

Abstract: Plant-specific EFFECTORS OF TRANSCRIPTION (ET) are characterised by a variable number of highly conserved ET repeats, which are involved in zinc and DNA binding. In addition, ETs share a GIY-YIG domain, involved in DNA nicking activity. It was hypothesised that ETs might act as epigenetic regulators. Here, methylome, transcriptome and phenotypic analyses were performed to investigate the role of ET factors and their involvement in DNA methylation in Arabidopsis thaliana. Comparative DNA methylation and transcriptome analyses in flowers and seedlings of et mutants revealed ET-specific differentially expressed genes and mostly independently characteristic, ET-specific differentially methylated regions. Loss of ET function results in pleiotropic developmental defects. The accumulation of cyclobutane pyrimidine dimers after ultraviolet stress in et mutants suggests an ET function in DNA repair.

RPW8/HR Repeats Predict NLR-dependent Hybrid Performance

Abstract: 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.

RST1 and RIPR connect the cytosolic RNA exosome to the Ski complex in Arabidopsis

Abstract: The RNA exosome is a key 3’-5’ exoribonuclease with an evolutionarily conserved structure and function. Its cytosolic functions require the co-factors SKI7 and the Ski complex. Here we demonstrate by co-purification experiments that the ARM repeat protein RESURRECTION1 (RST1) and RST1 INTERACTING PROTEIN (RIPR) connect the cytosolic Arabidopsis RNA exosome to the Ski complex. rst1 and ripr mutants accumulate RNA quality control siRNAs (rqc-siRNAs) produced by the post-transcriptional gene silencing (PTGS) machinery when mRNA degradation is compromised. The small RNA populations observed in rst1 and ripr mutants are also detected in mutants lacking the RRP45B/CER7 core exosome subunit. Thus, molecular and genetic evidence supports a physical and functional link between RST1, RIPR and the RNA exosome. Our data reveal the existence of additional cytosolic exosome co-factors besides the known SKI subunits. RST1 is not restricted to plants, as homologues with a similar domain architecture but unknown function exist in animals, including humans.

Finding genetic variants in plants without complete genomes

Abstract: Structural variants and presence/absence polymorphisms are common in plant genomes, yet they are routinely overlooked in genome-wide association studies (GWAS). Here, we expand the genetic variants detected in GWAS to include major deletions, insertions, and rearrangements. We first use raw sequencing data directly to derive short sequences, k-mers, that mark a broad range of polymorphisms independently of a reference genome. We then link k-mers associated with phenotypes to specific genomic regions. Using this approach, we re-analyzed 2,000 traits measured in Arabidopsis thaliana, tomato, and maize populations. Associations identified with k-mers recapitulate those found with single-nucleotide polymorphisms (SNPs), however, with stronger statistical support. Moreover, we identified new associations with structural variants and with regions missing from reference genomes. Our results demonstrate the power of performing GWAS before linking sequence reads to specific genomic regions, which allow detection of a wider range of genetic variants responsible for phenotypic variation.

Publisher Correction: Natural selection on the Arabidopsis thaliana genome in present and future climates.

Abstract: An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

RPW8/HR Repeats Control NLR Activation in A. thaliana

Abstract: In many plant species, conflicts between divergent elements of the immune system, especially nucleotide-binding oligomerization domain-like receptors (NLR), can lead to hybrid necrosis. Here, we report deleterious allele-specific interactions between an NLR and a non-NLR gene cluster, resulting in not one, but multiple hybrid necrosis cases in Arabidopsis thaliana. The NLR cluster is RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7), which can confer strain-specific resistance to oomycetes. The non-NLR cluster is RESISTANCE TO POWDERY MILDEW 8 (RPW8) / HOMOLOG OF RPW8 (HR), which can confer broad-spectrum resistance to both fungi and oomycetes. RPW8/HR proteins contain at the N-terminus a potential transmembrane domain, followed by a specific coiled-coil (CC) domain that is similar to a domain found in pore-forming toxins MLKL and HET-S from mammals and fungi. C-terminal to the CC domain is a variable number of 21- or 14-amino acid repeats, reminiscent of regulatory 21-amino acid repeats in fungal HET-S. The number of repeats in different RPW8/HR proteins along with the sequence of a short C-terminal tail predicts their ability to activate immunity in combination with specific RPP7 partners. Whether a larger or smaller number of repeats is more dangerous depends on the specific RPW8/HR autoimmune risk variant. Author Summary 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, RESISTANCE TO PERONOSPORA PARASITICA 7 (RPP7) and RESISTANCE TO POWDERY MILDEW 8 (RPW8)/HOMOLOG OF RPW8 (HR). The RPP7 locus encodes intracellular nucleotide binding site-leucine rich repeat (NLR) immune receptors that can confer strain-specific resistance to oomycetes, while the RPW8/HR locus encodes atypical resistance proteins, of which some 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, 21- or 14-amino acid long RPW8/HR repeats. We demonstrate that the number of RPW8/HR repeats and the short C-terminal tail correlate, in an allele-specific manner, with the severity of hybrid necrosis when these alleles are combined with RPP7 variants. We discuss these findings in light of sequence similarity between RPW8/HR and pore-forming toxins MLKL and HET-S from mammals and fungi.

Position-dependent effects of cytosine methylation on FWA expression in Arabidopsis thaliana

Abstract: Gene expression can be modulated by epigenetic modifications to chromatin, and variants of the same locus distinguished by fixed, heritable epigenetic differences are known as epialleles. DNA methylation at cytosines is a prominent epigenetic modification, particularly in plant genomes, that can modulate gene expression. There are several examples where epialleles are associated with differentially methylated regions that affect the expression of overlapping or close-by genes. However, there are also many differentially methylated regions that have not been assigned a biological function despite their proximity to genes. We investigated the positional importance of DNA methylation at the FWA (FLOWERING WAGENINGEN) locus in Arabidopsis thaliana, a paradigm for stable epialleles. We show that cytosine methylation can be established not only over the well-characterized SINE-derived repeat elements that overlap with the transcription start site, but also in more distal promoter regions. FWA silencing, however, is most effective when methylation covers the transcription start site.

Specific Regulation of TCP genes by miR319

Abstract: MicroRNAs (miRNAs) are major regulators of gene expression in multicellular organisms. Many of the evolutionarily conserved miRNAs in plants are encoded by small gene families. The miR159/miR319 family has six members of similar sequences sharing 17 nucleotides in Arabidopsis thaliana. The members of this miRNA family regulate genes encoding TCP (TEOSINTE BRANCHED1, CYCLOIDEA and PCF1/2) and MYB transcription factors. However, despite their sequence similarities, genetic evidence indicates that miR159 and miR319 fulfil different roles in vivo. Here, we confirm previous findings showing that TCP genes are not targeted by miR159. Thus, specific small sequence differences between the miRNAs allow for the specific regulation of TCP transcription factors by miR319 miRNAs.