Method development and application of Next Generation Sequencing in forward genetics

Forward genetic screens remain one of the main genetic tools to characterize gene functions in plants. Recent advances in Next Generation Sequencing (NGS) technology have greatly reduced the time required for mutant identification in forward genetic screening. The major advantage of NGS enabled mapping, known as mapping-by-sequencing, is the simultaneous marker identification and genotyping and identification of the genomic loci causing phenotypes. We have been among the first to show that mapping-by-sequencing can be performed even within the same genetic background using mutagen-induced changes as segregating markers. As a proof of this concept, we mapped a previously unknown suppressor of like heterochromatin protein1 (lhp1) mutant. We developed a computational pipeline for the same and integrated it into an existing mapping-by-sequencing pipeline called SHOREmap. Though mapping-by-sequencing is now being routinely used, less effort has been put in optimizing the experimental set-up. Therefore, we developed new computational pipeline called Pop-Seq simulator that can simulate different mapping populations and sequencing experiments. It simulates recombinant genomes by following empirical determined recombination frequency and landscape, which make simulations close to reality. Using Pop-Seq simulator we simulated different mapping-by-sequencing scenarios and created guidelines for mapping-by-sequencing experiments in Arabidopsis. Although mapping-by-sequencing has already become a standard method in Arabidopsis, the application in crops is hindered by the large genome sizes and the lack of complete reference genomes. Therefore, we have used the Pop-Seq simulator to extend our analysis on the experimental design of mapping- by-sequencing to two crop model species, rice and barley, in which next generation sequencing-based mapping becomes tangible reality. Besides, we have developed a reference-free method called NIKS (needle in the k-stack) that enables mapping-by- sequencing in species without pre-assembled reference sequence, gene annotation, or genetic map. NIKS directly compares genomes using k-mers from whole genome sequencing data to identify homozygous mutations and extend the sequence associated with mutation site by local de novo assembly. We have used ab initio gene structural prediction to annotate the effect of mutations, which led us to the identification of causal mutation. This method will facilitate mapping-by-sequencing in non-model species

Genome expansion of Arabis alpina linked with retrotransposition and reduced symmetric DNA methylation

This document is the Accepted Manuscript version, made available in accordance to Springer Nature Terms of reuse of archived manuscripts.Despite evolutionary conserved mechanisms to silence transposable element activity, there are drastic differences in the abundance of transposable elements even among closely related plant species. We conducted a de novo assembly for the 375 .Mb genome of the perennial model plant, Arabis alpina. Analysing this genome revealed long-lasting and recent transposable element activity predominately driven by Gypsy long terminal repeat retrotransposons, which extended the low-recombining pericentromeres and transformed large formerly euchromatic regions into repeat-rich pericentromeric regions. This reduced capacity for long terminal repeat retrotransposon silencing and removal in A. alpina co-occurs with unexpectedly low levels of DNA methylation. Most remarkably, the striking reduction of symmetrical CG and CHG methylation suggests weakened DNA methylation maintenance in A. alpina compared with Arabidopsis thaliana. Phylogenetic analyses indicate a highly dynamic evolution of some components of methylation maintenance machinery that might be related to the unique methylation in A. alpina.Peer reviewe

Kicking against the PRCs - a domesticated transposase antagonises silencing mediated by polycomb group proteins and is an accessory component of polycomb repressive complex 2

The Polycomb group (PcG) and trithorax group (trxG) genes play crucial roles in development by regulating expression of homeotic and other genes controlling cell fate. Both groups catalyse modifications of chromatin, particularly histone methylation, leading to epigenetic changes that affect gene activity. The trxG antagonizes the function of PcG genes by activating PcG target genes, and consequently trxG mutants suppress PcG mutant phenotypes. We previously identified the ANTAGONIST OF LIKE HETEROCHROMATIN PROTEIN1 (ALP1) gene as a genetic suppressor of mutants in the Arabidopsis PcG gene LIKE HETEROCHROMATIN PROTEIN1 (LHP1). Here, we show that ALP1 interacts genetically with several other PcG and trxG components and that it antagonizes PcG silencing. Transcriptional profiling reveals that when PcG activity is compromised numerous target genes are hyper-activated in seedlings and that in most cases this requires ALP1. Furthermore, when PcG activity is present ALP1 is needed for full activation of several floral homeotic genes that are repressed by the PcG. Strikingly, ALP1 does not encode a known chromatin protein but rather a protein related to PIF/Harbinger class transposases. Phylogenetic analysis indicates that ALP1 is broadly conserved in land plants and likely lost transposase activity and acquired a novel function during angiosperm evolution. Consistent with this, immunoprecipitation and mass spectrometry (IP-MS) show that ALP1 associates, in vivo, with core components of POLYCOMB REPRESSIVE COMPLEX 2 (PRC2), a widely conserved PcG protein complex which functions as a H3K27me3 histone methyltransferase. Furthermore, in reciprocal pulldowns using the histone methyltransferase CURLY LEAF (CLF), we identify not only ALP1 and the core PRC2 components but also plant-specific accessory components including EMBRYONIC FLOWER 1 (EMF1), a transcriptional repressor previously associated with PRC1-like complexes. Taken together our data suggest that ALP1 inhibits PcG silencing by blocking the interaction of the core PRC2 with accessory components that promote its HMTase activity or its role in inhibiting transcription. ALP1 is the first example of a domesticated transposase acquiring a novel function as a PcG component. The antagonistic interaction of a modified transposase with the PcG machinery is novel and may have arisen as a means for the cognate transposon to evade host surveillance or for the host to exploit features of the transposition machinery beneficial for epigenetic regulation of gene activity.Fil: Liang, Shih Chieh. University of Edinburgh; Reino UnidoFil: Hartwig, Ben. Max Planck Institute for Plant Breeding Research; AlemaniaFil: Perera, Pumi. University of Edinburgh; Reino UnidoFil: Mora Garcia, Santiago. Fundación Instituto Leloir; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Instituto de Investigaciones Bioquimicas de Buenos Aires; ArgentinaFil: de Leau, Erica. University of Edinburgh; Reino UnidoFil: Thornton, Harry. University of Edinburgh; Reino UnidoFil: Lima de Alves, Flavia. University of Edinburgh; Reino UnidoFil: Rapsilber, Juri. University of Edinburgh; Reino UnidoFil: Yang, Suxin. University of Edinburgh; Reino UnidoFil: James, Geo Velikkakam. Max Planck Institute for Plant Breeding Research; AlemaniaFil: Schneeberger, Korbinian. Max Planck Institute for Plant Breeding Research; AlemaniaFil: Finnegan, E. Jean. University of Edinburgh; Reino UnidoFil: Turck, Franziska. Max Planck Institute for Plant Breeding Research; AlemaniaFil: Goodrich, Justin. Mc Gill University; Canad

Fast Isogenic Mapping-by-Sequencing of Ethyl Methanesulfonate-Induced Mutant Bulks

Mapping-by-sequencing (or SHOREmapping) has revitalized the powerful concept of forward genetic screens in plants. However, as in conventional genetic mapping approaches, mapping-by-sequencing requires phenotyping of mapping populations established from crosses between two diverged accessions. In addition to the segregation of the focal phenotype, this introduces natural phenotypic variation, which can interfere with the recognition of quantitative phenotypes. Here, we demonstrate how mapping-by-sequencing and candidate gene identification can be performed within the same genetic background using only mutagen-induced changes as segregating markers. Using a previously unknown suppressor of mutants of like heterochromatin protein1 (lhp1), which in its functional form is involved in chromatin-mediated gene repression, we identified three closely linked ethyl methanesulfonate-induced changes as putative candidates. In order to assess allele frequency differences between such closely linked mutations, we introduced deep candidate resequencing using the new Ion Torrent Personal Genome Machine sequencing platform to our mutant identification pipeline and thereby reduced the number of causal candidate mutations to only one. Genetic analysis of two independent additional alleles confirmed that this mutation was causal for the suppression of lhp1

Chromosome-level assembly of Arabidopsis thaliana Ler reveals the extent of translocation and inversion polymorphisms

Resequencing or reference-based assemblies reveal large parts of the small-scale sequence variation. However, they typically fail to separate such local variation into colinear and rearranged variation, because they usually do not recover the complement of large-scale rearrangements, including transpositions and inversions. Besides the availability of hundreds of genomes of diverse Arabidopsis thaliana accessions, there is so far only one full-length assembled genome: the reference sequence. We have assembled 117 Mb of the A. thaliana Landsberg erecta (Ler) genome into five chromosome-equivalent sequences using a combination of short Illumina reads, long PacBio reads, and linkage information. Whole-genome comparison against the reference sequence revealed 564 transpositions and 47 inversions comprising ∼3.6 Mb, in addition to 4.1 Mb of nonreference sequence, mostly originating from duplications. Although rearranged regions are not different in local divergence from colinear regions, they are drastically depleted for meiotic recombination in heterozygotes. Using a 1.2-Mb inversion as an example, we show that such rearrangement-mediated reduction of meiotic recombination can lead to genetically isolated haplotypes in the worldwide population of A. thaliana Moreover, we found 105 single-copy genes, which were only present in the reference sequence or the Ler assembly, and 334 single-copy orthologs, which showed an additional copy in only one of the genomes. To our knowledge, this work gives first insights into the degree and type of variation, which will be revealed once complete assemblies will replace resequencing or other reference-dependent methods.Knowledge of the exact distribution of meiotic crossovers (COs) and gene conversions (GCs) is essential for understanding many aspects of population genetics and evolution, from haplotype structure and long-distance genetic linkage to the generation of new allelic variants of genes. To this end, we resequenced the four products of 13 meiotic tetrads along with 10 doubled haploids derived from Arabidopsis thaliana hybrids. GC detection through short reads has previously been confounded by genomic rearrangements. Rigid filtering for misaligned reads allowed GC identification at high accuracy and revealed an ∼80-kb transposition, which undergoes copy-number changes mediated by meiotic recombination. Non-crossover associated GCs were extremely rare most likely due to their short average length of ∼25-50 bp, which is significantly shorter than the length of CO-associated GCs. Overall, recombination preferentially targeted non-methylated nucleosome-free regions at gene promoters, which showed significant enrichment of two sequence motifs.This work was supported by Spanish Ministry of Economy and Competitiveness Centro de Excelencia Severo Ochoa 2013-2017 Grant SEV-2012-0208. L.Z. was supported by the International PhD scholarship program of La Caixa at CR

Arabidopsis thaliana DM2h (R8) within the Landsberg RPP1-like Resistance Locus Underlies Three Different Cases of EDS1-Conditioned Autoimmunity

Plants have a large panel of nucleotide-binding/leucine rich repeat (NLR) immune receptors which monitor host interference by diverse pathogen molecules (effectors) and trigger disease resistance pathways. NLR receptor systems are necessarily under tight control to mitigate the trade-off between induced defenses and growth. Hence, mis-regulated NLRs often cause autoimmunity associated with stunting and, in severe cases, necrosis. Nucleocytoplasmic ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) is indispensable for effector-triggered and autoimmune responses governed by a family of Toll-Interleukin1-Receptor-related NLR receptors (TNLs). EDS1 operates coincidently or immediately downstream of TNL activation to transcriptionally reprogram cells for defense. We show here that low levels of nuclear-enforced EDS1 are sufficient for pathogen resistance in Arabidopsis thaliana, without causing negative effects. Plants expressing higher nuclear EDS1 amounts have the genetic, phenotypic and transcriptional hallmarks of TNL autoimmunity. In a screen for genetic suppressors of nuclear EDS1 autoimmunity, we map multiple, independent mutations to one gene, DM2h, lying within the polymorphic DANGEROUS MIX2 cluster of TNL RPP1-like genes from A. thaliana accession Landsberg erecta (Ler). The DM2 locus is a known hotspot for deleterious epistatic interactions leading to immune-related incompatibilities between A. thaliana natural accessions. We find that DM2hLer underlies two further genetic incompatibilities involving the RPP1-likeLer locus and EDS1. We conclude that the DM2hLer TNL protein and nuclear EDS1 cooperate, directly or indirectly, to drive cells into an immune response at the expense of growth. A further conclusion is that regulating the available EDS1 nuclear pool is fundamental for maintaining homeostatic control of TNL immune pathways