Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Comprehensive analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88Linked to original sources
PubMed2006

The RNA polymerase III-dependent family of genes in hemiascomycetes: comparative RNomics, decoding strategies, transcription and evolutionary implications.

We present the first comprehensive analysis of RNA polymerase III (Pol III) transcribed genes in ten yeast genomes. This set includes all tRNA genes (tDNA) and genes coding for SNR6 (U6), SNR52, SCR1 and RPR1 RNA in the nine hemiascomycetes Saccharomyces cerevisiae, Saccharomyces castellii, Candida glabrata, Kluyveromyces waltii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and the archiascomycete Schizosaccharomyces pombe. We systematically analysed sequence specificities of tRNA genes, polymorphism, variability of introns, gene redundancy and gene clustering. Analysis of decoding strategies showed that yeasts close to S.cerevisiae use bacterial decoding rules to read the Leu CUN and Arg CGN codons, in contrast to all other known Eukaryotes. In D.hansenii and C.albicans, we identified a novel tDNA-Leu (AAG), reading the Leu CUU/CUC/CUA codons with an unusual G at position 32. A systematic 'p-distance tree' using the 60 variable positions of the tRNA molecule revealed that most tDNAs cluster into amino acid-specific sub-trees, suggesting that, within hemiascomycetes, orthologous tDNAs are more closely related than paralogs. We finally determined the bipartite A- and B-box sequences recognized by TFIIIC. These minimal sequences are nearly conserved throughout hemiascomycetes and were satisfactorily retrieved at appropriate locations in other Pol III genes.

Ascomycota↗
PubMed2006

PONGO: a web server for multiple predictions of all-alpha transmembrane proteins.

The annotation efforts of the BIOSAPIENS European Network of Excellence have generated several distributed annotation systems (DAS) with the aim of integrating Bioinformatics resources and annotating metazoan genomes (http://www.biosapiens.info). In this context, the PONGO DAS server (http://pongo.biocomp.unibo.it) provides the annotation on predictive basis for the all-alpha membrane proteins in the human genome, not only through DAS queries, but also directly using a simple web interface. In order to produce a more comprehensive analysis of the sequence at hand, this annotation is carried out with four selected and high scoring predictors: TMHMM2.0, MEMSAT, PRODIV and ENSEMBLE1.0. The stored and pre-computed predictions for the human proteins can be searched and displayed in a graphical view. However the web service allows the prediction of the topology of any kind of putative membrane proteins, regardless of the organism and more importantly with the same sequence profile for a given sequence when required. Here we present a new web server that incorporates the state-of-the-art topology predictors in a single framework, so that putative users can interactively compare and evaluate four predictions simultaneously for a given sequence. Together with the predicted topology, the server also displays a signal peptide prediction determined with SPEP. The PONGO web server is available at http://pongo.biocomp.unibo.it/pongo.

Humans↗
PubMed2006

CRSD: a comprehensive web server for composite regulatory signature discovery.

Transcription factors (TFs) and microRNAs play important roles in the regulation of human gene expression, and the study of their combinatory regulations of gene expression is a new research field. We constructed a comprehensive web server, the composite regulatory signature database (CRSD), that can be applied in investigating complex regulatory behaviors involving gene expression signatures (GESs), microRNA regulatory signatures (MRSs) and TF regulatory signatures (TRSs). Six well-known and large-scale databases, including the human UniGene, mature microRNAs, putative promoter, TRANSFAC, pathway and Gene Ontology (GO) databases, were integrated to provide the comprehensive analysis in CRSD. Two new genome-wide databases, of MRSs and TRSs, were also constructed and further integrated into CRSD. To accomplish the microarray data analysis at one go, several methods, including microarray data pretreatment, statistical and clustering analysis, iterative enrichment analysis and motif discovery, were closely integrated in the web server, which has not been the case in previous studies. Our implementation showed that the published literature could demonstrate the results of genome-wide enrichment analysis. We conclude that CRSD is a powerful and useful bioinformatic web server and may provide new insights into gene regulation networks. CRSD and the online tutorial are publicly available at http://biochip.nchu.edu.tw/crsd1/.

3' Untranslated Regions↗
PubMed2003

Specificity assessment from fractionation experiments (SAFE): a novel method to evaluate microarray probe specificity based on hybridisation stringencies.

The cDNA-chip technology is a highly versatile tool for the comprehensive analysis of gene expression at the transcript level. Although it has been applied successfully in expression profiling projects, there is an ongoing dispute concerning the quality of such expression data. The latter critically depends on the specificity of hybridisation. SAFE (specificity assessment from fractionation experiments) is a novel method to discriminate between non- specific cross-hybridisation and specific signals. We applied in situ fractionation of hybridised target on DNA-chips by means of repeated washes with increasing stringencies. Different fractions of hybridised target are washed off at defined stringencies and the collected fluorescence intensity data at each step comprise the fractionation curve. Based on characteristic features of the fractionation curve, unreliable data can be filtered and eliminated from subsequent analyses. The approach described here provides a novel experimental tool to identify probes that produce specific hybridisation signals in DNA-chip expression profiling approaches. The iterative use of the SAFE procedure will result in increasingly reliable sets of probes for microarray experiments and significantly improve the overall efficiency and reliability of RNA expression profiling data from DNA-chip experiments.

Animals↗
PubMed2006

Affinity selection of DNA-binding protein complexes using mRNA display.

Comprehensive analysis of DNA-protein interactions is important for mapping transcriptional regulatory networks on a genome-wide level. Here we present a new application of mRNA display for in vitro selection of DNA-binding protein heterodimeric complexes. Under improved selection conditions using a TPA-responsive element (TRE) as a bait DNA, known interactors c-fos and c-jun were simultaneously enriched about 100-fold from a model library (a 1:1:20 000 mixture of c-fos, c-jun and gst genes) after one round of selection. Furthermore, almost all kinds of the AP-1 family genes including c-jun, c-fos, junD, junB, atf2 and b-atf were successfully selected from an mRNA display library constructed from a mouse brain poly A(+) RNA after six rounds of selection. These results indicate that the mRNA display selection system can identify a variety of DNA-binding protein complexes in a single experiment. Since almost all transcription factors form heterooligomeric complexes to bind with their target DNA, this method should be most useful to search for DNA-binding transcription factor complexes.

Animals↗
PubMed2026

LAMBDA: a prophage detection benchmark for genomic language models.

Transformer-based genomic sequence models represent an emerging frontier in computational biology. Yet, their embeddings have not yet shown the same level of predictive power as natural and protein language models, highlighting a gap between current implementations and theoretical promise. Existing benchmarks for DNA language models primarily focus on classifying regulatory elements in eukaryotic genomes, leaving open the fundamental question of whether these models learn sequence-level features across whole genomes. We introduce LAMBDA, a benchmark designed to rigorously evaluate genome language model embeddings through phage-bacteria sequence discrimination across four categories of increasing complexity: probing tasks, fine-tuning assessments, diagnostic tests, and genome-wide prophage detection. Our comprehensive analysis of current genomic language models provides insight into the importance of training data selection relative to model size, the need for domain-specific training, and the capabilities and limitations of genomic language models for detecting prophage sequences. This benchmark represents a challenging genomic annotation task in the bacterial domain and addresses a key computational problem with direct relevance to microbiology and medicine.

Prophages↗
PubMed2002

Identification of AP-2-regulated genes by macroarray profiling of gene expression in human A375P melanoma.

Transcription factor AP-2 is a negative regulator of metastasis. Its expression is down regulated with progression of melanoma cells to metastasis. In this study, we performed macroarray profiling of gene expression of human A375P melanoma cells and their derivatives with overexpression of AP-2 and dominant-negative AP-2. Such comprehensive analysis lead to an identification of genes such as MMP-2, E-cadherin, melanoma adhesion molecule, early growth response 1, fibroblast growth factor 3, ubiquitin C, histone deacetylase 3 and integrin alpha 5,7, beta 3,5 as regulated by AP-2. Whereas some of these are known as AP-2-regulated genes, the others are not so far. Thus the study reports for the first time identification of new genes regulated by AP-2 that may be involved in metastasis of melanoma.

DNA-Binding Proteins↗
PubMed1997

Energy metabolism in exertional heat stroke with acute renal failure.

BACKGROUND: Heat stroke is the clinical syndrome produced when the body overheats. It can develop in the army and in healthy civilian populations who physically exert themselves in a hot and humid environment during the summer, and may result in a significant number of heat-related deaths. Since strenuous exercise is one of the major exacerbating and precipitating factors, the incidence of exertional heat stroke (ExHS) is high among military personnel undergoing military training. Furthermore, acute renal failure (ARF) may occur in 25% of patients with ExHS and it can cause metabolic alterations that affect amino acid, carbohydrate, and lipid metabolism. Adequate nutritional support is essential for the treatment of ARF. The most important determinant of nutrient requirement in ARF is the degree of hypercatabolism caused by disease associated with renal function impairment. Indirect calorimetry (IDCM) is the method by which metabolic rate is estimated from measurements of oxygen consumption and carbon dioxide production. It can also provide information about the type and rate of substrate utilization in vivo (protein, carbohydrate, and fat). METHOD: The present clinical study is a comprehensive analysis of metabolic changes which includes energy expenditure (EE) and substrate utilization in 10 patients with ExHS with ARF and 10 patients with exertional heat exhaustion (ExHE) by the use of IDCM. RESULTS: Serum urea nitrogen, creatinine, peak creatine phosphokinase levels and heart rate were significantly increased in ExHS patients during ARF stage. Serum albumin levels were significantly decreased in ExHS patients with ARF. Resting energy expenditure (REE) was increased in patients with ExHS induced ARF and was not correlated with body temperature (r = 0.421). The average increase in EE during ARF stage was about 24%. The respiratory quotient in patients with ExHS induced ARF was lower than that in normal subjects and also in patients with ExHE. Urea nitrogen appearance rate increased in patients with ExHS induced ARF and in patients with ExHE without ARF. The percentage of total REE derived from fat in ExHS induced ARF and ExHE increased, while in patients with ExHS induced ARF and ExHE, the percentages of total REE derived from carbohydrate and protein were lower than those in control subjects. CONCLUSIONS: The present results suggest that patients with exertional heat injury (both ExHS and ExHE) have hypermetabolism during the acute stage. Furthermore, patients with exertional heat-induced rhabdomyolysis and ARF have a moderately higher hypermetabolism than those without ARF during the acute stage. We believe that this mainly reflects a more pronounced reduction of the vital cell mass (muscle) in relation to body weight, and/or a compromised substrate oxidation in ExHS with ARF. Whether or not this subgroup of patients will require a higher energy/caloric support merits further investigation.

Acute Kidney Injury↗
PubMed2001

Evolution of alpha-amylases: architectural features and key residues in the stabilization of the (beta/alpha)(8) scaffold.

We provide a comprehensive analysis of the current enzymes with alpha-amylase activity (AAMYs) that belong to family 13 glycoside hydrolase (GH-13; 144 Archaea, Bacteria, and Eukaryota sequences from 87 different species). This study aims to further knowledge of the evolutionary molecular relationships among the sequences of their A and B domains with special emphasis on the correlation between what is observed in the structures and protein evolution. Multialignments for the A domain distinguish two clusters for sequences from Archaea organisms, eight for sequences from Bacteria organisms, and three for sequences from Eukaryota organisms. The clusters for Bacteria do not follow any strict taxonomic pathway; in fact, they are rather scattered. When we compared the A domains of sequences belonging to different kingdoms, we found that various pairs of clusters were significantly similar. Using either sequence similarity with crystallized structures or secondary-structure prediction methods, we identified in all AAMYs the eight putative beta-strands that constitute the beta-sheet in the TIM barrel of the A domain and studied the packing in its interior. We also discovered a "hidden homology" in the TIM barrel, an invariant Gly located upstream in the sequence before the conserved Asp in beta-strand 3. This Gly precedes an alpha-helix and is actively involved in capping its N-terminal end with a capping box. In all cases, a Schellman motif caps the C-terminal end of this helix.

Cluster Analysis↗
PubMed1998

Evolutionary relationships among the members of an ancient class of non-LTR retrotransposons found in the nematode Caenorhabditis elegans.

We took advantage of the massive amount of sequence information generated by the Caenorhabditis elegans genome project to perform a comprehensive analysis of a group of over 100 related sequences that has allowed us to describe two new C. elegans non-LTR retrotransposons. We named them Sam and Frodo. We also determined that several highly divergent subfamilies of both elements exist in C. elegans. It is likely that several master copies have been active at the same time in C. elegans, although only a few copies of both Sam and Frodo have characteristics that are compatible with them being active today. We discuss whether it is more appropriate under these circumstances to define only 2 elements corresponding to the most divergent groups of sequences or up to 16, considering each subfamily a different element. The C. elegans elements are related to other previously described non-LTR retrotransposons (CR1, found in different vertebrates; SR1, from the trematode Schistosoma; Q and T1, from the mosquito Anopheles). All of these elements, according to the analysis of their reverse transcriptases, form a monophyletic cluster that we call the "T1/CR1 subgroup." Elements of this subgroup are thus ancient components of the genome of animal species. However, we discuss the possibility that these elements may occasionally be horizontally transmitted.

Amino Acid Sequence↗
PubMed2000

Ty3/Gypsy retrotransposons: description of new Arabidopsis thaliana elements and evolutionary perspectives derived from comparative genomic data.

We performed a comprehensive analysis of the evolution of the Ty3/GYPSY: group of long-terminal-repeat retrotransposons (also known as METAVIRIDAE:). Exhaustive database searches allowed us to detect novel elements of this group. In particular, the Arabidopsis thaliana and Drosophila melanogaster genome sequencing projects have recently disclosed a large number of new Ty3/GYPSY: sequences. So far, elements of three different Ty3/GYPSY: lineages had been described for A. thaliana. Here, we describe six new lineages, which we have called Tit-for-tat1, Tit-for-tat2, Gimli, Gloin, Legolas, and Little Athila. We confirm that plant Ty3/GYPSY: elements form two main monophyletic groups. Moreover, our results suggest that at least four independent ancestral lineages existed before the monocot-dicot split, about 200 MYA. Twelve sequences from D. melanogaster that may correspond to new elements are also described. Some of these sequences are similar to those of OSVALDO: and Ulysses, two elements of the OSVALDO: clade that had never before been described for D. melanogaster. Comparative analyses of multiple organisms, some of them with completely sequenced genomes, show that the number of lineages of Ty3/GYPSY: elements is very variable. Thus, while only 1 lineage is present in Saccharomyces cerevisiae, at least 6 exist in Caenorhabditis elegans, at least 9 are present in the A. thaliana, and perhaps 20 are present in D. melanogaster. Finally, we suggest that the presence of a chromodomain-containing integrase, a feature of some closely related Ty3/GYPSY: elements of fungi, plants, and animals, may be used to define a new METAVIRIDAE: genus.

Amino Acid Sequence↗
PubMed2024

PSEUDO-RESPONSE REGULATOR 3b and transcription factor ABF3 modulate abscisic acid-dependent drought stress response in soybean.

The circadian system plays a pivotal role in facilitating the ability of crop plants to respond and adapt to fluctuations in their immediate environment effectively. Despite the increasing comprehension of PSEUDO-RESPONSE REGULATORs and their involvement in the regulation of diverse biological processes, including circadian rhythms, photoperiodic control of flowering, and responses to abiotic stress, the transcriptional networks associated with these factors in soybean (Glycine max (L.) Merr.) remain incompletely characterized. In this study, we provide empirical evidence highlighting the significance of GmPRR3b as a crucial mediator in regulating the circadian clock, drought stress response, and abscisic acid (ABA) signaling pathway in soybeans. A comprehensive analysis of DNA affinity purification sequencing and transcriptome data identified 795 putative target genes directly regulated by GmPRR3b. Among them, a total of 570 exhibited a significant correlation with the response to drought, and eight genes were involved in both the biosynthesis and signaling pathways of ABA. Notably, GmPRR3b played a pivotal role in the negative regulation of the drought response in soybeans by suppressing the expression of abscisic acid-responsive element-binding factor 3 (GmABF3). Additionally, the overexpression of GmABF3 exhibited an increased ability to tolerate drought conditions, and it also restored the hypersensitive phenotype of the GmPRR3b overexpressor. Consistently, studies on the manipulation of GmPRR3b gene expression and genome editing in plants revealed contrasting reactions to drought stress. The findings of our study collectively provide compelling evidence that emphasizes the significant contribution of the GmPRR3b-GmABF3 module in enhancing drought tolerance in soybean plants. Moreover, the transcriptional network of GmPRR3b provides valuable insights into the intricate interactions between this gene and the fundamental biological processes associated with plant adaptation to diverse environmental conditions.

Glycine max↗
PubMed1997

Amino acid substitutions preserve protein folding by conserving steric and hydrophobicity properties.

We present a comprehensive analysis of amino acid substitution patterns (sets of residues in a position of a multiple alignment) and conservation of physicochemical properties in alignments of protein sequences. Of the one million possible substitution patterns, only a few hundred account for the majority of aligned positions. Very similar distributions of substitution patterns are observed in all but one of the diverse databases of multiple alignments. In these substitution patterns we analyzed the conservation of 511 physicochemical and steric amino acid properties. Highest conservation was observed in those steric and transfer free energy-related properties that are crucial for folding. The best conserved steric properties include the minimal width of the side chains and their interactions with other residues. Among the hydrophobicity-related properties, charge and those properties that provide information on propensities to form secondary structures or side chain conformation, appear to be better conserved than pure hydrophobicity measures. Physicochemical sequence analysis based on the most conserved properties is expected to aid searching a protein sequence query against a database of multiple alignments, prediction of secondary and tertiary structures and protein engineering.

Amino Acids↗
PubMed2004

Phytoestrogens and their human metabolites show distinct agonistic and antagonistic properties on estrogen receptor alpha (ERalpha) and ERbeta in human cells.

Phytoestrogens exert pleiotropic effects on cellular signaling and show some beneficial effects on estrogen-dependent diseases. However, due to activation/inhibition of the estrogen receptors ERalpha or ERbeta, these compounds may induce or inhibit estrogen action and, therefore, have the potential to disrupt estrogen signaling. We performed a comprehensive analysis and potency comparison of phytoestrogens and their human metabolites for ER binding, induction/suppression of ERalpha and ERbeta transactivation, and coactivator recruitment in human cells. The soy-derived genistein, coumestrol, and equol displayed a preference for transactivation of ERbeta compared to ERalpha and were 10- to 100-fold less potent than diethylstilbestrol. In contrast, zearalenone was the most potent phytoestrogen tested and activated preferentially ERalpha. All other phytoestrogens tested, including resveratrol and human metabolites of daidzein and enterolactone, were weak ER agonists. Interestingly, the daidzein metabolites 3',4',7-isoflavone and 4',6,7-isoflavone were superagonists on ERalpha and ERbeta. All phytoestrogens tested showed reduced potencies to activate ERalpha and ERbeta compared to diethylstilbestrol on the estrogen-responsive C3 promoter compared to a consensus estrogen response element indicating a degree of promoter dependency. Zearalenone and resveratrol were antagonistic on both ERalpha and ERbeta at high doses. The phytoestrogens enhanced preferentially recruitment of GRIP1 to ERalpha similar to 17beta-estradiol. In contrast, for ERbeta no distinct preference for one coactivator (GRIP1 or SRC-1) was apparent and the overall coactivator association was less pronounced than for ERalpha. Due to their abundance and (anti)-estrogenic potencies, the soy-derived isoflavones, coumestrol, resveratrol, and zearalenone would appear to have the potential for effectively functioning as endocrine disruptors.

Cells, Cultured↗
PubMed2006

Whole-genome expression profiling defines the HrpL regulon of Pseudomonas syringae pv. tomato DC3000, allows de novo reconstruction of the Hrp cis clement, and identifies novel coregulated genes.

Pseudomonas syringae pv. tomato DC3000 is a model pathogen of tomato and Arabidopsis that uses a hypersensitive response and pathogenicity (Hrp) type III secretion system (T3SS) to deliver virulence effector proteins into host cells. Expression of the Hrp system and many effector genes is activated by the HrpL alternative sigma factor. Here, an open reading frame-specific whole-genome microarray was constructed for DC3000 and used to comprehensively identify genes that are differentially expressed in wild-type and deltahrpL strains. Among the genes whose differential regulation was statistically significant, 119 were upregulated and 76 were downregulated in the wild-type compared with the deltahrpL strain. Hierarchical clustering revealed a subset of eight genes that were upregulated particularly rapidly. Gibbs sampling of regions upstream of HrpL-activated operons revealed the Hrp promoter as the only identifiable regulatory motif and supported an iterative refinement involving real-time polymerase chain reaction testing of additional HrpL-activated genes and refinements in a hidden Markov model that can be used to predict Hrp promoters in P. syringae strains. This iterative bioinformatic-experimental approach to a comprehensive analysis of the HrpL regulon revealed a mix of genes controlled by HrpL, including those encoding most type III effectors, twin-arginine transport (TAT) substrates, other regulatory proteins, and proteins involved in the synthesis or metabolism of phytohormones, phytotoxins, and myo-inositol. This analysis provides an extensively verified, robust method for predicting Hrp promoters in P. syringae genomes, and it supports subsequent identification of effectors and other factors that likely are important to the host-specific virulence of P. syringae.

Bacterial Proteins↗
PubMed2006

Protein expression profiles in an endosymbiotic cyanobacterium revealed by a proteomic approach.

Molecular mechanisms behind adaptations in the cyanobacterium (Nostoc sp.) to a life in endosymbiosis with plants are still not clarified, nor are the interactions between the partners. To get further insights, the proteome of a Nostoc strain, freshly isolated from the symbiotic gland tissue of the angiosperm Gunnera manicata Linden, was analyzed and compared with the proteome of the same strain when free-living. Extracted proteins were separated by two-dimensional gel electrophoresis and were identified by matrix-assisted laser desorption/ionization-time of flight mass spectrometry combined with tandem mass spectrometry. Even when the higher percentage of differentiated cells (heterocysts) in symbiosis was compensated for, the majority of the proteins detected in the symbiotic cyanobacteria were present in the free-living counterpart, indicating that most cellular processes were common for both stages. However, differential expression profiling revealed a significant number of proteins to be down-regulated or missing in the symbiotic stage, while others were more abundant or only expressed in symbiosis. The differential protein expression was primarily connected to i) cell envelope-associated processes, including proteins involved in exopolysaccharide synthesis and surface and membrane associated proteins, ii) to changes in growth and metabolic activities (C and N), including upregulation of nitrogenase and proteins involved in the oxidative pentose phosphate pathway and downregulation of Calvin cycle enzymes, and iii) to the dark, microaerobic conditions offered inside the Gunnera gland cells, including changes in relative phycobiliprotein concentrations. This is the first comprehensive analysis of proteins in the symbiotic state.

Bacterial Proteins↗
PubMed2000

Stagonospora avenae secretes multiple enzymes that hydrolyze oat leaf saponins.

The phytopathogenic fungus Stagonospora avenae is able to infect oat leaves despite the presence of avenacoside saponins in the leaf tissue. In response to pathogen attack, avenacosides are converted into 26-desglucoavenacosides (26-DGAs), which possess antifungal activity. These molecules are comprised of a steroidal backbone linked to a branched sugar chain consisting of one alpha-L-rhamnose and two (avenacoside A) or three (avenacoside B) beta-D-glucose residues. Isolates of the fungus that are pathogenic to oats are capable of sequential hydrolysis of the sugar residues from the 26-DGAs. Degradation is initiated by removal of the L-rhamnose, which abolishes antifungal activity. The D-glucose residues are then hydrolyzed by beta-glucosidase activity. A comprehensive analysis of saponin-hydrolyzing activities was undertaken, and it was established that S. avenae isolate WAC1293 secretes three enzymes, one alpha-rhamnosidase and two beta-glucosidases, that carry out this hydrolysis. The major beta-glucosidase was purified and the gene encoding the enzyme cloned. The protein is similar to saponin-hydrolyzing enzymes produced by three other phytopathogenic fungi, Gaeumannomyces graminis, Septoria lycopersici, and Botrytis cinerea, and is a family 3 beta-glucosidase. The gene encoding the beta-glucosidase is expressed during infection of oat leaves but is not essential for pathogenicity.

Amino Acid Sequence↗
PubMed2002

Potato spindle tuber viroid strains of different pathogenicity induces and suppresses expression of common and unique genes in infected tomato.

Viroids are the smallest plant pathogens. These RNAs do not encode proteins and are not encapsidated, and yet they can replicate autonomously, move systemically, and cause diseases in infected plants. Notably, strains of a viroid with subtle differences in nucleotide sequences can cause dramatically different symptoms in infected plants. These features make viroids unique probes to investigate the role of a pathogenic RNA genome in triggering host responses. We conducted a comprehensive analysis of the differential gene expression patterns of tomato plants at various stages of infection by a mild and severe strain of Potato spindle tuber viroid (PSTVd). We also compared tomato gene expression altered by the PSTVd strains with that altered by Tobacco mosaic virus (TMV). Our analyses revealed that the two PSTVd strains altered expression of both common and unique tomato genes. These genes encode products involved in defense/stress response, cell wall structure, chloroplast function, protein metabolism, and other diverse functions. Five genes have unknown functions. Four genes are novel. The expression of some but not all of these genes was also altered by TMV infection. Our results indicate that viroids, although structurally simple, can trigger complex host responses. Further characterization of viroid-altered gene expression in a host plant should help understand viroid pathogenicity and, potentially, the mechanisms of RNA-mediated regulation of plant gene expression.

Gene Expression Regulation, Developmental↗