Publications

Bassel GW, Fung P, Chow TF, Foong JA, Provart NJ, Cutler SR

Plant Physiol. 2008 May;147(1):143-55

PubMed PMID: 18359847

Abstract

The transition from seed to seedling is mediated by germination, a complex process that starts with imbibition and completes with radicle emergence. To gain insight into the transcriptional program mediating germination, previous studies have compared the transcript profiles of dry, dormant, and germinating after-ripened Arabidopsis (Arabidopsis thaliana) seeds. While informative, these approaches did not distinguish the transcriptional responses due to imbibition, shifts in metabolism, or breaking of dormancy from those triggered by the initiation of germination. In this study, three mechanistically distinct small molecules that inhibit Arabidopsis seed germination (methotrexate, 2, 4-dinitrophenol, and cycloheximide) were identified using a small-molecule screen and used to probe the germination transcriptome. Germination-responsive transcripts were defined as those with significantly altered transcript abundance across all inhibitory treatments with respect to control germinating seeds, using data from ATH1 microarrays. This analysis identified numerous germination regulators as germination responsive, including the DELLA proteins GAI, RGA, and RGL3, the abscisic acid-insensitive proteins ABI4, ABI5, ABI8, and FRY1, and the gibberellin receptor GID1A. To help visualize these and other publicly available seed microarray data, we designed a seed mRNA expression browser using the electronic Fluorescent Pictograph platform. An overall decrease in gene expression and a 5-fold greater number of transcripts identified as statistically down-regulated in drug-inhibited seeds point to a role for mRNA degradation or turnover during seed germination. The genes identified in our study as responsive to germination define potential uncharacterized regulators of this process and provide a refined transcriptional signature for germinating Arabidopsis seeds.

Tan J, Pitout JD, Guttman DS

J. Clin. Microbiol. 2008 May;46(5):1870-2

PubMed PMID: 18337388

Abstract

We present an imipenem lysate metallo-beta-lactamase (MBL) functional assay. This assay eliminates false-positive results due to the bactericidal effects of EDTA, can be performed with inexpensive reagents available in most laboratories, and is as accurate as the MBL Etest. It is appropriate for both high-accuracy screens and laboratories in developing countries with limited resources.

Lewis JD, Abada W, Ma W, Guttman DS, Desveaux D

J. Bacteriol. 2008 Apr;190(8):2880-91

PubMed PMID: 18263728

Abstract

Pseudomonas syringae utilizes the type III secretion system to translocate effector proteins into plant cells, where they can contribute to the pathogen’s ability to infect and cause disease. Recognition of these effectors by resistance proteins induces defense responses that typically include a programmed cell death reaction called the hypersensitive response. The YopJ/HopZ family of type III effector proteins is a common family of effector proteins found in animal- and plant-pathogenic bacteria. The HopZ family in P. syringae includes HopZ1a(PsyA2), HopZ1b(PgyUnB647), HopZ1c(PmaE54326), HopZ2(Ppi895A) and HopZ3(PsyB728a). HopZ1a is predicted to be most similar to the ancestral hopZ allele and causes a hypersensitive response in multiple plant species, including Arabidopsis thaliana. Therefore, it has been proposed that host defense responses have driven the diversification of this effector family. In this study, we further characterized the hypersensitive response induced by HopZ1a and demonstrated that it is not dependent on known resistance genes. Further, we identified a novel virulence function for HopZ2 that requires the catalytic cysteine demonstrated to be required for protease activity. Sequence analysis of the HopZ family revealed the presence of a predicted myristoylation sequence in all members except HopZ3. We demonstrated that the myristoylation site is required for membrane localization of this effector family and contributes to the virulence and avirulence activities of HopZ2 and HopZ1a, respectively. This paper provides insight into the selective pressures driving virulence protein evolution by describing a detailed functional characterization of the diverse HopZ family of type III effectors with the model plant Arabidopsis.

Schreiber K, Ckurshumova W, Peek J, Desveaux D

Plant J. 2008 May;54(3):522-31

PubMed PMID: 18248597

Abstract

The study of plant pathogenesis and the development of effective treatments to protect plants from diseases could be greatly facilitated by a high-throughput pathosystem to evaluate small-molecule libraries for inhibitors of pathogen virulence. The interaction between the Gram-negative bacterium Pseudomonas syringae and Arabidopsis thaliana is a model for plant pathogenesis. However, a robust high-throughput assay to score the outcome of this interaction is currently lacking. We demonstrate that Arabidopsis seedlings incubated with P. syringae in liquid culture display a macroscopically visible ‘bleaching’ symptom within 5 days of infection. Bleaching is associated with a loss of chlorophyll from cotyledonary tissues, and is correlated with bacterial virulence. Gene-for-gene resistance is absent in the liquid environment, possibly because of the suppression of the hypersensitive response under these conditions. Importantly, bleaching can be prevented by treating seedlings with known inducers of plant defence, such as salicylic acid (SA) or a basal defence-inducing peptide of bacterial flagellin (flg22) prior to inoculation. Based on these observations, we have devised a high-throughput liquid assay using standard 96-well plates to investigate the P. syringae-Arabidopsis interaction. An initial screen of small molecules active on Arabidopsis revealed a family of sulfanilamide compounds that afford protection against the bleaching symptom. The most active compound, sulfamethoxazole, also reduced in planta bacterial growth when applied to mature soil-grown plants. The whole-organism liquid assay provides a novel approach to probe chemical libraries in a high-throughput manner for compounds that reduce bacterial virulence in plants.

Cutter AD

Genetics 2008 Mar;178(3):1661-72

PubMed PMID: 18245859

Abstract

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.

Cutter AD

Mol. Biol. Evol. 2008 Apr;25(4):778-86

PubMed PMID: 18234705

Abstract

Accurate inference of the dates of common ancestry among species forms a central problem in understanding the evolutionary history of organisms. Molecular estimates of divergence time rely on the molecular evolutionary prediction that neutral mutations and substitutions occur at the same constant rate in genomes of related species. This underlies the notion of a molecular clock. Most implementations of this idea depend on paleontological calibration to infer dates of common ancestry, but taxa with poor fossil records must rely on external, potentially inappropriate, calibration with distantly related species. The classic biological models Caenorhabditis and Drosophila are examples of such problem taxa. Here, I illustrate internal calibration in these groups with direct estimates of the mutation rate from contemporary populations that are corrected for interfering effects of selection on the assumption of neutrality of substitutions. Divergence times are inferred among 6 species each of Caenorhabditis and Drosophila, based on thousands of orthologous groups of genes. I propose that the 2 closest known species of Caenorhabditis shared a common ancestor <24 MYA (Caenorhabditis briggsae and Caenorhabditis sp. 5) and that Caenorhabditis elegans diverged from its closest known relatives <30 MYA, assuming that these species pass through at least 6 generations per year; these estimates are much more recent than reported previously with molecular clock calibrations from non-nematode phyla. Dates inferred for the common ancestor of Drosophila melanogaster and Drosophila simulans are roughly concordant with previous studies. These revised dates have important implications for rates of genome evolution and the origin of self-fertilization in Caenorhabditis.

McCann HC, Guttman DS

New Phytol. 2008;177(1):33-47

PubMed PMID: 18078471

Abstract

Many bacterial plant pathogens require the type III secretion system (T3SS) and its effector proteins (T3SEs) to invade and extract nutrients from their hosts successfully. While the molecular function of this system is being studied intensively, we know comparatively little about the evolutionary and ecological pressures governing its fate over time, and even less about the detailed mechanisms underlying and driving complex T3SS-mediated coevolutionary dynamics. In this review we summarize our current understanding of how host-pathogen interactions evolve, with a particular focus on the T3SS of bacterial plant pathogens. We explore the evolutionary origins of the T3SS relative to the closely related flagellar system, and investigate the evolutionary pressures on this secretion and translocation apparatus. We examine the evolutionary forces acting on T3SEs, and compare the support for vertical descent with modification of these virulence-associated systems (pathoadaptation) vs horizontal gene transfer. We address the evolutionary origins of T3SEs from the perspective of both the evolutionary mechanisms that generate new effectors, and the mobile elements that may be the source of novel genetic material. Finally, we propose a number of questions raised by these studies, which may serve to guide our thinking about these complex processes.

Stavrinides J, McCann HC, Guttman DS

Cell. Microbiol. 2008 Feb;10(2):285-92

PubMed PMID: 18034865

Abstract

Many bacterial pathogens require a type III secretion system (T3SS) and suite of type III secreted effectors (T3SEs) to successfully colonize their hosts, extract nutrients and consequently cause disease. T3SEs, in particular, are key components of the bacterial arsenal, as they function directly inside the host to disrupt or suppress critical components of the defence network. The development of host defence and surveillance systems imposes intense selective pressures on these bacterial virulence factors, resulting in a host-pathogen co-evolutionary arms race. This arms race leaves its genetic signature in the pattern and structure of natural genetic variation found in T3SEs, thereby permitting us to infer the specific evolutionary processes and pressures driving these interactions. In this review, we summarize our current knowledge of T3SS-mediated host-pathogen co-evolution. We examine the evolution of the T3SS and the T3SEs that traverse it, in both plant and animal pathosystems, and discuss the processes that maintain these important pathogenicity determinants within pathogen populations. We go on to examine the possible origins of T3SEs, the mechanisms that give rise to new T3SEs and the processes that underlie their evolution.

Lan H, Carson R, Provart NJ, Bonner AJ

BMC Bioinformatics 2007;8:358

PubMed PMID: 17888165

Abstract

BACKGROUND: Arabidopsis thaliana is the model species of current plant genomic research with a genome size of 125 Mb and approximately 28,000 genes. The function of half of these genes is currently unknown. The purpose of this study is to infer gene function in Arabidopsis using machine-learning algorithms applied to large-scale gene expression data sets, with the goal of identifying genes that are potentially involved in plant response to abiotic stress.

RESULTS: Using in house and publicly available data, we assembled a large set of gene expression measurements for A. thaliana. Using those genes of known function, we first evaluated and compared the ability of basic machine-learning algorithms to predict which genes respond to stress. Predictive accuracy was measured using ROC50 and precision curves derived through cross validation. To improve accuracy, we developed a method for combining these classifiers using a weighted-voting scheme. The combined classifier was then trained on genes of known function and applied to genes of unknown function, identifying genes that potentially respond to stress. Visual evidence corroborating the predictions was obtained using electronic Northern analysis. Three of the predicted genes were chosen for biological validation. Gene knockout experiments confirmed that all three are involved in a variety of stress responses. The biological analysis of one of these genes (At1g16850) is presented here, where it is shown to be necessary for the normal response to temperature and NaCl.

CONCLUSION: Supervised learning methods applied to large-scale gene expression measurements can be used to predict gene function. However, the ability of basic learning methods to predict stress response varies widely and depends heavily on how much dimensionality reduction is used. Our method of combining classifiers can improve the accuracy of such predictions – in this case, predictions of genes involved in stress response in plants – and it effectively chooses the appropriate amount of dimensionality reduction automatically. The method provides a useful means of identifying genes in A. thaliana that potentially respond to stress, and we expect it would be useful in other organisms and for other gene functions.

Urquhart W, Gunawardena AH, Moeder W, Ali R, Berkowitz GA, Yoshioka K

Plant Mol. Biol. 2007 Dec;65(6):747-61

PubMed PMID: 17885810

Abstract

The hypersensitive response (HR) involves programmed cell death (PCD) in response to pathogen infection. To investigate the pathogen resistance signaling pathway, we previously identified the Arabidopsis mutant cpr22, which displays constitutive activation of multiple defense responses including HR like cell death. The cpr22 mutation has been identified as a 3 kb deletion that fuses two cyclic nucleotide-gated ion channel (CNGC)-encoding genes, ATCNGC11 and ATCNGC12, to generate a novel chimeric gene, ATCNGC11/12. In this study, we conducted a characterization of cell death induced by transient expression of ATCNGC11/12 in Nicotiana benthamiana. Electron microscopic analysis of this cell death showed similar characteristics to PCD, such as plasma membrane shrinkage and vesicle formation. The hallmark of animal PCD, fragmentation of nuclear DNA, was also observed in ATCNGC11/12-induced cell death. The development of cell death was significantly suppressed by caspase-1 inhibitors, suggesting the involvement of caspases in this process. Recently, vacuolar processing enzyme (VPE) was isolated as the first plant caspase-like protein, which is involved in HR development. In VPE-silenced plants development of cell death induced by ATCNGC11/12 was much slower and weaker compared to control plants, suggesting the involvement of VPE as a caspase in ATCNGC11/12-induced cell death. Complementation analysis using a Ca2+ uptake deficient yeast mutant demonstrated that the ATCNGC11/12 channel is permeable to Ca2+. Additionally, calcium channel blockers such as GdCl3 inhibited ATCNGC11/12-induced HR formation, whereas potassium channel blockers did not. Taken together, these results indicate that the cell death that develops in the cpr22 mutant is indeed PCD and that the chimeric channel, ATCNGC11/12, is at the point of, or up-stream of the calcium signal necessary for the development of HR.