The Aeromonas salmonicida subsp. salmonicida exoproteome: determination of the complete repertoire of Type-Three Secretion System effectors and identification of other virulence factors

BACKGROUND Aeromonas salmonicida subsp. salmonicida, the etiologic agent of furunculosis, is a major pathogen of fisheries worldwide. Several virulence factors have been described, but the type-three secretion system (T3SS) is recognized as having a major effect on virulence by injecting effectors directly into fish cells. In this study we used high-throughput proteomics to display the differences between in vitro secretome of A. salmonicida wild-type (wt, hypervirulent, JF2267) and T3SS-deficient (isogenic ΔascV, extremely low-virulent, JF2747) strains in exponential and stationary phases of growth. RESULTS Results confirmed the secretion of effectors AopH, AexT, AopP and AopO via T3SS, and for the first time demonstrated the impact of T3SS in secretion of Ati2, AopN and ExsE that are known as effectors in other pathogens. Translocators, needle subunits, Ati1, and AscX were also secreted in supernatants (SNs) dependent on T3SS. AopH, Ati2, AexT, AopB and AopD were in the top seven most abundant excreted proteins. EF-G, EF-Tu, DnaK, HtpG, PNPase, PepN and MdeA were moderately secreted in wt SNs and predicted to be putative T3 effectors by bioinformatics. Pta and ASA_P5G088 were increased in wt SNs and T3-associated in other bacteria. Ten conserved cytoplasmic proteins were more abundant in wt SNs than in the ΔascV mutant, but without any clear association to a secretion system. T1-secreted proteins were predominantly found in wt SNs: OmpAI, OmpK40, DegQ, insulinase ASA_0716, hypothetical ASA_0852 and ASA_3619. Presence of T3SS components in pellets was clearly decreased by ascV deletion, while no impact was observed on T1- and T2SS. Our results demonstrated that the ΔascV mutant strain excreted well-described (VapA, AerA, AerB, GCAT, Pla1, PlaC, TagA, Ahe2, GbpA and enolase) and yet uncharacterized potential toxins, adhesins and enzymes as much as or even more than the wt strain. Other putative important virulence factors were not detected. CONCLUSIONS We demonstrated the whole in vitro secretome and T3SS repertoire of hypervirulent A. salmonicida. Several toxins, adhesins and enzymes that are not part of the T3SS secretome were secreted to a higher extent in the extremely low-virulent ΔascV mutant. All together, our results show the high importance of an intact T3SS to initiate the furunculosis and offer new information about the pathogenesis.

A METAGENOMIC STUDY OF THE TICK MIDGUT

A Metagenomic Study of the Tick Midgut Daniel Yuan, B.S. Supervisory Professor : Steven J. Norris, Ph.D. Southern tick–associated rash illness (STARI) or Master’s disease is a Lyme-like illness that occurs following bites by Amblyomma americanum, the lone-star tick. Clinical symptoms include a bull’s eye rash similar to the erythema migrans lesions of Lyme disease, as well as fever and joint pains. Lyme disease is caused by Borrelia burgdorferi and related spirochetes. However, B. burgdorferi has not been detected in STARI patients, or in ticks in the South Central U.S. The causative agent of STARI has not been identified, although it was once thought to be caused by another Borrelia species, Borrelia lonestari. Furthermore, while adult A. americanum have up to a 5.6% Borrelia lonestari infection rate, the prevalence of all Borrelia species in Texas ticks as a whole is not known. Previous studies indicate that 6%-30% of Northern Ixodes scapularis ticks are infected by Borrelia burgdorferi while only 10% of Northern A. americanum and I. scapularis ticks are infected by Borrelia species. The first specific aim of this project was to determine the bacterial community that inhabits the midgut of Texas and Northeastern ticks by using high throughput metagenomic sequencing to sequence bacterial 16S rDNA. Through the use of massively parallel 454 sequencing, we were able to individually sequence hundreds of thousands of 16S rDNA regions of the bacterial flora from 133 ticks from the New York, Missouri and Texas. The presence of previously confirmed endosymbionts, specifically the Rickettsia spp. and Coxiella spp., that are commonly found in ticks were confirmed, as well as some highly prevalent genera that were previously undocumented. Furthermore, multiple pathogenic genera sequences were often found in the same tick, suggesting the possibility of co-infection of multiple pathogenic species. The second specific aim was to use Borrelia specific primers to screen 344 individual ticks from Missouri, Texas and the Northeast to determine the prevalence of Borrelia species in ticks. To screen for Borrelia species, two housekeeping genes, uvrA and recG, were selected as well as the 16S-23S rDNA intergenic spacer. Ticks from Missouri, Texas and New York were screened. None of the Missouri or Texas ticks tested positive for Borrelia spp. The rate of I. scapularis infection by B.burgdorferi is dependent on tick feeding activity as well as reservoir availability. B. burgdorferi is endemic in the Northeast, sometimes reported as highly present in over 50% of all I. scapularis ticks. 11.6% of all New York ticks were positive for a species of Borrelia, however only 6.9% of all New York ticks were positive for B. burgdorferi. Despite being significantly lower than 50%, the results still fall in line with previous reports of about the prevalence of B. burgdorferi. 1.5% of all Texas ticks were positive for a Borrelia species, specifically B. lonestari. While this study was unable to identify the causative agent for STARI, 454 sequencing was able to provide a tremendous insight into the bacterial flora and possible pathogenic species of both the I. scapularis and the A. americanum tick.

NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

Angiotensin receptor agonistic autoantibody is highly prevalent in preeclampsia: correlation with disease severity.

Preeclampsia (PE), a syndrome affecting 5% of pregnancies, characterized by hypertension and proteinuria, is a leading cause of maternal and fetal morbidity and mortality. The condition is often accompanied by the presence of a circulating maternal autoantibody, the angiotensin II type I receptor agonistic autoantibody (AT(1)-AA). However, the prevalence of AT(1)-AA in PE remains unknown, and the correlation of AT(1)-AA titers with the severity of the disease remains undetermined. We used a sensitive and high-throughput luciferase bioassay to detect AT(1)-AA levels in the serum of 30 normal, 37 preeclamptic (10 mild and 27 severe), and 23 gestational hypertensive individuals. Here we report that AT(1)-AA is highly prevalent in PE ( approximately 95%). Next, by comparing the levels of AT(1)-AA among women with mild and severe PE, we found that the titer of AT(1)-AA is proportional to the severity of the disease. Intriguingly, among severe preeclamptic patients, we discovered that the titer of AT(1)-AA is significantly correlated with the clinical features of PE: systolic blood pressure (r=0.56), proteinuria (r=0.70), and soluble fms-like tyrosine kinase-1 level (r=0.71), respectively. Notably, only AT(1)-AA, and not soluble fms-like tyrosine kinase-1, levels are elevated in gestational hypertensive patients. These data serve as compelling clinical evidence that AT(1)-AA is highly prevalent in PE, and its titer is strongly correlated to the severity of the disease.

Intact flagellar motor of Borrelia burgdorferi revealed by cryo-electron tomography: evidence for stator ring curvature and rotor/C-ring assembly flexion.

The bacterial flagellar motor is a remarkable nanomachine that provides motility through flagellar rotation. Prior structural studies have revealed the stunning complexity of the purified rotor and C-ring assemblies from flagellar motors. In this study, we used high-throughput cryo-electron tomography and image analysis of intact Borrelia burgdorferi to produce a three-dimensional (3-D) model of the in situ flagellar motor without imposing rotational symmetry. Structural details of B. burgdorferi, including a layer of outer surface proteins, were clearly visible in the resulting 3-D reconstructions. By averaging the 3-D images of approximately 1,280 flagellar motors, a approximately 3.5-nm-resolution model of the stator and rotor structures was obtained. flgI transposon mutants lacked a torus-shaped structure attached to the flagellar rod, establishing the structural location of the spirochetal P ring. Treatment of intact organisms with the nonionic detergent NP-40 resulted in dissolution of the outermost portion of the motor structure and the C ring, providing insight into the in situ arrangement of the stator and rotor structures. Structural elements associated with the stator followed the curvature of the cytoplasmic membrane. The rotor and the C ring also exhibited angular flexion, resulting in a slight narrowing of both structures in the direction perpendicular to the cell axis. These results indicate an inherent flexibility in the rotor-stator interaction. The FliG switching and energizing component likely provides much of the flexibility needed to maintain the interaction between the curved stator and the relatively symmetrical rotor/C-ring assembly during flagellar rotation.

Complete Genome Sequence of Treponema paraluiscuniculi, Strain Cuniculi A: The Loss of Infectivity to Humans Is Associated with Genome Decay.

Treponema paraluiscuniculi is the causative agent of rabbit venereal spirochetosis. It is not infectious to humans, although its genome structure is very closely related to other pathogenic Treponema species including Treponema pallidum subspecies pallidum, the etiological agent of syphilis. In this study, the genome sequence of Treponema paraluiscuniculi, strain Cuniculi A, was determined by a combination of several high-throughput sequencing strategies. Whereas the overall size (1,133,390 bp), arrangement, and gene content of the Cuniculi A genome closely resembled those of the T. pallidum genome, the T. paraluiscuniculi genome contained a markedly higher number of pseudogenes and gene fragments (51). In addition to pseudogenes, 33 divergent genes were also found in the T. paraluiscuniculi genome. A set of 32 (out of 84) affected genes encoded proteins of known or predicted function in the Nichols genome. These proteins included virulence factors, gene regulators and components of DNA repair and recombination. The majority (52 or 61.9%) of the Cuniculi A pseudogenes and divergent genes were of unknown function. Our results indicate that T. paraluiscuniculi has evolved from a T. pallidum-like ancestor and adapted to a specialized host-associated niche (rabbits) during loss of infectivity to humans. The genes that are inactivated or altered in T. paraluiscuniculi are candidates for virulence factors important in the infectivity and pathogenesis of T. pallidum subspecies.

Determinants of Acidobacteria activity inferred from the relative abundances of 16S rRNA transcripts in German grassland and forest soils

16S rRNA genes and transcripts of Acidobacteria were investigated in 57 grassland and forest soils of three different geographic regions. Acidobacteria contributed 9-31% of bacterial 16S rRNA genes whereas the relative abundances of the respective transcripts were 4-16%. The specific cellular 16S rRNA content (determined as molar ratio of rRNA:rRNA genes) ranged between 3 and 80, indicating a low in situ growth rate. Correlations with flagellate numbers, vascular plant diversity and soil respiration suggest that biotic interactions are important determinants of Acidobacteria 16S rRNA transcript abundances in soils. While the phylogenetic composition of Acidobacteria differed significantly between grassland and forest soils, high throughput denaturing gradient gel electrophoresis and terminal restriction fragment length polymorphism fingerprinting detected 16S rRNA transcripts of most phylotypes in situ. Partial least squares regression suggested that chemical soil conditions such as pH, total nitrogen, C:N ratio, ammonia concentrations and total phosphorus affect the composition of this active fraction of Acidobacteria. Transcript abundance for individual Acidobacteria phylotypes was found to correlate with particular physicochemical (pH, temperature, nitrogen or phosphorus) and, most notably, biological parameters (respiration rates, abundances of ciliates or amoebae, vascular plant diversity), providing culture-independent evidence for a distinct niche specialization of different Acidobacteria even from the same subdivision.

Monitoring of alcohol markers by capillary electrophoresis.

Work dealing with the monitoring of alcohol markers by CE performed during the past two decades led to the development of assays for carbohydrate-deficient transferrin (CDT), ethyl sulfate, ethyl glucuronide, and phosphatidylethanol in body fluids and first attempts for the detection of the urinary 5-hydroxytryptophol/5-hydroxyindoleacetic acid ratio and stable hemoglobin acetaldehyde adducts. Most notably are assays for CDT that have been commercialized and are being used in many laboratories under routine conditions. This paper provides insight into the development, specifications, and use of the currently known CE-based assays suitable to detect alcohol markers. The achievements reached so far indicate that CE is an attractive technology for monitoring alcohol markers. This is particularly seen with the CDT assays that do not require an elaborate sample pretreatment and thus could be fully automated for high-throughput analyses on multicapillary instruments.

SNPs for parentage testing and traceability in globally diverse breeds of sheep

DNA-based parentage determination accelerates genetic improvement in sheep by increasing pedigree accuracy. Single nucleotide polymorphism (SNP) markers can be used for determining parentage and to provide unique molecular identifiers for tracing sheep products to their source. However, the utility of a particular "parentage SNP" varies by breed depending on its minor allele frequency (MAF) and its sequence context. Our aims were to identify parentage SNPs with exceptional qualities for use in globally diverse breeds and to develop a subset for use in North American sheep. Starting with genotypes from 2,915 sheep and 74 breed groups provided by the International Sheep Genomics Consortium (ISGC), we analyzed 47,693 autosomal SNPs by multiple criteria and selected 163 with desirable properties for parentage testing. On average, each of the 163 SNPs was highly informative (MAF≥0.3) in 48±5 breed groups. Nearby polymorphisms that could otherwise confound genetic testing were identified by whole genome and Sanger sequencing of 166 sheep from 54 breed groups. A genetic test with 109 of the 163 parentage SNPs was developed for matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. The scoring rates and accuracies for these 109 SNPs were greater than 99% in a panel of North American sheep. In a blinded set of 96 families (sire, dam, and non-identical twin lambs), each parent of every lamb was identified without using the other parent's genotype. In 74 ISGC breed groups, the median estimates for probability of a coincidental match between two animals (PI), and the fraction of potential adults excluded from parentage (PE) were 1.1×10(-39) and 0.999987, respectively, for the 109 SNPs combined. The availability of a well-characterized set of 163 parentage SNPs facilitates the development of high-throughput genetic technologies for implementing accurate and economical parentage testing and traceability in many of the world's sheep breeds.

Identification of animal Pasteurellaceae by MALDI-TOF mass spectrometry

Species of the family Pasteurellaceae play an important role as primary or opportunistic, predominantly respiratory, pathogens in domestic and wild animals. Some of them cause severe disease with high economic losses in commercial animal husbandry. Hence, rapid and accurate differentiation of Pasteurellaceae is important and signifies a particular challenge to diagnostic laboratories. Identification and differentiation of Pasteurellaceae is mostly done using phenotypic tests or genetic identification based on sequence similarity of housekeeping genes, such as the rrs gene encoding the 16S ribosomal RNA (16S rRNA). Both approaches are time consuming, laborious, and costly, therefore often delaying the final diagnosis of disease or epidemics. Matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry represents an alternative rapid and reliable method for the differentiation of most members of the family Pasteurellaceae. It is able to differentiate within a few minutes the currently known 18 genera and most of the over 60 species and subspecies of Pasteurellaceae including many members encountered in veterinary diagnostic laboratories. A few closely related species and subspecies that cannot be discriminated by MALDI-TOF are easily identified further by complementary simple tests, such as hemolysis done simultaneously or routinely during pathogen isolation.

A genomic perspective on a new bacterial genus and species from the Alcaligenaceae family, Basilea psittacipulmonis

BACKGROUND A novel Gram-negative, non-haemolytic, non-motile, rod-shaped bacterium was discovered in the lungs of a dead parakeet (Melopsittacus undulatus) that was kept in captivity in a petshop in Basel, Switzerland. The organism is described with a chemotaxonomic profile and the nearly complete genome sequence obtained through the assembly of short sequence reads. RESULTS Genome sequence analysis and characterization of respiratory quinones, fatty acids, polar lipids, and biochemical phenotype is presented here. Comparison of gene sequences revealed that the most similar species is Pelistega europaea, with BLAST identities of only 93% to the 16S rDNA gene, 76% identity to the rpoB gene, and a similar GC content (~43%) as the organism isolated from the parakeet, DSM 24701 (40%). The closest full genome sequences are those of Bordetella spp. and Taylorella spp. High-throughput sequencing reads from the Illumina-Solexa platform were assembled with the Edena de novo assembler to form 195 contigs comprising the ~2 Mb genome. Genome annotation with RAST, construction of phylogenetic trees with the 16S rDNA (rrs) gene sequence and the rpoB gene, and phylogenetic placement using other highly conserved marker genes with ML Tree all suggest that the bacterial species belongs to the Alcaligenaceae family. Analysis of samples from cages with healthy parakeets suggested that the newly discovered bacterial species is not widespread in parakeet living quarters. CONCLUSIONS Classification of this organism in the current taxonomy system requires the formation of a new genus and species. We designate the new genus Basilea and the new species psittacipulmonis. The type strain of Basilea psittacipulmonis is DSM 24701 (= CIP 110308 T, 16S rDNA gene sequence Genbank accession number JX412111 and GI 406042063).

Monogenic Human Skin Disorders

Human genodermatoses represent a broad and partly confusing spectrum of countless rare diseases with confluent and overlapping phenotypes often impeding a precise diagnosis in an affected individual. High-throughput sequencing techniques have expedited the identification of novel genes and have dramatically simplified the establishment of genetic diagnoses in such heterogeneous disorders. The precise genetic diagnosis of a skin disorder is crucial for the appropriate counselling of patients and their relatives regarding the course of the disease, prognosis and recurrence risks. Understanding the underlying pathophysiology is a prerequisite to understanding the disease and developing specific, targeted or individualized therapeutic approaches. We aimed to create a comprehensive overview of human genodermatoses and their respective genetic aetiology known to date. We hope this may represent a useful tool in guiding dermatologists towards genetic diagnoses, providing patients with individual knowledge on the respective disorder and applying novel research findings to clinical practice.

Normal platelet activation profile in patients with peripheral arterial disease on aspirin.

BACKGROUND Peripheral arterial disease (PAD) is a progressive vascular disease associated with a high risk of cardiovascular morbidity and death. Antithrombotic prevention is usually applied by prescribing the antiplatelet agent aspirin. However, in patients with PAD aspirin fails to provide protection against myocardial infarction and death, only reducing the risk of ischemic stroke. Platelets may play a role in disease development, but this has not been tested by proper mechanistic studies. In the present study, we performed a systematic evaluation of platelet reactivity in whole blood from patients with PAD using two high-throughput assays, i.e. multi-agonist testing of platelet activation by flow cytometry and multi-parameter testing of thrombus formation on spotted microarrays. METHODS Blood was obtained from 40 patients (38 on aspirin) with PAD in majority class IIa/IIb and from 40 age-matched control subjects. Whole-blood flow cytometry and multiparameter thrombus formation under high-shear flow conditions were determined using recently developed and validated assays. RESULTS Flow cytometry of whole blood samples from aspirin-treated patients demonstrated unchanged high platelet responsiveness towards ADP, slightly elevated responsiveness after glycoprotein VI stimulation, and decreased responsiveness after PAR1 thrombin receptor stimulation, compared to the control subjects. Most parameters of thrombus formation under flow were similarly high for the patient and control groups. However, in vitro aspirin treatment caused a marked reduction in thrombus formation, especially on collagen surfaces. When compared per subject, markers of ADP- and collagen-induced integrin activation (flow cytometry) strongly correlated with parameters of collagen-dependent thrombus formation under flow, indicative of a common, subject-dependent regulation of both processes. CONCLUSION Despite of the use of aspirin, most platelet activation properties were in the normal range in whole-blood from class II PAD patients. These data underline the need for more effective antithrombotic pharmacoprotection in PAD.

Metabolomics reveals herbivore-induced metabolites of resistance and susceptibility in maize leaves and roots

Plants respond to herbivory by reprogramming their metabolism. Most research in this context has focused on locally induced compounds that function as toxins or feeding deterrents. We developed an ultra-high-pressure liquid chromatography time-of-flight mass spectrometry (UHPLC-TOF-MS)-based metabolomics approach to evaluate local and systemic herbivore-induced changes in maize leaves, sap, roots and root exudates without any prior assumptions about their function. Thirty-two differentially regulated compounds were identified from Spodoptera littoralis-infested maize seedlings and isolated for structure assignment by microflow nuclear magnetic resonance (CapNMR). Nine compounds were quantified by a high throughput direct nano-infusion tandem mass spectrometry/mass spectrometry (MS/MS) method. Leaf infestation led to a marked local increase of 1,3-benzoxazin-4-ones, phospholipids, N-hydroxycinnamoyltyramines, azealic acid and tryptophan. Only few changes were found in the root metabolome, but 1,3-benzoxazin-4-ones increased in the vascular sap and root exudates. The role of N-hydroxycinnamoyltyramines in plant–herbivore interactions is unknown, and we therefore tested the effect of the dominating p-coumaroyltyramine on S. littoralis. Unexpectedly, p-coumaroyltyramine was metabolized by the larvae and increased larval growth, possibly by providing additional nitrogen to the insect. Taken together, this study illustrates that herbivore attack leads to the induction of metabolites that can have contrasting effects on herbivore resistance in the leaves and roots.

Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

Background Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest. Methods Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea. Results We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics. Conclusions The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.