Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

Molecular cophylogenetic relationships between European bats and their ectoparasitic mites (Acari, Spinturnicidae).

Cospeciation between host-parasite species is generally thought to result in mirror-image congruent phylogenies. Incongruence can be explained by mechanisms such as host switching, duplication, failure to speciate and sorting events. To investigate the level of association in the host-parasite relationship between Spinturnicid mites and their bat hosts, we constructed the phylogenetic tree of the genus Spinturnix (Acari, Mesostigmata) and compared it to the host phylogeny. We sequenced 938bp of the mitochondrial 16S rDNA and Cytochrome Oxydase subunit I (COI) genes among eleven morphospecies of Spinturnix collected on 20 European Vespertilionid and Rhinolophid bat species. Phylogenetic reconstruction of hosts and parasites showed statistical evidence for cospeciation and suggested that their evolutionary history involved also failure to speciate events and host switches. The latter seem to be mainly promoted by similar roosting habits of the host. As currently understood, host associations of Spinturnicid mites likely results from a complex interaction between the phylogenetic history of the host and the behaviour and the ecology of both parasite and host.

On homology searches by protein Blast and the characterization of the age of genes

Background: It has been shown in a variety of organisms, including mammals, that genes that appeared recently in evolution, for example orphan genes, evolve faster than older genes. Low functional constraints at the time of origin of novel genes may explain these results. However, this observation has been recently attributed to an artifact caused by the inability of Blast to detect the fastest genes in different eukaryotic genomes. Distinguishing between these two possible explanations would be of great importance for any studies dealing with the taxon distribution of proteins and the origin of novel genes. Results: Here we used simulations of protein sequences to examine the capacity of Blast to detect proteins of diverse evolutionary rates in the different species of an eukaryotic phylogenetic tree that included metazoans, fungi and plants. We simulated the evolution of protein genes with the same evolutionary rates than those observed in functional mammalian genes and with among-site rate heterogeneity. Under these conditions, we found that only a very small percentage of simulated ancestral eukaryotic proteins was affected by the Blast artifact. We show that the good detectability of Blast is due to the heterogeneity of protein evolutionary rates at different sites, since only a small conserved motif in a sequence suffices to detect its homologues. Our results indicate that Blast, at least when applied within eukaryotes, only misses homologues of extremely fast-evolving sequences, which are rare in the mammalian genome, as well as sequences evolving homogeneously or pseudogenes.Conclusion: Although great care should be exercised in the recognition of remote homologues, most functional mammalian genes can be detected in eukaryotic genomes by Blast. That is, the majority of functional mammalian genes are not as fast as for not being detected in other metazoans, fungi or plants, if they had been present in these organisms. Thus, the correlation previously found between age and rate seems not to be due to a pure Blast artifact, at least for mammals. This may have important implications to understand the mechanisms by which novel genes originate.

Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal

Background: Myotragus balearicus was an endemic bovid from the Balearic Islands (Western Mediterranean) that became extinct around 6,000-4,000 years ago. The Myotragus evolutionary lineage became isolated in the islands most probably at the end of the Messinian crisis, when the desiccation of the Mediterranean ended, in a geological date established at 5.35 Mya. Thus, the sequences of Myotragus could be very valuable for calibrating the mammalian mitochondrial DNA clock and, in particular, the tree of the Caprinae subfamily, to which Myotragus belongs. Results: We have retrieved the complete mitochondrial cytochrome b gene (1,143 base pairs), plus fragments of the mitochondrial 12S gene and the nuclear 28S rDNA multi-copy gene from a well preserved Myotragus subfossil bone. The best resolved phylogenetic trees, obtained with the cytochrome b gene, placed Myotragus in a position basal to the Ovis group. Using the calibration provided by the isolation of Balearic Islands, we calculated that the initial radiation of caprines can be dated at 6.2 ± 0.4 Mya. In addition, alpine and southern chamois, considered until recently the same species, split around 1.6 ± 0.3 Mya, indicating that the two chamois species have been separated much longer than previously thought. Conclusion: Since there are almost no extant endemic mammals in Mediterranean islands, the sequence of the extinct Balearic endemic Myotragus has been crucial for allowing us to use the Messinian crisis calibration point for dating the caprines phylogenetic tree.

PhylomeDB v4: zooming into the plurality of evolutionary histories of a genome

Phylogenetic trees representing the evolutionary relationships of homologous genes are the entry point for many evolutionary analyses. For instance, the use of a phylogenetic tree can aid in the inference of orthology and paralogy relationships, and in the detection of relevant evolutionary events such as gene family expansions and contractions, horizontal gene transfer, recombination or incomplete lineage sorting. Similarly, given the plurality of evolutionary histories among genes encoded in a given genome, there is a need for the combined analysis of genome-wide collections of phylogenetic trees (phylomes). Here, we introduce a new release of PhylomeDB (http://phylomedb.org), a public repository of phylomes. Currently, PhylomeDB hosts 120 public phylomes, comprising >1.5 million maximum likelihood trees and multiple sequence alignments. In the current release, phylogenetic trees are annotated with taxonomic, protein-domain arrangement, functional and evolutionary information. PhylomeDB is also a major source for phylogeny-based predictions of orthology and paralogy, covering >10 million proteins across 1059 sequenced species. Here we describe newly implemented PhylomeDB features, and discuss a benchmark of the orthology predictions provided by the database, the impact of proteome updates and the use of the phylome approach in the analysis of newly sequenced genomes and transcriptomes.

Nuclear DNA content variation in life history phases of the Bonnemaisoniaceae (Rhodophyta)

Nuclear DNA content in gametophytes and sporophytes or the prostrate phases of the following species of Bonnemaisoniaceae (Asparagopsis armata, Asparagopsis taxiformis, Bonnemaisonia asparagoides, Bonnemaisonia clavata and Bonnemaisonia hamifera) were estimated by image analysis and static microspectrophotometry using the DNA-localizing fluorochrome DAPI (4′, 6-diamidino-2-phenylindole, dilactate) and the chicken erythrocytes standard. These estimates expand on the Kew database of DNA nuclear content. DNA content values for 1C nuclei in the gametophytes (spermatia and vegetative cells) range from 0.5 pg to 0.8 pg, and for 2C nuclei in the sporophytes or the prostrate phases range from 1.15-1.7 pg. Although only the 2C and 4C values were observed in the sporophyte or the prostrate phase, in the vegetative cells of the gametophyte the values oscillated from 1C to 4C, showing the possible start of endopolyploidy. The results confirm the alternation of nuclear phases in these Bonnemaisoniaceae species, in those that have tetrasporogenesis, as well as those that have somatic meiosis. The availability of a consensus phylogenetic tree for Bonnemaisoniaceae has opened the way to determine evolutionary trends in DNA contents. Both the estimated genome sizes and the published chromosome numbers for Bonnemaisoniaceae suggest a narrow range of values consistent with the conservation of an ancestral genome.

Divergent evolution and purifying selection of the flaA gene sequences in Aeromonas

BACKGROUND: The bacterial flagellum is the most important organelle of motility in bacteria and plays a key role in many bacterial lifestyles, including virulence. The flagellum also provides a paradigm of how hierarchical gene regulation, intricate protein-protein interactions and controlled protein secretion can result in the assembly of a complex multi-protein structure tightly orchestrated in time and space. As if to stress its importance, plants and animals produce receptors specifically dedicated to the recognition of flagella. Aside from motility, the flagellum also moonlights as an adhesion and has been adapted by humans as a tool for peptide display. Flagellar sequence variation constitutes a marker with widespread potential uses for studies of population genetics and phylogeny of bacterial species. RESULTS: We sequenced the complete flagellin gene (flaA) in 18 different species and subspecies of Aeromonas. Sequences ranged in size from 870 (A. allosaccharophila) to 921 nucleotides (A. popoffii). The multiple alignment displayed 924 sites, 66 of which presented alignment gaps. The phylogenetic tree revealed the existence of two groups of species exhibiting different FlaA flagellins (FlaA1 and FlaA2). Maximum likelihood models of codon substitution were used to analyze flaA sequences. Likelihood ratio tests suggested a low variation in selective pressure among lineages, with an omega ratio of less than 1 indicating the presence of purifying selection in almost all cases. Only one site under potential diversifying selection was identified (isoleucine in position 179). However, 17 amino acid positions were inferred as sites that are likely to be under positive selection using the branch-site model. Ancestral reconstruction revealed that these 17 amino acids were among the amino acid changes detected in the ancestral sequence. CONCLUSION: The models applied to our set of sequences allowed us to determine the possible evolutionary pathway followed by the flaA gene in Aeromonas, suggesting that this gene have probably been evolving independently in the two groups of Aeromonas species since the divergence of a distant common ancestor after one or several episodes of positive selection. REVIEWERS: This article was reviewed by Alexey Kondrashov, John Logsdon and Olivier Tenaillon (nominated by Laurence D Hurst).

Decoupled evolution of floral traits and climatic preferences in a clade of Neotropical Gesneriaceae.

BACKGROUND: Major factors influencing the phenotypic diversity of a lineage can be recognized by characterizing the extent and mode of trait evolution between related species. Here, we compared the evolutionary dynamics of traits associated with floral morphology and climatic preferences in a clade composed of the genera Codonanthopsis, Codonanthe and Nematanthus (Gesneriaceae). To test the mode and specific components that lead to phenotypic diversity in this group, we performed a Bayesian phylogenetic analysis of combined nuclear and plastid DNA sequences and modeled the evolution of quantitative traits related to flower shape and size and to climatic preferences. We propose an alternative approach to display graphically the complex dynamics of trait evolution along a phylogenetic tree using a wide range of evolutionary scenarios. RESULTS: Our results demonstrated heterogeneous trait evolution. Floral shapes displaced into separate regimes selected by the different pollinator types (hummingbirds versus insects), while floral size underwent a clade-specific evolution. Rates of evolution were higher for the clade that is hummingbird pollinated and experienced flower resupination, compared with species pollinated by bees, suggesting a relevant role of plant-pollinator interactions in lowland rainforest. The evolution of temperature preferences is best explained by a model with distinct selective regimes between the Brazilian Atlantic Forest and the other biomes, whereas differentiation along the precipitation axis was characterized by higher rates, compared with temperature, and no regime or clade-specific patterns. CONCLUSIONS: Our study shows different selective regimes and clade-specific patterns in the evolution of morphological and climatic components during the diversification of Neotropical species. Our new graphical visualization tool allows the representation of trait trajectories under parameter-rich models, thus contributing to a better understanding of complex evolutionary dynamics.

Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal

Background: Myotragus balearicus was an endemic bovid from the Balearic Islands (Western Mediterranean) that became extinct around 6,000-4,000 years ago. The Myotragus evolutionary lineage became isolated in the islands most probably at the end of the Messinian crisis, when the desiccation of the Mediterranean ended, in a geological date established at 5.35 Mya. Thus, the sequences of Myotragus could be very valuable for calibrating the mammalian mitochondrial DNA clock and, in particular, the tree of the Caprinae subfamily, to which Myotragus belongs. Results: We have retrieved the complete mitochondrial cytochrome b gene (1,143 base pairs), plus fragments of the mitochondrial 12S gene and the nuclear 28S rDNA multi-copy gene from a well preserved Myotragus subfossil bone. The best resolved phylogenetic trees, obtained with the cytochrome b gene, placed Myotragus in a position basal to the Ovis group. Using the calibration provided by the isolation of Balearic Islands, we calculated that the initial radiation of caprines can be dated at 6.2 ± 0.4 Mya. In addition, alpine and southern chamois, considered until recently the same species, split around 1.6 ± 0.3 Mya, indicating that the two chamois species have been separated much longer than previously thought. Conclusion: Since there are almost no extant endemic mammals in Mediterranean islands, the sequence of the extinct Balearic endemic Myotragus has been crucial for allowing us to use the Messinian crisis calibration point for dating the caprines phylogenetic tree.

HIV-1 Transmission During Recent Infection and During Treatment Interruptions as Major Drivers of New Infections in the Swiss HIV Cohort Study.

BACKGROUND: Reducing the fraction of transmissions during recent human immunodeficiency virus (HIV) infection is essential for the population-level success of "treatment as prevention". METHODS: A phylogenetic tree was constructed with 19 604 Swiss sequences and 90 994 non-Swiss background sequences. Swiss transmission pairs were identified using 104 combinations of genetic distance (1%-2.5%) and bootstrap (50%-100%) thresholds, to examine the effect of those criteria. Monophyletic pairs were classified as recent or chronic transmission based on the time interval between estimated seroconversion dates. Logistic regression with adjustment for clinical and demographic characteristics was used to identify risk factors associated with transmission during recent or chronic infection. FINDINGS: Seroconversion dates were estimated for 4079 patients on the phylogeny, and comprised between 71 (distance, 1%; bootstrap, 100%) to 378 transmission pairs (distance, 2.5%; bootstrap, 50%). We found that 43.7% (range, 41%-56%) of the transmissions occurred during the first year of infection. Stricter phylogenetic definition of transmission pairs was associated with higher recent-phase transmission fraction. Chronic-phase viral load area under the curve (adjusted odds ratio, 3; 95% confidence interval, 1.64-5.48) and time to antiretroviral therapy (ART) start (adjusted odds ratio 1.4/y; 1.11-1.77) were associated with chronic-phase transmission as opposed to recent transmission. Importantly, at least 14% of the chronic-phase transmission events occurred after the transmitter had interrupted ART. CONCLUSIONS: We demonstrate a high fraction of transmission during recent HIV infection but also chronic transmissions after interruption of ART in Switzerland. Both represent key issues for treatment as prevention and underline the importance of early diagnosis and of early and continuous treatment.

Diarrhea outbreaks in suckling piglets due to rotavirus group C single and mixed (rotavirus groups A and B) infections

Porcine group A rotavirus (PoRVA) is a major cause of neonatal diarrhea in suckling and recently weaned piglets worldwide. The involvement of non-group A rotavirus in cases of neonatal diarrhea in piglets are sporadic. In Brazil there are no reports of the porcine rotavirus group C (PoRVC) as etiologic agent of the diarrhea outbreaks in piglets. The aim of this study was to describe the identification of rotavirus group C in single and in mixed infection with rotavirus groups A and B in three neonatal diarrhea outbreaks in suckling (<21-day-old) piglets, with 70% to 80% and 20% to 25% of morbidity and lethality rates, respectively, in three pig herds located in the state of Santa Catarina, Brazil. The diagnosis of PoRV in the diarrheic fecal samples was performed using polyacrylamide gel electrophoresis (PAGE) to identify the presence of porcine rotavirus groups A, B (PoRVB), and C, and by RT-PCR (PoRVA and PoRVC) and semi-nested (SN)-PCR (PoRVB) to partially amplify the VP4 (VP8*)-VP7, NSP2, and VP6 genes of PoRVA, PoRVB, and PoRVC, respectively. One RT-PCR (PoRVA and PoRVC) and SN-PCR (PoRVB) product of each group of rotavirus of each diarrhea outbreak was submitted to nucleotide (nt) sequence analysis. Based on the PAGE technique, 4 (25%) and 1 (6.25%) of the 16 diarrheic fecal samples evaluated in the first outbreak presented PoRVA and PoRVC electropherotype, respectively, and 11 (68.75%) were negative. In the second outbreak, 3 (42.85%) of the 7 fecal samples evaluated presented PoRVA electropherotype, and in 3 (42.85%) and in 1 (14.3%) fecal samples were detected inconclusive and negative results, respectively. Three (30%) of the 10 fecal samples of the third outbreak presented PoRVC electropherotype; 5 (50%) and 2 (20%) samples showed negative and inconclusive results, respectively. Based on the RT-PCR and SN-PCR assays in the first neonatal diarrhea outbreak, PoRVC was detected in 13 (81.2%) of the 16 diarrheic fecal samples evaluated. PoRVC single infection was identified in 4 (25%) of these samples and mixed infections with PoRVA and PoRVB in 9 (56.2%) fecal samples. All of the seven diarrheic fecal samples evaluated from the second neonatal diarrhea outbreak were positive for PoRVC, whereas its mixed infection with other PoRV groups was detected in 4 (57.2%) samples. In the third outbreak, PoRVC in single infection was detected in all of the 10 diarrheic fecal samples analyzed. In the nt sequence analysis, the PoRVA strains of the first and second outbreaks demonstrated higher nt identity with G4P[6] and G9P[23] genotypes, respectively. The PoRVB strains (first and second outbreaks) and the PoRVC strains (first, second, and third outbreaks) showed higher nt identity and clustered in the phylogenetic tree with PoRVB and PoRVC strains that belong to the N4 and I1 genotypes, respectively. This is the first description in Brazil of the involvement of PoRVC in the etiology of diarrhea outbreaks in suckling piglets. The results of this study demonstrated that PoRVC, in both single and mixed infections, is an important enteropathogen involved in neonatal diarrhea outbreaks in piglets and that the use of more sensitive diagnostic techniques allows the identification of mixed infections involving two or even three groups of PoRV, which may be more common than previously reported.

Monitoring and molecular characterization of group Drotavirus in Brazilian poultry farms

Abstract: Rotaviruses are etiological agents of diarrhea both in humans and in several animal species. Data on avian Group D rotaviruses (RVD) are scarce, especially in Brazil. We detected RVD in 4 pools of intestinal contents of broilers, layer and broiler breeders out of a total of 111 pools from 8 Brazilian states, representing an occurrence of 3.6%, by a specific RVD RT-PCR targeting the VP6 gene. Phylogenetic tree confirmed that the Brazilian strains belong to group D and 3 of the sequences were identical in terms of amino acid whereas one showed 99.5% identity with the others. The sequences described in this study are similar to other sequences previously detected in Brazil, confirming the conserved nature of the VP6 protein.

Unraveling the functional divergence of membrane-bound pyrophosphatases

Inorganic pyrophosphatases (PPases) are enzymes that hydrolyze pyrophosphate (PPi)which is produced as a byproduct in many important growth related processes e.g. in the biosynthesis of DNA, proteins and lipids. PPases can be either soluble or membranebound. Membrane-bound PPases (mPPases) are ion transporters that couple the energy released during PPi hydrolysis to Na+ or H+ transport. When I started the project, only three Na+-transporting mPPases were known to exist. In this study, I aimed to confirm if Na+-transport is a common function of mPPases. Furthermore, the amino acid residues responsible for determining the transporter specificity were unknown. I constructed a phylogenetic tree for mPPases and selected the representative bacterial and archaeal mPPases to be investigated. I expressed different prokaryotic mPPases in Escherichia coli, isolated these as inverted membrane vesicles and characterized their functions. In the first project I identified four new Na+-PPases, two K+-dependent H+-PPases and one K+-independent mPPase. The residues determining the transporter specificity were identified by site-directed mutagenesis. I showed that the conserved glutamate residues are important for specificity, though are not the only residues that influence it. This research clarified the ion transport specificities throughout the mPPase phylogenetic tree, and revealed that Na+ transport is a widespread function of mPPases. In addition, it became clear that the transporter specificity can be predicted from the amino acid sequence in combination with a phylogenetic analysis. In the second project, I identified a novel class of mPPases, which is capable of transporting both Na+ and H+ ions and is mainly found in bacteria of the human gastrointestinal tract. The physiological role of these novel enzymes may be to help the bacteria survive in the demanding conditions of the host. In the third project, I characterized the Chlorobium limicola Na+-PPase and found that this and related mPPases are able to transport H+ ions at subphysiological Na+ concentrations. In addition, the H+-transport activity was shown to be a common function of all studied Na+-PPases at low Na+ concentrations. I observed that mutating gate-lysine to asparagine eliminated the H+ but not the Na+ ion transport function, indicating the important role of the residue in the transport of H+. In the fourth project, I characterized the unknown and evolutionary divergent mPPase clade of the phylogenetic tree. The enzymes belonging to this clade are able to transport H+ ions and, based on their sequence, were expected to be K+- and Na+-independent. The sequences of membrane-bound PPase are usually highly conserved, but the enzymes belonging to this clade are more divergent and usually contain 100−150 extra amino acid residues compared to other known mPPases. Despite the vast sequence differences, these mPPases have the full set of important residues and, surprisingly, are regulated by Na+ and K+ ions. These enzymes are mainly of bacterial origin.

Genotyping hepatitis C virus from hemodialysis patients in Central Brazil by line probe assay and sequence analysis

The present study examined the distribution of hepatitis C virus (HCV) genotypes and subtypes in a hemodialysis population in Goiás State, Central Brazil, and evaluated the efficiency of two genotyping methods: line probe assay (LiPA) based on the 5' noncoding region and nucleotide sequencing of the nonstructural 5B (NS5B) region of the genome. A total of 1095 sera were tested for HCV RNA by RT-nested PCR of the 5' noncoding region. The LiPA assay was able to genotype all 131 HCV RNA-positive samples. Genotypes 1 (92.4%) and 3 (7.6%) were found. Subtype 1a (65.7%) was the most prevalent, followed by subtypes 1b (26.7%) and 3a (7.6%). Direct nucleotide sequencing of 340 bp from the NS5B region was performed in 106 samples. The phylogenetic tree showed that 98 sequences (92.4%) were classified as genotype 1, subtypes 1a (72.6%) and 1b (19.8%), and 8 sequences (7.6%) as subtype 3a. The two genotyping methods gave concordant results within HCV genotypes and subtypes in 100 and 96.2% of cases, respectively. Only four samples presented discrepant results, with LiPA not distinguishing subtypes 1a and 1b. Therefore, HCV genotype 1 (subtype 1a) is predominant in hemodialysis patients in Central Brazil. By using sequence analysis of the NS5B region as a reference standard method for HCV genotyping, we found that LiPA was efficient at the genotype level, although some discrepant results were observed at the subtype level (sensitivity of 96.1% for subtype 1a and 95.2% for subtype 1b). Thus, analysis of the NS5B region permitted better discrimination between HCV subtypes, as required in epidemiological investigations.

VAP-1 as drug target in inflammation : structural features determining ligand interactions

Vascular adhesion protein-1 (VAP-1), which belongs to the copper amine oxidases (CAOs), is a validated drug target in inflammatory diseases. Inhibition of VAP-1 blocks the leukocyte trafficking to sites of inflammation and alleviates inflammatory reactions. In this study, a novel set of potent pyridazinone inhibitors is presented together with their X-ray structure complexes with VAP-1. The crystal structure of serum VAP-1 (sVAP-1) revealed an imidazole binding site in the active site channel and, analogously, the pyridazinone inhibitors were designed to bind into the channel. This is the first time human VAP-1 has been crystallized with a reversible inhibitor and the structures reveal detailed information of the binding mode on the atomic level. Similarly to some earlier studied inhibitors of human VAP-1, the designed pyridazinone inhibitors bind rodent VAP-1 with a lower affinity than human VAP-1. Therefore, we made homology models of rodent VAP-1 and compared human and rodent enzymes to determine differences that might affect the inhibitor binding. The comparison of the crystal structures of the human VAP-1 and the mouse VAP-1 homology model revealed key differences important for the species specific binding properties. In general, the channel in mouse VAP-1 is more narrow and polar than the channel in human VAP-1, which is wider and more hydrophobic. The differences are located in the channel leading to the active site, as well as, in the entrance to the active site channel. The information obtained from these studies is of great importance for the development and design of drugs blocking the activity of human VAP-1, as rodents are often used for in vivo testing of candidate drugs. In order to gain more insight into the selective binding properties of the different CAOs in one species a comprehensive evolutionary study of mammalian CAOs was performed. We found that CAOs can be classified into sub-families according to the residues X1 and X2 of the Thr/Ser-X1-X2-Asn-Tyr-Asp active site motif. In the phylogenetic tree, CAOs group into diamine oxidase, retina specific amine oxidase and VAP-1/serum amine oxidase clades based on the residue in the position X2. We also found that VAP-1 and SAO can be further differentiated based on the residue in the position X1. This is the first large-scale comparison of CAO sequences, which explains some of the reasons for the unique substrate specificities within the CAO family.