- Research article
- Open Access
Comparative evolutionary genomics of the HADH2 gene encoding Aβ-binding alcohol dehydrogenase/17β-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10)
© Marques et al; licensee BioMed Central Ltd. 2006
- Received: 02 April 2006
- Accepted: 09 August 2006
- Published: 09 August 2006
The Aβ-binding alcohol dehydrogenase/17β-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10) is an enzyme involved in pivotal metabolic processes and in the mitochondrial dysfunction seen in the Alzheimer's disease. Here we use comparative genomic analyses to study the evolution of the HADH2 gene encoding ABAD/HSD10 across several eukaryotic species.
Both vertebrate and nematode HADH2 genes showed a six-exon/five-intron organization while those of the insects had a reduced and varied number of exons (two to three). Eutherian mammal HADH2 genes revealed some highly conserved noncoding regions, which may indicate the presence of functional elements, namely in the upstream region about 1 kb of the transcription start site and in the first part of intron 1. These regions were also conserved between Tetraodon and Fugu fishes. We identified a conserved alternative splicing event between human and dog, which have a nine amino acid deletion, causing the removal of the strand βF. This strand is one of the seven strands that compose the core β-sheet of the Rossman fold dinucleotide-binding motif characteristic of the short chain dehydrogenase/reductase (SDR) family members. However, the fact that the substrate binding cleft residues are retained and the existence of a shared variant between human and dog suggest that it might be functional. Molecular adaptation analyses across eutherian mammal orthologues revealed the existence of sites under positive selection, some of which being localized in the substrate-binding cleft and in the insertion 1 region on loop D (an important region for the Aβ-binding to the enzyme). Interestingly, a higher than expected number of nonsynonymous substitutions were observed between human/chimpanzee and orangutan, with six out of the seven amino acid replacements being under molecular adaptation (including three in loop D and one in the substrate binding loop).
Our study revealed that HADH2 genes maintained a reasonable conserved organization across a large evolutionary distance. The conserved noncoding regions identified among mammals and between pufferfishes, the evidence of an alternative splicing variant conserved between human and dog, and the detection of positive selection across eutherian mammals, may be of importance for further research on ABAD/HSD10 function and its implication in the Alzheimer's disease.
- Eutherian Mammal
- Alternative Splice Event
- Alternative Splice Variant
- Substrate Binding Cleft
The enzyme Aβ-binding alcohol dehydrogenase/17β-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10) belongs to short chain dehydrogenase/reductase (SDR) family and its distinct features include the capacity to bind amyloid-beta peptide (Aβ)  and the ability to use a broad array of substrates, encompassing 3-hydroxyacyl-CoA derivatives, steroids, alcohols, and β-hydroxybutyrate [2–6]. ABAD/HSD10 has function in the mitochondria and is expressed in several tissues, including brain, liver and gonads . The broad expression pattern together with the multiple substrate specificities enables the enzyme to participate in several metabolic processes (reviewed in ), namely, the oxidation of fatty acids and branched-chain amino acids [6, 8], sex steroids metabolism in gonads [9, 10] and oxidation of steroid modulators of GABAA receptors in brain .
ABAD/HSD10 was found to bind Aβ in a yeast-two hybrid screen against human brain and HeLa cDNA libraries . Subsequently, various studies [12, 13] provided evidence that ABAD/HSD10 can mediate the cytotoxic effects of Aβ in the mitochondrial compartment, thus contributing to the mitochondrial dysfunction seen in Alzheimer's disease (AD) (reviewed in ). The structure determination of human ABAD/HSD10 in complex with Aβ  revealed that the enzyme's loop D is the binding site of Aβ, and that binding of Aβ to ABAD/HSD10 leads to distortion of the ABAD/HSD10 structure, possibly inhibiting its enzymatic activity. Transgenic mice overexpressing ABAD/HSD10 in an Aβ-rich environment displayed neuronal oxidative stress, cell death and accelerated decline in spatial learning and memory [12, 13]. In addition, ABAD/HSD10 expression was reported to be enhanced in brains from patients with AD when compared with brains from non-demented age-matched controls [12, 13].
Although ABAD/HSD10 is an enzyme with function in the mitochondria it is encoded by a nuclear gene termed HADH2 (official name). The human HADH2 gene, mapped at chromosome Xp11.2, consists of six exons and five introns  and encodes a protein with 261 amino acids. The three-dimensional structure of ABAD/HSD10 has been determined with great resolution [2, 15]. The enzyme has a homotetrameric structure, with each subunit containing a Rossman fold dinucleotide-binding motif, composed of a core β-sheet of seven parallel strands flanked by six α-helices, which is involved in the interaction with the cofactor (NAD). The C-terminal portion of the enzyme is involved in substrate binding and harbors the Ser/Lys/Tyr catalytic triad characteristic of the SDR family members .
Here we performed a thorough comparative genomic analysis of HADH2 genes from 21 organisms, including several species of mammals, amphibians, fishes, insects and nematodes. We provide insights into the evolution of such HADH2 genes, namely its genomic organization, patterns of sequence conservation, alternative splicing variants and phylogenetic relationships. Moreover, we provide evidence suggesting signatures of positive selection across eutherian mammal ABAD/HSD10 proteins, which may be of importance for more applied biomedical research on enzyme function and its implication in the Alzheimer's disease.
Cross-species comparison of HADH2 gene organization
Accsession numbers, chromosome location and GC content of HADH2 genes
Gene length (Kb)
human (Homo sapiens)
chromosome X (p11.2)
chimpanzee (Pan troglodytes)
orangutan (Pongo Pygmaeus)
rhesus monkey (Macaca mulatta)
rat (Rattus norvegicus)
mouse (Mus musculus)
dog (Canis familiaris)b
cat (Felis catus)c
cow (Bos Taurus)
pig (Sus scrofa)d
opossum (Monodelphis domestica)e
western clawed frog (Xenopus tropicalis)f
african clawed frog (Xenopus laevis)d
zebrafish (Danio rerio)
Fugu (Takifugu rubripes)c
Tetraodon (Tetraodon nigroviridis)
fruitfly (Drosophila melanogaster)
mosquito (Anopheles gambiae)c
honeybee (Apis mellifera)
Caenorhabditis briggsae c
The gene overall GC content was similar among the eutherian mammals (48–55%), but substantially heterogeneous among fishes (33–57%), insects (18–60%) and nematodes (34–43%). The GC content of third-codon position (GC3) showed a high heterogeneity across HADH2 genes, while GC content in the first (GC1) and second-codon (GC2) positions were more homogeneous. Previously, it was reported that in most vertebrate genomes, GC3 levels are correlated with the GC level of the isochore region containing the gene [17, 18]. Accordingly, the vertebrate HADH2 genes with high GC3 values should probably be localized in GC-rich isochores, which also preferentially allow the accumulation of SINE repeats. Conversely, the comparatively low GC3 values observed in rodents, amphibians and zebrafish HADH2 genes should probably reflect a shift toward GC-poorer regions. Interestingly, honeybee HADH2 gene was remarkably AT-biased, which was related both to the very low GC3 content and the extreme predominance of AT-rich repeats in its introns (Additional file 1).
Eutherian mammals HADH2 genes were localized in the X chromosome (Table 1), confirming the known high synteny conservation of eutherian X chromosomes . Tetraodon HADH2 gene was localized on chromosome 9, which was previously reported to show a high synteny with human Xp11.2 region . Zebrafish HADH2 gene was, however, localized on chromosome 23, which has a very low synteny either with human chromosome X or with Tetraodon chromosome 9 .
Patterns of nucleotide sequence conservation among HADH2 genes
The comparisons among fish (Figure 5B) revealed that conservation between pufferfishes and zebrafish sequences was largely limited to coding regions and to a microRNA (family let-7) positioned upstream of the HADH2 gene, reflecting the deep divergence of the Euteleostei (110–160 million years; Myr) . Despite the much shorter evolutionary distance separating Fugu and Tetraodon (18–30 Myr) , the intron identity between the two pufferfishes (67.7%) was lower than the intron identity between human and each of the other non-primate eutherian mammals (> 70 Myr divergence) . This might be related with the reported high neutral nucleotide substitution rate between Fugu and Tetraodon, which was shown to be greater than that between human and mouse . C. elegans and C. briggsae HADH2 genes were substantially conserved in coding regions but not in the introns (Figure 5C), which is supported by previous inferences of the reduced conservation of noncoding sequences between the genomes of these nematodes, due to their large divergence (~100 Myr) . Insect HADH2 genes were conserved only in HADH2 coding regions. The divergent organization of insect HADH2 genes and the extremely high content of AT-rich repeats in the honeybee introns ruled out any conservation among insect introns or upstream regions.
High conservation of noncoding sequences across multiple species may indicate the presence of functional elements. The Multipipmaker analysis (Figure 5) showed the existence of gap-free alignments in noncoding regions with more than 70% identity, conserved across all human/eutherian mammals and to a lesser extent between Tetraodon and Fugu. Within the HADH2 upstream region, the portion at about 1 kb of the transcription start site showed an increase in conservation, with the alignment from the positions 897 bp to 1239 bp (human sequence) exhibiting > 73% identity for human/eutherian mammals. The upstream regions of Tetraodon and Fugu also displayed highly conserved alignments, in particular one, positioned very close to the miRNA, with 88% identity along 126 bp. Considering intronic regions, the first part of intron 1 contained the largest (130 bp) gap-free alignment conserved (72–73% with the human sequence) between human, carnivores and cow, which also overlapped with some shorter human/rodent gap-free alignments with ≥70% identity. Interestingly, the largest conserved gap-free alignment (72% identity along 89 bp) between Tetraodon and Fugu introns was also in the beginning of pufferfishes intron 1. Among eutherian mammals, regions of high conservation were also detected in the other introns. In the smaller intron 4, the last 50 bp alignment showed an identity above 75% for all human/eutherian mammal comparisons. The intron 5, which is involved in an HADH2 alternative splicing event similar in human and dog, showed an identity above 70% in the last 69 bp for all human/eutherian mammal comparisons. The 3'UTR region, which is frequently highly conserved among orthologue genes due to its involvement in pos-transcriptional control, displayed an identity above 70% in the human/eutherian mammal comparisons.
Patterns of ABAD/HSD10 amino acid sequence conservation
Adaptative selection on eutherian mammal ABAD/HSD10 proteins
Results of the gene level approach (PAML) applied in eutherian mammal HADH2 genes
p0 = 0.824, (p1 = 0.176)
ω0 = 0.035 (ω1 = 1)
M1a vs M2a
p0 = 0.824, p1 = 0.114, (ps = 0.062)
ω0 = 0.035 (ω1 = 1) ws = 1
p = 0.091, q = 0.371
p0 = 0.919, (p1 = 0.106), q = 0.728
M7 vs M8
ps = 0.08098, ωs = 1
One ratio (M0)
ω = 0.176
M0 vs free-ratio
ω0 = 1.302 (foreg.)
M0 vs two-ratio
ω1 = 0.167 (back.)
Results of the protein-level approach (TreeSAAP) showing the amino acid sites and physicochemical amino acid properties influenced by positive-destabilizing selection among eutherian mammal ABAD/HSD10 proteins.
Amino acid site
Amino acid site
Substrate binding-cleft sites
Pα, α m , pH i
N s , B r, μ
B l , R F , α n , R a , H p , pK', H, N s , B r , μ
K 0 , M w , V 0 , M w , H t
Pβ, B l , K 0 , R a
K 0 , C a , M w , V 0
pH i , F, P, α c , P t
Sites in insertion 1 region of loop D
N s , B r , E l , H p , E t
B l , C a , H, Mw, V 0 , F
P c , pK', F, P t
Pα, α m
Pα, E sm
N s , Pβ, B r , H p , E t
α n , R a , H p
C, R a
P r , P, α m
P c , pK', F, P t
N s , Pβ, B l , R F , R a , H p , H t
α n , R a , H p , α c
Pα, E sm
K 0 , M w , V 0
N s , B r , H p , E t , Pβ
E sm , Pα
pK', F, P r , α c
E sm , Pα
C a , M w , V 0
The evidence of positive selection, suggested both by the gene-level and the protein-level approaches in the lineage separating human/chimpanzee from orangutan reflects both the particular higher nonsynonymous than synonymous substitutions in that lineage and also the fact of most of the amino acid changes were nonconservative, thus likely to result in modifications of the physicochemical amino acid properties (Figure 10). Human/chimpanzee and orangutan differed in seven and five (three in case of chimpanzee) nonsynonymous and synonymous substitutions, respectively. By contrast, the other eutherian mammal branches revealed a considerably greater number of synonymous than nonsynonymous substitutions (results not shown). Among primates, human and chimpanzee proteins were identical while a higher number of amino acid differences was observed between human (and chimpanzee) and orangutan (seven) relatively to the orangutan and macaca (three) (Figure 10). This suggests that different rates of protein evolution might have occurred in the evolution of the great ape ABAD/HSD10 proteins. Of the six positively selected sites in the branch separating human/chimpanzee from orangutan, three were localized in loop D (98, 106 and 108) and one in substrate-binding loop (214). It is also interesting to note that in two of the six positive selected sites (106 and 108), the changes occurring in the human/chimpanzee sequences were not seen in the other eutherian ABAD/HSD10 proteins, and in three sites (80, 119 and 214), amino acid changes occurred exclusively in the human/chimpanzee sequences with the other eutherian sequences containing the same amino acid.
Here, we present a comprehensive comparative genomic analysis of the HADH1 gene encoding ABAD/17β-HSD10 across 21 species. HADH2 gene revealed a substantial conserved organization across a large evolutionary distance: vertebrate and nematode HADH2 genes showed a six-exon/five-intron gene organization, but insects showed a reduced and varied number of exons (two to three). In general, HADH2 genes were less than 3.20 Kb long (Table 1; Figure 1), with the exception of zebrafish gene which was much larger (9.38 Kb). At nucleotide level, a notable characteristic of HADH2 genes was the extensive variation in GC3 content.
The reduced conservation in noncoding regions between human and the non-eutherian vertebrate orthologues, reinforce previous conclusions that in general, few noncoding functional sequences remain conserved for large evolutionary distances . Several studies [34–36] have shown that the conservation in non-coding regions between human and distant vertebrate species (e.g., fishes) is restricted to a subset of genes that are involved in pivotal biological processes such as development and transcription regulation. These genes display a high density of conserved elements in their introns and intergenic regions, which is related with the need to preserve complex and crucial regulatory mechanisms in basic vertebrate development [36, 37]. A recent broad study  about intron conserved elements between human and various vertebrates, including chicken and fishes, reported that multispecies conserved noncoding sequences distribution is not uniform across human introns. Indeed, the longer introns of the genes involved in development and transcription regulation showed a tendency to accumulate conserved sequences, while the majority of relatively short introns (< 9 Kb) displayed none or few conserved elements . In view of the limited conservation between human and distantly vertebrates, it has been assumed that for the majority of human genes, comparisons between multiple, moderately related species might represent a better strategy to search for potential regulatory elements , although cautious is required as the high degree of similarity might also reflect low substitutions rates of evolution. It is possible that some of the HADH2 conserved noncoding sequences, identified in the comparisons between human and the five moderately related eutherian mammals might indeed represent regulatory elements, and thus be good candidates for functional experimental studies. For instance, the significant sequence similarity throughout the upstream regions, namely in the region about 1 kb of the transcription start site, may suggest the presence of regulatory elements, likely involved in transcription. The first part of intron 1 was found to be highly conserved among eutherian mammals and interestingly, it was the highest conserved intronic region between Tetraodon and Fugu. Given that intron-associated regulatory elements on genes tend to be localized preferentially in intron 1 , such pattern of conservation in eutherian mammals and in pufferfishes genes suggests that HADH2 intron 1 may potentially contain a regulatory element.
Previous studies [41, 42] reported an increase in both exon and intron conservation in the regions flanking conserved alternative splice sites. The conserved alternative splicing event between human and dog HADH2 genes raises the possibility that some of the conserved intronic regions, particularly in the introns flanking exon 5 (intron 4 and 5), might reflect the presence of splicing regulatory elements subject to purifying selection. In addition, the high conservation found in the exons 5 and 6 likely reflect their importance in coding for residues of the substrate binding cleft, but it might also be related with the presence of conserved exonic splicing regulatory elements. The conservation of human and dog HADH2 alternative donor splice sites in other eutherian mammals, suggests they can potentially also express an identical alternative splicing variant. However, as previously noted [43, 44], conservation of a splice site is not enough to predict the existence of a variant. Thus, the evaluation of identical alternative splice variants in the other mammals, including the one identified solely in dog, will help to further elucidate the potential functional importance of the conservation of intronic elements in eutherian HADH2 genes.
Curiously, the two identified alternative splicing events interfere with the enzyme C-terminal, an important functional region. In the dog variant 3 (169 amino acid), the last 92 amino acids are replaced by seven new ones. As this leads to the loss of two residues of the catalytic triad and of the substrate binding loop residues, likely, the dog variant 3 is non-functional. The production of non-functional transcripts is not uncommon, with many genes using that as a mechanism to control the mRNA expression levels [45, 46]. In human and dog variant 2, the loss of strand βF, which is one of the seven strands that compose the core β-sheet of the ABAD/HSD10 Rossmann fold dinucleotide-binding motif, may have both structural and functional consequences, as this strand is adjacent to the region containing the substrate binding loop and is also involved in the substrate binding . However, the fact that the substrate binding cleft residues are retained and the existence of a shared variant between human and dog suggest that it might be functional. The assessment of the alternative splice variants functionality, as well as the mechanisms regulating splice sites selection, and of the regional expression levels, will be fundamental to determine the implications of the alternative splicing in the physiological and pathophysiological functions of ABAD/HSD10. Of special importance is how ABAD/HSD10 alternative splice variants behave in an Aβ-rich environment? In this respect, since the human variant 2 has a normal N-terminal region, which is responsible for the interaction with Aβ , it may still retain the ability to interact with Aβ. Moreover, the human variant 2 is supported by an mRNA sequence (BC008708) derived from neuroblastoma cell lines, suggesting that it may be expressed in the human brain.
At protein level, ABAD/HSD10 showed a high similarity across the different species, namely in regions comprehending the substrate binding cleft and the subunit association, clearly suggesting that ABAD/HSD10 maintained a substantial structural and functional conservation across a very large evolutionary distance. Indeed, previous studies have demonstrated that the rat and fruitfly orthologues exhibit enzymatic activities similar to those of the human enzyme [2, 47], suggesting that ABAD/HSD10 might have important functions both in vertebrates and invertebrates. However, despite being recognized that the broad substrate specificity of ABAD/HSD10 enables it to participate in several metabolic pathways, the physiological properties of the enzyme are not yet completely understood. In mammals, ABAD/HSD10 was suggested to have an important role in metabolism of sex steroid hormones . It was reported to be expressed in the Leydig cells of testes from various mammals and the differentiation-dependent expression of ABAD/HSD10 in rodent testes suggested that this enzyme might contribute to protecting Leydig cells from the effects of estrogens . The great importance of fruitfly ABAD/HSD10 (termed scully) was demonstrated by the mutational inactivation of the enzyme, which induced a lethal phenotype during embryonic and pupal development, with mutants displaying non-functional gonads, lipid accumulation and aberrant mitochondria . Human ABAD/HSD10 deficiency causes a disorder in which the isoleucine degradation is impaired . Patients with this deficiency show severe neurological abnormalities, including psychomotor retardation and progressive infantile neurodegeneration. A beneficial role in the cellular response to metabolic stress was attributed to the ABAD/HSD10 enzyme, due to its ability to utilize the ketone body β-hydroxybutyrate as a substrate . ABAD/HSD10 may also be important in the stabilization of mitochondrial function  and in the maintenance of normal functions of GABAergic neurons .
Despite the high degree of conservation among eutherian ABAD/HSD10 proteins, suggesting strong purifying selection pressures, we investigated signatures of positive selection by using both a gene-level and protein-level approaches. The failure of the gene-level approach in providing significant evidence of molecular adaptation across eutherian mammal orthologues, probably reflect the known lack of power of the used LRTs in detect positive selection when divergence between sequences in the data set is low . By contrast, the significant evidence of positive selection provided by the protein level approach indicates its ability to identify molecular adaptation even when proteins are highly conserved . An interesting finding, was the detection of positive selection in important functional regions, in particular in the lineage separating human/chimpanzee from orangutan, which, contrasting with other eutherian mammal branches, accumulated a higher number of nonsynonymous than synonymous substitutions. Of the seven residues differing between human and orangutan, the protein-level approach detected six to be under molecular adaptation, four of which localized in particularly important functional regions. Specifically, three sites were localized in the loop D (site 98 in a region belonging to the substrate binding cleft and sites 106 and 108 in insertion 1 region) and a fourth site (214) was localized in the substrate binding loop. The previous sites 98 and 108, plus three additional sites (95, 102, 103) localized in loop D and site 202 belonging to the substrate binding loop were also identified to be under positive selection in other eutherian mammal branches. The potential functional meaning of the positive selection in eutherian mammal ABAD/HSD10 proteins, particularly in the regions belonging to the substrate binding cleft is, however, difficult to ascertain given the enzyme multiple substrate specificities and participation in various metabolic pathways. As stated behind, ABAD/HSD10 seems to have maintained a considerable structural and functional conservation across a very large evolutionary distance. However, it is important to say that although the identification of positively selected amino acid sites do not necessarily prove that such amino acid replacements modify the protein function, their occurrence in functional important sites provide a strong evidence for further functional experimental analyses. In this respect, particularly appealing is the evidence of positive selection in loop D, given its involvement in Aβ binding. Recently, the determination of the crystal structure of ABAD/HSD10 bound to Aβ and mutational studies on loop D furnished strong evidence for this loop functions as the binding site for Aβ .
The sequencing of various genomes provided the opportunity to study the molecular evolution of the HADH2 gene encoding ABAD/HSD10. Our study revealed that HADH2 genes maintained a very similar organization and substantial conservation at amino acid level over more than one billion years. The identification of a conserved alternative splicing event between human and dog and highly conserved noncoding regions among eutherian mammals may provide a framework for further investigation of HADH2 gene regulation. The evidence of positive selection across eutherian mammal ABAD/HSD10 proteins may be of importance for more applied biomedical research on the enzyme function and its implication in the Alzheimer's disease.
Database search of HADH2 gene sequences
The sequences of the HADH2 gene encoding ABAD/HSD10 protein (synonymous names include SCHAD, ERAB, MHBD, and scully for the Drosophila orthologue) were retrieved from Ensembl , NCBI  and TIGR  databases (Table 1). Gene sequences were identified either by TBLASTN searches within the various species genome sequence projects using known ABAD/HSD10 amino acid sequences as queries or manually reconstructed from whole genome shotgun (WGS) traces through MEGABLAST searches. For some species, the HADH2 genes retrieved were only partially sequenced (see figure 1). Additionally, BLAST searches were performed against species specific mRNA reference sequence (RefSeq) databases at NCBI to detect alternative splicing variants.
The amount and composition of repetitive elements was investigated using RepeatMasker , CENSOR  and Tandem Repeats Finder v.3.01 . After removing repetitive motifs, Pip software analysis  was used to align and identify patterns of sequence conservation across vertebrates, fishes and nematodes HADH2 genes. The 2 kb region upstream the HADH2 genes was included in the analysis, excepting for cat and orangutan genes, where only a smaller portion was available (1 kb for cat and 200 bp for orangutan). Protein sequences were aligned with CLUSTALW . Sliding window percent amino acid identity analyses (excluding indels) were conducted using Swaap 1.0.2 .
HADH2 sequences were investigated for variation in base composition (or compositional bias), mutational saturation, and gene conversion, which are events known to disturb phylogenetic reconstructions. The chi-square test of homogeneity implemented in TREE-PUZZLE v5.2  was used to evaluate variation in base composition at each codon position. GENECONV v1.81  was employed, using the default settings, to detect recombination/gene conversion events in the data set. To test for mutational saturation, we plotted the number of transitions and transversions from first, second, and third-codon positions against the pairwise genetic distances. The SYM+G+I model was identified with Modeltest v3.06  as the best evolutionary model fitting the data. Transitions and transversions for all pairwise sequence comparisons were calculated using MEGA v3.1 , whereas genetic distances were calculated in PAUP v.4.0 b10 . In the absence of mutational saturation, genetic distances and nucleotide substitutions give a linear relationship. Conversely, in case of nucleotide saturation, genetic distances are larger than substitutions [66, 67].
The phylogenetic relationships among HADH2 sequences from different species were determined using Maximum-likelihood (ML) and Bayesian methods, implemented in PAUP v.4.0 b10  and MRBAYES v3.1 , respectively. The ML tree was reconstructed through a heuristic search with ten random additions of taxa and tree bisection-reconnection (TBR) branch swapping algorithm. Bootstrap support (BS) values were estimated with 100 replicates. In the Bayesian analysis, four markov chains were run for 500,000 generations with burn-in values of 2,500 generations and trees being sampled every 100 generations. Bayesian posterior probabilities (BPP) were used to evaluate branch support. Both trees were rooted using C. elegans and C. briggsae sequences as outgroups.
Detection of positive selection
Positive selection analyses were restricted to eutherian mammals to avoid violations in the evolutionary assumptions, i.e. absence of nucleotide saturation and base compositional bias, which requires closely related sequences (see Table 1 and Figure 6). We used two strategies to identify positive selection: (i) a gene-level approach based on the ratio (ω) of nonsynonymous (d N ) to synonymous (d S ) substitutions rate (i.e., ω = dN/dS), and (ii) a protein-level approach which evaluates the physicochemical importance of amino acid changes on the protein structure. The unrooted eutherian mammals ML tree was used in the analyses.
The gene-level approach implemented in PAML v3.14  uses likelihood ratio tests (LRT) to compare two nested models, a model that does not account for sites with ω > 1 (null model) and a model that does (positive selection model) . We used two LRTs based on site specific models, which compare the null models M1a and M7 against the alternative models (positive selection models) M2a and M8, respectively. The posterior probability of a site being under positive selection was obtained using the Bayes Empirical Bayes (BEB) method implemented in PAML . We also constructed two LRTs based on branch models, the first compares one ratio model (M0, assumes the same ω ratio for all branches) with the free-ratios model (allows an independent ω ratio for each branch) and the second compares model M0 with the two-ratio model (assumes a ω ratio for foreground branch different from that of background branch). The two-ratio model was applied in a specific branch (see results for details).
The protein-level approach implemented in TreeSAAP  measures the selective influences of 31 physicochemical properties across a phylogenetic tree following McClellan and McCracken method . The program uses a gradient of categories to classify each property change from conservative to radical, and calculates a z-score which indicates the direction of selection. In our analysis, we were interested in detecting positive-destabilizing selection as this results in radical structural or functional shifts in local regions of the protein, thus, being unambiguously correlated with molecular adaptation. An amino acid property is said to be affected by positive-destabilizing selection when the frequency of changes in radical magnitude categories exceeds the frequency(s) expected by chance, as indicated by positive z-scores.
A. Marques was supported by a PhD grant (SFRH/BD/19228/2004) from Fundação para a Ciência e a Tecnologia. Comments made by three anonymous referees improved a previous version of this manuscript.
- Yan SD, Fu J, Soto C, Chen X, Zhu H, Al-Mohanna F, Collision K, Zhu A, Stern E, Saído T, Tohyama M, Ogawa S, Roher A, Stern D: An intracellular protein that binds amyloid-beta peptide and mediates neurotoxicity in Alzheimer's disease. Nature. 1997, 389: 689-695. 10.1038/39522.PubMedView ArticleGoogle Scholar
- Powell AJ, Read JA, Banfield MJ, Gunn-Moore F, S Yan D, Lustbader J, Stern AR, Stern DM, Brady RL: Recognition of structurally diverse substrates by type II 3-hydroxyacyl-CoA dehydrogenase (HADH II)/amyloid-bbinding alcohol dehydrogenase (ABAD). J Mol Biol. 2000, 303: 311-327. 10.1006/jmbi.2000.4139.PubMedView ArticleGoogle Scholar
- Furuta S, Kobaysashi A, Miyazawa S, Hashimoto T: Cloning and expression of cDNA for a newly identified isoenzyme of bovine liver 3-hydroxyacyl-CoA dehydrogenase and its import into mitochondria. Biochim Biophys Acta. 1997, 1350: 317-324.PubMedView ArticleGoogle Scholar
- He XY, Mwerz G, Mehta P, Schulz H, Yang SY: Human brain short chain L-3-hydroxyacyl coenzymeA dehydrogenase is a single domain multifunctional enzyme. J Biol Chem. 1999, 274: 15014-15019. 10.1074/jbc.274.21.15014.PubMedView ArticleGoogle Scholar
- He XY, Merz G, Yang YZ, Mehta P, Schulz H, Yang SY: Characterization and localization of human type10 17b-hydroxysteroid dehydrogenase. Eur J Biochem. 2001, 268: 4899-4907. 10.1046/j.0014-2956.2001.02421.2421.x.PubMedView ArticleGoogle Scholar
- Ofman R, Ruiter JP, Feenstra M, Duran M, Poll-The BT, Zschocke J, Ensenauer R, Lehnert W, Sass JO, Sperl W, Wanders RJ: 2-Methyl-3-hydroxybutyryl-CoA dehydrogenase deficiency is caused by mutations in the HADH2 gene. Am J Hum Genet. 2003, 72: 1300-1307. 10.1086/375116.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang SY, He XY, Schulz H: Multiple functions of type 10 17beta-hydroxysteroid dehydrogenase. Trends Endocrinol Metab. 2005, 16: 167-75. 10.1016/j.tem.2005.03.006.PubMedView ArticleGoogle Scholar
- He XY, Schulz H, Yang SY: A human brain L-3-hydroxyacyl-coenzyme A dehydrogenase is identical to an amyloid beta-peptide-binding protein involved in Alzheimer's disease. J Biol Chem. 1998, 273: 10741-10746. 10.1074/jbc.273.17.10741.PubMedView ArticleGoogle Scholar
- He XY, Yang YZ, Peehl DM, Lauderdale A, Schulz H, Yang SY: Oxidative 3alpha-hydroxysteroid dehydrogenase activity of human type 10 17beta-hydroxysteroid dehydrogenase. J Steroid Biochem Mol Biol. 2003, 87: 191-198. 10.1016/j.jsbmb.2003.07.007.PubMedView ArticleGoogle Scholar
- Ivell R, Balvers M, Anand RJ, Paust HJ, Mckinnell C, Sharpe R: Differentiation-dependent expression of 17beta-hydroxysteroid dehydrogenase, type 10, in the rodent testis: effect of aging in Leydig cells. Endocrinology. 2003, 144: 3130-7. 10.1210/en.2002-0082.PubMedView ArticleGoogle Scholar
- He XY, Wegiel J, Yang SY: Intracellular oxidation of allopregnanolone by human brain type 10 17beta-hydroxysteroid dehydrogenase. Brain Res. 2005, 1040: 29-35. 10.1016/j.brainres.2005.01.022.PubMedView ArticleGoogle Scholar
- Lustbader JW, Cirilli M, Lin C, Xu HW, Takuma K, Wang N, Caspersen C, Chen X, Pollak S, Chaney M, Trinchese F, Liu S, Gunn-Moore F, Lue LF, Walker DG, Kuppusamy P, Zewier ZL, Arancio O, Stern D, Yan SS, Wu H: ABAD directly links Aβ to mitochondrial toxicity in Alzheimer's disease. Science. 2004, 304: 448-452. 10.1126/science.1091230.PubMedView ArticleGoogle Scholar
- Takuma K, Yao J, Xu H, Chen X, Luddy J, Trillat AC, Stern DM, Arancio O, Yan SS: ABAD enhances Abeta-induced cell stress via mitochondrial dysfunction. FASEB J. 2005, 19: 597-8.PubMedGoogle Scholar
- Yan SD, Stern DM: Mitochondrial dysfunction and Alzheimer's disease: role of amyloid-beta peptide alcohol dehydrogenase (ABAD). Int J Exp Pathol. 2005, 86: 161-71. 10.1111/j.0959-9673.2005.00427.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Kissinger CR, Rejto PA, Pelletier LA, Thomson JA, Showalter RE, Abreo MA, Agree CS, Margosiak S, Meng JJ, Aust RM, Vanderpool D, Li B, Tempczyk-Russell A, Villafranca JE: Crystal structure of human ABAD/HSD10 with a bound inhibitor: implications for design of Alzheimer's disease therapeutics. J Mol Biol. 2004, 342: 943-52. 10.1016/j.jmb.2004.07.071.PubMedView ArticleGoogle Scholar
- Jornvall H, Persson B, Krook M, Atrian S, Gonzalez-Duarte R, Jeffery J, Ghosh D: Short-chain dehydrogenases/reductases (SDR). Biochemistry. 1995, 34: 6003-13. 10.1021/bi00018a001.PubMedView ArticleGoogle Scholar
- Clay O, Caccio S, Zoubakn S, Mouchiroud D, Bernardi G: Human coding and non-coding DNA: compositional correlations. Mol Phylogenet Evol. 1996, 5: 2-12. 10.1006/mpev.1996.0002.PubMedView ArticleGoogle Scholar
- D'Onofrio G, Bernardi G: A universal compositional correlation among codon positions. Gene. 1992, 110: 81-88. 10.1016/0378-1119(92)90447-W.PubMedView ArticleGoogle Scholar
- Ohno S: Sex chromosomes and sex-linked genes. 1967, Berlin: Springer-VerlagView ArticleGoogle Scholar
- Kohn M, Kehrer-Sawatzki H, Vogel W, Graves JA, Hameister H: Wide genome comparisons reveal the origins of the human X chromosome. Trends Genet. 2004, 20: 598-603. 10.1016/j.tig.2004.09.008.PubMedView ArticleGoogle Scholar
- Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS: The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res. 2005, 15: 1307-14. 10.1101/gr.4134305.PubMedPubMed CentralView ArticleGoogle Scholar
- Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28: 4364-4375. 10.1093/nar/28.21.4364.PubMedPubMed CentralView ArticleGoogle Scholar
- Mount SM: Genomic Sequence, Splicing, and Gene Annotation. Am J Hum Genet. 2000, 67: 788-792. 10.1086/303098.PubMedPubMed CentralView ArticleGoogle Scholar
- Burset M, Seledtsov IA, Solovyev VV: SpliceDB: database of canonical and non-canonical mammalian splice sites. Nucleic Acids Res. 2001, 29: 255-9. 10.1093/nar/29.1.255.PubMedPubMed CentralView ArticleGoogle Scholar
- Stamm S, Zhang MQ, Marr TG, Helfman DM: A sequence compilation and comparison of exons that are alternatively spliced in neurons. Nucleic Acids Res. 1994, 22: 1515-26.PubMedPubMed CentralView ArticleGoogle Scholar
- Thanaraj TA, Clark F: Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions. Nucleic Acids Res. 2001, 29: 2581-2593. 10.1093/nar/29.12.2581.PubMedPubMed CentralView ArticleGoogle Scholar
- Murphy WJ, Eizirik E, O'Brien SJ, Madsen O, Scally M, Douady CJ, Teeling E, Ryder OA, Stanhope MJ, de Jong WW, Springer MS: Resolution of the early placental mammal radiation using Bayesian phylogenetics. Science. 2001, 294: 2348-2351. 10.1126/science.1067179.PubMedView ArticleGoogle Scholar
- Wu CI, Li WH: Evidence for higher rates of nucleotide substitution in rodents than in man. Proc Natl Acad Sci. 1985, 82: 1741-1745. 10.1073/pnas.82.6.1741.PubMedPubMed CentralView ArticleGoogle Scholar
- Nelson JS: Fishes of the world. 1994, New York, WileyGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.PubMedView ArticleGoogle Scholar
- Springer MS, Murphy WJ, Eizirik E, O'Brien SJ: Placental mammal diversification and the Cretaceous-Tertiary boundary. Proc Natl Acad Sci. 2003, 100: 1056-61. 10.1073/pnas.0334222100.PubMedPubMed CentralView ArticleGoogle Scholar
- Stein LD, Bao Z, Blasiar D, Blumenthal T, Brent MR, Chen N, Chinwalla A, Clarke L, Clee C, Coghlan A: The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics. PLoS Biol. 2003, 1: e45-10.1371/journal.pbio.0000045.PubMedPubMed CentralView ArticleGoogle Scholar
- Thomas JW, Touchman JW, Blakesley RW, Bouffard GG, Beckstrom-Sternberg SM, Margulies EH, Blanchette M, Siepel AC, Thomas PJ, McDowell JC: Comparative analyses of multi-species sequences from targeted genomic regions. Nature. 2003, 424: 788-793. 10.1038/nature01858.PubMedView ArticleGoogle Scholar
- Bagheri-Fam S, Ferraz C, Demaille J, Scherer G, Pfeifer D: Comparative genomics of the SOX9 region in human and Fugu rubripes: conservation of short regulatory sequence elements within large intergenic regions. Genomics. 2001, 78: 73-82. 10.1006/geno.2001.6648.PubMedView ArticleGoogle Scholar
- Goode DK, Snell P, Elgar G: Comparative analysis of vertebrate Shh genes identifies novel conserved non-coding sequence. Mamm Genome. 2003, 14: 192-201. 10.1007/s00335-002-3052-z.PubMedView ArticleGoogle Scholar
- Woolfe A, Goodson M, Goode DK, Snell P, McEwen GK, Vavouri T, Smith SF, North P, Callaway H, Kelly K, Walter K, Abnizova I, Gilks W, Edwards YJK, Cooke JE, Elgar G: Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol. 2005, 3: e7-10.1371/journal.pbio.0030007.PubMedPubMed CentralView ArticleGoogle Scholar
- Ovcharenko I, Loots GG, Nobrega MA, Hardison RC, Miller W, Stubbs L: Evolution and functional classification of vertebrate gene deserts. Genome Res. 2005, 15: 137-145. 10.1101/gr.3015505.PubMedPubMed CentralView ArticleGoogle Scholar
- Sironi M, Menozzi G, Comi GP, Cagliani R, Bresolin N, Pozzoli U: Analysis of intronic conserved elements indicates that functional complexity might represent a major source of negative selection on non-coding sequences. Hum Mol Genet. 2005, 14: 2533-2546. 10.1093/hmg/ddi257.PubMedView ArticleGoogle Scholar
- Boffelli D, Nobrega MA, Rubin EM: Comparative genomics at the vertebrate extremes. Nat Ver Genet. 2004, 5: 456-65. 10.1038/nrg1350.View ArticleGoogle Scholar
- Majewski J, Ott J: Distribution and characterization of regulatory elements in the human genome. Genome Res. 2002, 12: 1827-1836. 10.1101/gr.606402.PubMedPubMed CentralView ArticleGoogle Scholar
- Baek G, Green P: Sequence conservation, relative isoform frequencies, and nonsense-mediated decay in evolutionarily conserved alternative splicing. Proc Natl Acad Sci. 2005, 102: 12813-1288. 10.1073/pnas.0506139102.PubMedPubMed CentralView ArticleGoogle Scholar
- Sorek R, Ast G: Intronic sequences flanking alternatively spliced exons are conserved between human and mouse. Genome Res. 2003, 13: 1631-7. 10.1101/gr.1208803.PubMedPubMed CentralView ArticleGoogle Scholar
- Fulop C, Cs-Szabo G, Glant TT: Species-specific alternative splicing of the epidermal growth factor-like domain 1 of cartilage aggrecan. Biochem J. 1996, 319: 935-940.PubMedPubMed CentralView ArticleGoogle Scholar
- Laverdiere M, Beaudoin J, Lavigueur A: Species-specific regulation of alternative splicing in the C-terminal region of the p53 tumor suppressor gene. Nucl Acids Res. 2000, 28: 1489-1497. 10.1093/nar/28.6.1489.PubMedPubMed CentralView ArticleGoogle Scholar
- Bingham PM, Chou TB, Mims I, Zachar Z: On/off regulation of gene expression at the level of splicing. Trends Genet. 1988, 4: 134-138. 10.1016/0168-9525(88)90136-9.PubMedView ArticleGoogle Scholar
- Lareau LF, Green RE, Bhatnagar RS, Brenner SE: The evolving roles of alternative splicing. Curr Opin Struct Biol. 2004, 14: 273-82. 10.1016/j.sbi.2004.05.002.PubMedView ArticleGoogle Scholar
- Shafqat N, Marschall HU, Filling C, Nordling E, Wu XQ, Bjork L, Thyberg J, Martensson E, Salim S, Jornvall H, Oppermann U: Expanded substrate screenings of human and Drosophila type 10 17beta-hydroxysteroid dehydrogenases (HSDs) reveal multiple specificities in bile acid and steroid hormone metabolism: characterization of multifunctional 3alpha/7alpha/7beta/17beta/20beta/21-HSD. Biochem J. 2003, 376: 49-60. 10.1042/BJ20030877.PubMedPubMed CentralView ArticleGoogle Scholar
- Torroja L, Ortuno-Sahagunet D, Ferrus A, Hammer B, Barbas JA: scully, an essential gene of Drosophila, is homologous to mammalian mitochondrial type II L-3-hydroxyacyl-CoA dehydrogenase/amyloid-b peptide-binding protein. J Cell Biol. 1998, 141: 1009-1017. 10.1083/jcb.141.4.1009.PubMedPubMed CentralView ArticleGoogle Scholar
- Yan SD, Zhu Y, Stern ED, Hwang YC, Hori O, Ogawa S, Frosch MP, Connolly ES, McTaggert R, Pinsky DJ, Clarke S, Stern DM, Ramasamy R: Amyloid beta-peptide-binding alcohol dehydrogenase is a component of the cellular response to nutricional stress. J Biol Chem. 2000, 275: 27100-27109.PubMedGoogle Scholar
- Anisimova M, Bielawski JP, Yang Z: Accuracy and power of bayes prediction of amino acid sites under positive selection. Mol Biol Evol. 2002, 19: 950-958.PubMedView ArticleGoogle Scholar
- McClellan DA, Palfreyman EJ, Smith MJ, Moss JL, Christensen RG, Sailsbery JK: Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins. Mol Biol Evol. 2005, 22: 437-455. 10.1093/molbev/msi028.PubMedView ArticleGoogle Scholar
- Ensembl. [http://www.ensembl.org]
- NCBI. [http://www.ncbi.nlm.nih.gov]
- TIGR. [http://www.tigr.org]
- RepeatMasker. [http://repeatmasker.org]
- CENSOR. [http://www.girinst.org]
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acid Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.PubMedPubMed CentralView ArticleGoogle Scholar
- Schwartz S, Zhang Z, Frazer KA, Smit A, Riemer C, Bouck J, Gibbs R, Hardison R, Miller W: PipMaker – a web server for aligning two genomic DNA sequences. Genome Res. 2000, 10 (4): 577-586. 10.1101/gr.10.4.577.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: ClustalW: improving the sensitivity of progressive multiple sequence alignment through sequence weighing, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.PubMedPubMed CentralView ArticleGoogle Scholar
- Pride DT, Blaser MJ: Concerted evolution between duplicated genetic elements in Helicobacter pylori. J Mol Bio. 2002, 316: 629-42. 10.1006/jmbi.2001.5311.View ArticleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.PubMedView ArticleGoogle Scholar
- Sawyer SA: GENECONV: a computer package for the statistical detection of gene conversion. Distributed by the author. 1999, University of Washington in St. Louis, Department of Mathematics, [http://www.math.wustl.edu/~sawyer]Google Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView ArticleGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated Software for Molecular Evolutionary Genetics Analysis and Sequence Alignment. Briefings in Bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.PubMedView ArticleGoogle Scholar
- Swofford DL: PAUP *: phylogenetic analysis using parsimony (*and other methods). Version 4.0b10. 2002, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Graybeal A: Evaluating the phylogenetic utility of genes: search for genes informative about deep divergences among vertebrates. Syst Biol. 1994, 43: 174-193. 10.2307/2413460.View ArticleGoogle Scholar
- Irwin DM, Kocher TD, Wilson AC: Evolution of the cytochrome b gene of mammals. J Mol Evol. 1991, 32: 128-144.PubMedView ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F: MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13: 555-556.PubMedGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen A-MK: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 431-449.PubMedPubMed CentralGoogle Scholar
- Yang Z, Wong WSW, Nielsen R: Bayes Empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.PubMedView ArticleGoogle Scholar
- Woolley S, Johnson J, Smith MJ, Crandall KA, McClellan DA: TreeSAAP: selection on amino acid properties using phylogenetic trees. Bioinformatics. 2003, 22: 671-2. 10.1093/bioinformatics/btg043.View ArticleGoogle Scholar
- McClellan DA, McCracken KG: Estimating the influence of selection on the variable amino acid sites of the cytochrome b protein functional domain. Mol Biol Evol. 2001, 18: 917-925.PubMedView ArticleGoogle Scholar
- Humphrey W, Dalke A, Schulten K: VMD – Visual Molecular Dynamics. J Molec Graphics. 1996, 14: 33-38. 10.1016/0263-7855(96)00018-5.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.