Arabidopsis-Based Dual-Layered Biological Network Analysis Elucidates Fully Modulated Pathways Related to Sugarcane Resistance on Biotrophic Pathogen Infection.

We assembled a dual-layered biological network to study the roles of resistance gene analogs (RGAs) in the resistance of sugarcane to infection by the biotrophic fungus causing smut disease. Based on sugarcane-Arabidopsis orthology, the modeling used metabolic and protein-protein interaction (PPI) data from Arabidopsis thaliana (from Kyoto Encyclopedia of Genes and Genomes (KEGG) and BioGRID databases) and plant resistance curated knowledge for Viridiplantae obtained through text mining of the UniProt/SwissProt database. With the network, we integrated functional annotations and transcriptome data from two sugarcane genotypes that differ significantly in resistance to smut and applied a series of analyses to compare the transcriptomes and understand both signal perception and transduction in plant resistance. We show that the smut-resistant sugarcane has a larger arsenal of RGAs encompassing transcriptionally modulated subnetworks with other resistance elements, reaching hub proteins of primary metabolism. This approach may benefit molecular breeders in search of markers associated with quantitative resistance to diseases in non-model systems.

Characterization of genes responsive to osmotic and oxidative stresses of the sugarcane bacterial pathogen Leifsonia xyli subsp. xyli.

Leifsonia xyli subsp. xyli (Lxx) colonizes the xylem vessels of sugarcane, a plant niche where microorganisms are highly exposed to oxidative and osmotic stresses. This study performed an in silico analysis of the genome of Lxx and characterized 16 genes related to the detoxification of oxidative species (peroxidases, O2- dismutases, and methionine reductases) and to the production and transport of osmolytes and analyzed their expression in vitro after 30, 60, and 120 min of exposure to H2O2 or PEG. The PAGE activity of superoxide dismutase (Mn-SOD as confirmed by inhibition tests) and of catalase (CAT) and the accumulation of trehalose were also assessed. Exposure to H2O2 increased the expression of most oxidative-responsive genes and decreased the expression of those related to osmotic responses, whereas the opposite occurred after exposure to PEG. The isoform profiles of CAT and Mn-SOD shifted in response to H2O2 but not to PEG and Lxx cells accumulated more trehalose over time after exposure to PEG compared with non-exposed cells. The experimental results validated the in silico analysis and indicated that this obligate endophytic parasite has multiple and functional mechanisms to combat the stresses imposed by its host.

Genome survey of resistance gene analogs in sugarcane: genomic features and differential expression of the innate immune system from a smut-resistant genotype.

BACKGROUND: Resistance genes composing the two-layer immune system of plants are thought as important markers for breeding pathogen-resistant crops. Many have been the attempts to establish relationships between the genomic content of Resistance Gene Analogs (RGAs) of modern sugarcane cultivars to its degrees of resistance to diseases such as smut. However, due to the highly polyploid and heterozygous nature of sugarcane genome, large scale RGA predictions is challenging. RESULTS: We predicted, searched for orthologs, and investigated the genomic features of RGAs within a recently released sugarcane elite cultivar genome, alongside the genomes of sorghum, one sugarcane ancestor (Saccharum spontaneum), and a collection of de novo transcripts generated for six modern cultivars. In addition, transcriptomes from two sugarcane genotypes were obtained to investigate the roles of RGAs differentially expressed (RGADE) in their distinct degrees of resistance to smut. Sugarcane references lack RGAs from the TNL class (Toll-Interleukin receptor (TIR) domain associated to nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains) and harbor elevated content of membrane-associated RGAs. Up to 39% of RGAs were organized in clusters, and 40% of those clusters shared synteny. Basically, 79% of predicted NBS-encoding genes are located in a few chromosomes. S. spontaneum chromosome 5 harbors most RGADE orthologs responsive to smut in modern sugarcane. Resistant sugarcane had an increased number of RGAs differentially expressed from both classes of RLK (receptor-like kinase) and RLP (receptor-like protein) as compared to the smut-susceptible. Tandem duplications have largely contributed to the expansion of both RGA clusters and the predicted clades of RGADEs. CONCLUSIONS: Most of smut-responsive RGAs in modern sugarcane were potentially originated in chromosome 5 of the ancestral S. spontaneum genotype. Smut resistant and susceptible genotypes of sugarcane have a distinct pattern of RGADE. TM-LRR (transmembrane domains followed by LRR) family was the most responsive to the early moment of pathogen infection in the resistant genotype, suggesting the relevance of an innate immune system. This work can help to outline strategies for further understanding of allele and paralog expression of RGAs in sugarcane, and the results should help to develop a more applied procedure for the selection of resistant plants in sugarcane.

Genome-wide alternative splicing landscapes modulated by biotrophic sugarcane smut pathogen.

Alternative splicing (AS) promotes transcriptome and proteome diversity during growth, development, and stress responses in eukaryotes. Genome-wide studies of AS in sugarcane (Saccharum spp.) are lacking, mainly due to the absence of a high-quality sequenced reference genome, sugarcane's large, complex genome, and the variable chromosome numbers and polyploidy of sugarcane cultivars. Here, we analyzed changes in the sugarcane isoform-level transcriptome and AS landscape during infection with the smut fungus (Sporisorium scitamineum) using a hybrid approach involving Sorghum bicolor reference-based and Trinity de novo mapping tools. In total, this analysis detected 16,039 and 15,379 transcripts (≥2 FPKM) at 5 and 200 days after infection, respectively. A conservative estimate of isoform-level expression suggested that approximately 5,000 (14%) sugarcane genes undergo AS. Differential expression analysis of the alternatively spliced genes in healthy and smut-infected sugarcane revealed 896 AS events modulated at different stages of infection. Gene family and gene ontology functional enrichment analysis of the differentially spliced genes revealed overrepresentation of functional categories related to the cell wall, defense, and redox homeostasis pathways. Our study provides novel insight into the AS landscape of sugarcane during smut disease interactions.

Corrigendum: Revisiting Meiosis in Sugarcane: Chromosomal Irregularities and the Prevalence of Bivalent Configurations.

[This corrects the article DOI: 10.3389/fgene.2018.00213.].

Revisiting Meiosis in Sugarcane: Chromosomal Irregularities and the Prevalence of Bivalent Configurations.

Traditional sugarcane cultivars (Saccharum officinarum) proved highly susceptible to diseases, and this led breeders to progress to interspecific crosses resulting in disease resistance. A backcrossing program to S. officinarum was then required to boost sucrose content. Clonal selection across generations and incorporation of other germplasm into cultivated backgrounds established the (narrow) genetic base of modern cultivars (Saccharum spp.), which have a man-made genome. The genome complexity has inspired several molecular studies that have elucidated aspects of sugarcane genome constitution, architecture, and cytogenetics. However, there is a critical shortage of information on chromosome behavior throughout meiosis in modern cultivars. In this study, we examined the microsporogenesis of a contemporary variety, providing a detailed analysis of the meiotic process and chromosome association at diakinesis, using FISH with centromeric probes. Chromosomal abnormalities were documented by examining high quality preparations of pollen mother cells (700 in total). Approximately 70% of the cells showed abnormalities, such as metaphase chromosomes not lined up at the plate, lagging chromosomes and chromosomal bridges, and tetrad cells with micronuclei. Some dyads with asynchronous behavior were also observed. Due to the hybrid composition of the sugarcane genome, we suggest that bivalent incomplete pairing may occur in the first prophase leading to univalency. The presence of rod bivalents showing the lagging tendency is consistent with a reduction in chiasma frequency. Finally, the presence of chromatin bridges indicates the indirect occurrence of chromosomal inversions, although chromosome fragments were not clearly recognized. Possible reasons for such meiotic abnormalities and the large prevalence of bivalent formation are discussed.

Comparative Genomics of Smut Pathogens: Insights From Orphans and Positively Selected Genes Into Host Specialization.

Host specialization is a key evolutionary process for the diversification and emergence of new pathogens. However, the molecular determinants of host range are poorly understood. Smut fungi are biotrophic pathogens that have distinct and narrow host ranges based on largely unknown genetic determinants. Hence, we aimed to expand comparative genomics analyses of smut fungi by including more species infecting different hosts and to define orphans and positively selected genes to gain further insights into the genetics basis of host specialization. We analyzed nine lineages of smut fungi isolated from eight crop and non-crop hosts: maize, barley, sugarcane, wheat, oats, Zizania latifolia (Manchurian rice), Echinochloa colona (a wild grass), and Persicaria sp. (a wild dicot plant). We assembled two new genomes: Ustilago hordei (strain Uhor01) isolated from oats and U. tritici (strain CBS 119.19) isolated from wheat. The smut genomes were of small sizes, ranging from 18.38 to 24.63 Mb. U. hordei species experienced genome expansions due to the proliferation of transposable elements and the amount of these elements varied among the two strains. Phylogenetic analysis confirmed that Ustilago is not a monophyletic genus and, furthermore, detected misclassification of the U. tritici specimen. The comparison between smut pathogens of crop and non-crop hosts did not reveal distinct signatures, suggesting that host domestication did not play a dominant role in shaping the evolution of smuts. We found that host specialization in smut fungi likely has a complex genetic basis: different functional categories were enriched in orphans and lineage-specific selected genes. The diversification and gain/loss of effector genes are probably the most important determinants of host specificity.

Metabolome Dynamics of Smutted Sugarcane Reveals Mechanisms Involved in Disease Progression and Whip Emission.

Sugarcane smut disease, caused by the biotrophic fungus Sporisorium scitamineum, is characterized by the development of a whip-like structure from the plant meristem. The disease causes negative effects on sucrose accumulation, fiber content and juice quality. The aim of this study was to exam whether the transcriptomic changes already described during the infection of sugarcane by S. scitamineum result in changes at the metabolomic level. To address this question, an analysis was conducted during the initial stage of the interaction and through disease progression in a susceptible sugarcane genotype. GC-TOF-MS allowed the identification of 73 primary metabolites. A set of these compounds was quantitatively altered at each analyzed point as compared with healthy plants. The results revealed that energetic pathways and amino acid pools were affected throughout the interaction. Raffinose levels increased shortly after infection but decreased remarkably after whip emission. Changes related to cell wall biosynthesis were characteristic of disease progression and suggested a loosening of its structure to allow whip growth. Lignin biosynthesis related to whip formation may rely on Tyr metabolism through the overexpression of a bifunctional PTAL. The altered levels of Met residues along with overexpression of SAM synthetase and ACC synthase genes suggested a role for ethylene in whip emission. Moreover, unique secondary metabolites antifungal-related were identified using LC-ESI-MS approach, which may have potential biomarker applications. Lastly, a putative toxin was the most important fungal metabolite identified whose role during infection remains to be established.

Sugarcane smut: shedding light on the development of the whip-shaped sorus.

Background and Aims: Sugarcane smut is caused by the fungus Sporisorium scitamineum (Ustilaginales/Ustilaginomycotina/Basidiomycota), which is responsible for losses in sugarcane production worldwide. Infected plants show a profound metabolic modification resulting in the development of a whip-shaped structure (sorus) composed of a mixture of plant tissues and fungal hyphae. Within this structure, ustilospores develop and disseminate the disease. Despite the importance of this disease, a detailed histopathological analysis of the plant-pathogen interaction is lacking. Methods: The whip-shaped sorus was investigated using light microscopy, scanning and transmission electron microscopy, histochemical tests and epifluorescence microscopy coupled with deconvolution. Key Results: Sorus growth is mediated by intercalary meristem activity at the base of the sorus, where the fungus causes partial host cell wall degradation and formation of intercellular spaces. Sporogenesis in S. scitamineum is thallic, with ustilospore initials in intercalary or terminal positions, and mostly restricted to the base of the sorus. Ustilospore maturation is centrifugal in relation to the ground parenchyma and occurs throughout the sorus median region. At the apex of the sorus, the fungus produces sterile cells and promotes host cell detachment. Hyphae are present throughout the central axis of the sorus (columella). The plant cell produces callose around the intracellular hyphae as well as inside the papillae at the infection site. Conclusions: The ontogeny of the whip-shaped sorus suggests that the fungus can cause the acropetal growth in the intercalary meristem. The sporogenesis of S. scitamineum was described in detail, demonstrating that the spores are formed exclusively at the base of the whip. Light was also shed on the nature of the sterile cells. The presence of the fungus alters the host cell wall composition, promotes its degradation and causes the release of some peripheral cells of the sorus. Finally, callose was observed around fungal hyphae in infected cells, suggesting that deposition of callose by the host may act as a structural response to fungal infection.

Functional analysis of oxidative burst in sugarcane smut-resistant and -susceptible genotypes.

MAIN CONCLUSION: Smut pathogen induced an early modulation of the production and scavenging of reactive oxygen species during defence responses in resistant sugarcane that coincided with the developmental stages of fungal growth. Sporisorium scitamineum is the causal agent of sugarcane smut disease. In this study, we characterized sugarcane reactive oxygen species (ROS) metabolism in response to the pathogen in smut-resistant and -susceptible genotypes. Sporisorium scitamineum teliospore germination and appressorium formation coincided with H2O2 accumulation in resistant plants. The superoxide dismutase (SOD) activity was not responsive in any of the genotypes; however, a higher number of isoenzymes were detected in resistant plants. In addition, related to resistance were lipid peroxidation, a decrease in catalase (CAT), and an increase in glutathione S-transferase (GST) activities and an earlier transcript accumulation of ROS marker genes (CAT3, CATA, CATB, GST31, GSTt3, and peroxidase 5-like). Furthermore, based on proteomic data, we suggested that the source of the increased hydrogen peroxide (H2O2) may be due to a protein of the class III peroxidase, which was inhibited in the susceptible genotype. H2O2 is sensed and probably transduced through overlapping systems related to ascorbate-glutathione and thioredoxin to influence signalling pathways, as revealed by the presence of thioredoxin h-type, ascorbate peroxidase, and guanine nucleotide-binding proteins in the infected resistant plants. Altogether, our data depicted the balance of the oxidative burst and antioxidant enzyme activity in the outcome of this interaction.

Molecular variability and genetic relationship among Brazilian strains of the sugarcane smut fungus.

Sporisorium scitamineum is the fungus that causes sugarcane smut disease. Despite of the importance of sugarcane for Brazilian agribusiness and the persistence of the pathogen in most cropping areas, genetic variation studies are still missing for Brazilian isolates. In this study, sets of isolates were analyzed using two molecular markers (AFLP and telRFLP) and ITS sequencing. Twenty-two whips were collected from symptomatic plants in cultivated sugarcane fields of Brazil. A total of 41 haploid strains of compatible mating types were selected from individual teliospores and used for molecular genetic analyses. telRFLP and ITS analyses were expanded to six Argentine isolates, where the sugarcane smut was first recorded in America. Genetic relationship among strains suggests the human-mediated dispersal of S. scitamineum within the Brazilian territory and between the two neighboring countries. Two genetically distinct groups were defined by the combined analysis of AFLP and telRFLP. The opposite mating-type strains derived from single teliospores were clustered together into these main groups, but had not always identical haplotypes. telRFLP markers analyzed over two generations of selfing and controlled outcrossing confirmed the potential for emergence of new variants and occurrence of recombination, which are relevant events for evolution of virulence and environmental adaptation.

RNAseq Transcriptional Profiling following Whip Development in Sugarcane Smut Disease.

Sugarcane smut disease is caused by the biotrophic fungus Sporisorium scitamineum. The disease is characterized by the development of a whip-like structure from the primary meristems, where billions of teliospores are produced. Sugarcane smut also causes tillering and low sucrose and high fiber contents, reducing cane productivity. We investigated the biological events contributing to disease symptoms in a smut intermediate-resistant sugarcane genotype by examining the transcriptional profiles (RNAseq) shortly after inoculating the plants and immediately after whip emission. The overall picture of disease progression suggests that premature transcriptional reprogramming of the shoot meristem functions continues until the emergence of the whip. The guidance of this altered pattern is potentially primarily related to auxin mobilization in addition to the involvement of other hormonal imbalances. The consequences associated with whip emission are the modulation of typical meristematic functions toward reproductive organ differentiation, requiring strong changes in carbon partitioning and energy production. These changes include the overexpression of genes coding for invertases and trehalose-6P synthase, as well as other enzymes from key metabolic pathways, such as from lignin biosynthesis. This is the first report describing changes in the transcriptional profiles following whip development, providing a hypothetical model and candidate genes to further study sugarcane smut disease progression.

Complete Genome Sequence of Sporisorium scitamineum and Biotrophic Interaction Transcriptome with Sugarcane.

Sporisorium scitamineum is a biotrophic fungus responsible for the sugarcane smut, a worldwide spread disease. This study provides the complete sequence of individual chromosomes of S. scitamineum from telomere to telomere achieved by a combination of PacBio long reads and Illumina short reads sequence data, as well as a draft sequence of a second fungal strain. Comparative analysis to previous available sequences of another strain detected few polymorphisms among the three genomes. The novel complete sequence described herein allowed us to identify and annotate extended subtelomeric regions, repetitive elements and the mitochondrial DNA sequence. The genome comprises 19,979,571 bases, 6,677 genes encoding proteins, 111 tRNAs and 3 assembled copies of rDNA, out of our estimated number of copies as 130. Chromosomal reorganizations were detected when comparing to sequences of S. reilianum, the closest smut relative, potentially influenced by repeats of transposable elements. Repetitive elements may have also directed the linkage of the two mating-type loci. The fungal transcriptome profiling from in vitro and from interaction with sugarcane at two time points (early infection and whip emergence) revealed that 13.5% of the genes were differentially expressed in planta and particular to each developmental stage. Among them are plant cell wall degrading enzymes, proteases, lipases, chitin modification and lignin degradation enzymes, sugar transporters and transcriptional factors. The fungus also modulates transcription of genes related to surviving against reactive oxygen species and other toxic metabolites produced by the plant. Previously described effectors in smut/plant interactions were detected but some new candidates are proposed. Ten genomic islands harboring some of the candidate genes unique to S. scitamineum were expressed only in planta. RNAseq data was also used to reassure gene predictions.

Leptospira interrogans serovar copenhageni harbors two lexA genes involved in SOS response.

Bacteria activate a regulatory network in response to the challenges imposed by DNA damage to genetic material, known as the SOS response. This system is regulated by the RecA recombinase and by the transcriptional repressor lexA. Leptospira interrogans is a pathogen capable of surviving in the environment for weeks, being exposed to a great variety of stress agents and yet retaining its ability to infect the host. This study aims to investigate the behavior of L. interrogans serovar Copenhageni after the stress induced by DNA damage. We show that L. interrogans serovar Copenhageni genome contains two genes encoding putative LexA proteins (lexA1 and lexA2) one of them being potentially acquired by lateral gene transfer. Both genes are induced after DNA damage, but the steady state levels of both LexA proteins drop, probably due to auto-proteolytic activity triggered in this condition. In addition, seven other genes were up-regulated following UV-C irradiation, recA, recN, dinP, and four genes encoding hypothetical proteins. This set of genes is potentially regulated by LexA1, as it showed binding to their promoter regions. All these regions contain degenerated sequences in relation to the previously described SOS box, TTTGN 5CAAA. On the other hand, LexA2 was able to bind to the palindrome TTGTAN10TACAA, found in its own promoter region, but not in the others. Therefore, the L. interrogans serovar Copenhageni SOS regulon may be even more complex, as a result of LexA1 and LexA2 binding to divergent motifs. New possibilities for DNA damage response in Leptospira are expected, with potential influence in other biological responses such as virulence.

Xylella fastidiosa comparative genomic database is an information resource to explore the annotation, genomic features, and biology of different strains.

The Xylella fastidiosa comparative genomic database is a scientific resource with the aim to provide a user-friendly interface for accessing high-quality manually curated genomic annotation and comparative sequence analysis, as well as for identifying and mapping prophage-like elements, a marked feature of Xylella genomes. Here we describe a database and tools for exploring the biology of this important plant pathogen. The hallmarks of this database are the high quality genomic annotation, the functional and comparative genomic analysis and the identification and mapping of prophage-like elements. It is available from web site http://www.xylella.lncc.br.

Xylella fastidiosa comparative genomic database is an information resource to explore the annotation, genomic features, and biology of different strains

The Xylella fastidiosa comparative genomic database is a scientific resource with the aim to provide a user-friendly interface for accessing high-quality manually curated genomic annotation and comparative sequence analysis, as well as for identifying and mapping prophage-like elements, a marked feature of Xylella genomes. Here we describe a database and tools for exploring the biology of this important plant pathogen. The hallmarks of this database are the high quality genomic annotation, the functional and comparative genomic analysis and the identification and mapping of prophage-like elements. It is available from web site http://www.xylella.lncc.br.

The Cryphonectria parasitica mitochondrial rns gene: plasmid-like elements, introns and homing endonucleases.

The mt-rns gene of Cryphonectria parasitica is 9872bp long and includes two group I and two group II introns. An analysis of intronic protein-encoding sequences revealed that LAGLIDADG ORFs, which usually are associated with group I introns, were transferred at least twice into group II introns. A plasmid-like mitochondrial element (plME) that appears in high amounts in previously mutagen-induced mit1 and mit2 hypovirulent mutants of the Ep155 standard virulent strain of C. parasitica was found to be derived from a short region of the mt-rns gene, including the exon 1 and most of the first intron. The plME is a 4.2-kb circular, multimeric DNA and an autonomously-replicating mtDNA fragment. Although sexual transmission experiments indicate that the plME does not directly cause hypovirulence, its emergence is one manifestation of the many complex molecular and genetic events that appear to underlie this phenotype.

Characterization of new IS elements and studies of their dispersion in two subspecies of Leifsonia xyli.

BACKGROUND: Leifsonia xyli is a xylem-inhabiting bacterial species comprised of two subspecies: L. xyli subsp. xyli (Lxx) and L. xyli subsp. cynodontis (Lxc). Lxx is the causal agent of ratoon stunting disease in sugarcane commercial fields and Lxc colonizes the xylem of several grasses causing either mild or no symptoms of disease. The completely sequenced genome of Lxx provided insights into its biology and pathogenicity. Since IS elements are largely reported as an important source of bacterial genome diversification and nothing is known about their role in chromosome architecture of L. xyli, a comparative analysis of Lxc and Lxx elements was performed. RESULTS: Sample sequencing of Lxc genome and comparative analysis with Lxx complete DNA sequence revealed a variable number of IS transposable elements acting upon genomic diversity. A detailed characterization of Lxc IS elements and a comparative review with IS elements of Lxx are presented. Each genome showed a unique set of elements although related to same IS families when considering features such as similarity among transposases, inverted and direct repeats, and element size. Most of the Lxc and Lxx IS families assigned were reported to maintain transposition at low levels using translation regulatory mechanisms, consistent with our in silico analysis. Some of the IS elements were found associated with rearrangements and specific regions of each genome. Differences were also found in the effect of IS elements upon insertion, although none of the elements were preferentially associated with gene disruption. A survey of transposases among genomes of Actinobacteria showed no correlation between phylogenetic relatedness and distribution of IS families. By using Southern hybridization, we suggested that diversification of Lxc isolates is also mediated by insertion sequences in probably recent events. CONCLUSION: Collectively our data indicate that transposable elements are involved in genome diversification of Lxc and Lxx. The IS elements were probably acquired after the divergence of the two subspecies and are associated with genome organization and gene contents. In addition to enhancing understanding of IS element dynamics in general, these data will contribute to our ongoing comparative analyses aimed at understanding the biological differences of the Lxc and Lxx.

Temporal and spatial expression of the myostatin gene during chicken embryo development.

Myostatin is a potent growth and differentiation factor involved in skeletal muscle tissue formation in vertebrates. In the present study, temporal and spatial expression patterns of myostatin transcripts were investigated in chicken embryos. Myostatin mRNA was detected by RT-PCR analysis in embryos collected immediately after oviposition (stage HH1) and persisted until the fifth day of incubation (stage HH26). Whole-mount in situ hybridization revealed myostatin to be expressed in the ventral myotomal region of mature somites, thus confirming the importance of myostatin in skeletal muscle tissue formation during avian embryogenesis. A smaller myostatin transcript was also identified. This transcript appears to have resulted from an alternative splicing event from common GT-AG processing sites. Analysis of the amino acid sequence generated from this alternative transcript confirmed the presence of a truncated protein that lacks the C terminal region, including the cysteine domains characteristic of the TGF-beta super family. The temporal and spatial patterns of myostatin expression presented in this study agree with the proposed role of myostatin as modulator of muscle cell proliferation.

Comparative bioinformatic analysis of genes expressed in common bean (Phaseolus vulgaris L.) seedlings.

To rapidly and cost-effectively generate gene expression data, we developed an annotated unigene database of common bean (Phaseolus vulgaris L.). In this study, 3 cDNA libraries were constructed from the bean breeding line SEL1308, 1 from young leaf and 2 from seedlings inoculated or not inoculated with the fungal pathogen Colletotrichum lindemuthianum (Sacc. & Magnus) Briosi & Cavara, which causes anthracnose in common bean. To this date, 5255 single-pass sequences have been included in the database after selection based on sequence quality. These ESTs were trimmed and clustered using the computer programs Phred and CAP3 to form a unigene collection of 3126 unique sequences. Within clusters, 318 single nucleotide polymorphisms (SNPs) and 68 insertions-deletions (indels) were found, indicating the presence of paralogous gene families in our database. Each unigene sequence was analyzed for possible function using their similarity to known genes represented in the GenBank database and classified into 14 categories. Only 314 unigenes showed significant similarities to Phaseolus genomic sequences and P. vulgaris ESTs, which indicates that 90% (2818 unigenes) of our database represent newly discovered common bean genes. In addition, 12% (387 unigenes) were shown to be specific to common bean. This study represents a first step towards the discovery of novel genes in beans and a valuable source of molecular markers for expressed gene tagging and mapping.