Graham J. Etherington

Genome analyses of an aggressive and invasive lineage of the Irish potato famine pathogen.

Pest and pathogen losses jeopardise global food security and ever since the 19th century Irish famine, potato late blight has exemplified this threat. The causal oomycete pathogen, Phytophthora infestans, undergoes major population shifts in agricultural systems via the successive emergence and migration of asexual lineages. The phenotypic and genotypic bases of these selective sweeps are largely unknown but management strategies need to adapt to reflect the changing pathogen population. Here, we used molecular markers to document the emergence of a lineage, termed 13_A2, in the European P. infestans population, and its rapid displacement of other lineages to exceed 75% of the pathogen population across Great Britain in less than three years. We show that isolates of the 13_A2 lineage are among the most aggressive on cultivated potatoes, outcompete other aggressive lineages in the field, and overcome previously effective forms of plant host resistance. Genome analyses of a 13_A2 isolate revealed extensive genetic and expression polymorphisms particularly in effector genes. Copy number variations, gene gains and losses, amino-acid replacements and changes in expression patterns of disease effector genes within the 13_A2 isolate likely contribute to enhanced virulence and aggressiveness to drive this population displacement. Importantly, 13_A2 isolates carry intact and in planta induced Avrblb1, Avrblb2 and Avrvnt1 effector genes that trigger resistance in potato lines carrying the corresponding R immune receptor genes Rpi-blb1, Rpi-blb2, and Rpi-vnt1.1. These findings point towards a strategy for deploying genetic resistance to mitigate the impact of the 13_A2 lineage and illustrate how pathogen population monitoring, combined with genome analysis, informs the management of devastating disease epidemics.

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

BackgroundThe 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.ResultsAs 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.ConclusionsBy 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.

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

Summary 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.

Recombination Analysis Tool (RAT): a program for the high-throughput detection of recombination

MOTIVATION Recombination can be a prevailing drive in shaping genome evolution. RAT (Recombination Analysis Tool) is a Java-based tool for investigating recombination events in any number of aligned sequences (protein or DNA) of any length (short viral sequences to full genomes). It is an uncomplicated and intuitive application and allows the user to view only the regions of sequence alignments they are interested in. RESULTS RAT was applied to viral sequences. Its utility was demonstrated through the detection of a known recombinant of HIV and a detailed analysis of Noroviruses, the most common cause of viral gastroenteritis in humans. AVAILABILITY RAT, along with a users guide, is freely available from http://jic-bioinfo.bbsrc.ac.uk/bioinformatics-research/staff/graham_etherington/RAT.htm.

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

BackgroundThe 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.ResultUsing 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.ConclusionRenSeq 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.

From pathogen genomes to host plant processes: the power of plant parasitic oomycetes.

Recent pathogenomic research on plant parasitic oomycete effector function and plant host responses has resulted in major conceptual advances in plant pathology, which has been possible thanks to the availability of genome sequences.

Transcriptome analysis of the sulfate deficiency response in the marine microalga Emiliania huxleyi

The response to sulfate deficiency of plants and freshwater green algae has been extensively analysed by system biology approaches. By contrast, seawater sulfate concentration is high and very little is known about the sulfur metabolism of marine organisms. Here, we used a combination of metabolite analysis and transcriptomics to analyse the response of the marine microalga Emiliania huxleyi as it acclimated to sulfate limitation. Lowering sulfate availability in artificial seawater from 25 to 5 mM resulted in significant reduction in growth and intracellular concentrations of dimethylsulfoniopropionate and glutathione. Sulfate-limited E. huxleyi cells showed increased sulfate uptake but sulfate reduction to sulfite did not seem to be regulated. Sulfate limitation in E. huxleyi affected expression of 1718 genes. The vast majority of these genes were upregulated, including genes involved in carbohydrate and lipid metabolism, and genes involved in the general stress response. The acclimation response of E. huxleyi to sulfate deficiency shows several similarities to the well-described responses of Arabidopsis and Chlamydomonas, but also has many unique features. This dataset shows that even though E. huxleyi is adapted to constitutively high sulfate concentration, it retains the ability to re-program its gene expression in response to reduced sulfate availability.

Investigation of the potential for triterpene synthesis in rice through genome mining and metabolic engineering

The first committed step in sterol biosynthesis in plants involves the cyclization of 2,3-oxidosqualene by the oxidosqualene cyclase (OSC) enzyme cycloartenol synthase. 2,3-Oxidosqualene is also a precursor for triterpene synthesis. Antimicrobial triterpenes are common in dicots, but seldom found in monocots, with the notable exception of oat. Here, through genome mining and metabolic engineering, we investigate the potential for triterpene synthesis in rice. The first two steps in the oat triterpene pathway are catalysed by a divergent OSC (AsbAS1) and a cytochrome P450 (CYP51). The genes for these enzymes form part of a metabolic gene cluster. To investigate the origins of triterpene synthesis in monocots, we analysed systematically the OSC and CYP51 gene families in rice. We also engineered rice for elevated triterpene content. We discovered a total of 12 OSC and 12 CYP51 genes in rice and uncovered key events in the evolution of triterpene synthesis. We further showed that the expression of AsbAS1 in rice leads to the accumulation of the simple triterpene, β-amyrin. These findings provide new insights into the evolution of triterpene synthesis in monocots and open up opportunities for metabolic engineering for disease resistance in rice and other cereals.

EXPRSS: an Illumina based high-throughput expression-profiling method to reveal transcriptional dynamics

BackgroundNext Generation Sequencing technologies have facilitated differential gene expression analysis through RNA-seq and Tag-seq methods. RNA-seq has biases associated with transcript lengths, lacks uniform coverage of regions in mRNA and requires 10–20 times more reads than a typical Tag-seq. Most existing Tag-seq methods either have biases or not high throughput due to use of restriction enzymes or enzymatic manipulation of 5’ ends of mRNA or use of RNA ligations.ResultsWe have developed EXpression Profiling through Randomly Sheared cDNA tag Sequencing (EXPRSS) that employs acoustic waves to randomly shear cDNA and generate sequence tags at a relatively defined position (~150-200 bp) from the 3′ end of each mRNA. Implementation of the method was verified through comparative analysis of expression data generated from EXPRSS, NlaIII-DGE and Affymetrix microarray and through qPCR quantification of selected genes. EXPRSS is a strand specific and restriction enzyme independent tag sequencing method that does not require cDNA length-based data transformations. EXPRSS is highly reproducible, is high-throughput and it also reveals alternative polyadenylation and polyadenylated antisense transcripts. It is cost-effective using barcoded multiplexing, avoids the biases of existing SAGE and derivative methods and can reveal polyadenylation position from paired-end sequencing.ConclusionsEXPRSS Tag-seq provides sensitive and reliable gene expression data and enables high-throughput expression profiling with relatively simple downstream analysis.

Arabidopsis late blight: Infection of a nonhost plant by Albugo laibachii enables full colonization by Phytophthora infestans

The oomycete pathogen Phytophthora infestans causes potato late blight, and as a potato and tomato specialist pathogen, is seemingly poorly adapted to infect plants outside the Solanaceae. Here, we report the unexpected finding that P. infestans can infect Arabidopsis thaliana when another oomycete pathogen, Albugo laibachii, has colonized the host plant. The behaviour and speed of P. infestans infection in Arabidopsis pre‐infected with A. laibachii resemble P. infestans infection of susceptible potato plants. Transcriptional profiling of P. infestans genes during infection revealed a significant overlap in the sets of secreted‐protein genes that are induced in P. infestans upon colonization of potato and susceptible Arabidopsis, suggesting major similarities in P. infestans gene expression dynamics on the two plant species. Furthermore, we found haustoria of A. laibachii and P. infestans within the same Arabidopsis cells. This Arabidopsis—A. laibachii—P. infestans tripartite interaction opens up various possibilities to dissect the molecular mechanisms of P. infestans infection and the processes occurring in co‐infected Arabidopsis cells.