TreeDomViewer: a tool for the visualization of phylogeny and protein domain structure.

Phylogenetic analysis and examination of protein domains allow accurate genome annotation and are invaluable to study proteins and protein complex evolution. However, two sequences can be homologous without sharing statistically significant amino acid or nucleotide identity, presenting a challenging bioinformatics problem. We present TreeDomViewer, a visualization tool available as a web-based interface that combines phylogenetic tree description, multiple sequence alignment and InterProScan data of sequences and generates a phylogenetic tree projecting the corresponding protein domain information onto the multiple sequence alignment. Thereby it makes use of existing domain prediction tools such as InterProScan. TreeDomViewer adopts an evolutionary perspective on how domain structure of two or more sequences can be aligned and compared, to subsequently infer the function of an unknown homolog. This provides insight into the function assignment of, in terms of amino acid substitution, very divergent but yet closely related family members. Our tool produces an interactive scalar vector graphics image that provides orthological relationship and domain content of proteins of interest at one glance. In addition, PDF, JPEG or PNG formatted output is also provided. These features make TreeDomViewer a valuable addition to the annotation pipeline of unknown genes or gene products. TreeDomViewer is available at http://www.bioinformatics.nl/tools/treedom/.

Tracing the evolution of a large protein complex in the eukaryotes, NADH:ubiquinone oxidoreductase (Complex I).

The increasing availability of sequenced genomes enables the reconstruction of the evolutionary history of large protein complexes. Here, we trace the evolution of NADH:ubiquinone oxidoreductase (Complex I), which has increased in size, by so-called supernumary subunits, from 14 subunits in the bacteria to 30 in the plants and algae, 37 in the fungi and 46 in the mammals. Using a combination of pair-wise and profile-based sequence comparisons at the levels of proteins and the DNA of the sequenced eukaryotic genomes, combined with phylogenetic analyses to establish orthology relationships, we were able to (1) trace the origin of six of the supernumerary subunits to the alpha-proteobacterial ancestor of the mitochondria, (2) detect previously unidentified homology relations between subunits from fungi and mammals, (3) detect previously unidentified subunits in the genomes of several species and (4) document several cases of gene duplications among supernumerary subunits in the eukaryotes. One of these, a duplication of N7BM (B17.2), is particularly interesting as it has been lost from genomes that have also lost Complex I proteins, making it a candidate for a Complex I interacting protein. A parsimonious reconstruction of eukaryotic Complex I evolution shows an initial increase in size that predates the separation of plants, fungi and metazoa, followed by a gradual adding and incidental losses of subunits in the various evolutionary lineages. This evolutionary scenario is in contrast to that for Complex I in the prokaryotes, for which the combination of several separate, and previously independently functioning modules into a single complex has been proposed.

Rickettsia phylogenomics: unwinding the intricacies of obligate intracellular life.

BACKGROUND: Completed genome sequences are rapidly increasing for Rickettsia, obligate intracellular alpha-proteobacteria responsible for various human diseases, including epidemic typhus and Rocky Mountain spotted fever. In light of phylogeny, the establishment of orthologous groups (OGs) of open reading frames (ORFs) will distinguish the core rickettsial genes and other group specific genes (class 1 OGs or C1OGs) from those distributed indiscriminately throughout the rickettsial tree (class 2 OG or C2OGs). METHODOLOGY/PRINCIPAL FINDINGS: We present 1823 representative (no gene duplications) and 259 non-representative (at least one gene duplication) rickettsial OGs. While the highly reductive (approximately 1.2 MB) Rickettsia genomes range in predicted ORFs from 872 to 1512, a core of 752 OGs was identified, depicting the essential Rickettsia genes. Unsurprisingly, this core lacks many metabolic genes, reflecting the dependence on host resources for growth and survival. Additionally, we bolster our recent reclassification of Rickettsia by identifying OGs that define the AG (ancestral group), TG (typhus group), TRG (transitional group), and SFG (spotted fever group) rickettsiae. OGs for insect-associated species, tick-associated species and species that harbor plasmids were also predicted. Through superimposition of all OGs over robust phylogeny estimation, we discern between C1OGs and C2OGs, the latter depicting genes either decaying from the conserved C1OGs or acquired laterally. Finally, scrutiny of non-representative OGs revealed high levels of split genes versus gene duplications, with both phenomena confounding gene orthology assignment. Interestingly, non-representative OGs, as well as OGs comprised of several gene families typically involved in microbial pathogenicity and/or the acquisition of virulence factors, fall predominantly within C2OG distributions. CONCLUSION/SIGNIFICANCE: Collectively, we determined the relative conservation and distribution of 14354 predicted ORFs from 10 rickettsial genomes across robust phylogeny estimation. The data, available at PATRIC (PathoSystems Resource Integration Center), provide novel information for unwinding the intricacies associated with Rickettsia pathogenesis, expanding the range of potential diagnostic, vaccine and therapeutic targets.

An NB-LRR protein required for HR signalling mediated by both extra- and intracellular resistance proteins.

Tomato (Solanum lycopersicum) Cf resistance genes confer hypersensitive response (HR)-associated resistance to strains of the pathogenic fungus Cladosporium fulvum that express the matching avirulence (Avr) gene. Previously, we identified an Avr4-responsive tomato (ART) gene that is required for Cf-4/Avr4-induced HR in Nicotiana benthamiana as demonstrated by virus-induced gene silencing (VIGS). The gene encodes a CC-NB-LRR type resistance (R) protein analogue that we have designated NRC1 (NB-LRR protein required for HR-associated cell death 1). Here we describe that knock-down of NRC1 in tomato not only affects the Cf-4/Avr4-induced HR but also compromises Cf-4-mediated resistance to C. fulvum. In addition, VIGS using NRC1 in N. benthamiana revealed that this protein is also required for the HR induced by the R proteins Cf-9, LeEix, Pto, Rx and Mi. Transient expression of NRC1(D481V), which encodes a constitutively active NRC1 mutant protein, triggers an elicitor-independent HR. Subsequently, we transiently expressed this auto-activating protein in N. benthamiana silenced for genes known to be involved in HR signalling, thereby allowing NRC1 to be positioned in an HR signalling pathway. We found that NRC1 requires RAR1 and SGT1 to be functional, whereas it does not require NDR1 and EDS1. As the Cf-4 protein requires EDS1 for its function, we hypothesize that NRC1 functions downstream of EDS1. We also found that NRC1 acts upstream of a MAP kinase pathway. We conclude that Cf-mediated resistance signalling requires a downstream NB-LRR protein that also functions in cell death signalling pathways triggered by other R proteins.

Changes in gene expression during male meiosis in Petunia hybrida.

We analyzed changes in gene expression during male meiosis in Petunia by combining the meiotic staging of pollen mother cells from a single anther with cDNA-AFLP transcript profiling of mRNA from the synchronously developing sister anthers. The transcript profiling experiments focused on the identification of genes with a modulated expression profile during meiosis, while premeiotic archesporial cells and postmeiotic microspores served as a reference. About 8000 transcript tags, estimated at 30% of the total transcriptome, were generated, of which around 6% exhibited a modulated gene expression pattern at meiosis. Cluster analysis revealed a transcriptional cascade that coincides with the initiation and progression through all stages of the two meiotic divisions. Fragments that exhibited high expression specifically during meiosis I were characterized further by sequencing; 90 out of the 293 sequenced fragments showed homology with known genes, belonging to a wide range of gene classes, including previously characterized meiotic genes. In-situ hybridization experiments were performed to determine the spatial expression pattern for five selected transcript tags. Its concurrence with cDNA-AFLP transcript profiles indicates that this is an excellent approach to study genes involved in specialized processes such as meiosis. Our data set provides the potential to unravel unique meiotic genes that are as yet elusive to reverse genetics approaches.