Epigenomic Diversity in a Global Collection of Arabidopsis thaliana Accessions.

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

Genome and time-of-day transcriptome of Wolffia australiana link morphological minimization with gene loss and less growth control.

Rootless plants in the genus Wolffia are some of the fastest growing known plants on Earth. Wolffia have a reduced body plan, primarily multiplying through a budding type of asexual reproduction. Here, we generated draft reference genomes for Wolffia australiana (Benth.) Hartog & Plas, which has the smallest genome size in the genus at 357 Mb and has a reduced set of predicted protein-coding genes at about 15,000. Comparison between multiple high-quality draft genome sequences from W. australiana clones confirmed loss of several hundred genes that are highly conserved among flowering plants, including genes involved in root developmental and light signaling pathways. Wolffia has also lost most of the conserved nucleotide-binding leucine-rich repeat (NLR) genes that are known to be involved in innate immunity, as well as those involved in terpene biosynthesis, while having a significant overrepresentation of genes in the sphingolipid pathways that may signify an alternative defense system. Diurnal expression analysis revealed that only 13% of Wolffia genes are expressed in a time-of-day (TOD) fashion, which is less than the typical ∼40% found in several model plants under the same condition. In contrast to the model plants Arabidopsis and rice, many of the pathways associated with multicellular and developmental processes are not under TOD control in W. australiana, where genes that cycle the conditions tested predominantly have carbon processing and chloroplast-related functions. The Wolffia genome and TOD expression data set thus provide insight into the interplay between a streamlined plant body plan and optimized growth.

The complex architecture and epigenomic impact of plant T-DNA insertions.

The bacterium Agrobacterium tumefaciens has been the workhorse in plant genome engineering. Customized replacement of native tumor-inducing (Ti) plasmid elements enabled insertion of a sequence of interest called Transfer-DNA (T-DNA) into any plant genome. Although these transfer mechanisms are well understood, detailed understanding of structure and epigenomic status of insertion events was limited by current technologies. Here we applied two single-molecule technologies and analyzed Arabidopsis thaliana lines from three widely used T-DNA insertion collections (SALK, SAIL and WISC). Optical maps for four randomly selected T-DNA lines revealed between one and seven insertions/rearrangements, and the length of individual insertions from 27 to 236 kilobases. De novo nanopore sequencing-based assemblies for two segregating lines partially resolved T-DNA structures and revealed multiple translocations and exchange of chromosome arm ends. For the current TAIR10 reference genome, nanopore contigs corrected 83% of non-centromeric misassemblies. The unprecedented contiguous nucleotide-level resolution enabled an in-depth study of the epigenome at T-DNA insertion sites. SALK_059379 line T-DNA insertions were enriched for 24nt small interfering RNAs (siRNA) and dense cytosine DNA methylation, resulting in transgene silencing via the RNA-directed DNA methylation pathway. In contrast, SAIL_232 line T-DNA insertions are predominantly targeted by 21/22nt siRNAs, with DNA methylation and silencing limited to a reporter, but not the resistance gene. Additionally, we profiled the H3K4me3, H3K27me3 and H2A.Z chromatin environments around T-DNA insertions using ChIP-seq in SALK_059379, SAIL_232 and five additional T-DNA lines. We discovered various effect s ranging from complete loss of chromatin marks to the de novo incorporation of H2A.Z and trimethylation of H3K4 and H3K27 around the T-DNA integration sites. This study provides new insights into the structural impact of inserting foreign fragments into plant genomes and demonstrates the utility of state-of-the-art long-range sequencing technologies to rapidly identify unanticipated genomic changes.

Albugo candida race diversity, ploidy and host-associated microbes revealed using DNA sequence capture on diseased plants in the field.

Physiological races of the oomycete Albugo candida are biotrophic pathogens of diverse plant species, primarily the Brassicaceae, and cause infections that suppress host immunity to other pathogens. However, A. candida race diversity and the consequences of host immunosuppression are poorly understood in the field. We report a method that enables sequencing of DNA of plant pathogens and plant-associated microbes directly from field samples (Pathogen Enrichment Sequencing: PenSeq). We apply this method to explore race diversity in A. candida and to detect A. candida-associated microbes in the field (91 A. candida-infected plants). We show with unprecedented resolution that each host plant species supports colonization by one of 17 distinct phylogenetic lineages, each with an unique repertoire of effector candidate alleles. These data reveal the crucial role of sexual and asexual reproduction, polyploidy and host domestication in A. candida specialization on distinct plant species. Our bait design also enabled phylogenetic assignment of DNA sequences from bacteria and fungi from plants in the field. This paper shows that targeted sequencing has a great potential for the study of pathogen populations while they are colonizing their hosts. This method could be applied to other microbes, especially to those that cannot be cultured.

Identification and rapid mapping of a gene conferring broad-spectrum late blight resistance in the diploid potato species Solanum verrucosum through DNA capture technologies.

KEY MESSAGE: A broad-spectrum late blight disease-resistance gene from Solanum verrucosum has been mapped to potato chromosome 9. The gene is distinct from previously identified-resistance genes. We have identified and characterised a broad-spectrum resistance to Phytophthora infestans from the wild Mexican species Solanum verrucosum. Diagnostic resistance gene enrichment (dRenSeq) revealed that the resistance is not conferred by previously identified nucleotide-binding, leucine-rich repeat genes. Utilising the sequenced potato genome as a reference, two complementary enrichment strategies that target resistance genes (RenSeq) and single/low-copy number genes (Generic-mapping enrichment Sequencing; GenSeq), respectively, were deployed for the rapid, SNP-based mapping of the resistance through bulked-segregant analysis. Both approaches independently positioned the resistance, referred to as Rpi-ver1, to the distal end of potato chromosome 9. Stringent post-enrichment read filtering identified a total of 64 informative SNPs that corresponded to the expected ratio for significant polymorphisms in the parents as well as the bulks. Of these, 61 SNPs are located on potato chromosome 9 and reside within 27 individual genes, which in the sequenced potato clone DM locate to positions 45.9 to 60.9 Mb. RenSeq- and GenSeq-derived SNPs within the target region were converted into allele-specific PCR-based KASP markers and further defined the position of the resistance to a 4.3 Mb interval at the bottom end of chromosome 9 between positions 52.62-56.98 Mb.

NLR-parser: rapid annotation of plant NLR complements.

MOTIVATION: The repetitive nature of plant disease resistance genes encoding for nucleotide-binding leucine-rich repeat (NLR) proteins hampers their prediction with standard gene annotation software. Motif alignment and search tool (MAST) has previously been reported as a tool to support annotation of NLR-encoding genes. However, the decision if a motif combination represents an NLR protein was entirely manual. RESULTS: The NLR-parser pipeline is designed to use the MAST output from six-frame translated amino acid sequences and filters for predefined biologically curated motif compositions. Input reads can be derived from, for example, raw long-read sequencing data or contigs and scaffolds coming from plant genome projects. The output is a tab-separated file with information on start and frame of the first NLR specific motif, whether the identified sequence is a TNL or CNL, potentially full or fragmented. In addition, the output of the NB-ARC domain sequence can directly be used for phylogenetic analyses. In comparison to other prediction software, the highly complex NB-ARC domain is described in detail using several individual motifs.

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB-LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations.

RenSeq is a NB-LRR (nucleotide binding-site leucine-rich repeat) gene-targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB-LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB-LRRs and can be accessed through a genome browser that we provide. We compared these NB-LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum 'Heinz 1706' extended the NB-LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co-segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi-ber2) and S. ruiz-ceballosii (Rpi-rzc1), we were able to apply RenSeq successfully to identify markers that co-segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy-to-adapt Galaxy pipelines.

Defining the full tomato NB-LRR resistance gene repertoire using genomic and cDNA RenSeq.

BACKGROUND: The availability of draft crop plant genomes allows the prediction of the full complement of genes that encode NB-LRR resistance gene homologs, enabling a more targeted breeding for disease resistance. Recently, we developed the RenSeq method to reannotate the full NB-LRR gene complement in potato and to identify novel sequences that were not picked up by the automated gene prediction software. Here, we established RenSeq on the reference genome of tomato (Solanum lycopersicum) Heinz 1706, using 260 previously identified NB-LRR genes in an updated Solanaceae RenSeq bait library. RESULT: Using 250-bp MiSeq reads after RenSeq on genomic DNA of Heinz 1706, we identified 105 novel NB-LRR sequences. Reannotation included the splitting of gene models, combination of partial genes to a longer sequence and closing of assembly gaps. Within the draft S. pimpinellifolium LA1589 genome, RenSeq enabled the annotation of 355 NB-LRR genes. The majority of these are however fragmented, with 5'- and 3'-end located on the edges of separate contigs. Phylogenetic analyses show a high conservation of all NB-LRR classes between Heinz 1706, LA1589 and the potato clone DM, suggesting that all sub-families were already present in the last common ancestor. A phylogenetic comparison to the Arabidopsis thaliana NB-LRR complement verifies the high conservation of the more ancient CCRPW8-type NB-LRRs. Use of RenSeq on cDNA from uninfected and late blight-infected tomato leaves allows the avoidance of sequence analysis of non-expressed paralogues. CONCLUSION: RenSeq is a promising method to facilitate analysis of plant resistance gene complements. The reannotated tomato NB-LRR complements, phylogenetic relationships and chromosomal locations provided in this paper will provide breeders and scientists with a useful tool to identify novel disease resistance traits. cDNA RenSeq enables for the first time next-gen sequencing approaches targeted to this very low-expressed gene family without the need for normalization.

Identification and localisation of the NB-LRR gene family within the potato genome.

BACKGROUND: The potato genome sequence derived from the Solanum tuberosum Group Phureja clone DM1-3 516 R44 provides unparalleled insight into the genome composition and organisation of this important crop. A key class of genes that comprises the vast majority of plant resistance (R) genes contains a nucleotide-binding and leucine-rich repeat domain, and is collectively known as NB-LRRs. RESULTS: As part of an effort to accelerate the process of functional R gene isolation, we performed an amino acid motif based search of the annotated potato genome and identified 438 NB-LRR type genes among the ~39,000 potato gene models. Of the predicted genes, 77 contain an N-terminal toll/interleukin 1 receptor (TIR)-like domain, and 107 of the remaining 361 non-TIR genes contain an N-terminal coiled-coil (CC) domain. Physical map positions were established for 370 predicted NB-LRR genes across all 12 potato chromosomes. The majority of NB-LRRs are physically organised within 63 identified clusters, of which 50 are homogeneous in that they contain NB-LRRs derived from a recent common ancestor. CONCLUSIONS: By establishing the phylogenetic and positional relationship of potato NB-LRRs, our analysis offers significant insight into the evolution of potato R genes. Furthermore, the data provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from Solanum species.

A complex resistance locus in Solanum americanum recognizes a conserved Phytophthora effector.

Late blight caused by Phytophthora infestans greatly constrains potato production. Many Resistance (R) genes were cloned from wild Solanum species and/or introduced into potato cultivars by breeding. However, individual R genes have been overcome by P. infestans evolution; durable resistance remains elusive. We positionally cloned a new R gene, Rpi-amr1, from Solanum americanum, that encodes an NRC helper-dependent CC-NLR protein. Rpi-amr1 confers resistance in potato to all 19 P. infestans isolates tested. Using association genomics and long-read RenSeq, we defined eight additional Rpi-amr1 alleles from different S. americanum and related species. Despite only ~90% identity between Rpi-amr1 proteins, all confer late blight resistance but differentially recognize Avramr1 orthologues and paralogues. We propose that Rpi-amr1 gene family diversity assists detection of diverse paralogues and alleles of the recognized effector, facilitating durable resistance against P. infestans.

Identification of Avramr1 from Phytophthora infestans using long read and cDNA pathogen-enrichment sequencing (PenSeq).

Potato late blight, caused by the oomycete pathogen Phytophthora infestans, significantly hampers potato production. Recently, a new Resistance to Phytophthora infestans (Rpi) gene, Rpi-amr1, was cloned from a wild Solanum species, Solanum americanum. Identification of the corresponding recognized effector (Avirulence or Avr) genes from P. infestans is key to elucidating their naturally occurring sequence variation, which in turn informs the potential durability of the cognate late blight resistance. To identify the P. infestans effector recognized by Rpi-amr1, we screened available RXLR effector libraries and used long read and cDNA pathogen-enrichment sequencing (PenSeq) on four P. infestans isolates to explore the untested effectors. Using single-molecule real-time sequencing (SMRT) and cDNA PenSeq, we identified 47 highly expressed effectors from P. infestans, including PITG_07569, which triggers a highly specific cell death response when transiently coexpressed with Rpi-amr1 in Nicotiana benthamiana, suggesting that PITG_07569 is Avramr1. Here we demonstrate that long read and cDNA PenSeq enables the identification of full-length RXLR effector families and their expression profile. This study has revealed key insights into the evolution and polymorphism of a complex RXLR effector family that is associated with the recognition by Rpi-amr1.

Epigenetic silencing of a multifunctional plant stress regulator.

The central regulator of the ethylene (ET) signaling pathway, which controls a plethora of developmental programs and responses to environmental cues in plants, is ETHYLENE-INSENSITIVE2 (EIN2). Here we identify a chromatin-dependent regulatory mechanism at EIN2 requiring two genes: ETHYLENE-INSENSITIVE6 (EIN6), which is a H3K27me3 demethylase also known as RELATIVE OF EARLY FLOWERING6 (REF6), and EIN6 ENHANCER (EEN), the Arabidopsis homolog of the yeast INO80 chromatin remodeling complex subunit IES6 (INO EIGHTY SUBUNIT). Strikingly, EIN6 (REF6) and the INO80 complex redundantly control the level and the localization of the repressive histone modification H3K27me3 and the histone variant H2A.Z at the 5' untranslated region (5'UTR) intron of EIN2. Concomitant loss of EIN6 (REF6) and the INO80 complex shifts the chromatin landscape at EIN2 to a repressive state causing a dramatic reduction of EIN2 expression. These results uncover a unique type of chromatin regulation which safeguards the expression of an essential multifunctional plant stress regulator.

High contiguity Arabidopsis thaliana genome assembly with a single nanopore flow cell.

The handheld Oxford Nanopore MinION sequencer generates ultra-long reads with minimal cost and time requirements, which makes sequencing genomes at the bench feasible. Here, we sequence the gold standard Arabidopsis thaliana genome (KBS-Mac-74 accession) on the bench with the MinION sequencer, and assemble the genome using typical consumer computing hardware (4 Cores, 16 Gb RAM) into chromosome arms (62 contigs with an N50 length of 12.3 Mb). We validate the contiguity and quality of the assembly with two independent single-molecule technologies, Bionano optical genome maps and Pacific Biosciences Sequel sequencing. The new A. thaliana KBS-Mac-74 genome enables resolution of a quantitative trait locus that had previously been recalcitrant to a Sanger-based BAC sequencing approach. In summary, we demonstrate that even when the purpose is to understand complex structural variation at a single region of the genome, complete genome assembly is becoming the simplest way to achieve this goal.

Accelerated cloning of a potato late blight-resistance gene using RenSeq and SMRT sequencing.

Global yields of potato and tomato crops have fallen owing to potato late blight disease, which is caused by Phytophthora infestans. Although most commercial potato varieties are susceptible to blight, many wild potato relatives show variation for resistance and are therefore a potential source of Resistance to P. infestans (Rpi) genes. Resistance breeding has exploited Rpi genes from closely related tuber-bearing potato relatives, but is laborious and slow. Here we report that the wild, diploid non-tuber-bearing Solanum americanum harbors multiple Rpi genes. We combine resistance (R) gene sequence capture (RenSeq) with single-molecule real-time (SMRT) sequencing (SMRT RenSeq) to clone Rpi-amr3i. This technology should enable de novo assembly of complete nucleotide-binding, leucine-rich repeat receptor (NLR) genes, their regulatory elements and complex multi-NLR loci from uncharacterized germplasm. SMRT RenSeq can be applied to rapidly clone multiple R genes for engineering pathogen-resistant crops.

Targeted capture and sequencing of gene-sized DNA molecules.

Targeted capture provides an efficient and sensitive means for sequencing specific genomic regions in a high-throughput manner. To date, this method has mostly been used to capture exons from the genome (the exome) using short insert libraries and short-read sequencing technology, enabling the identification of genetic variants or new members of large gene families. Sequencing larger molecules results in the capture of whole genes, including intronic and intergenic sequences that are typically more polymorphic and allow the resolution of the gene structure of homologous genes, which are often clustered together on the chromosome. Here, we describe an improved method for the capture and single-molecule sequencing of DNA molecules as large as 7 kb by means of size selection and optimized PCR conditions. Our approach can be used to capture, sequence, and distinguish between similar members of the NB-LRR gene family-key genes in plant immune systems.

Hypersensitive response to Potato virus Y in potato cultivar Sárpo Mira is conferred by the Ny-Smira gene located on the long arm of chromosome IX.

Potato virus Y (PVY, Potyvirus) is the fifth most important plant virus worldwide in terms of economic and scientific impact. It infects members of the family Solanaceae and causes losses in potato, tomato, tobacco, pepper and petunia production. In potato and its wild relatives, two types of resistance genes against PVY have been identified. While Ry genes confer symptomless extreme resistance, Ny genes cause a hypersensitive response visible as local necrosis that may also be able to prevent the virus from spreading under certain environmental conditions. The potato cultivar Sárpo Mira originates from Hungary and is highly resistant to PVY, although the source of this resistance remains unknown. We show that cv. Sárpo Mira reacts with a hypersensitive response leading to necrosis after PVYNTN infection in detached leaf, whole plant and grafting assays. The hypersensitivity to PVYNTN segregated amongst 140 individuals of tetraploid progeny of cvs. Sárpo Mira × Maris Piper in a 1:1 ratio, indicating that it was conferred by a single, dominant gene in simplex. Moreover, we identified five DNA markers linked to this trait and located the underlying locus (Ny-Smira) to the long arm of potato chromosome IX. This position corresponds to the location of the Rychc and Ny-1 genes for PVY resistance. A simple PCR marker, located 1 cM from the Ny-Smira gene, can be recommended for selection of PVY-resistant progeny of cv. Sárpo Mira.

Genomic DNA library preparation for resistance gene enrichment and sequencing (RenSeq) in plants.

Enrichment of genomic DNA for genes of interest prior to high-throughput sequencing offers an efficient and cost-effective approach to reduce genome complexity. Target enrichment typically yields higher read-depth for selected genes and is therefore suitable for determination of sequence polymorphisms and enables multiplexing of samples. Target enrichment also provides a means to annotate specific gene families within the sequenced organisms without the requirements for gene models. Here we describe enrichment procedures for NB-LRR-type plant resistance genes that can, for example, be used to establish the NB-LRR gene complements of individual plants and to map resistances more rapidly using a bulked segregant analysis.

Elevating crop disease resistance with cloned genes.

Essentially all plant species exhibit heritable genetic variation for resistance to a variety of plant diseases caused by fungi, bacteria, oomycetes or viruses. Disease losses in crop monocultures are already significant, and would be greater but for applications of disease-controlling agrichemicals. For sustainable intensification of crop production, we argue that disease control should as far as possible be achieved using genetics rather than using costly recurrent chemical sprays. The latter imply CO2 emissions from diesel fuel and potential soil compaction from tractor journeys. Great progress has been made in the past 25 years in our understanding of the molecular basis of plant disease resistance mechanisms, and of how pathogens circumvent them. These insights can inform more sophisticated approaches to elevating disease resistance in crops that help us tip the evolutionary balance in favour of the crop and away from the pathogen. We illustrate this theme with an account of a genetically modified (GM) blight-resistant potato trial in Norwich, using the Rpi-vnt1.1 gene isolated from a wild relative of potato, Solanum venturii, and introduced by GM methods into the potato variety Desiree.