Comparative Genomic Analysis of Bacterial Data in BV-BRC: An Example Exploring Antimicrobial Resistance.

As genomic and related data continue to expand, research biologists are often hampered by the computational hurdles required to analyze their data. The National Institute of Allergy and Infectious Diseases (NIAID) established the Bioinformatics Resource Centers (BRC) to assist researchers with their analysis of genome sequence and other omics-related data. Recently, the PAThosystems Resource Integration Center (PATRIC), the Influenza Research Database (IRD), and the Virus Pathogen Database and Analysis Resource (ViPR) BRCs merged to form the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) at https://www.bv-brc.org/ . The combined BV-BRC leverages the functionality of the original resources for bacterial and viral research communities with a unified data model, enhanced web-based visualization and analysis tools, and bioinformatics services. Here we demonstrate how antimicrobial resistance data can be analyzed in the new resource.

An automated system for interrogating the evolution of microbial endosymbiosis.

Inter-kingdom endosymbiotic interactions between bacteria and eukaryotic cells are critical to human health and disease. However, the molecular mechanisms that drive the emergence of endosymbiosis remain obscure. Here, we describe the development of a microfluidic system, named SEER (S̲ystem for the E̲volution of E̲ndosymbiotic R̲elationships), that automates the evolutionary selection of bacteria with enhanced intracellular survival and persistence within host cells, hallmarks of endosymbiosis. Using this system, we show that a laboratory strain of Escherichia coli that initially possessed limited abilities to survive within host cells, when subjected to SEER selection, rapidly evolved to display a 55-fold enhancement in intracellular survival. Notably, molecular dissection of the evolved strains revealed that a single-point mutation in a flexible loop of CpxR, a gene regulator that controls bacterial stress responses, substantially contributed to this intracellular survival. Taken together, these results establish SEER as the first microfluidic system for investigating the evolution of endosymbiosis, show the importance of CpxR in endosymbiosis, and set the stage for evolving bespoke inter-kingdom endosymbiotic systems with novel or emergent properties.

Identification of new protein-protein and protein-DNA interactions linked with wood formation in Populus trichocarpa.

Cellular processes, such as signal transduction and cell wall deposition, are organized by macromolecule interactions. Experimentally determined protein-protein interactions (PPIs) and protein-DNA interactions (PDIs) relevant to woody plant development are sparse. To begin to develop a Populus trichocarpa Torr. & A. Gray wood interactome, we applied the yeast-two-hybrid (Y2H) assay in different ways to enable the discovery of novel PPIs and connected networks. We first cloned open reading frames (ORFs) for 361 genes markedly upregulated in secondary xylem compared with secondary phloem and performed a binary Y2H screen with these proteins. By screening a xylem cDNA library for interactors of a subset of these proteins and then recapitulating the process by using a subset of the interactors as baits, we ultimately identified 165 PPIs involving 162 different ORFs. Thirty-eight transcription factors (TFs) included in our collection of P. trichocarpa wood ORFs were used in a Y1H screen for binding to promoter regions of three genes involved in lignin biosynthesis resulting in 40 PDIs involving 20 different TFs. The network incorporating both the PPIs and PDIs included 14 connected subnetworks, with the largest having 132 members. Protein-protein interactions and PDIs validated previous reports and also identified new candidate wood formation proteins and modules through their interactions with proteins and promoters known to be involved in secondary cell wall synthesis. Selected examples are discussed including a PPI between Mps one binder (MOB1) and a mitogen-activated protein kinase kinase kinase kinase (M4K) that was further characterized by assays confirming the PPI as well as its effect on subcellular localization. Mapping of published transcriptomic data showing developmentally detailed expression patterns across a secondary stem onto the network supported that the PPIs and PDIs are relevant to wood formation, and also illustrated that wood-associated interactions involve gene products that are not upregulated in secondary xylem.

Characterization and RNA-seq analysis of underperformer, an activation-tagged potato mutant.

The potato cv. Bintje and a Bintje activation-tagged mutant, underperformer (up) were compared. Mutant up plants grown in vitro were dwarf, with abundant axillary shoot growth, greater tuber yield, altered tuber traits and early senescence compared to wild type. Under in vivo conditions, the dwarf and early senescence phenotypes of the mutant remained, but the up plants exhibited a lower tuber yield and fewer axillary shoots compared to wild type. Southern blot analyses indicated a single T-DNA insertion in the mutant, located on chromosome 10. Initial PCR-based gene expression studies indicated transcriptional activation/repression of several genes in the mutant flanking the insertion. The gene immediately flanking the right border of the T-DNA insertion, which encoded an uncharacterized Broad complex, Tramtrac, Bric-a-brac; also known as Pox virus and Zinc finger (BTB/POZ) domain-containing protein (StBTB/POZ1) containing an Armadillo repeat region, was up-regulated in the mutant. Global gene expression comparisons between Bintje and up using RNA-seq on leaves from 60 day-old plants revealed a dataset of over 1,600 differentially expressed genes. Gene expression analyses suggested a variety of biological processes and pathways were modified in the mutant, including carbohydrate and lipid metabolism, cell division and cell cycle activity, biotic and abiotic stress responses, and proteolysis.

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence.

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States.

UNLABELLED: Coronaviruses are known to infect humans and other animals and cause respiratory and gastrointestinal diseases. Here we report the emergence of porcine epidemic diarrhea virus (PEDV) in the United States and determination of its origin, evolution, and genotypes based on temporal and geographical evidence. Histological lesions in small intestine sections of affected pigs and the complete genomic sequences of three emergent strains of PEDV isolated from outbreaks in Minnesota and Iowa were characterized. Genetic and phylogenetic analyses of the three U.S. strains revealed a close relationship with Chinese PEDV strains and their likely Chinese origin. The U.S. PEDV strains underwent evolutionary divergence, which can be classified into two sublineages. The three emergent U.S. strains are most closely related to a strain isolated in 2012 from Anhui Province in China, which might be the result of multiple recombination events between different genetic lineages or sublineages of PEDV. Molecular clock analysis of the divergent time based on the complete genomic sequences is consistent with the actual time difference, approximately 2 to 3 years, of the PED outbreaks between China (December 2010) and the United States (May 2013). The finding that the emergent U.S. PEDV strains share unique genetic features at the 5'-untranslated region with a bat coronavirus provided further support of the evolutionary origin of PEDV from bats and potential cross-species transmission. The data from this study have important implications for understanding the ongoing PEDV outbreaks in the United States and will guide future efforts to develop effective preventive and control measures against PEDV. IMPORTANCE: The sudden emergence of porcine epidemic diarrhea virus (PEDV), a coronavirus, for the first time in the United States causes significant economic and public health concerns. Since its recognition in May 2013, PEDV has rapidly spread across the United States, resulting in high mortality in piglets in more than 17 States now. The ongoing outbreaks of Middle East respiratory syndrome coronavirus in humans from countries in or near the Arabian Peninsula and the historical deadly nature of the 2002 outbreaks of severe acute respiratory syndrome coronavirus create further anxiety over the emergence of PEDV in the United States due to the lack of scientific information about the origin and evolution of this emerging coronavirus. Here we report the detailed genetic characterization, origin, and evolution of emergent PEDV strains in the United States. The results provide much needed information to devise effective preventive and control strategies against PEDV in the United States.

An active ac/ds transposon system for activation tagging in tomato cultivar m82 using clonal propagation.

Tomato (Solanum lycopersicum) is a model organism for Solanaceae in both molecular and agronomic research. This project utilized Agrobacterium tumefaciens transformation and the transposon-tagging construct Activator (Ac)/Dissociator (Ds)-ATag-Bar_gosGFP to produce activation-tagged and knockout mutants in the processing tomato cultivar M82. The construct carried hygromycin resistance (hyg), green fluorescent protein (GFP), and the transposase (TPase) of maize (Zea mays) Activator major transcript X054214.1 on the stable Ac element, along with a 35S enhancer tetramer and glufosinate herbicide resistance (BAR) on the mobile Ds-ATag element. An in vitro propagation strategy was used to produce a population of 25 T0 plants from a single transformed plant regenerated in tissue culture. A T1 population of 11,000 selfed and cv M82 backcrossed progeny was produced from the functional T0 line. This population was screened using glufosinate herbicide, hygromycin leaf painting, and multiplex polymerase chain reaction (PCR). Insertion sites of transposed Ds-ATag elements were identified through thermal asymmetric interlaced PCR, and resulting product sequences were aligned to the recently published tomato genome. A population of 509 independent, Ds-only transposant lines spanning all 12 tomato chromosomes has been developed. Insertion site analysis demonstrated that more than 80% of these lines harbored Ds insertions conducive to activation tagging. The capacity of the Ds-ATag element to alter transcription was verified by quantitative real-time reverse transcription-PCR in two mutant lines. The transposon-tagged lines have been immortalized in seed stocks and can be accessed through an online database, providing a unique resource for tomato breeding and analysis of gene function in the background of a commercial tomato cultivar.

Comparative genomics of early-diverging Brucella strains reveals a novel lipopolysaccharide biosynthesis pathway.

UNLABELLED: Brucella species are Gram-negative bacteria that infect mammals. Recently, two unusual strains (Brucella inopinata BO1T and B. inopinata-like BO2) have been isolated from human patients, and their similarity to some atypical brucellae isolated from Australian native rodent species was noted. Here we present a phylogenomic analysis of the draft genome sequences of BO1T and BO2 and of the Australian rodent strains 83-13 and NF2653 that shows that they form two groups well separated from the other sequenced Brucella spp. Several important differences were noted. Both BO1T and BO2 did not agglutinate significantly when live or inactivated cells were exposed to monospecific A and M antisera against O-side chain sugars composed of N-formyl-perosamine. While BO1T maintained the genes required to synthesize a typical Brucella O-antigen, BO2 lacked many of these genes but still produced a smooth LPS (lipopolysaccharide). Most missing genes were found in the wbk region involved in O-antigen synthesis in classic smooth Brucella spp. In their place, BO2 carries four genes that other bacteria use for making a rhamnose-based O-antigen. Electrophoretic, immunoblot, and chemical analyses showed that BO2 carries an antigenically different O-antigen made of repeating hexose-rich oligosaccharide units that made the LPS water-soluble, which contrasts with the homopolymeric O-antigen of other smooth brucellae that have a phenol-soluble LPS. The results demonstrate the existence of a group of early-diverging brucellae with traits that depart significantly from those of the Brucella species described thus far. IMPORTANCE: This report examines differences between genomes from four new Brucella strains and those from the classic Brucella spp. Our results show that the four new strains are outliers with respect to the previously known Brucella strains and yet are part of the genus, forming two new clades. The analysis revealed important information about the evolution and survival mechanisms of Brucella species, helping reshape our knowledge of this important zoonotic pathogen. One discovery of special importance is that one of the strains, BO2, produces an O-antigen distinct from any that has been seen in any other Brucella isolates to date.

Transposon tagging in diploid strawberry.

Fragaria vesca was transformed with a transposon tagging construct harbouring amino terminally deleted maize transposase and EGFP (Ac element), NPTII, CaMV 35S promoter (P35S) driving transposase and mannopine synthase promoter (Pmas) driving EGFP (Ds element). Of 180 primary transgenics, 48 were potential launch pads, 72 were multiple insertions or chimaeras, and 60 exhibited somatic transposition. T1 progeny of 32 putative launch pads were screened by multiplex PCR for transposition. Evidence of germ-line transposition occurred in 13 putative launch pads; however, the transposition frequency was too low in three for efficient recovery of transposants. The transposition frequency in the remaining launch pads ranged from 16% to 40%. After self-pollination of the T0 launch pads, putative transposants in the T1 generation were identified by multiplex PCR. Sequencing of hiTAIL-PCR products derived from nested primers within the Ds end sequences (either P35S at the left border or the inverted repeat at the right border) of T1 plants revealed transposition of the Ds element to distant sites in the strawberry genome. From more than 2400 T1 plants screened, 103 unique transposants have been identified, among which 17 were somatic transpositions observed in the T0 generation. Ds insertion sites were dispersed among various gene elements [exons (15%), introns (23%), promoters (30%), 3' UTRs (17%) as well as intergenically (15%)]. Three-primer (one on either side of the Ds insertion and one within the Ds T-DNA) PCR could be used to identify homozygous T2 transposon-tagged plants. The mutant collection has been catalogued in an on-line database.

Comparative genomics of early-diverging Brucella strains reveals a novel lipopolysaccharide biosynthesis pathway.

UNLABELLED: Brucella species are Gram-negative bacteria that infect mammals. Recently, two unusual strains (Brucella inopinata BO1(T) and B. inopinata-like BO2) have been isolated from human patients, and their similarity to some atypical brucellae isolated from Australian native rodent species was noted. Here we present a phylogenomic analysis of the draft genome sequences of BO1(T) and BO2 and of the Australian rodent strains 83-13 and NF2653 that shows that they form two groups well separated from the other sequenced Brucella spp. Several important differences were noted. Both BO1(T) and BO2 did not agglutinate significantly when live or inactivated cells were exposed to monospecific A and M antisera against O-side chain sugars composed of N-formyl-perosamine. While BO1(T) maintained the genes required to synthesize a typical Brucella O-antigen, BO2 lacked many of these genes but still produced a smooth LPS (lipopolysaccharide). Most missing genes were found in the wbk region involved in O-antigen synthesis in classic smooth Brucella spp. In their place, BO2 carries four genes that other bacteria use for making a rhamnose-based O-antigen. Electrophoretic, immunoblot, and chemical analyses showed that BO2 carries an antigenically different O-antigen made of repeating hexose-rich oligosaccharide units that made the LPS water-soluble, which contrasts with the homopolymeric O-antigen of other smooth brucellae that have a phenol-soluble LPS. The results demonstrate the existence of a group of early-diverging brucellae with traits that depart significantly from those of the Brucella species described thus far. IMPORTANCE: This report examines differences between genomes from four new Brucella strains and those from the classic Brucella spp. Our results show that the four new strains are outliers with respect to the previously known Brucella strains and yet are part of the genus, forming two new clades. The analysis revealed important information about the evolution and survival mechanisms of Brucella species, helping reshape our knowledge of this important zoonotic pathogen. One discovery of special importance is that one of the strains, BO2, produces an O-antigen distinct from any that has been seen in any other Brucella isolates to date.

Phylogeny of gammaproteobacteria.

The phylogeny of the large bacterial class Gammaproteobacteria has been difficult to resolve. Here we apply a telescoping multiprotein approach to the problem for 104 diverse gammaproteobacterial genomes, based on a set of 356 protein families for the whole class and even larger sets for each of four cohesive subregions of the tree. Although the deepest divergences were resistant to full resolution, some surprising patterns were strongly supported. A representative of the Acidithiobacillales routinely appeared among the outgroup members, suggesting that in conflict with rRNA-based phylogenies this order does not belong to Gammaproteobacteria; instead, it (and, independently, "Mariprofundus") diverged after the establishment of the Alphaproteobacteria yet before the betaproteobacteria/gammaproteobacteria split. None of the orders Alteromonadales, Pseudomonadales, or Oceanospirillales were monophyletic; we obtained strong support for clades that contain some but exclude other members of all three orders. Extreme amino acid bias in the highly A+T-rich genome of Candidatus Carsonella prevented its reliable placement within Gammaproteobacteria, and high bias caused artifacts that limited the resolution of the relationships of other insect endosymbionts, which appear to have had multiple origins, although the unbiased genome of the endosymbiont Sodalis acted as an attractor for them. Instability was observed for the root of the Enterobacteriales, with nearly equal subsets of the protein families favoring one or the other of two alternative root positions; the nematode symbiont Photorhabdus was identified as a disruptor whose omission helped stabilize the Enterobacteriales root.

Genome sequences of three agrobacterium biovars help elucidate the evolution of multichromosome genomes in bacteria.

The family Rhizobiaceae contains plant-associated bacteria with critical roles in ecology and agriculture. Within this family, many Rhizobium and Sinorhizobium strains are nitrogen-fixing plant mutualists, while many strains designated as Agrobacterium are plant pathogens. These contrasting lifestyles are primarily dependent on the transmissible plasmids each strain harbors. Members of the Rhizobiaceae also have diverse genome architectures that include single chromosomes, multiple chromosomes, and plasmids of various sizes. Agrobacterium strains have been divided into three biovars, based on physiological and biochemical properties. The genome of a biovar I strain, A. tumefaciens C58, has been previously sequenced. In this study, the genomes of the biovar II strain A. radiobacter K84, a commercially available biological control strain that inhibits certain pathogenic agrobacteria, and the biovar III strain A. vitis S4, a narrow-host-range strain that infects grapes and invokes a hypersensitive response on nonhost plants, were fully sequenced and annotated. Comparison with other sequenced members of the Alphaproteobacteria provides new data on the evolution of multipartite bacterial genomes. Primary chromosomes show extensive conservation of both gene content and order. In contrast, secondary chromosomes share smaller percentages of genes, and conserved gene order is restricted to short blocks. We propose that secondary chromosomes originated from an ancestral plasmid to which genes have been transferred from a progenitor primary chromosome. Similar patterns are observed in select Beta- and Gammaproteobacteria species. Together, these results define the evolution of chromosome architecture and gene content among the Rhizobiaceae and support a generalized mechanism for second-chromosome formation among bacteria.

A robust species tree for the alphaproteobacteria.

The branching order and coherence of the alphaproteobacterial orders have not been well established, and not all studies have agreed that mitochondria arose from within the Rickettsiales. A species tree for 72 alphaproteobacteria was produced from a concatenation of alignments for 104 well-behaved protein families. Coherence was upheld for four of the five orders with current standing that were represented here by more than one species. However, the family Hyphomonadaceae was split from the other Rhodobacterales, forming an expanded group with Caulobacterales that also included Parvularcula. The three earliest-branching alphaproteobacterial orders were the Rickettsiales, followed by the Rhodospirillales and then the Sphingomonadales. The principal uncertainty is whether the expanded Caulobacterales group is more closely associated with the Rhodobacterales or the Rhizobiales. The mitochondrial branch was placed within the Rickettsiales as a sister to the combined Anaplasmataceae and Rickettsiaceae, all subtended by the Pelagibacter branch. Pelagibacter genes will serve as useful additions to the bacterial outgroup in future evolutionary studies of mitochondrial genes, including those that have transferred to the eukaryotic nucleus.

GeneTrees: a phylogenomics resource for prokaryotes.

The GeneTrees phylogenomics system pursues comparative genomic analyses from the perspective of gene phylogenies for individual genes. The GeneTrees project has the goal of providing detailed evolutionary models for all protein-coding gene components of the fully sequenced genomes. Currently, a database of alignments and trees for all protein sequences for 325 fully sequenced and annotated prokaryote genomes is available. The prokaryote database contains 890,000 protein sequences organized into over 100,000 alignments, each described by a phylogenetic tree. An original homology group discovery tool assembles sets of related proteins from all versus all pairwise alignments. Multiple alignments for each homology group are stored and subjected to phylogenetic tree inference. A graphical web interface provides visual exploration of the GeneTrees database. Homology groups can be queried by sequence identifiers or annotation terms. Genomes can be browsed visually on a gene map of each chromosome or plasmid. Phylogenetic trees with support values are displayed in conjunction with the associated sequence alignment. A variety of classes of information can be selected to label the tree tips to aid in visual evaluation of annotation and gene function. This web interface is available at http://genetrees.vbi.vt.edu.

Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis.

Draft genome sequences have been determined for the soybean pathogen Phytophthora sojae and the sudden oak death pathogen Phytophthora ramorum. Oömycetes such as these Phytophthora species share the kingdom Stramenopila with photosynthetic algae such as diatoms, and the presence of many Phytophthora genes of probable phototroph origin supports a photosynthetic ancestry for the stramenopiles. Comparison of the two species' genomes reveals a rapid expansion and diversification of many protein families associated with plant infection such as hydrolases, ABC transporters, protein toxins, proteinase inhibitors, and, in particular, a superfamily of 700 proteins with similarity to known oömycete avirulence genes.

The xylem and phloem transcriptomes from secondary tissues of the Arabidopsis root-hypocotyl.

The growth of secondary xylem and phloem depends on the division of cells in the vascular cambium and results in an increase in the diameter of the root and stem. Very little is known about the genetic mechanisms that control cambial activity and the differentiation of secondary xylem and phloem cell types. To begin to identify new genes required for vascular cell differentiation and function, we performed genome-wide expression profiling of xylem and phloem-cambium isolated from the root-hypocotyl of Arabidopsis (Arabidopsis thaliana). Gene expression in the remaining nonvascular tissue was also profiled. From these transcript profiles, we assembled three sets of genes with expression significantly biased toward xylem, phloem-cambium, or nonvascular tissue. We also assembled three two-tissue sets of genes with expression significantly biased toward xylem/phloem-cambium, xylem/nonvascular, or phloem-cambium/nonvascular tissues. Localizations predicted by transcript profiles were supported by results from promoter-reporter and reverse transcription-polymerase chain reaction experiments with nine xylem- or phloem-cambium-biased genes. An analysis of the members of the phloem-cambium gene set suggested that some genes involved in regulating primary meristems are also regulators of the cambium. Secondary phloem was implicated in the synthesis of auxin, glucosinolates, cytokinin, and gibberellic acid. Transcript profiles also supported the importance of class III HD ZIP and KANADI transcription factors as regulators of radial patterning during secondary growth, and identified several members of the G2-like, NAC, AP2, MADS, and MYB transcription factor families that may play roles as regulators of xylem or phloem cell differentiation and activity.

PIML: the Pathogen Information Markup Language.

MOTIVATION: A vast amount of information about human, animal and plant pathogens has been acquired, stored and displayed in varied formats through different resources, both electronically and otherwise. However, there is no community standard format for organizing this information or agreement on machine-readable format(s) for data exchange, thereby hampering interoperation efforts across information systems harboring such infectious disease data. RESULTS: The Pathogen Information Markup Language (PIML) is a free, open, XML-based format for representing pathogen information. XSLT-based visual presentations of valid PIML documents were developed and can be accessed through the PathInfo website or as part of the interoperable web services federation known as ToolBus/PathPort. Currently, detailed PIML documents are available for 21 pathogens deemed of high priority with regard to public health and national biological defense. A dynamic query system allows simple queries as well as comparisons among these pathogens. Continuing efforts are being taken to include other groups' supporting PIML and to develop more PIML documents. AVAILABILITY: All the PIML-related information is accessible from http://www.vbi.vt.edu/pathport/pathinfo/

The S8 serine, C1A cysteine and A1 aspartic protease families in Arabidopsis.

The Arabidopsis thaliana genome has over 550 protease sequences representing all five catalytic types: serine, cysteine, aspartic acid, metallo and threonine (MEROPS peptidase database, http://merops.sanger.ac.uk/), which probably reflect a wide variety of as yet unidentified functions performed by plant proteases. Recent indications that the 26S proteasome, a T1 family-threonine protease, is a regulator of light and hormone responsive signal transduction highlight the potential of proteases to participate in many aspects of plant growth and development. Recent discoveries that proteases are required for stomatal distribution, embryo development and disease resistance point to wider roles for four additional multigene families that include some of the most frequently studied (yet poorly understood) plant proteases: the subtilisin-like, serine proteases (family S8), the papain-like, cysteine proteases (family C1A), the pepsin-like, aspartic proteases (family A1) and the plant matrixin, metalloproteases (family M10A). In this report, 54 subtilisin-like, 30 papain-like and 59 pepsin-like proteases from Arabidopsis, are compared with S8, C1A and A1 proteases known from other plant species at the functional, phylogenetic and gene structure levels. Examples of structural conservation between S8, C1A and A1 genes from rice, barley, tomato and soybean and those from Arabidopsis are noted, indicating that some common, essential plant protease roles were established before the divergence of monocots and eudicots. Numerous examples of tandem duplications of protease genes and evidence for a variety of restricted expression patterns suggest that a high degree of specialization exists among proteases within each family. We propose that comprehensive analysis of the functions of these genes in Arabidopsis will firmly establish serine, cysteine and aspartic proteases as regulators and effectors of a wide range of plant processes.