Cycles of satellite and transposon evolution in Arabidopsis centromeres.

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

vcf2gwas: Python API for comprehensive GWAS analysis using GEMMA.

MOTIVATION: Genome-wide association study (GWAS) requires a researcher to perform a multitude of different actions during analysis. From editing and formatting genotype and phenotype information to running the analysis software to summarizing and visualizing the results. A typical GWAS workflow poses a significant challenge of utilizing the command-line, manual text-editing and requiring knowledge of one or more programming/scripting languages, especially for newcomers. RESULTS: vcf2gwas is a package that provides a convenient pipeline to perform all of the steps of a traditional GWAS workflow by reducing it to a single command-line input of a Variant Call Format file and a phenotype data file. In addition, all the required software is installed with the package. vcf2gwas also implements several useful features enhancing the reproducibility of GWAS analysis. AVAILABILITY AND IMPLEMENTATION: The source code of vcf2gwas is available under the GNU General Public License. The package can be easily installed using conda. Installation instructions and a manual including tutorials can be accessed on the package website at https://github.com/frankvogt/vcf2gwas. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multiple Sources of Introduction of North American Arabidopsis thaliana from across Eurasia.

Large-scale movement of organisms across their habitable range, or migration, is an important evolutionary process that can shape genetic diversity and influence the adaptive spread of alleles. Although human migrations have been studied in great detail with modern and ancient genomes, recent anthropogenic influence on reducing the biogeographical constraints on the migration of nonnative species has presented opportunities in several study systems to ask the questions about how repeated introductions shape genetic diversity in the introduced range. We present an extensive overview of population structure of North American Arabidopsis thaliana by studying a set of 500 whole-genome sequenced and over 2,800 RAD-seq genotyped individuals in the context of global diversity represented by Afro-Eurasian genomes. We use methods based on haplotype and rare-allele sharing as well as phylogenetic modeling to identify likely sources of introductions of extant N. American A. thaliana from the native range in Africa and Eurasia. We find evidence of admixture among the introduced lineages having increased haplotype diversity and reduced mutational load. We also detect signals of selection in immune-system-related genes that may impart qualitative disease resistance to pathogens of bacterial and oomycete origin. We conclude that multiple introductions to a nonnative range can rapidly enhance the adaptive potential of a colonizing species by increasing haplotypic diversity through admixture. Our results lay the foundation for further investigations into the functional significance of admixture.

Host-associated microbe PCR (hamPCR) enables convenient measurement of both microbial load and community composition.

The ratio of microbial population size relative to the amount of host tissue, or 'microbial load', is a fundamental metric of colonization and infection, but it cannot be directly deduced from microbial amplicon data such as 16S rRNA gene counts. Because existing methods to determine load, such as serial dilution plating, quantitative PCR, and whole metagenome sequencing add substantial cost and/or experimental burden, they are only rarely paired with amplicon sequencing. We introduce host-associated microbe PCR (hamPCR), a robust strategy to both quantify microbial load and describe interkingdom microbial community composition in a single amplicon library. We demonstrate its accuracy across multiple study systems, including nematodes and major crops, and further present a cost-saving technique to reduce host overrepresentation in the library prior to sequencing. Because hamPCR provides an accessible experimental solution to the well-known limitations and statistical challenges of compositional data, it has far-reaching potential in culture-independent microbiology.

What natural variation can teach us about resistance durability.

Breeding a crop variety to be resistant to a pathogen usually takes years. This is problematic because pathogens, with short generation times and fluid genomes, adapt quickly to overcome resistance. The triumph of the pathogen is not inevitable, however, as there are numerous examples of durable resistance, particularly in wild plants. Which factors then contribute to such resistance stability over millennia? We review current knowledge of wild and agricultural pathosystems, detailing the importance of genetic, species and spatial heterogeneity in the prevention of pathogen outbreaks. We also highlight challenges associated with increasing resistance diversity in crops, both in light of pathogen (co-)evolution and breeding practices. Historically it has been difficult to incorporate heterogeneity into agriculture due to reduced efficiency in harvesting. Recent advances implementing computer vision and automation in agricultural production may improve our ability to harvest mixed genotype and mixed species plantings, thereby increasing resistance durability.

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

Microorganisms from all domains of life establish associations with plants. Although some harm the plant, others antagonize pathogens or prime the plant immune system, support the acquisition of nutrients, tune plant hormone levels, or perform additional services. Most culture-independent plant microbiome research has focused on amplicon sequencing of the 16S rRNA gene and/or the internal transcribed spacer (ITS) of rRNA genomic loci, which show the relative abundance of the microbes to each other. Here, we describe shotgun sequencing of 275 wild Arabidopsis thaliana leaf microbiomes from southwest Germany, with additional bacterial 16S and eukaryotic ITS1 rRNA amplicon data from 176 of these samples. Shotgun data, which unlike the amplicon data capture the ratio of microbe to plant DNA, enable scaling of microbial read abundances to reflect the microbial load on the host. In a more cost-effective hybrid strategy, we show they also allow a similar scaling of amplicon data to overcome compositionality problems. Our wild plants were dominated by bacterial sequences, with eukaryotes contributing only a minority of reads. 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. Critically, considering absolute microbial load led to fundamentally different conclusions about microbiome assembly and the interaction networks among major taxa.

Role of lysine residues of the Magnaporthe oryzae effector AvrPiz-t in effector- and PAMP-triggered immunity.

Magnaporthe oryzae is an important fungal pathogen of both rice and wheat. However, how M. oryzae effectors modulate plant immunity is not fully understood. Previous studies have shown that the M. oryzae effector AvrPiz-t targets the host ubiquitin-proteasome system to manipulate plant defence. In return, two rice ubiquitin E3 ligases, APIP6 and APIP10, ubiquitinate AvrPiz-t for degradation. To determine how lysine residues contribute to the stability and function of AvrPiz-t, we generated double (K1,2R-AvrPiz-t), triple (K1,2,3R-AvrPiz-t) and lysine-free (LF-AvrPiz-t) mutants by mutating lysines into arginines in AvrPiz-t. LF-AvrPiz-t showed the highest protein accumulation when transiently expressed in rice protoplasts. When co-expressed with APIP10 in Nicotiana benthamiana, LF-AvrPiz-t was more stable than AvrPiz-t and was less able to degrade APIP10. The avirulence of LF-AvrPiz-t on Piz-t:HA plants was less than that of AvrPiz-t, which led to resistance reduction and lower accumulation of the Piz-t:HA protein after inoculation with the LF-AvrPiz-t-carrying isolate. Chitin- and flg22-induced production of reactive oxygen species (ROS) was higher in LF-AvrPiz-t than in AvrPiz-t transgenic plants. In addition, LF-AvrPiz-t transgenic plants were less susceptible than AvrPiz-t transgenic plants to a virulent isolate. Furthermore, both AvrPiz-t and LF-AvrPiz-t interacted with OsRac1, but the suppression of OsRac1-mediated ROS generation by LF-AvrPiz-t was significantly lower than that by AvrPiz-t. Together, these results suggest that the lysine residues of AvrPiz-t are required for its avirulence and virulence functions in rice.

The Monocot-Specific Receptor-like Kinase SDS2 Controls Cell Death and Immunity in Rice.

Programmed cell death (PCD) plays critical roles in plant immunity but must be regulated to prevent excessive damage. The E3 ubiquitin ligase SPL11 negatively regulates PCD and immunity in plants. We show that SPL11 cell-death suppressor 2 (SDS2), an S-domain receptor-like kinase, positively regulates PCD and immunity in rice by engaging and regulating SPL11 and related kinases controlling defense responses. An sds2 mutant shows reduced immune responses and enhanced susceptibility to the blast fungus Magnaporthe oryzae. Conversely, SDS2 over-expression induces constitutive PCD accompanied by elevated immune responses and enhanced resistance to M. oryzae. SDS2 interacts with and phosphorylates SPL11, which in turn ubiquitinates SDS2, leading to its degradation. In addition, SDS2 interacts with related receptor-like cytoplasmic kinases, OsRLCK118/176, that positively regulate immunity by phosphorylating the NADPH oxidase OsRbohB to stimulate ROS production. Thus, a plasma membrane-resident protein complex consisting of SDS2, SPL11, and OsRLCK118/176 controls PCD and immunity in rice.

OsWRKY67 Plays a Positive Role in Basal and XA21-Mediated Resistance in Rice.

WRKY proteins play important roles in transcriptional reprogramming in plants in response to various stresses including pathogen attack. In this study, we functionally characterized a rice WRKY gene, OsWRKY67, whose expression is upregulated against pathogen challenges. Activation of OsWRKY67 by T-DNA tagging significantly improved the resistance against two rice pathogens, Magnaporthe oryzae and Xanthomonas oryzae pv. oryzae (Xoo). Reactive oxygen species (ROS) rapidly accumulated in OsWRKY67 activation mutant lines in response to elicitor treatment, compared with the controls. Overexpression of OsWRKY67 in rice confirmed enhanced disease resistance, but led to a restriction of plant growth in transgenic lines with high levels of OsWRKY67 protein. OsWRKY67 RNAi lines significantly reduced resistance to M. oryzae and Xoo isolates tested, and abolished XA21-mediated resistance, implying the possibility of broad-spectrum resistance from OsWRKY67. Transcriptional activity and subcellular localization assays indicated that OsWRKY67 is present in the nucleus where it functions as a transcriptional activator. Quantitative PCR revealed that the pathogenesis-related genes, PR1a, PR1b, PR4, PR10a, and PR10b, are upregulated in OsWRKY67 overexpression lines. Therefore, these results suggest that OsWRKY67 positively regulates basal and XA21-mediated resistance, and is a promising candidate for genetic improvement of disease resistance in rice.

Quantification of hydrogen peroxide in plant tissues using Amplex Red.

Reactive oxygen species (ROS) are by-products of photosynthesis and respiration in plant tissues. Abiotic and biotic stressors also induce the production and temporary accumulation of ROS in plants, including hydrogen peroxide (H2O2), whereby they can act as secondary messengers/chemical mediators in plant defense signaling and lead to programmed cell death. H2O2 acts as a hub for critical information flow in plants. Despite such key roles in fundamental cellular processes, reliable determination of H2O2 levels in plant tissues is hard to achieve. We optimized an Amplex Red-based quantitation method for H2O2 estimation from plant tissue lysate. The standard limit of detection and quantitation was determined as 6 and 18picomol respectively. In this study we also quantified constitutive and/or induced levels of H2O2 in three model plants, Pinus nigra (Austrian pine), Oryza sativa (rice), and Arabidopsis thaliana. Overall, assay sensitivity was in the nmolg-1 FW range. Commonly used additives for H2O2 extraction such as activated charcoal, ammonium sulfate, perchloric acid, polyvinylpolypyrrolidone, and trichloroacetic acid either degraded H2O2 directly or interfered with the Amplex Red assay. Finally, We measured stability of Amplex Red working solution over one month of storage at -80°C and found it to be significantly stable over time. With appropriate modifications, this optimized method should be applicable to any plant tissue.

Temporal patterns of damage and decay kinetics of DNA retrieved from plant herbarium specimens.

Herbaria archive a record of changes of worldwide plant biodiversity harbouring millions of specimens that contain DNA suitable for genome sequencing. To profit from this resource, it is fundamental to understand in detail the process of DNA degradation in herbarium specimens. We investigated patterns of DNA fragmentation and nucleotide misincorporation by analysing 86 herbarium samples spanning the last 300 years using Illumina shotgun sequencing. We found an exponential decay relationship between DNA fragmentation and time, and estimated a per nucleotide fragmentation rate of 1.66 × 10(-4) per year, which is six times faster than the rate estimated for ancient bones. Additionally, we found that strand breaks occur specially before purines, and that depurination-driven DNA breakage occurs constantly through time and can to a great extent explain decreasing fragment length over time. Similar to what has been found analysing ancient DNA from bones, we found a strong correlation between the deamination-driven accumulation of cytosine to thymine substitutions and time, which reinforces the importance of substitution patterns to authenticate the ancient/historical nature of DNA fragments. Accurate estimations of DNA degradation through time will allow informed decisions about laboratory and computational procedures to take advantage of the vast collection of worldwide herbarium specimens.

The E3 Ligase APIP10 Connects the Effector AvrPiz-t to the NLR Receptor Piz-t in Rice.

Although nucleotide-binding domain, leucine-rich repeat (NLR) proteins are the major immune receptors in plants, the mechanism that controls their activation and immune signaling remains elusive. Here, we report that the avirulence effector AvrPiz-t from Magnaporthe oryzae targets the rice E3 ligase APIP10 for degradation, but that APIP10, in return, ubiquitinates AvrPiz-t and thereby causes its degradation. Silencing of APIP10 in the non-Piz-t background compromises the basal defense against M. oryzae. Conversely, silencing of APIP10 in the Piz-t background causes cell death, significant accumulation of Piz-t, and enhanced resistance to M. oryzae, suggesting that APIP10 is a negative regulator of Piz-t. We show that APIP10 promotes degradation of Piz-t via the 26S proteasome system. Furthermore, we demonstrate that AvrPiz-t stabilizes Piz-t during M. oryzae infection. Together, our results show that APIP10 is a novel E3 ligase that functionally connects the fungal effector AvrPiz-t to its NLR receptor Piz-t in rice.

Identification and characterization of suppressor mutants of spl11- mediated cell death in rice.

Lesion mimic mutants have been used to dissect programmed cell death (PCD) and defense-related pathways in plants. The rice lesion-mimic mutant spl11 exhibits race nonspecific resistance to the bacterial pathogen Xanthomonas oryzae pv. oryzae and the fungal pathogen Magnaporthe oryzae. Spl11 encodes an E3 ubiquitin ligase and is a negative regulator of PCD in rice. To study the regulation of Spl11-mediated PCD, we performed a genetic screen and identified three spl11 cell-death suppressor (sds) mutants. These suppressors were characterized for their resistance to X. oryzae pv. oryzae and M. oryzae and for their expression of defense-related genes. The suppression of the cell-death phenotypes was generally correlated with reduced expression of defense-related genes. When rice was challenged with avirulent isolates of M. oryzae, the disease phenotype was unaffected in the sds mutants, indicating that the suppression might be Spl11-mediated pathway specific and may only be involved in basal defense. In addition, we mapped one of the suppressor mutations to a 140-kb interval on the long arm of rice chromosome 1. Identification and characterization of these sds mutants should facilitate our efforts to elucidate the Spl11-mediated PCD pathway.

The Magnaporthe oryzae effector AvrPiz-t targets the RING E3 ubiquitin ligase APIP6 to suppress pathogen-associated molecular pattern-triggered immunity in rice.

Although the functions of a few effector proteins produced by bacterial and oomycete plant pathogens have been elucidated in recent years, information for the vast majority of pathogen effectors is still lacking, particularly for those of plant-pathogenic fungi. Here, we show that the avirulence effector AvrPiz-t from the rice blast fungus Magnaporthe oryzae preferentially accumulates in the specialized structure called the biotrophic interfacial complex and is then translocated into rice (Oryza sativa) cells. Ectopic expression of AvrPiz-t in transgenic rice suppresses the flg22- and chitin-induced generation of reactive oxygen species (ROS) and enhances susceptibility to M. oryzae, indicating that AvrPiz-t functions to suppress pathogen-associated molecular pattern (PAMP)-triggered immunity in rice. Interaction assays show that AvrPiz-t suppresses the ubiquitin ligase activity of the rice RING E3 ubiquitin ligase APIP6 and that, in return, APIP6 ubiquitinates AvrPiz-t in vitro. Interestingly, agroinfection assays reveal that AvrPiz-t and AvrPiz-t Interacting Protein 6 (APIP6) are both degraded when coexpressed in Nicotiana benthamiana. Silencing of APIP6 in transgenic rice leads to a significant reduction of flg22-induced ROS generation, suppression of defense-related gene expression, and enhanced susceptibility of rice plants to M. oryzae. Taken together, our results reveal a mechanism in which a fungal effector targets the host ubiquitin proteasome system for the suppression of PAMP-triggered immunity in plants.