RegulonDB v8.0: omics data sets, evolutionary conservation, regulatory phrases, cross-validated gold standards and more.

This article summarizes our progress with RegulonDB (http://regulondb.ccg.unam.mx/) during the past 2 years. We have kept up-to-date the knowledge from the published literature regarding transcriptional regulation in Escherichia coli K-12. We have maintained and expanded our curation efforts to improve the breadth and quality of the encoded experimental knowledge, and we have implemented criteria for the quality of our computational predictions. Regulatory phrases now provide high-level descriptions of regulatory regions. We expanded the assignment of quality to various sources of evidence, particularly for knowledge generated through high-throughput (HT) technology. Based on our analysis of most relevant methods, we defined rules for determining the quality of evidence when multiple independent sources support an entry. With this latest release of RegulonDB, we present a new highly reliable larger collection of transcription start sites, a result of our experimental HT genome-wide efforts. These improvements, together with several novel enhancements (the tracks display, uploading format and curational guidelines), address the challenges of incorporating HT-generated knowledge into RegulonDB. Information on the evolutionary conservation of regulatory elements is also available now. Altogether, RegulonDB version 8.0 is a much better home for integrating knowledge on gene regulation from the sources of information currently available.

Exopolyphosphatase of Pseudomonas aeruginosa is essential for the production of virulence factors, and its expression is controlled by NtrC and PhoB acting at two interspaced promoters.

The exopolyphosphatase (Ppx) of Pseudomonas aeruginosa is encoded by the PA5241 gene (ppx). Ppx catalyses the hydrolysis of inorganic polyphosphates to orthophosphate (Pi). In the present work, we identified and characterized the promoter region of ppx and its regulation under environmental stress conditions. The role of Ppx in the production of several virulence factors was demonstrated through studies performed on a ppx null mutant. We found that ppx is under the control of two interspaced promoters, dually regulated by nitrogen and phosphate limitation. Under nitrogen-limiting conditions, its expression was controlled from a σ(54)-dependent promoter activated by the response regulator NtrC. However, under Pi limitation, the expression was controlled from a σ(70) promoter, activated by PhoB. Results obtained from the ppx null mutant demonstrated that Ppx is involved in the production of virulence factors associated with both acute infection (e.g. motility-promoting factors, blue/green pigment production, C6-C12 quorum-sensing homoserine lactones) and chronic infection (e.g. rhamnolipids, biofilm formation). Molecular and physiological approaches used in this study indicated that P. aeruginosa maintains consistently proper levels of Ppx regardless of environmental conditions. The precise control of ppx expression appeared to be essential for the survival of P. aeruginosa and the occurrence of either acute or chronic infection in the host.

Homorepeat variability within the human population.

Genetic variation within populations plays a crucial role in driving evolution. Unlike the average protein sequence, the evolution of homorepeats can be influenced by DNA replication slippage, when DNA polymerases either add or skip repeats of nucleotides. While there are some diseases known to be caused by abnormal changes in the length of amino acid homorepeats, naturally occurring variations in homorepeat length remain relatively unexplored. In our study, we examined the variation in amino acid homorepeat length of human individuals by analyzing 125 748 exomes, as well as 15 708 whole genomes. Our analyses revealed significant variability in homorepeat length across the human population, indicating that these motifs are prone to mutations at higher rates than non repeat sequences. We focused our study on glutamine homorepeats, also known as polyQ sequences, and found that shorter polyQ sequences tend to exhibit greater length variation, while longer ones primarily undergo deletions. Notably, polyQ sequencesthat are more conserved across primates tend to show less variation within the human population, indicating stronger selective pressure to maintain their length. Overall, our results demonstrate that there is large natural variation in the length of homorepeats within the human population, with no apparent impact on observable traits.

RegulonDB version 7.0: transcriptional regulation of Escherichia coli K-12 integrated within genetic sensory response units (Gensor Units).

RegulonDB (http://regulondb.ccg.unam.mx/) is the primary reference database of the best-known regulatory network of any free-living organism, that of Escherichia coli K-12. The major conceptual change since 3 years ago is an expanded biological context so that transcriptional regulation is now part of a unit that initiates with the signal and continues with the signal transduction to the core of regulation, modifying expression of the affected target genes responsible for the response. We call these genetic sensory response units, or Gensor Units. We have initiated their high-level curation, with graphic maps and superreactions with links to other databases. Additional connectivity uses expandable submaps. RegulonDB has summaries for every transcription factor (TF) and TF-binding sites with internal symmetry. Several DNA-binding motifs and their sizes have been redefined and relocated. In addition to data from the literature, we have incorporated our own information on transcription start sites (TSSs) and transcriptional units (TUs), obtained by using high-throughput whole-genome sequencing technologies. A new portable drawing tool for genomic features is also now available, as well as new ways to download the data, including web services, files for several relational database manager systems and text files including BioPAX format.

Genetic changes during a laboratory adaptive evolution process that allowed fast growth in glucose to an Escherichia coli strain lacking the major glucose transport system.

BACKGROUND: Escherichia coli strains lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS), which is the major bacterial component involved in glucose transport and its phosphorylation, accumulate high amounts of phosphoenolpyruvate that can be diverted to the synthesis of commercially relevant products. However, these strains grow slowly in glucose as sole carbon source due to its inefficient transport and metabolism. Strain PB12, with 400% increased growth rate, was isolated after a 120 hours adaptive laboratory evolution process for the selection of faster growing derivatives in glucose. Analysis of the genetic changes that occurred in the PB12 strain that lacks PTS will allow a better understanding of the basis of its growth adaptation and, therefore, in the design of improved metabolic engineering strategies for enhancing carbon diversion into the aromatic pathways. RESULTS: Whole genome analyses using two different sequencing methodologies: the Roche NimbleGen Inc. comparative genome sequencing technique, and high throughput sequencing with Illumina Inc. GAIIx, allowed the identification of the genetic changes that occurred in the PB12 strain. Both methods detected 23 non-synonymous and 22 synonymous point mutations. Several non-synonymous mutations mapped in regulatory genes (arcB, barA, rpoD, rna) and in other putative regulatory loci (yjjU, rssA and ypdA). In addition, a chromosomal deletion of 10,328 bp was detected that removed 12 genes, among them, the rppH, mutH and galR genes. Characterization of some of these mutated and deleted genes with their functions and possible functions, are presented. CONCLUSIONS: The deletion of the contiguous rppH, mutH and galR genes that occurred simultaneously, is apparently the main reason for the faster growth of the evolved PB12 strain. In support of this interpretation is the fact that inactivation of the rppH gene in the parental PB11 strain substantially increased its growth rate, very likely by increasing glycolytic mRNA genes stability. Furthermore, galR inactivation allowed glucose transport by GalP into the cell. The deletion of mutH in an already stressed strain that lacks PTS is apparently responsible for the very high mutation rate observed.

Correction: Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

[This corrects the article DOI: 10.1371/journal.pone.0249773.].

High-throughput transcriptome sequencing and comparative analysis of Escherichia coli and Schizosaccharomyces pombe in respiratory and fermentative growth.

In spite of increased complexity in eukaryotes compared to prokaryotes, several basic metabolic and regulatory processes are conserved. Here we explored analogies in the eubacteria Escherichia coli and the unicellular fission yeast Schizosaccharomyces pombe transcriptomes under two carbon sources: 2% glucose; or a mix of 2% glycerol and 0.2% sodium acetate using the same growth media and growth phase. Overall, twelve RNA-seq libraries were constructed. A total of 593 and 860 genes were detected as differentially expressed for E. coli and S. pombe, respectively, with a log2 of the Fold Change ≥ 1 and False Discovery Rate ≤ 0.05. In aerobic glycolysis, most of the expressed genes were associated with cell proliferation in both organisms, including amino acid metabolism and glycolysis. In contrast in glycerol/acetate condition, genes related to flagellar assembly and membrane proteins were differentially expressed such as the general transcription factors fliA, flhD, flhC, and flagellum assembly genes were detected in E. coli, whereas in S. pombe genes for hexose transporters, integral membrane proteins, galactose metabolism, and ncRNAs related to cellular stress were overexpressed. In general, our study shows that a conserved "foraging behavior" response is observed in these eukaryotic and eubacterial organisms in gluconeogenic carbon sources.

Specialization of the Reiterated Copies of the Heterodimeric Integration Host Factor Genes in Geobacter sulfurreducens.

Integration host factor (IHF) is a widely distributed small heterodimeric protein member of the bacterial Nucleoid-Associated Proteins (NAPs), implicated in multiple DNA regulatory processes. IHF recognizes a specific DNA sequence and induces a large bend of the nucleic acid. IHF function has been mainly linked with the regulation of RpoN-dependent promoters, where IHF commonly recognizes a DNA sequence between the enhancer-binding region and the promoter, facilitating a close contact between the upstream bound activator and the promoter bound, RNA polymerase. In most proteobacteria, the genes encoding IHF subunits (ihfA and ihfB) are found in a single copy. However, in some Deltaproteobacteria, like Geobacter sulfurreducens, those genes are duplicated. To date, the functionality of IHF reiterated encoding genes is unknown. In this work, we achieved the functional characterization of the ihfA-1, ihfA-2, ihfB-1, and ihfB-2 from G. sulfurreducens. Unlike the ΔihfA-2 or ΔihfB-1 strains, single gene deletion in ihfA-1 or ihfB-2, provokes an impairment in fumarate and Fe(III) citrate reduction. Accordingly, sqRT-PCR experiments showed that ihfA-1 and ihfB-2 were expressed at higher levels than ihfA-2 and ihfB-1. In addition, RNA-Seq analysis of the ΔihfA-1 and ΔihfB-2 strains revealed a total of 89 and 122 differentially expressed genes, respectively. Furthermore, transcriptional changes in 25 genes were shared in both mutant strains. Among these genes, we confirmed the upregulation of the pilA-repressor, GSU1771, and downregulation of the triheme-cytochrome (pgcA) and the aconitate hydratase (acnA) genes by RT-qPCR. EMSA experiments also demonstrated the direct binding of IHF to the upstream promoter regions of GSU1771, pgcA and acnA. PilA changes in ΔihfA-1 and ΔihfB-2 strains were also verified by immunoblotting. Additionally, heme-staining of subcellular fractions in ΔihfA-1 and ΔihfB-2 strains revealed a remarkable deficit of c-type cytochromes. Overall, our data indicate that at least during fumarate and Fe(III) citrate reduction, the functional IHF regulator is likely assembled by the products of ihfA-1 and ihfB-2. Also, a role of IHF controlling expression of multiple genes (other than RpoN-dependent) affects G. sulfurreducens physiology and extracellular electron transfer.

Transcriptome Analysis Reveals Cr(VI) Adaptation Mechanisms in Klebsiella sp. Strain AqSCr.

Klebsiella sp. strain AqSCr, isolated from Cr(VI)-polluted groundwater, reduces Cr(VI) both aerobically and anaerobically and resists up 34 mM Cr(VI); this resistance is independent of the ChrA efflux transporter. In this study, we report the whole genome sequence and the transcriptional profile by RNA-Seq of strain AqSCr under Cr(VI)-adapted conditions and found 255 upregulated and 240 downregulated genes compared to controls without Cr(VI) supplementation. Genes differentially transcribed were mostly associated with oxidative stress response, DNA repair and replication, sulfur starvation response, envelope-osmotic stress response, fatty acid (FA) metabolism, ribosomal subunits, and energy metabolism. Among them, genes not previously associated with chromium resistance, for example, cybB, encoding a putative superoxide oxidase (SOO), gltA2, encoding an alternative citrate synthase, and des, encoding a FA desaturase, were upregulated. The sodA gene encoding a manganese superoxide dismutase was upregulated in the presence of Cr(VI), whereas sodB encoding an iron superoxide dismutase was downregulated. Cr(VI) resistance mechanisms in strain AqSCr seem to be orchestrated by the alternative sigma factors fecl, rpoE, and rpoS (all of them upregulated). Membrane lipid analysis of the Cr(IV)-adapted strain showed a lower proportion of unsaturated lipids with respect to the control, which we hypothesized could result from unsaturated lipid peroxidation followed by degradation, together with de novo synthesis mediated by the upregulated FA desaturase-encoding gene, des. This report helps to elucidate both Cr(VI) toxicity targets and global bacterial response to Cr(VI).

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

RegulonDB (version 6.0): gene regulation model of Escherichia coli K-12 beyond transcription, active (experimental) annotated promoters and Textpresso navigation.

RegulonDB (http://regulondb.ccg.unam.mx/) is the primary reference database offering curated knowledge of the transcriptional regulatory network of Escherichia coli K12, currently the best-known electronically encoded database of the genetic regulatory network of any free-living organism. This paper summarizes the improvements, new biology and new features available in version 6.0. Curation of original literature is, from now on, up to date for every new release. All the objects are supported by their corresponding evidences, now classified as strong or weak. Transcription factors are classified by origin of their effectors and by gene ontology class. We have now computational predictions for sigma(54) and five different promoter types of the sigma(70) family, as well as their corresponding -10 and -35 boxes. In addition to those curated from the literature, we added about 300 experimentally mapped promoters coming from our own high-throughput mapping efforts. RegulonDB v.6.0 now expands beyond transcription initiation, including RNA regulatory elements, specifically riboswitches, attenuators and small RNAs, with their known associated targets. The data can be accessed through overviews of correlations about gene regulation. RegulonDB associated original literature, together with more than 4000 curation notes, can now be searched with the Textpresso text mining engine.

Selection for unequal densities of sigma70 promoter-like signals in different regions of large bacterial genomes.

The evolutionary processes operating in the DNA regions that participate in the regulation of gene expression are poorly understood. In Escherichia coli, we have established a sequence pattern that distinguishes regulatory from nonregulatory regions. The density of promoter-like sequences, that could be recognizable by RNA polymerase and may function as potential promoters, is high within regulatory regions, in contrast to coding regions and regions located between convergently transcribed genes. Moreover, functional promoter sites identified experimentally are often found in the subregions of highest density of promoter-like signals, even when individual sites with higher binding affinity for RNA polymerase exist elsewhere within the regulatory region. In order to see the generality of this pattern, we have analyzed 43 additional genomes belonging to most established bacterial phyla. Differential densities between regulatory and nonregulatory regions are detectable in most of the analyzed genomes, with the exception of those that have evolved toward extreme genome reduction. Thus, presence of this pattern follows that of genes and other genomic features that require weak selection to be effective in order to persist. On this basis, we suggest that the loss of differential densities in the reduced genomes of host-restricted pathogens and symbionts is an outcome of the process of genome degradation resulting from the decreased efficiency of purifying selection in highly structured small populations. This implies that the differential distribution of promoter-like signals between regulatory and nonregulatory regions detected in large bacterial genomes confers a significant, although small, fitness advantage. This study paves the way for further identification of the specific types of selective constraints that affect the organization of regulatory regions and the overall distribution of promoter-like signals through more detailed comparative analyses among closely related bacterial genomes.

Marine Archaeon Methanosarcina acetivorans Enhances Polyphosphate Metabolism Under Persistent Cadmium Stress.

Phosphate metabolism was studied to determine whether polyphosphate (polyP) pools play a role in the enhanced resistance against Cd2+ and metal-removal capacity of Cd2+-preadapted (CdPA) Methanosarcina acetivorans. Polyphosphate kinase (PPK), exopolyphosphatase (PPX) and phosphate transporter transcript levels and their activities increased in CdPA cells compared to control (Cnt) cells. K+ inhibited recombinant Ma-PPK and activated Ma-PPX, whereas divalent cations activated both enzymes. Metal-binding polyP and thiol-containing molecule contents, Cd2+-removal, and biofilm synthesis were significantly higher in CdPA cells >Cnt cells plus a single addition of Cd2+>Cnt cells. Also, CdPA cells showed a higher number of cadmium, sulfur, and phosphorus enriched-acidocalcisomes than control cells. Biochemical and physiological phenotype exhibited by CdPA cells returned to that of Cnt cells when cultured without Cd2+. Furthermore, no differences in the sequenced genomes upstream and downstream of the genes involved in Cd2+ resistance were found between CdPA and Cnt cells, suggesting phenotype loss rather than genome mutations induced by chronic Cd2+-exposure. Instead, a metabolic adaptation induced by Cd2+ stress was apparent. The dynamic ability of M. acetivorans to change its metabolism, depending on the environmental conditions, may be advantageous to remove cadmium in nature and biodigesters.

Inferring modules of functionally interacting proteins using the Bond Energy Algorithm.

BACKGROUND: Non-homology based methods such as phylogenetic profiles are effective for predicting functional relationships between proteins with no considerable sequence or structure similarity. Those methods rely heavily on traditional similarity metrics defined on pairs of phylogenetic patterns. Proteins do not exclusively interact in pairs as the final biological function of a protein in the cellular context is often hold by a group of proteins. In order to accurately infer modules of functionally interacting proteins, the consideration of not only direct but also indirect relationships is required. In this paper, we used the Bond Energy Algorithm (BEA) to predict functionally related groups of proteins. With BEA we create clusters of phylogenetic profiles based on the associations of the surrounding elements of the analyzed data using a metric that considers linked relationships among elements in the data set. RESULTS: Using phylogenetic profiles obtained from the Cluster of Orthologous Groups of Proteins (COG) database, we conducted a series of clustering experiments using BEA to predict (upper level) relationships between profiles. We evaluated our results by comparing with COG's functional categories, And even more, with the experimentally determined functional relationships between proteins provided by the DIP and ECOCYC databases. Our results demonstrate that BEA is capable of predicting meaningful modules of functionally related proteins. BEA outperforms traditionally used clustering methods, such as k-means and hierarchical clustering by predicting functional relationships between proteins with higher accuracy. CONCLUSION: This study shows that the linked relationships of phylogenetic profiles obtained by BEA is useful for detecting functional associations between profiles and extending functional modules not found by traditional methods. BEA is capable of detecting relationship among phylogenetic patterns by linking them through a common element shared in a group. Additionally, we discuss how the proposed method may become more powerful if other criteria to classify different levels of protein functional interactions, as gene neighborhood or protein fusion information, is provided.

Systematic discovery of analogous enzymes in thiamin biosynthesis.

In all genome-sequencing projects completed to date, a considerable number of 'gaps' have been found in the biochemical pathways of the respective species. In many instances, missing enzymes are displaced by analogs, functionally equivalent proteins that have evolved independently and lack sequence and structural similarity. Here we fill such gaps by analyzing anticorrelating occurrences of genes across species. Our approach, applied to the thiamin biosynthesis pathway comprising approximately 15 catalytic steps, predicts seven instances in which known enzymes have been displaced by analogous proteins. So far we have verified four predictions by genetic complementation, including three proteins for which there was no previous experimental evidence of a role in the thiamin biosynthesis pathway. For one hypothetical protein, biochemical characterization confirmed the predicted thiamin phosphate synthase (ThiE) activity. The results demonstrate the ability of our computational approach to predict specific functions without taking into account sequence similarity.

Acetate biostimulation as an effective treatment for cleaning up alkaline soil highly contaminated with Cr(VI).

Stimulation of microbial reduction of Cr(VI) to the less toxic and less soluble Cr(III) through electron donor addition has been regarded as a promising approach for the remediation of chromium-contaminated soil and groundwater sites. However, each site presents different challenges; local physicochemical characteristics and indigenous microbial communities influence the effectiveness of the biostimulation processes. Here, we show microcosm assays stimulation of microbial reduction of Cr(VI) in highly alkaline and saline soil samples from a long-term contaminated site in Guanajuato, Mexico. Acetate was effective promoting anaerobic microbial reduction of 15 mM of Cr(VI) in 25 days accompanied by an increase in pH from 9 to 10. Our analyses showed the presence of Halomonas, Herbaspirillum, Nesterenkonia/Arthrobacter, and Bacillus species in the soil sample collected. Moreover, from biostimulated soil samples, it was possible to isolate Halomonas spp. strains able to grow at 32 mM of Cr(VI). Additionally, we found that polluted groundwater has bacterial species different to those found in soil samples with the ability to resist and reduce chromate using acetate and yeast extract as electron donors.

The unphosphorylated form of the PilR two-component system regulates pilA gene expression in Geobacter sulfurreducens.

In Geobacter sulfurreducens, metal reduction and generation of bioelectricity require the participation of several elements, and among them, the type IV pili has an essential role. The pilus is composed of multiple PilA monomers. Expression of pilA gene depends mainly on the σ54 factor and the response regulator protein PilR. In this work, we characterized the role of the PilS-PilR two-component system in the regulation of the pilA gene expression. Experimental evidence indicates that PilS is autophosphorylated at the His-334 residue, which in turn is transferred to the conserved Asp-53 in PilR. Contrary to other PilS-PilR systems, substitution D53N in PilR resulted in higher activation of the pilA gene. By using a pilA::luxCDABE fusion with different promoter fragments and in vitro DNA-binding assays, we demonstrated the existence of multiple functional PilR binding sites. A regulatory model in which the non-phosphorylated PilR protein directs activation of pilA expression by binding to two sites in the promoter region of this gene is presented.

Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing.

Understanding the genetic structure of Native American populations is important to clarify their diversity, demographic history, and to identify genetic factors relevant for biomedical traits. Here, we show a demographic history reconstruction from 12 Native American whole genomes belonging to six distinct ethnic groups representing the three main described genetic clusters of Mexico (Northern, Southern, and Maya). Effective population size estimates of all Native American groups remained below 2,000 individuals for up to 10,000 years ago. The proportion of missense variants predicted as damaging is higher for undescribed (~ 30%) than for previously reported variants (~ 15%). Several variants previously associated with biological traits are highly frequent in the Native American genomes. These findings suggest that the demographic and adaptive processes that occurred in these groups shaped their genetic architecture and could have implications in biological processes of the Native Americans and Mestizos of today.

Insights on the evolution of trehalose biosynthesis.

BACKGROUND: The compatible solute trehalose is a non-reducing disaccharide, which accumulates upon heat, cold or osmotic stress. It was commonly accepted that trehalose is only present in extremophiles or cryptobiotic organisms. However, in recent years it has been shown that although higher plants do not accumulate trehalose at significant levels they have actively transcribed genes encoding the corresponding biosynthetic enzymes. RESULTS: In this study we show that trehalose biosynthesis ability is present in eubacteria, archaea, plants, fungi and animals. In bacteria there are five different biosynthetic routes, whereas in fungi, plants and animals there is only one. We present phylogenetic analyses of the trehalose-6-phosphate synthase (TPS) and trehalose-phosphatase (TPP) domains and show that there is a close evolutionary relationship between these domains in proteins from diverse organisms. In bacteria TPS and TPP genes are clustered, whereas in eukaryotes these domains are fused in a single protein. CONCLUSION: We have demonstrated that trehalose biosynthesis pathways are widely distributed in nature. Interestingly, several eubacterial species have multiple pathways, while eukaryotes have only the TPS/TPP pathway. Vertebrates lack trehalose biosynthetic capacity but can catabolise it. TPS and TPP domains have evolved mainly in parallel and it is likely that they have experienced several instances of gene duplication and lateral gene transfer.

The genome project of Taenia solium.

We have constituted a consortium of key laboratories at the National Autonomous University of Mexico to carry out a genomic project for Taenia solium. This project will provide powerful resources for the study of taeniasis/cysticercosis, and, in conjunction with the Echinococcus granulosus and Echinococcus multilocularis genome project of expressed sequence tags (ESTs), will mark the advent of genomics for cestode parasites. Our project is planned in two consecutive stages. The first stage is being carried out to determine some basic parameters of the T. solium genome. Afterwards, we will evaluate the best strategy for the second stage, a full blown genome project. We have estimated the T. solium genome size by two different approaches: cytofluorometry on isolated cyton nuclei, as well as a probabilistic calculation based on approximately 2000 sequenced genomic clones, approximately 3000 ESTs, resulting in size estimates of 270 and 251 Mb, respectively. In terms of sequencing, our goal for the first stage is to characterize several thousand EST's (from adult worm and cysticerci cDNA libraries) and genomic clones. Results obtained so far from about 16,000 sequenced ESTs from the adult stage, show that only about 40% of the T. solium coding sequences have a previously sequenced homologue. Many of the best hits are found with mammalian genes, especially with humans. However, 1.5% of the hits lack homologues in humans, making these genes immediate candidates for investigation on pharmaco-therapy, diagnostics and vaccination. Most T. solium ESTs are related to gene regulation, and signal transduction. Other important functions are housekeeping, metabolism, cell division, cytoskeleton, proteases, vacuolar transport, hormone response, and extracellular matrix activities. Preliminary results also suggest that the genome of T. solium is not highly repetitive.