Tackling critical parameters in metazoan meta-barcoding experiments: a preliminary study based on coxI DNA barcode.

Nowadays DNA meta-barcoding is a powerful instrument capable of quickly discovering the biodiversity of an environmental sample by integrating the DNA barcoding approach with High Throughput Sequencing technologies. It mainly consists of the parallel reading of informative genomic fragment/s able to discriminate living entities. Although this approach has been widely studied, it still needs optimization in some necessary steps requested in its advanced accomplishment. A fundamental element concerns the standardization of bioinformatic analyses pipelines. The aim of the present study was to underline a number of critical parameters of laboratory material preparation and taxonomic assignment pipelines in DNA meta-barcoding experiments using the cytochrome oxidase subunit-I (coxI) barcode region, known as a suitable molecular marker for animal species identification. We compared nine taxonomic assignment pipelines, including a custom in-house method, based on Hidden Markov Models. Moreover, we evaluated the potential influence of universal primers amplification bias in qPCR, as well as the correlation between GC content with taxonomic assignment results. The pipelines were tested on a community of known terrestrial invertebrates collected by pitfall traps from a chestnut forest in Italy. Although the present analysis was not exhaustive and needs additional investigation, our results suggest some potential improvements in laboratory material preparation and the introduction of additional parameters in taxonomic assignment pipelines. These include the correct setup of OTU clustering threshold, the calibration of GC content affecting sequencing quality and taxonomic classification, as well as the evaluation of PCR primers amplification bias on the final biodiversity pattern. Thus, careful attention and further validation/optimization of the above-mentioned variables would be required in a DNA meta-barcoding experimental routine.

Towards barcode markers in Fungi: an intron map of Ascomycota mitochondria.

BACKGROUND: A standardized and cost-effective molecular identification system is now an urgent need for Fungi owing to their wide involvement in human life quality. In particular the potential use of mitochondrial DNA species markers has been taken in account. Unfortunately, a serious difficulty in the PCR and bioinformatic surveys is due to the presence of mobile introns in almost all the fungal mitochondrial genes. The aim of this work is to verify the incidence of this phenomenon in Ascomycota, testing, at the same time, a new bioinformatic tool for extracting and managing sequence databases annotations, in order to identify the mitochondrial gene regions where introns are missing so as to propose them as species markers. METHODS: The general trend towards a large occurrence of introns in the mitochondrial genome of Fungi has been confirmed in Ascomycota by an extensive bioinformatic analysis, performed on all the entries concerning 11 mitochondrial protein coding genes and 2 mitochondrial rRNA (ribosomal RNA) specifying genes, belonging to this phylum, available in public nucleotide sequence databases. A new query approach has been developed to retrieve effectively introns information included in these entries. RESULTS: After comparing the new query-based approach with a blast-based procedure, with the aim of designing a faithful Ascomycota mitochondrial intron map, the first method appeared clearly the most accurate. Within this map, despite the large pervasiveness of introns, it is possible to distinguish specific regions comprised in several genes, including the full NADH dehydrogenase subunit 6 (ND6) gene, which could be considered as barcode candidates for Ascomycota due to their paucity of introns and to their length, above 400 bp, comparable to the lower end size of the length range of barcodes successfully used in animals. CONCLUSION: The development of the new query system described here would answer the pressing requirement to improve drastically the bioinformatics support to the DNA Barcode Initiative. The large scale investigation of Ascomycota mitochondrial introns performed through this tool, allowing to exclude the introns-rich sequences from the barcode candidates exploration, could be the first step towards a mitochondrial barcoding strategy for these organisms, similar to the standard approach employed in metazoans.

The RHNumtS compilation: features and bioinformatics approaches to locate and quantify Human NumtS.

BACKGROUND: To a greater or lesser extent, eukaryotic nuclear genomes contain fragments of their mitochondrial genome counterpart, deriving from the random insertion of damaged mtDNA fragments. NumtS (Nuclear mt Sequences) are not equally abundant in all species, and are redundant and polymorphic in terms of copy number. In population and clinical genetics, it is important to have a complete overview of NumtS quantity and location. Searching PubMed for NumtS or Mitochondrial pseudo-genes yields hundreds of papers reporting Human NumtS compilations produced by in silico or wet-lab approaches. A comparison of published compilations clearly shows significant discrepancies among data, due both to unwise application of Bioinformatics methods and to a not yet correctly assembled nuclear genome. To optimize quantification and location of NumtS, we produced a consensus compilation of Human NumtS by applying various bioinformatics approaches. RESULTS: Location and quantification of NumtS may be achieved by applying database similarity searching methods: we have applied various methods such as Blastn, MegaBlast and BLAT, changing both parameters and database; the results were compared, further analysed and checked against the already published compilations, thus producing the Reference Human Numt Sequences (RHNumtS) compilation. The resulting NumtS total 190. CONCLUSION: The RHNumtS compilation represents a highly reliable reference basis, which may allow designing a lab protocol to test the actual existence of each NumtS. Here we report preliminary results based on PCR amplification and sequencing on 41 NumtS selected from RHNumtS among those with lower score. In parallel, we are currently designing the RHNumtS database structure for implementation in the HmtDB resource. In the future, the same database will host NumtS compilations from other organisms, but these will be generated only when the nuclear genome of a specific organism has reached a high-quality level of assembly.

Genome duplication and gene-family evolution: the case of three OXPHOS gene families.

DNA duplication is one of the main forces acting on the evolution of organisms because it creates the raw genetic material that natural selection can subsequently modify. Duplicated regions are mainly due to "errors" in different phases of meiosis, but DNA transposable elements and reverse transcription also contribute to amplify and move the genomic material to different genomic locations. As a result, redundancy affects genomes to variable degrees: from the single gene to the whole genome (WGD). Gene families are clusters of genes created by duplication and their size reflects the number of duplicated genes, called paralogs, in each species. The aim of this review is to describe the state of the art in the identification and analysis of gene families in eukaryotes, with specific attention to those generated by ancient large scale events in vertebrates (WGD or large segmental duplications). As a case study, we report our work on the evolution of gene families encoding subunits of the five OXPHOS (oxidative phosphorylation) complexes, fundamental and highly conserved in all respiring cells. Although OXPHOS gene families are smaller than the general trend in nuclear gene families, some exceptions are observed, such as three gene families with at least two paralogs in vertebrates. These gene families encode cytochrome c (Cyt c, the electron shuttle protein between complex III and IV), Lipid Binding Protein (LBP, the channel protein of complex V which transfers protons through the inner mitochondrial membrane) and the MLRQ subunit (MLRQ, a supernumerary subunit of the large complex I, with unknown function). We provide a two-step approach, based on structural genomic data, to demonstrate that these gene families should have arisen through WGD (or large segmental duplication) events at the origin of vertebrates and, only afterwards, underwent species-specific events of further gene duplications and loss. In summary, this review reflects the need to apply genome comparative approaches, deriving from both "classical" molecular phylogenetic analysis and "new" genome map analysis, to successfully define the complex evolutionary relations between gene family members which, in turn, are essential to obtain any other comparative phylogenetic or functional results.

Variability of the mitochondrial genome in mammals at the inter-species/intra-species boundary.

Genomic variations represent the molecular basis of the biodiversity of living organisms on which selection operates to generate evolution. In eukaryotes, genomic variability can be experienced in both nuclear and organellar, i.e. mitochondrial and plastid (where present), genomes, which can follow completely different evolution pathways, as revealed by comparative genomics analyses. In Metazoa, for which a substantial number of complete genome sequences are available (nuclear, but mainly mitochondrial), we are just starting to grasp the selective pressures operating on some basic features of the genome as a whole. In this brief review, we discuss the variability of the mitochondrial metazoan genome, with particular reference to mitochondrial DNA in mammals. In light of the recent assumption that a small segment of mitochondrial DNA may be used, particularly in Metazoa, as a species marker, some data on mitochondrial gene variability at the inter-species/intra-species boundary are reported. Intra-species variability has been evaluated in four mammalian species, Homo sapiens, Bos taurus, Sus scrofa and Canis familiaris, whereas the relationship between intra- and inter-species variability has been investigated in Bos taurus and Bos indicus.

p53FamTaG: a database resource of human p53, p63 and p73 direct target genes combining in silico prediction and microarray data.

BACKGROUND: The p53 gene family consists of the three genes p53, p63 and p73, which have polyhedral non-overlapping functions in pivotal cellular processes such as DNA synthesis and repair, growth arrest, apoptosis, genome stability, angiogenesis, development and differentiation. These genes encode sequence-specific nuclear transcription factors that recognise the same responsive element (RE) in their target genes. Their inactivation or aberrant expression may determine tumour progression or developmental disease. The discovery of several protein isoforms with antagonistic roles, which are produced by the expression of different promoters and alternative splicing, widened the complexity of the scenario of the transcriptional network of the p53 family members. Therefore, the identification of the genes transactivated by p53 family members is crucial to understand the specific role for each gene in cell cycle regulation. We have combined a genome-wide computational search of p53 family REs and microarray analysis to identify new direct target genes. The huge amount of biological data produced has generated a critical need for bioinformatic tools able to manage and integrate such data and facilitate their retrieval and analysis. DESCRIPTION: We have developed the p53FamTaG database (p53 FAMily TArget Genes), a modular relational database, which contains p53 family direct target genes selected in the human genome searching for the presence of the REs and the expression profile of these target genes obtained by microarray experiments. p53FamTaG database also contains annotations of publicly available databases and links to other experimental data. The genome-wide computational search of the REs was performed using PatSearch, a pattern-matching program implemented in the DNAfan tool. These data were integrated with the microarray results we produced from the overexpression of different isoforms of p53, p63 and p73 stably transfected in isogenic cell lines, allowing the comparative study of the transcriptional activity of all the proteins in the same cellular background.p53FamTaG database is available free at http://www2.ba.itb.cnr.it/p53FamTaG/ CONCLUSION: p53FamTaG represents a unique integrated resource of human direct p53 family target genes that is extensively annotated and provides the users with an efficient query/retrieval system which displays the results of our microarray experiments and allows the export of RE sequences. The database was developed for supporting and integrating high-throughput in silico and experimental analyses and represents an important reference source of knowledge for research groups involved in the field of oncogenesis, apoptosis and cell cycle regulation.

Molecular evolution of the neurotrophin family members and their Trk receptors.

Neurotrophins are structurally related proteins regulating brain development and function. Molecular evolution studies of neurotrophins and their receptors are essential for understanding the mechanisms underlying the coevolution processes of these gene families and how they correlate with the increased complexity of the vertebrate nervous system. In order to improve our current knowledge of the molecular evolution of neurotrophins and receptors, we have collected all information available in the literature and analyzed the genome database for each of them. Statistical analysis of aminoacid and nucleotide sequences of the neurotrophin and Trk family genes was applied to both complete genes and mature sequences, and different phylogenetic methods were used to compare aminoacid and nucleotide sequences variability among the different species. All collected data favor a model in which several rounds of genome duplications might have facilitated the generation of the many different neurotrophins and the acquisition of specific different functions correlated with the increased complexity of the vertebrate nervous system during evolution. We report findings that refine the structure of the evolutionary trees for neurotrophins and Trk receptors families, indicate different rates of evolution for each member of the two families, and newly demonstrate that the NGF-like genes found in Fowlpox and Canarypox viruses are closely related to reptile NGF.

Phylogenetic analyses of complete mitochondrial genome sequences suggest a basal divergence of the enigmatic rodent Anomalurus.

BACKGROUND: Phylogenetic relationships between Lagomorpha, Rodentia and Primates and their allies (Euarchontoglires) have long been debated. While it is now generally agreed that Rodentia constitutes a monophyletic sister-group of Lagomorpha and that this clade (Glires) is sister to Primates and Dermoptera, higher-level relationships within Rodentia remain contentious. RESULTS: We have sequenced and performed extensive evolutionary analyses on the mitochondrial genome of the scaly-tailed flying squirrel Anomalurus sp., an enigmatic rodent whose phylogenetic affinities have been obscure and extensively debated. Our phylogenetic analyses of the coding regions of available complete mitochondrial genome sequences from Euarchontoglires suggest that Anomalurus is a sister taxon to the Hystricognathi, and that this clade represents the most basal divergence among sampled Rodentia. Bayesian dating methods incorporating a relaxed molecular clock provide divergence-time estimates which are consistently in agreement with the fossil record and which indicate a rapid radiation within Glires around 60 million years ago. CONCLUSION: Taken together, the data presented provide a working hypothesis as to the phylogenetic placement of Anomalurus, underline the utility of mitochondrial sequences in the resolution of even relatively deep divergences and go some way to explaining the difficulty of conclusively resolving higher-level relationships within Glires with available data and methodologies.

The MitoDrome database annotates and compares the OXPHOS nuclear genes of Drosophila melanogaster, Drosophila pseudoobscura and Anopheles gambiae.

The oxidative phosphorylation (OXPHOS) is the primary energy-producing process of all aerobic organisms and the only cellular function under the dual control of both the mitochondrial and the nuclear genomes. Functional characterization and evolutionary study of the OXPHOS system is of great importance for the understanding of many as yet unclear aspects of nucleus-mitochondrion genomic co-evolution and co-regulation gene networks. The MitoDrome database is a web-based database which provides genomic annotations about nuclear genes of Drosophila melanogaster encoding for mitochondrial proteins. Recently, MitoDrome has included a new section annotating genomic information about OXPHOS genes in Drosophila pseudoobscura and Anopheles gambiae and their comparative analysis with their Drosophila melanogaster and human counterparts. The introduction of this new comparative annotation section into MitoDrome is expected to be a useful resource for both functional and structural genomics related to the OXPHOS system.

Metazoan OXPHOS gene families: evolutionary forces at the level of mitochondrial and nuclear genomes.

Mitochondrial and nuclear DNAs contribute to encode the whole mitochondrial protein complement. The two genomes possess highly divergent features and properties, but the forces influencing their evolution, even if different, require strong coordination. The gene content of mitochondrial genome in all Metazoa is in a frozen state with only few exceptions and thus mitochondrial genome plasticity especially concerns some molecular features, i.e. base composition, codon usage, evolutionary rates. In contrast the high plasticity of nuclear genomes is particularly evident at the macroscopic level, since its redundancy represents the main feature able to introduce genetic material for evolutionary innovations. In this context, genes involved in oxidative phosphorylation (OXPHOS) represent a classical example of the different evolutionary behaviour of mitochondrial and nuclear genomes. The simple DNA sequence of Cytochrome c oxidase I (encoded by the mitochondrial genome) seems to be able to distinguish intra- and inter-species relations between organisms (DNA Barcode). Some OXPHOS subunits (cytochrome c, subunit c of ATP synthase and MLRQ) are encoded by several nuclear duplicated genes which still represent the trace of an ancient segmental/genome duplication event at the origin of vertebrates.

The fatty acid synthase gene is a conserved p53 family target from worm to human.

The discovery that the p53 family consists of three members (p53, p63 and p73) in vertebrates and of a single homolog in invertebrates has raised the challenge of understanding the functions of the ancestor and how they have evolved and differentiated within the duplicated genes in vertebrates. Here, we report that the fatty acid synthase (FAS) gene, encoding for a key enzyme involved in the biogenesis of membrane lipids in rapidly proliferating cells, is a conserved target of the p53 family throughout the evolution. We show that CEP-1, the C. elegans p53 homolog, is able to bind the two p53 family responsive elements (REs) identified in the worm fasn-1 gene. Moreover, we demonstrate that fasn-1 expression is modulated by CEP-1 in vivo, by comparing wild-type and CEP-1 knockout worms. In human, luciferase and chromatin immunoprecipitation assays demonstrate that TAp73alpha and DeltaNp63alpha, but not p53, TAp73beta and TAp63alpha bind the two p53 REs of the human FASN gene. We show that the ectopic expression of TAp73alpha and DeltaNp63alpha leads to an increase of FASN mRNA levels, while their silencing produces a decrease of FASN expression. Furthermore, we present data showing a correlation between DeltaNp63alpha and FASN expression in cellular proliferation. Of relevant importance is that fasn-1 is the first CEP-1 direct target gene identified so far in C. elegans and our results suggest a new CEP-1 role in cellular proliferation and development, besides the one already described in apoptosis of germ cells. These data confirm the hypothesis that the ancestral functions of the single invertebrate gene may have been spread out among the three vertebrate members, each of them have acquired specific role in cell cycle regulation.

Evolution of ATP synthase subunit c and cytochrome c gene families in selected Metazoan classes.

To investigate the integrated evolution of mitochondrial and nuclear genomes in the eukaryotic cell, we have focused our attention on OXPHOS (oxidative phosphorylation) gene families which encode proteins involved in the main mitochondrial function. The present study reports the phylogenetic analysis of two OXPHOS gene families: ATP synthase subunit c (or lipid binding protein, LBP) and Cytochrome c (Cytc). Both gene families possess a higher expansion trend than the typically low duplication rate of OXPHOS genes in Metazoa, but follow a completely different evolutionary history, especially in mammals. LBP is represented by three well conserved isoforms in all mammals (P1, P2, P3): only P3 possesses a clearly conserved isoform in all Vertebrates, P1 and P2 were already present before the bird-mammal divergence and there are preliminary evidence from the in silico analysis that P1, the most evolutionary divergent isoform, is poorly expressed and not regulated by NRF1. In contrast, Cytc family presents at least two duplicated genes in all the analysed Vertebrates, is subject to a high expansion trend, especially of processed pseudogenes in mammals, and some events of gain and loss of function can be supposed.

MitoRes: a resource of nuclear-encoded mitochondrial genes and their products in Metazoa.

BACKGROUND: Mitochondria are sub-cellular organelles that have a central role in energy production and in other metabolic pathways of all eukaryotic respiring cells. In the last few years, with more and more genomes being sequenced, a huge amount of data has been generated providing an unprecedented opportunity to use the comparative analysis approach in studies of evolution and functional genomics with the aim of shedding light on molecular mechanisms regulating mitochondrial biogenesis and metabolism. In this context, the problem of the optimal extraction of representative datasets of genomic and proteomic data assumes a crucial importance. Specialised resources for nuclear-encoded mitochondria-related proteins already exist; however, no mitochondrial database is currently available with the same features of MitoRes, which is an update of the MitoNuc database extensively modified in its structure, data sources and graphical interface. It contains data on nuclear-encoded mitochondria-related products for any metazoan species for which this type of data is available and also provides comprehensive sequence datasets (gene, transcript and protein) as well as useful tools for their extraction and export. DESCRIPTION: MitoRes http://www2.ba.itb.cnr.it/MitoRes/ consolidates information from publicly external sources and automatically annotates them into a relational database. Additionally, it also clusters proteins on the basis of their sequence similarity and interconnects them with genomic data. The search engine and sequence management tools allow the query/retrieval of the database content and the extraction and export of sequences (gene, transcript, protein) and related sub-sequences (intron, exon, UTR, CDS, signal peptide and gene flanking regions) ready to be used for in silico analysis. CONCLUSION: The tool we describe here has been developed to support lab scientists and bioinformaticians alike in the characterization of molecular features and evolution of mitochondrial targeting sequences. The way it provides for the retrieval and extraction of sequences allows the user to overcome the obstacles encountered in the integrative use of different bioinformatic resources and the completeness of the sequence collection allows intra- and interspecies comparison at different biological levels (gene, transcript and protein).

Evolution of nuclearly encoded mitochondrial genes in Metazoa.

All Metazoan nuclear genomes underwent a continuous process of both complete and partial genetic material gain and loss. The forces modulating these events are also subject to the strict interaction between nuclear and mitochondrial (mt) genome. In this context we investigate the evolution of nuclear genes encoding proteins which target the mitochondrion, with a particular attention to genes involved in oxidative phosphorylation (OXPHOS), one of the most ancient and conserved functions. To examine thoroughly the evolutionary strategies that preserve OXPHOS and coordinate the two cellular genomes, a comparative analysis has been carried out for 78 OXPHOS gene families in several Metazoa (insects, tunicates, fishes and mammals). We demonstrate that the duplication rate of OXPHOS genes increases passing from invertebrates to vertebrates consistently with the total increase in genome size, but all species are prone to negatively select OXPHOS duplicates compared to the general trend of nuclear gene families. These results are consistent with the 'balance hypothesis' and, at least in insects, the expression of duplicate genes is low and strongly testis-biased.

UTRdb and UTRsite: a collection of sequences and regulatory motifs of the untranslated regions of eukaryotic mRNAs.

The 5' and 3' untranslated regions of eukaryotic mRNAs play crucial roles in the post-transcriptional regulation of gene expression through the modulation of nucleo-cytoplasmic mRNA transport, translation efficiency, subcellular localization and message stability. UTRdb is a curated database of 5' and 3' untranslated sequences of eukaryotic mRNAs, derived from several sources of primary data. Experimentally validated functional motifs are annotated (and also collated as the UTRsite database) and cross-links to genomic and protein data are provided. The integration of UTRdb with genomic and protein data has allowed the implementation of a powerful retrieval resource for the selection and extraction of UTR subsets based on their genomic coordinates and/or features of the protein encoded by the relevant mRNA (e.g. GO term, PFAM domain, etc.). All internet resources implemented for retrieval and functional analysis of 5' and 3' untranslated regions of eukaryotic mRNAs are accessible at http://www.ba.itb.cnr.it/UTR/.

Pseudogenes in metazoa: origin and features.

The complete genome sequences with their annotations are a considerable resource in biology, particularly in understanding the global structure of the genetic material at the molecular level. The reason why some eukaryotic genomes contain large quantities of apparently unnecessary DNA, namely pseudogenes, while others seem to invest in more efficient thinning processes or are equipped with protection systems against parasitic elements still remains a mystery. Several genome-wide surveys have been undertaken to identify pseudogenes in the completely sequenced genome, bringing to light some differences both in their amount and distribution. Since pseudogenes are important resources in evolutionary and comparative genomics - as 'molecular fossils' - in this paper, a survey on the origins, features, abundance and localisation of the different pseudogenes is reported. As an example of genes producing processed pseudogenes, some experimental data obtained in the authors' laboratories from the study of a nuclear gene coding for the mitochondrial transcription factor A (mtTFA), a key regulator of mitochondrial biogenesis, are also reported.

DNAfan: a software tool for automated extraction and analysis of user-defined sequence regions.

SUMMARY: DNAfan (DNA Feature ANalyzer) is a tool combining sequence-filtering and pattern searching. DNAfan automatically extracts user-defined sets of sequence fragments from large sequence sets. Fragments are defined by annotated gene feature keys and co- or non-occurring patterns within the feature or close to it. A gene feature parser and a pattern-based filter tool localizes and extracts the specific subset of sequences. The selected sequence data can subsequently be retrieved for analyses or further processed with DNAfan to find the occurrence of specific patterns or structural motifs. DNAfan is a powerful tool for pattern analysis. Its filter features restricts the pattern search to a well-defined set of sequences, allowing drastic reduction in false positive hits. AVAILABILITY: http://bighost.ba.itb.cnr.it:8080/Framework.

The evolution of the adenine nucleotide translocase family.

Homologous genes are grouped into families whose evolution may be different in the various organisms. For the variety of the processes and the well-known mechanism of gene gain and gene loss, which takes place in genome evolution, we deal in comparative analyses with a "one-to-many" or a "many-to-many" relationship between homologous genes going from invertebrates to vertebrates. In this scenario, it is important to understand how gene function has been preserved and in addition the innovations originated in a given lineage or species. The phylogenetic relations between gene family members and their molecular clock behavior may be very helpful to elucidate their functional fates in various organisms. This in turn can direct laboratory experiments and practical applications. In order to track the evolutionary history of the ANT gene family, we have collected and analyzed 46 sequences from fungi to mammals. Phylogenetic analyses have been performed on nucleotide and amino acidic sequences which have produced basically the same results. We observe the presence of multiple isoforms both in lower and higher eukaryotic species, thus a "many-to-many" correspondence between genes. The molecular phylogeny of ANT genes, reported in the present study, allows to date the time of divergence of ANT isoforms in various lineages. Furthermore, the logo analysis has been carried out to characterize the conservation features of ANT proteins particularly in their three similar domains originated by duplication.

Comparative genomics: the evolutionary history of the Bcl-2 family.

The Bcl-2 family is a huge family composed of various members, occurring in all animals, which are key regulators of apoptosis, the cell death program critical for cell survival and development, tissue homeostasis, and protection against pathogens. The members of the Bcl-2 family can be divided into pro-apoptotic and anti-apoptotic proteins. A delicate balance between these members exists in each cell and the regulations of these two groups of proteins determines whether the cell survives or undergoes apoptosis. Bcl-2 family proteins are characterized by distinct domains. All members possess at least one of the four motifs known as Bcl-2 homology domains (BH1 to BH4). Most pro-survival members which can inhibit apoptosis facing a wide variety of cytotoxic insults, contain at least BH1 and BH2 domains; those most similar to Bcl-2 have all four BH domains. All the pro-apoptosis family members possess BH3 domain which is the central domain. For the first time, a global phylogenetic analysis of all Bcl-2 family members is presented here. We have analyzed the genes known so far that have a different composition of the functional domains BH1, BH2, BH3 and BH4. The analyses were performed both on complete sequences (124 sites analyzed) and on single domains. We present the results obtained using both approaches. We have also analyzed the amino acid profile and the degree of conservation of the BH3 domains of pro- and anti-apoptotic proteins. The results of our phylogenetic analyses show that a clear-cut clustering into pro- and anti-apoptotic products, reproducible with different evolutionary methods, could also be obtained by analyzing restricted areas such as the BH1 and BH2 domains. It is noteworthy that even when the analysis is performed only on the BH3 domain, we have two clear-cut clusters. The evolutionary analysis of gene family members is a valuable tool to predict their functions and guide experimental assays to validate predictions. Once the functions of all the components are known, it will be possible to study the process in a holistic way.