TY - JOUR AU - Shi, Peng AB - Abstract Auditory detection is essential for survival and reproduction of vertebrates, yet the genetic changes underlying the evolution and diversity of hearing are poorly documented. Recent discoveries concerning prestin, which is responsible for cochlear amplification by electromotility, provide an opportunity to redress this situation. We identify prestin genes from the genomes of 14 vertebrates, including three fishes, one amphibian, one lizard, one bird, and eight mammals. An evolutionary analysis of these sequences and 34 previously known prestin genes reveals for the first time that this hearing gene was under positive selection in the most recent common ancestor (MRCA) of tetrapods. This discovery might document the genetic basis of enhanced high sound sensibility in tetrapods. An investigation of the adaptive gain and evolution of electromotility, an important evolutionary innovation for the highest hearing ability of mammals, detects evidence for positive selections on the MRCA of mammals, therians, and placentals, respectively. It is suggested that electromotility determined by prestin might initially appear in the MRCA of mammals, and its functional improvements might occur in the MRCA of therian and placental mammals. Our patch clamp experiments further support this hypothesis, revealing the functional divergence of voltage-dependent nonlinear capacitance of prestin from platypus, opossum, and gerbil. Moreover, structure-based cdocking analyses detect positively selected amino acids in the MRCA of placental mammals that are key residues in sulfate anion transport. This study provides new insights into the adaptation and functional diversity of hearing sensitivity in vertebrates by evolutionary and functional analysis of the hearing gene prestin. prestin, adaptation, selection, electromotility, tetrapods, NLC Introduction The vertebrate auditory system is responsible for the detection of sound signals in external environments, which is essential for survival and reproduction. The ability to hear is typically characterized by parameters that quantify frequency selectivity and sensitivity. These parameters generally include the high-frequency hearing limit (Fmax) (Masterton et al. 1969). Previous behavioral studies reveal that Fmax varies enormously across different vertebrate lineages (Fay 1996). For example, most fishes can detect sounds ranging from ∼0.44 to ∼0.74 kHz (Fay 1988; Popper and Fay 1997), although a few species are able to detect ultrasonic sounds (Popper 2000; Mann et al. 2001). In most species of frogs, Fmax values are about three to five times higher than those in most fishes, but ∼5-fold lower than those of most of birds (Fay 1988). In general, mammals detect frequencies above 20 kHz, and this is an important characteristic distinguishing them from other tetrapods (Fay 1988; Manley 1990, 2000). High-frequency hearing detection tends to increase phylogenetically from fishes to amphibians to mammals. Certainly, a few species in each lineage have higher frequency detection than their cohorts (Popper 2000; Mann et al. 2001; Feng et al. 2006). This trend is also documented in anatomical and electrophysiological studies (Heffner and Heffner 1991; Fay 1996; Coffin et al. 2004), yet little is known about the underlying genetic variation and evolutionary trajectories of the gene(s) involved in high-frequency detection. Prestin, a membrane motor protein driving the elongation of outer hair cells (OHCs), is a major determinant of cochlear amplification for hearing-frequency hearing (Brownell et al. 1985; Belyantseva et al. 2000; Zheng et al. 2000; Brownell et al. 2001). Prestin, also named SLC26A5, is a member of the solute carrier 26 (SLC26) gene family. It encodes anion exchangers capable of transporting a wide variety of monovalent and divalent anions (Sindic et al. 2007; Detro-Dassen et al. 2008). Some 10–12 transmembrane domains linked by intra- and extracellular loops in prestin are flanked by cytoplasmic N- and C-termini (Navaratnam et al. 2005; Zheng et al. 2005). Prestin, a motor protein essential for electromotility (Zheng et al. 2000; Dallos and Fakler 2002; Liberman et al. 2002), can be measured by its robust voltage-dependent nonlinear capacitance (NLC) (Santos-Sacchi 1991). In mammals, NLC and high-frequency hearing increase to a greater extent than in other vertebrates without electromotility (Fay 1988; Manley 1990, 2000; Tan et al. 2011). In addition, prestin-knockout mice show significantly reduced high-frequency sensitivity (Liberman et al. 2002; Wu et al. 2004). Mutations in this gene lead to nonsyndromicauditory loss in humans (Liu et al. 2003). Phylogenetic analysis of prestin sequence data unite echolocating bats that have the ability of high-frequency sound detection (Li et al. 2008), although these bats are not monophyletic in the species tree (Teeling et al. 2005). Surprisingly, echolocating whales and bats also cluster together in phylogenetic analyses of amino acid sequences (Li et al. 2010; Liu, Cotton, et al. 2010). The number of amino acid replacements is positively corrected with an increase in ability to detect high frequencies in whales (Liu, Rossiter, et al. 2010). Taken together, these studies strongly suggest that prestin is an important factor in high-frequency detection. Herein, we examine the evolutionary dynamics and selection pressure on prestin genes in vertebrate lineages. We focus on detecting selection pressure in the initial emergence of land vertebrates because these animals tend to detect higher frequencies of sounds than do fishes. This is followed by an investigation of the unique gain of somatic electromotility in mammalian OHCs, the basis for cochlear amplification (Brownell et al. 1985; Ashmore 1987). This amplification system gives mammals the greatest ability to detect the highest frequencies among vertebrates. Our analyses of 48 vertebrate prestin genes detect multiple episodes of adaptive evolution among vertebrates, suggesting that its function might change as organisms evolve. Furthermore, our results from whole-cell patch clamp functional experiments on platypus, opossum, and gerbil prestin support this hypothesis. This study provides insights into the adaptation and functional diversity of high frequency hearing in vertebrates. Materials and Methods Identification of prestin Genes and Collection of Fmax Data in Mammals Thirty-four prestin sequences obtained from GenBank included two birds, two fishes, and 30 mammals. In addition, we searched for prestin sequences in Ensembl (http://www.ensembl.org/) and NCBI (http://www.ncbi.nlm.nih.gov/) from 14 vertebrate genomes that have high genome coverage (≥6×). The taxa included three fishes (medaka, Oryzias latipes; stickleback, Gasterosteus aculeatus, and fugu, Takifugu rubripes), one amphibian (frog, Xenopus tropicalis), one nonavian reptile (lizard, Anolis carolinensis), one bird (turkey, Meleagris gallopavo), and eight mammals (chimpanzee, Pan troglodytes; gorilla, Gorilla gorilla; orangutan, Pongo pygmaeus; common marmoset, Callithrix jacchus; elephant, Loxodonta africana; guinea pig, Cavia porcellus; panda, Ailuropoda melanoleuca; and horse, Equus caballus) (supplementary table 1, Supplementary Material online). We used our previous pipeline for identifying genes with multiple exons (Yang, Shi, et al. 2005;,Liu et al. 2011) as briefly follows: TBlastN was employed to search prestin sequences in the genome databases using previously known prestin protein sequences, and known protein sequences of prestins and the best hit genomic sequences were used to conduct a protein–DNA comparison using Wise2 (http://www.ebi.ac.uk/Tools/Wise2/index.html), which provided the exon/intron structures and the full-length protein sequences and cDNA sequences of the putative genes. We did not search low-coverage mammalian genomes (about 2× coverage) from each major vertebrate lineage because sequence mining of multiexon genes was not feasible without chromosomal assemblies. Two analyses were performed in order to exclude false prestin homologies. First, we blasted the putative genes in GenBank to ensure the best hits were known prestin genes. Second, we constructed a tree using the neighbor-joining method (Saitou and Nei 1987) with protein Poisson distances (Nei and Kumar 2000). We employed SLC26A6, the closest related gene to prestin in the SLC26 gene family (Franchini and Elgoyhen 2006) as the outgroup to root the tree (supplementary fig. 1, Supplementary Material online). We collected Fmax values, which measured high-frequency hearing and defined the highest frequency audible at 60 dB sound pressure level (SPL) (Fay 1996), from mammals in order to detect differences among species. Values were collected from two monotremes, the platypus (Gates et al. 1974) and echidna (Mills and Shepherd 2001), three marsupials, the opossum (Reimer 1995), quolls (Dasyursus; Aitkin 1995) and brushtail possum (Trichosurus; Aitkin 1995), and 23 placentals, including the elephant (Heffner and Heffner 1982), guinea pig (Heffner et al. 1971), horseshoe bat (Long 1977), rat (Gourevitch 1965), bottlenose dolphin (Herman and Arbeit 1973), mouse (Ehret 1976), chinchilla (Clack 1966), macaque (Clack 1966), cat (Nienhuys and Clark 1979), human (Wier et al. 1977), squirrel monkey (Green 1975), yellow baboon (Hienz et al. 1982), lemur (Mitchell et al. 1971), bushbaby (Heffner et al. 1969b), tree shrew (Heffner et al. 1969a), dog (Heffner 1983), ferret (Kelly et al. 1986), kangaroo rat (Heffner and Masterton 1980), gerbil (Ryan 1976), rabbit (Heffner and Masterton 1980), little brown bat (Dalland 1965), common harbor seal (Mohl 1968), and frog-eating bat (Ryan et al. 1983). Evolutionary Analyses The prestin sequences were initially aligned using ClustalW (Chenna et al. 2003) followed by manual adjustments. Pairwise comparisons of the numbers of synonymous substitutions per synonymous site (dS) and nonsynonymous substitutions per nonsynonymous site (dN) were estimated by the modified Nei–Gojobori method (Zhang et al. 1998) in MEGA5 (Tamura et al. 2011). Tree-based selection tests were calculated by the branch-site likelihood method implemented in PAML4 (Yang, Wong, et al. 2005; Zhang et al. 2005) because this test is conservative and robust against violations of various model assumptions, unlike the previous version of the test (Yang and Nielsen 2002), which has a relatively high rate of false positives (Zhang 2004). This method compared two models. One model defined four classes of sites in terms of ω, the ratio of dN to dS. Codons conserved throughout the tree were assigned to site class 0 (estimated 0< ω0 <1). Site class 1 contained codons that were neutral throughout the tree (estimated ω1 = 1). Site class 2a included codons conserved on the background versus foreground branches (0 < ω0 < 1). These were assumed to fall under positive selection (ω2 > 1) on the foreground branches. Finally, site class 2b was comprised of codons assumed to be neutral (ω1 = 1) on the background branches but positive selection on the foreground branch. The corresponding null model fixed ω2 = 1, which differed from the alternative model. A likelihood ratio test was used to compare the two models. When the likelihood of the alternative model was significantly higher than that of the null model, it was assumed to indicate positive selection on the foreground branches. Where tests indicated positive selection, we recorded the sites under selection according to high posterior probabilities (≥0.95) following Bayes empirical Bayes (BEB) prediction (Yang, Wong, et al. 2005). Gene Synthesis, Cell Culture, and Transient Transfection The entire coding regions of platypus and opossum were synthesized and cloned into the expression vector pcDNA3.1(−) (Invitrogen Inc.). Correct orientation and reading frames were verified by sequencing analysis. In addition, the expression vector of gerbil's prestin was gifted from Dallos' lab (Northwestern University). HEK293 cells were grown in 35-mm dishes containing Dulbecco's modified Eagle's medium supplemented with 10% bovine calf serum (Invitrogen Inc.). When cell confluence reached roughly 50–60% of the surface area of the dishes, cotransfection of the expression vectors of the prestins and pEGFP-N1 (GFP) were accomplished using lipofectamine 2000 transfection reagent (Invitrogen Inc.). We used a ratio of 1 μg (GFP)/3 μg (prestin), added to 10-μl lipofectamine. The pEGFP-N1 plasmid generated a cytoplasmic EGFP protein as an independent marker for successful transfection of cells. After 24- to 48-h incubation, the successfully transfected cells were used for NLC measurements. Electrophysiological Experiments for NLC Measurements NLC was measured using whole-cell patch-clamp recordings that were performed by HEKA EPC 10 USB (HEKA Instruments Inc.) at room temperature (22–26 °C). Electrodes were pulled from borosilicate glass with resistances of 2.5–4 MΩ and filled with the internal solution containing: 140 mM CsCl, 2 mM MgCl2, 10 mM EGTA, and 10 mM HEPES. The cells were bathed during the recordings in an external solution containing 120 mM NaCl, 20 mM TEA-Cl, 2 mM CoCl, 2 mM MgCl2, 10 mM HEPES, and 5 mM glucose. Both solutions were adjusted to pH 7.2. Osmolarities of the internal and external solutions were adjusted to 300 and 320 mOsml−1 with glucose, respectively. Voltage-dependent capacitance was measured using the stair-step voltage protocol to obtain the parameters of charge movement (Huang and Santos-Sacchi 1993). Voltage was stepped from −140 to 100 mV in 10 mV increments of 10 ms each. The capacitive currents were sampled at 100 kHz and low-pass filtered at 5 kHz using PatchMaster software (HEKA Instruments Inc.). For each voltage, the measured membrane capacitance (Cm) was plotted as a function of membrane voltage (Vm) and fitted with the derivative of a two-state Boltzmann function:   Cm=Qmaxαexp[α(Vm−V1/2)](1+exp[−α(Vm−V1/2)])2+Clin, where Qmax was the maximum charge transfer, V1/2 was the voltage at which the maximum charge was equally distributed across the membrane, Clin was the linear capacitance, and α = ze/kT was the slope factor of the voltage dependence of charge transfer where k was Boltzmann's constant, T was absolute temperature, z was valence, and e was electron charge. The Clin was proportional to the surface area of the membrane (cell size). To compare the magnitude of NLC obtained from different cells with different levels of prestin expression as a function of cell size, we normalized the NLC by the linear capacitance of the cells. Because differences in Qmax could have been caused by cell size, the charge movement was normalized to Clin. This quantity, designated as charge density, had units of fC/pF. Predicting Prestin Protein Tertiary Structure and Structural Superposition The tertiary structure of prestin was predicted online by the profile–profile matching algorithms implemented in Phyre (http://www.sbg.bio.ic.ac.uk/∼phyre/html/casp8.html). The quality of the predicted proteins was estimated by E value, where an E value < 0.001 corresponded to an estimated precision >95% (Kelley and Sternberg 2009). Structural superposition was performed using the CE software package (Shindyalov and Bourne 1998). Root mean square deviation (RMSD) values between two protein structures were used to measure the degree of structural similarity for superposition. RMSD values <4 Å were assumed to indicate similar overall structures (Tung and Yang 2007). Binding Affinity of Prestin Protein and Anions The binding affinity between prestin and anions (Cl−, HCO3−, and SO42−) was assessed using CDOCKER (Wu et al. 2003) in Accelrys Discovery Studio 2.1 (Accelrys Software Inc.), a Charmm-based molecular docking tool for analyzing receptor-ligand interactions. The binding affinity was quantified by CDOCKER interaction energy (CIE), the interaction energy between the proteins and their binding ligands. Generally, CIE < 0 denoted the attraction of proteins and ligands; the lower the value the higher the attraction. In contrast, CIE > 0 indicated repulsion between ligands and proteins, that is, ligands could not bind to acceptor proteins; the higher the CIE value the stronger the repulsive force. Results and Discussion Vertebrate prestin Genes Herein, a total of 48 prestin sequences representing all major vertebrate lineages except for turtles and crocodilians were subjected to analyses (fig. 1). Alignment of the amino acid sequences (supplementary fig. 2, Supplementary Material online) revealed that the N- and C-termini contained the highest levels of variability. We computed Poisson-corrected evolutionary distances for each region in each major lineage to further quantify regional sequence diversity. Gene regions were assigned according to the prestin membrane topology proposed by several lines of functional evidence (Ludwig et al. 2001; Matsuda et al. 2004; Deak et al. 2005; Rajagopalan et al. 2006). The mean sequence distance in the N-terminal region (NR) (fig. 2) was significantly higher than that of the full-length sequence in both teleosts (P = 0.0047, Fisher's Exact Test) and tetrapods (P = 0.0027). The same was true for the C-terminus (teleosts: P = 0.0005; tetrapods: P = 0.0001). In contrast, the sequences of the transmembrane region (TR) in both lineages were more conserved than those of other regions. These comparisons revealed that NRs and C-terminal regions (CRs) evolved faster than TRs at the protein sequence level in both major vertebrate lineages. This result taken together with a previous result in mammals (Okoruwa et al. 2008) suggested high levels of variations in NRs and CRs might comprise a general pattern in vertebrate prestin genes. FIG.1. View large Download slide Phylogeny of vertebrates considered herein (Murphy et al. 2004). Black blocks mark the species in which the prestin genes are newly identified in this study. The species of bats and whales have been collapsed for illustration and detailed species information is presented in supplementary table 1 (Supplementary Material online). Positive selection tests are indicated by square labels. FIG.1. View large Download slide Phylogeny of vertebrates considered herein (Murphy et al. 2004). Black blocks mark the species in which the prestin genes are newly identified in this study. The species of bats and whales have been collapsed for illustration and detailed species information is presented in supplementary table 1 (Supplementary Material online). Positive selection tests are indicated by square labels. FIG. 2. View largeDownload slide Poisson distance of different regions of the prestin gene in teleosts and tetrapods. The dashed lines show the average Poisson distance for the full length of the protein in teleosts (light gray) and in tetrapods (dark). An asterisk indicates the significant difference of Poisson distance between different regions and the full-length protein. FIG. 2. View largeDownload slide Poisson distance of different regions of the prestin gene in teleosts and tetrapods. The dashed lines show the average Poisson distance for the full length of the protein in teleosts (light gray) and in tetrapods (dark). An asterisk indicates the significant difference of Poisson distance between different regions and the full-length protein. Adaptive Evolution of prestin Genes during the Emergence of Tetrapods To further understand the evolutionary dynamics and selective pressure on prestin genes between fishes and tetrapods, we compared the mean nonsynonymous nucleotide substitution distances between teleosts and tetrapods in NR, CR, TRs, extracellular regions (ERs), and intracellular regions (IRs) using the modified Nei–Gojobori method (Zhang et al. 1998). As presented in figure 3, dN values were higher in NRs and CRs than in ER, IR, and TR. This suggested either that positive selection had acted to favor amino acid replacements in the terminal regions or that some of the amino acids in these regions enjoyed fewer functional constraints. FIG. 3. View largeDownload slide Pairwise synonymous (dS) (filled circles) and nonsynonymous (dN) nucleotide distances (open triangles), and dN/dS ratio (filled columns) for different regions of the prestin gene between teleosts and tetrapods. FIG. 3. View largeDownload slide Pairwise synonymous (dS) (filled circles) and nonsynonymous (dN) nucleotide distances (open triangles), and dN/dS ratio (filled columns) for different regions of the prestin gene between teleosts and tetrapods. To distinguish between positive selection and the absence of functional constraints, we compared dN and dS, because the former was expected to exceed the latter in cases of positive selection. We estimated the mean synonymous distance for the five regions between teleosts and tetrapods and found that the mean dN values were much lower than dS values in ER, IR, and TR. In contrast, the dN/dS ratio was larger than 1.0 in NR, although the mean dN value was not significantly higher than dS(fig. 3). The absence of significance between dN and dS might have owed to the saturation of synonymous and nonsynonymous substitutions (Tanaka and Nei 1989). Positive selection between teleosts and tetrapods was difficult to detect by pairwise comparisons at this high level of sequence divergence (0.5 < dS < 1) because subsequent substitutions may have hidden the signal of positive selection, especially if the emergence of tetrapods occurred in a relatively short evolutionary time span (Zhang et al. 1998; Shi et al. 2003). To rectify this methodological limitation, we used an improved branch-site likelihood method (Zhang et al. 2005) to look for selection signal while restricting our signal test to the ancestral branch of tetrapods. This method more reliably and powerfully detected positive selection at this level of sequence divergence and saturation by computational simulation (Anisimova et al. 2001, 2002). We assigned the ancestral branch of tetrapods as the foreground branch and all others as background branches (fig. 1). This approach obtained a significantly higher likelihood for the alternative model than that of the null model (P < 0.001, χ2 test), suggesting positive selection on prestin in the ancestor of tetrapods. Conservative BEB (Yang, Wong, et al. 2005; Zhang et al. 2005) identified 20 sites under positive selection and with posterior probabilities > 0.95 (table 1). In addition, almost half of these sites were located in NRs and CRs. This could have explained the higher dN values found in these two regions. Table 1. Detection of Positive Selection in the Different Lineages of Prestin Orthologous Genes.     a Twice the difference between the log likelihood of the alternative model and that of the null model. The modified model A with ω2 fixed at 1 is the null model. The modified model A is used as the alternative model. b Multiple testing corrections are performed. c ω values are the nonsynonymous/synonymous rate ratios. p0 is the proportion of codons that have ω0 in all branches, p1 is the proportion of codons that have ω1 = 1, in all branches, p2a is the proportion of codons that have ω0 in the background branches but ω2 in the foreground branches, and p2b is the proportion of codons that have ω1 in the background branches but ω2 in the foreground branches. Note that as long as ω2 significantly exceeds 1 (as indicated by the likelihood ratio test), its exact value has little biological meaning due to the large estimation error. d Sites with the Bayes empirical Bayes posterior probabilities higher than 95% are shown. The sites are indexed by the amino acids at the site in the gerbil prestin. The sites in different regions are shown as followings: italic in NR and CR, single underlined in TR, double underlined in ER, and boxed in IR. View Large Functional and morphological audition innovations in tetrapods are believed to be adaptations for processing airborne sound (Fritzsch 1991). For example, when the ancestral tetrapods moved onto land, the spiracular pouch in fishes transformed into a tympanic middle ear, the hyomandibular bone transformed into the stapes (Gaupp 1898, 1913; Werner 1960; Thomson 1966; Lombard and Bolt 1988), and the basilar papilla formed as a unique sensory adaptation for airborne sound detection (Retzius 1881, 1884). Consistent with these morphological innovations, the rate of amino acid replacements in prestin genes, which has been shown to be highly related to high-frequency sensitivity, dramatically changed at the same time. The positive selection tests indicated an accelerated rate of amino acid changes in the most recent common ancestor (MRCA) of tetrapods, suggesting that prestin genes might have been involved in the functional shift from low-frequency hearing in fishes to higher frequency audition in tetrapods. To our knowledge, this analysis provided the first evidence that positive selection on prestin genes might have enhanced the high-frequency hearing of tetrapods. Multiple Episodic Adaptive Events on Mammalian prestin Genes Sound amplification via electromotility gives mammals a superior ability to hear high-frequency sound (Fay 1988; Manley 1990, 2000). This represents a major evolutionary advancement. Therefore, it is important to evaluate the evolutionary tempo and mode of prestin in mammals because the protein is responsible for somatic electromotility (Brownell et al. 1985; Ashmore 1987; Zheng et al. 2000). Using placental and marsupial prestin sequences, Franchini and Elgoyhen (2006) report positive selection signals in the MRCA of mammals. However, because mammals consist of three major clades of monotremes, marsupials, and placentals, the absence of monotreme prestin sequences in their study precludes understandings of the evolution of this gene and its function in mammals. The available platypus genome (Warren et al. 2008) and a clone of platypus prestin (Okoruwa et al. 2008) provide an opportunity to reexamine selective pressures on prestin genes for the MRCA of all mammals. The test assigns the ancestral branch, that leading to all mammals, to the foreground and all other amniotes to background branches (fig. 1). The results (table 1) show a significantly higher likelihood of the alternative model than that of the null model (after multiple testing correction P < 0.05, χ2 test), indicating that the MRCA of mammals experienced positive selection. This result supports the conclusion that the origin of electromotility in mammals likely happened in the MRCA of all mammals (Franchini and Elgoyhen 2006). Furthermore, analyses of all available audiogram data from monotremes, marsupials, and placentals reveals that the average upper hearing limit in monotremes (12 kHz) is significantly lower than that of therians (38 kHz, P < 0.01, t-test ) and placentals (61.9 kHz, P < 0.01). These observations lead to the hypothesis that prestin experienced additional adaptive selection for detecting high-frequency sound shortly after the origin of electromotility during the evolution of mammals. To test this hypothesis, we examined selection pressures on prestin genes on the ancestral branches that lead to therians and placentals respectively. The two branches were separately assigned as foreground branches and all other amniotes as background branches (fig. 1). Significant signals of positive selection were detected on the ancestral branches of therians and placentals, respectively (table 1). We also identified 23, 8, and 6 positively selected sites in the MRCA of mammals, therians, and placentals, respectively (table 1). The inferred positively selected amino acid substitutions on these three branches might have provided evidence for functional divergence. Functional Variations of prestin in Platypus, Opossum, and Gerbil If the ancestral branches of mammals, therians, and placentals are indeed under positive selection as described above, we would expect that functional changes of prestin might be observed in monotreme, marsupial, and placental mammals. Voltage-dependent NLC is one of the unique characteristics of prestin, and it is often used to measure prestin function (Santos-Sacchi 1991). NLC exhibits a bell-shaped dependence on membrane potential, and it can be fitted with the first derivative of a two-state Boltzmann function (Santos-Sacchi 1991, Oliver et al. 2001). Consequently, we measure the NLC of prestin genes from the platypus, opossum, and gerbil, the representative monotreme, marsupial, and placental mammals, respectively. As shown in figure 4, functional variations occur in three mammalian prestin genes. For platypus prestin, the parameters of NLC from 12 cells after fitting two-state Boltzmann function are as follows: Qmax/Clin = 15 ± 2.3fC/pF, V1/2 = −38.4 ±5.3 mV, 1/α = 40.57 ± 2.4 mV (mean ± S.E.). These values are consistent with previous reports (Tan et al. 2011). Compared with the platypus' prestin, 1/α value of opossum normalized from 25 cells is significantly lower (P = 0.01, Student's t-test, fig. 4A), suggesting changes in the reactivity of prestin to the membrane charge transfer. Other parameters are also well fitted a two-state Boltzmann function curve, with following values of Qmax/Clin = 9.7 ± 1.6 fC/pF, V1/2 = −30.1 ± 3.7 mV, and α = 58.98 ± 4.3 mV. FIG. 4. View largeDownload slide Voltage-dependent membrane capacitance (NLC) of HEK cells transiently transfected with prestin orthologs of platypus, opossum, and gerbil, respectively. The capacitance–voltage plots are fitted with two-state Boltzmann function. (A) 1/α values of prestin in three species are analyzed: platypus, 40.57 ± 2.4 mV (n = 12); opossum, 58.98 ± 4.3 mV (n = 25); and gerbil, 35.82 ± 2.7 mV (n = 20). (B) Comparison of V1/2 for three mammalian species: platypus, −38.4 ± 5.3 mV (n = 12); opossum, −30.1 ± 3.7 mV (n = 25); and gerbil, −67.9 ± 4.2 mV (n = 20). (C) Charge density for all three species is as follows: platypus, 15 ± 2.3 fC/pF (n = 12); opossum, 9.7 ± 1.6 fC/pF (n = 25); and gerbil, 19.7 ± 2.7 fC/pF (n = 20). All values are mean ± S.E. *P < 0.05, **P < 0.01. FIG. 4. View largeDownload slide Voltage-dependent membrane capacitance (NLC) of HEK cells transiently transfected with prestin orthologs of platypus, opossum, and gerbil, respectively. The capacitance–voltage plots are fitted with two-state Boltzmann function. (A) 1/α values of prestin in three species are analyzed: platypus, 40.57 ± 2.4 mV (n = 12); opossum, 58.98 ± 4.3 mV (n = 25); and gerbil, 35.82 ± 2.7 mV (n = 20). (B) Comparison of V1/2 for three mammalian species: platypus, −38.4 ± 5.3 mV (n = 12); opossum, −30.1 ± 3.7 mV (n = 25); and gerbil, −67.9 ± 4.2 mV (n = 20). (C) Charge density for all three species is as follows: platypus, 15 ± 2.3 fC/pF (n = 12); opossum, 9.7 ± 1.6 fC/pF (n = 25); and gerbil, 19.7 ± 2.7 fC/pF (n = 20). All values are mean ± S.E. *P < 0.05, **P < 0.01. In comparison with platypus and opossum prestin, the peak voltage of NLC (V1/2) of gerbil prestin is significantly shifted toward the hyperpolarizing direction (P < 0.01, Student's t-test, fig. 4B), with a value of −67.9 mv (n = 20). Furthermore, the curve-fitting parameters of charge density (Qmax/Clin = 19.7 ± 2.7 fC/pF) and 1/α value (35.82 ± 2.7 mV) also differ significantly from that of opossum prestin (P < 0.01, Student's t-test, fig. 4C). In addition to functional variation among the three mammalian lineages, functional changes occur between nonmammalian vertebrates and mammals and between fish and tetrapods. For example, whereas the prestin gene of all three mammals exhibits a robust bell-shaped voltage-dependent NLC, those of the zebrafish and chicken do not. Furthermore, the magnitude of NLC in chicken prestin is considerably larger than that of zebrafish (Tan et al. 2011). The generation of robust bell-shaped NLC in mammalian prestins from platypus and opossum to gerbil, as well as the functional improvement of NLC in nonmammalian vertebrates (e.g. zebrafish and chicken), might owe to multiple positive selection events on the ancestral branches of tetrapods, all mammals, therians, and placentals. This possibility requires that positively selected sites involve functional changes and site-directed mutagenesis studies provide strong support for this. For example, the NLC experiment of chimera gerbil prestin, constructed by exchanging 225M, a positively selected site on the ancestral branch of therians, into the corresponding site of gerbil prestin reveals functional changes of NLC and the motility of prestin-expressing cells (Kumano et al. 2009). Another site 260C on the same branch and site 196C on the ancestral branch of placentals also play an important role in functional changes of prestin (Rajagopalan et al. 2006; Kumano et al. 2009; McGuire et al. 2010). The mutagenesis of 415C, a positively selected site in the MRCA of mammals, can significantly decrease the magnitudes of NLC, suggesting 415C is required for the increase of NLC in mammals (McGuire et al. 2010). Amino acid 192C is inferred to be positively selected on the ancestral branch of tetrapods, and it plays an important role in charge movement of prestin (McGuire et al. 2010). All of these sites are in the list of positively selected sites (table 1). Thus, our sequence analysis is consistent with experimental results, and it may help discover more key functional sites of prestin, which is associated with the gain and subsequent evolution of NLC and electromotility. More importantly, our evolutionary analysis of prestin genes, taken together with behavioral and functional analyses, suggests that the gene experienced at least three adaptive selection events in mammals alone. The functional improvement of prestin might be a very complex stepwise process in mammals. Positive Selection in the MRCA of Placental Mammals for Transporting Sulfate In addition to changing NLC of electromotility, the function of prestin differs substantially among vertebrates with respect to its ability to transport sulfate. For example, in the zebrafish and chicken, prestin acts as an electrogenic antiporter exchanging SO42− for Cl− with a 1:1 stoichiometry (Schaechinger and Oliver 2007). This plesiomorphic function is not known to occur in mammals, including gerbils (Oliver et al. 2001). The vertebrate lineage that experienced this functional change remains to be identified. Regardless, functional data suggest that this transformation might be associated with structural conformation changes of prestin (Schaechinger and Oliver 2007). The 3D structure of prestin facilitates an evaluation of whether or not structural conformation changes are involved in functional sulfate transport. Whereas the 3D structure of the C-terminus is known (Pasqualetto et al. 2010), it remains unknown for TRs. Because these are important functional domains for anion transport (Bai et al. 2009; McGuire et al. 2010), our understanding of the changes in sulfate transport mechanisms are limited. Upon using Phyre to predict the 3D structure of TRs of prestin, the best hit of the gerbil's prestin is a chloride channel (PDB ID = 1ots) (Dutzler et al. 2003; E value = 3.1 × 10−5). All other vertebrate prestins hit the same model (1ots) with a predicted high accuracy (≥95%). Significant E values are found to range from 10−3 to 10−5, although the similarity between the template and query sequences is not so high (alignments in supplementary fig. 3, Supplementary Material online). Technological limitations of 3D modeling and the relatively few available crystallographic structures of membrane proteins require us to evaluate the reliability of the predicted 3D structure of prestin. The predicted structure agrees with the key properties of prestin. First, the predicted structure has 12-transmembrane helixes, and this is consistent with the secondary topology demonstrated by most functional and modeling assays (Oliver et al. 2001; Deak et al. 2005; Rajagopalan et al. 2006). Second, both template and prestin are anion channels that function to conduct Cl− across cell membranes in all vertebrates (Oliver et al. 2001; Dutzler et al. 2003; Schaechinger and Oliver 2007). Third, the accuracy of our 3D model can be validated by comparisons of the predicted functions by molecular docking and functional assays. Whereas a positive CIE value implies little or no affinity between prestin and anions, a negative value suggests that prestin can drive anions freely through the membrane. If our predicted prestin structure holds true, then all of the prestins in nonmammals and mammals should have a high affinity for Cl− and HCO3− and negative CIE values. Here, Cl− and HCO3− should be freely transported by prestins in all vertebrates, as evidenced by functional experiments on the zebrafish, chicken, and gerbil (Oliver et al. 2001; Schaechinger and Oliver 2007). As expected, the CIEs between prestins and Cl−/HCO3− in the zebrafish, chicken, and gerbil are all negative (data not shown). For SO42−, functional assays (Oliver et al. 2001; Schaechinger and Oliver 2007) predict that the zebrafish and chicken prestins will have high-binding affinities and the gerbil's prestin should have a lower affinity. Consistent with the functional assays, CIE values of prestin-SO42− in the zebrafish and chicken are −10.968 and −10.068, respectively, suggesting that their prestins can easily bind and transport SO42−. In contrast, the CIE value of prestin-SO42− in the gerbil is 5.19, indicating that this mammal cannot as easily transport SO42− (fig. 5A). Therefore, our modeled structure for prestin is congruent with the key aspects of transporting SO42−, and Cl−/HCO3−. The structure appears to provide a reliable means for determining the functional changes for permeability to SO42− and for inferring positively selected amino acid sites associated with changes in permeability in mammals. FIG. 5. View largeDownload slide (A) Distribution of CIE values between prestins and SO42− in representative vertebrate species. Asterisks indicate genes that provide functional evidence. (B) CIE values before and after artificial mutations of positively selected sites in the opossum and gerbil. FIG. 5. View largeDownload slide (A) Distribution of CIE values between prestins and SO42− in representative vertebrate species. Asterisks indicate genes that provide functional evidence. (B) CIE values before and after artificial mutations of positively selected sites in the opossum and gerbil. CIE values of prestin-SO42− can be used to evaluate when functional change might have occurred in the vertebrates. The values for the zebrafish, frog, lizard, chicken, platypus, and opossum are negative, ranging from −11.26 to −10.939. In contrast, placental mammals have positive CIE values, ranging from 4.36 to 5.951 (fig. 5A). Thus, the functional change in sulfate transport most likely has its origin in the MRCA of placental mammals and the positively selected amino acids in the MRCA of these animals might be associated with changes in the ability to transport sulfate. We tested whether the functional change depended on the structural conformation caused by positive selection or not by performing structural analyses after artificially exchanging positively selected sites identified on the branch of placental MRCA. First, we created two chimera prestins: chimera opossum prestin (CP1) and chimera gerbil prestin (CP2). CP1 was constructed by inserting the positively selected sites of placental mammals into the corresponding sites of opossum prestin. CP2 was established by exchanging the positively selected sites of gerbil prestin with the corresponding sites of opossum prestin. Second, the 3D structures of CP1 and CP2 were modeled and reliability of the modeling was validated. Both chimera prestins also hit the same model (1ots) with highly predicted accuracy and significant E values (supplementary table 2 and alignments in supplementary fig. 3, Supplementary Material online). When inferred positively selected sites in the prestin sequence of opossum were replaced by those from the gerbil prestin, the CIE value changed from −9.85 to 43.89 (fig. 5B), suggesting these positively selected amino acids affected the ability to transport sulfate. This might have owed to positively selected amino acids changing the local structural conformation of the pore region of the channel. Superposition of the opossum prestin and CP1 structures detected a mismatch between one region in opossum prestin (256–260) and the corresponding part in the CP1 (256–260). The side chain of LYS256 in the CP1 projected into the pore, and this might have blocked the entryway of the anion channel for the sulfate (fig. 6A). The ability of the gerbil prestin to transport sulfate was fully rescued by replacing positively selected sites with the amino acids of opossum prestin (fig. 5B). Superposition of the structures of the gerbil prestin and CP2 showed that the helix turn region (431–436) in the CP2 was similar to that of the opossum prestin (fig. 6B). It might have swung out of the channel's pore to allow sulfate penetration and, hence, obtained a negative CIE value (fig. 5B). FIG. 6. View largeDownload slide Structural superposition of the prestins. (A) Structural superposition of the opossum prestin (green color) and its chimera CP1 (tan color). Amino acids under positive selection are marked by a red ball and stick. The “▴” denotes positively selected amino acids in the chimera. The side chain of the LYS256 in opossum prestin (purple ball and stick) and the corresponding residue (also LYS256) in CP1 (blue ball and stick) are shown. (B) Structural superposition of the gerbil prestin (yellow color) and the chimera CP2 (cyan color). Only those amino acids under positive selection are marked by red ball and stick. The“▴” denotes positively selected amino acids in the CP2. The helix turn of gerbil prestin (431–436) and the CP2 (431–436) are colored by purple and blue, respectively. Note that the ▴ PRO240 is not labeled because this site in the structure prediction of CP2 was not available. FIG. 6. View largeDownload slide Structural superposition of the prestins. (A) Structural superposition of the opossum prestin (green color) and its chimera CP1 (tan color). Amino acids under positive selection are marked by a red ball and stick. The “▴” denotes positively selected amino acids in the chimera. The side chain of the LYS256 in opossum prestin (purple ball and stick) and the corresponding residue (also LYS256) in CP1 (blue ball and stick) are shown. (B) Structural superposition of the gerbil prestin (yellow color) and the chimera CP2 (cyan color). Only those amino acids under positive selection are marked by red ball and stick. The“▴” denotes positively selected amino acids in the CP2. The helix turn of gerbil prestin (431–436) and the CP2 (431–436) are colored by purple and blue, respectively. Note that the ▴ PRO240 is not labeled because this site in the structure prediction of CP2 was not available. Results from the cdocking analysis and artificial mutagenesis of positively selected sites on the ancestral branch of placental mammals support the above findings to some extent. However, exactly how the positively selected residues influence the changes of SO42− transport ability remains to be detailed. One possible scenario is that the positively selected residues are located on the pore of the channel, and they directly bind to SO42−and Cl−. Another scenario is that the positively selected sites surround the pore and indirectly bind the anions; this would change the conformation and, thus, result in defective sulfate transport. The absence of functional data precludes the unambiguous selection of one scenario. Regardless, our cdocking results more strongly support the latter possibility. The loss of SO42− transport ability seems to have happened in the MRCA of placental mammals. The positively selected amino acids in this MRCA appear to be related to changes in sulfate transport ability. The prestin of the MRCA of placental mammals, as well as that in the majority of, if not all, placental mammals, has a novel anion transport function that might further enhance the ability to detect high frequencies. Interestingly, our audiogram analysis supports this conjecture. Placental mammals generally possess a superior ability to detect higher frequencies (average upper hearing limit 61.9 kHz) than do marsupials (average upper hearing limit 38 kHz). Conclusion Our evolutionary analysis of prestin genes from 48 vertebrates provides evidence for multiple instances of positive selection and functional divergence events during vertebrate evolution. Prestin appears to have undergone positive selection during the emergence of tetrapods and, for the first time, adapted hearing for a terrestrial lifestyle. Moreover, our analyses indicate three independent adaptive events in the evolution of mammalian prestin genes. The first is predicted to have occurred in the MRCA of mammals and this possibly resulted from the gain of NLC and electromotility. The second adaptive event seems to have occurred in the MRCA of therian mammals and this might be related with the functional improvement of electromotility. The third appeared in the MRCA of placentals, which is associated with the ability of prestin to transport solutes, further enhances high-frequency detection. Our functional experiments support sequentially functional enhancements of prestin in monotremes, marsupials, and placentals, respectively. Combined with other results on bats and whales (Li et al. 2008, 2010; Liu, Cotton, et al. 2010; Liu, Rossiter, et al. 2010), these findings suggest that prestin genes underwent at least six positive selection events during the evolution of vertebrates. This discovery represents an unusually detailed understanding of how adaptation leads to functional diversity for the perception of high-frequency sound. Acknowledgments We thank Professor Peter Dallos (Northwestern University, Illinois, USA) for the expression clone of the gerbil prestin. We also thank Professor Chen Zhang for consultation regarding electrophysiological techniques. Valuable comments were made by the members of the Shi lab. This work was supported by a start-up fund of “Hundreds-Talent Program” from Chinese Academy of Sciences and by grants from Key Project from National Natural Science Foundation of China (30930015) to P.S. and by a Visiting Professorship for Senior International Scientists and Discovery Grant A3148 from the Natural Sciences and Engineering Research Council (Canada) to R.W.M. References Aitkin L.  The auditory neurobiology of marsupials: a review,  Hear Res. ,  1995, vol.  82 (pg.  257- 266) Google Scholar CrossRef Search ADS PubMed  Anisimova M,  Bielawski JP,  Yang Z.  Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution,  Mol Biol Evol. ,  2001, vol.  18 (pg.  1585- 1592) Google Scholar CrossRef Search ADS PubMed  Anisimova M,  Bielawski JP,  Yang Z.  Accuracy and power of Bayes prediction of amino acid sites under positive selection,  Mol Biol Evol. ,  2002, vol.  19 (pg.  950- 958) Google Scholar CrossRef Search ADS PubMed  Ashmore JF.  A fast motile response in guinea-pig outer hair cells: the cellular basis of the cochlear amplifier,  J Physiol. ,  1987, vol.  388 (pg.  323- 347) Google Scholar CrossRef Search ADS PubMed  Bai JP,  Surguchev A,  Montoya S,  Aronson PS,  Santos-Sacchi J,  Navaratnam D.  Prestin's anion transport and voltage-sensing capabilities are independent,  Biophys J. ,  2009, vol.  96 (pg.  3179- 3186) Google Scholar CrossRef Search ADS PubMed  Belyantseva IA,  Adler HJ,  Curi R,  Frolenkov GI,  Kachar B.  Expression and localization of prestin and the sugar transporter GLUT-5 during development of electromotility in cochlear outer hair cells,  J Neurosci. ,  2000, vol.  20 pg.  RC116  Google Scholar PubMed  Brownell WE,  Bader CR,  Bertrand D,  de Ribaupierre Y.  Evoked mechanical responses of isolated cochlear outer hair cells,  Science ,  1985, vol.  227 (pg.  194- 196) Google Scholar CrossRef Search ADS PubMed  Brownell WE,  Spector AA,  Raphael RM,  Popel AS.  Micro- and nanomechanics of the cochlear outer hair cell,  Annu Rev Biomed Eng. ,  2001, vol.  3 (pg.  169- 194) Google Scholar CrossRef Search ADS PubMed  Chenna R,  Sugawara H,  Koike T,  Lopez R,  Gibson TJ,  Higgins DG,  Thompson JD.  Multiple sequence alignment with the Clustal series of programs,  Nucleic Acids Res. ,  2003, vol.  31 (pg.  3497- 3500) Google Scholar CrossRef Search ADS PubMed  Clack TD.  Effect of signal duration on the auditory sensitivity of humans and monkeys (Macaca mulatta),  J Acoust Soc Am. ,  1966, vol.  40 (pg.  1140- 1146) Google Scholar CrossRef Search ADS PubMed  Coffin A,  Kelley M,  Manley GA,  Popper AN.  Manley GA,  Fay RR,  Popper AN.  Evolution of sensory hair cells,  Evolution of the vertebrate auditory system ,  2004 New York Springer-Verlag(pg.  55- 94) Dalland JI.  Hearing sensitivity in bats,  Science ,  1965, vol.  150 (pg.  1185- 1186) Google Scholar CrossRef Search ADS PubMed  Dallos P,  Fakler B.  Prestin, a new type of motor protein,  Nat Rev Mol Cell Biol. ,  2002, vol.  3 (pg.  104- 111) Google Scholar CrossRef Search ADS PubMed  Deak L,  Zheng J,  Orem A,  Du GG,  Aguinaga S,  Matsuda K,  Dallos P.  Effects of cyclic nucleotides on the function of prestin,  J Physiol. ,  2005, vol.  563 (pg.  483- 496) Google Scholar CrossRef Search ADS PubMed  Detro-Dassen S,  Schanzler M,  Lauks H,  Martin I,  zuBerstenhorst SM,  Nothmann D,  Torres-Salazar D,  Hidalgo P,  Schmalzing G,  Fahlke C.  Conserved dimeric subunit stoichiometry of SLC26 multifunctional anion exchangers,  J Biol Chem. ,  2008, vol.  283 (pg.  4177- 4188) Google Scholar CrossRef Search ADS PubMed  Dutzler R,  Campbell EB,  MacKinnon R.  Gating the selectivity filter in ClC chloride channels,  Science ,  2003, vol.  300 (pg.  108- 112) Google Scholar CrossRef Search ADS PubMed  Ehret G.  Critical bands and filter characteristics in the ear of the housemouse (Mus musculus),  Biol Cybern. ,  1976, vol.  24 (pg.  35- 42) Google Scholar CrossRef Search ADS PubMed  Fay R. ,  Hearing in vertebrates: a psychophysics databook ,  1988 Winnetka (IL) Hill-Fay Associates Fay R.  Webster DB,  Fay RR,  Popper AN.  Structure and function in sound discrimination among vertebrates,  The evolutionary biology of hearing ,  1996 New York Springer-Verlag(pg.  246- 247) Feng AS,  Narins PM,  Xu CH,  Lin WY,  Yu ZL,  Qiu Q,  Xu ZM,  Shen JX.  Ultrasonic communication in frogs,  Nature ,  2006, vol.  440 (pg.  333- 336) Google Scholar CrossRef Search ADS PubMed  Franchini LF,  Elgoyhen AB.  Adaptive evolution in mammalian proteins involved in cochlear outer hair cell electromotility,  Mol Phylogenet Evol. ,  2006, vol.  41 (pg.  622- 635) Google Scholar CrossRef Search ADS PubMed  Fritzsch B.  Webster DB,  Fay RR,  Popper AN.  The water-to-land transition: evolution of the tetrapod basilar papilla, middle ear, and auditory nuclei,  The evolutionary biology of hearing ,  1991 New York Springer-Verlagpg.  18  Gates GR,  Saunders JC,  Bock GR,  Aitkin LM,  Elliott MA.  Peripheral auditory function in the platypus, Ornithorhynchus anatinus,  J Acoust Soc Am. ,  1974, vol.  56 (pg.  152- 156) Google Scholar CrossRef Search ADS PubMed  Gaupp E.  Ontogenese und phylogenese des schalleitenden apparates bei den wirbeltieren,  Ergeb Anat Entwickl. ,  1898, vol.  8 (pg.  990- 1149) Gaupp E.  Die reichertsche theorie (Hammer, Amboss und Kieferfrage),  Arch Anat Physiol Anat Abt Suppl. ,  1913(pg.  1- 416) Gourevitch G.  Auditory masking in the rat,  J Acoust Soc Am. ,  1965, vol.  37 (pg.  439- 443) Google Scholar CrossRef Search ADS PubMed  Green S.  Auditory sensitivity and equal loudness in the squirrel monkey (Saimiris ciureus),  J Exp Anal Behav. ,  1975, vol.  23 (pg.  255- 264) Google Scholar CrossRef Search ADS PubMed  Heffner H,  Masterton B.  Hearing in glires: domestic rabbit, cotton rat, feral house mouse, and kangaroo rat,  J Acoust Soc Am. ,  1980, vol.  68 (pg.  1584- 1599) Google Scholar CrossRef Search ADS   Heffner HE.  Hearing in large and small dogs: absolute thresholds and size of the tympanic membrane,  Behav Neurosci. ,  1983, vol.  97 (pg.  310- 318) Google Scholar CrossRef Search ADS   Heffner HE,  Ravizza R,  Masterton B.  Hearing in primitive mammals: III, tree shrew (Tupaia glis),  J Aud Res. ,  1969a, vol.  9 (pg.  12- 18) Heffner HE,  Ravizza RJ,  Masterton B.  Hearing in primitive mammals: IV, bushbaby (Galago senegalensis),  J Aud Res. ,  1969b, vol.  9 (pg.  19- 23) Heffner R,  Heffner H,  Masterton B.  Behavioral measurements of absolute and frequency-difference thresholds in guinea pig,  J Acoust Soc Am. ,  1971, vol.  49 (pg.  1888- 1895) Google Scholar CrossRef Search ADS PubMed  Heffner RS,  Heffner HE.  Hearing in the elephant (Elephas maximus): absolute sensitivity, frequency discrimination, and sound localization,  J Comp Physiol Psychol. ,  1982, vol.  96 (pg.  926- 944) Google Scholar CrossRef Search ADS PubMed  Heffner RS,  Heffner HE.  Webster DB,  Fay RR,  Popper AN.  Evolution of sound localization in mammals,  The evolutionary biology of hearing ,  1991 New York Springer-Verlag(pg.  691- 711) Herman LM,  Arbeit WR.  Stimulus control and auditory discrimination learning sets in the bottlenose dolphin,  J Exp Anal Behav. ,  1973, vol.  19 (pg.  379- 394) Google Scholar CrossRef Search ADS PubMed  Hienz RD,  Turkkan JS,  Harris AH.  Pure tone thresholds in the yellow baboon (Papio cynocephalus),  Hear Res. ,  1982, vol.  8 (pg.  71- 75) Google Scholar CrossRef Search ADS PubMed  Huang G,  Santos-Sacchi J.  Mapping the distribution of the outer hair cell motility voltage sensor by electrical amputation,  Biophys J. ,  1993, vol.  65 (pg.  2228- 2236) Google Scholar CrossRef Search ADS PubMed  Kelley LA,  Sternberg MJ.  Protein structure prediction on the Web: a case study using the Phyre server,  Nat Protoc. ,  2009, vol.  4 (pg.  363- 371) Google Scholar CrossRef Search ADS PubMed  Kelly JB,  Kavanagh GL,  Dalton JC.  Hearing in the ferret (Mustela putorius): thresholds for pure tone detection,  Hear Res. ,  1986, vol.  24 (pg.  269- 275) Google Scholar CrossRef Search ADS PubMed  Kumano S,  Tan X,  He DZ,  Iida K,  Murakoshi M,  Wada H.  Mutation-induced reinforcement of prestin-expressing cells,  Biochem Biophys Res Commun. ,  2009, vol.  389 (pg.  569- 574) Google Scholar CrossRef Search ADS PubMed  Li G,  Wang J,  Rossiter SJ,  Jones G,  Cotton JA,  Zhang S.  The hearing gene Prestin reunites echolocating bats,  Proc Natl Acad Sci. U S A. ,  2008, vol.  105 (pg.  13959- 13964) Google Scholar CrossRef Search ADS PubMed  Li Y,  Liu Z,  Shi P,  Zhang J.  The hearing gene Prestin unites echolocating bats and whales,  Curr Biol. ,  2010, vol.  20 (pg.  R55- R56) Google Scholar CrossRef Search ADS PubMed  Liberman MC,  Gao J,  He DZ,  Wu X,  Jia S,  Zuo J.  Prestin is required for electromotility of the outer hair cell and for the cochlear amplifier,  Nature ,  2002, vol.  419 (pg.  300- 304) Google Scholar CrossRef Search ADS PubMed  Liu XZ,  Ouyang XM,  Xia XJ, et al.  (17 co-authors) Prestin, a cochlear motor protein, is defective in non-syndromic hearing loss,  Hum Mol Genet. ,  2003, vol.  12 (pg.  1155- 1162) Google Scholar CrossRef Search ADS PubMed  Liu Y,  Cotton JA,  Shen B,  Han X,  Rossiter SJ,  Zhang S.  Convergent sequence evolution between echolocating bats and dolphins,  Curr Biol. ,  2010, vol.  20 (pg.  R53- R54) Google Scholar CrossRef Search ADS PubMed  Liu Y,  Rossiter SJ,  Han X,  Cotton JA,  Zhang S.  Cetaceans on a molecular fast track to ultrasonic hearing,  Curr Biol. ,  2010, vol.  20 (pg.  1834- 1839) Google Scholar CrossRef Search ADS PubMed  Liu Z,  Li S,  Wang W,  Xu D,  Murphy RW,  Shi P.  Parallel evolution of KCNQ4 in echolocating bats,  PLoS One. ,  2011, vol.  6 pg.  e26618  Google Scholar CrossRef Search ADS PubMed  Lombard RE,  Bolt JR.  Fritzsch B,  Ryan M,  Wilczynski W,  Hetherington T,  Walkowiak W.  The evolution of the stapes in Paleozoic tetrapods,  The evolution of the amphibian auditory system ,  1988 New York Wiley and Sons(pg.  37- 67) Long GR.  Masked auditory thresholds from the bat, Rhinolophus ferrumequinum,  J Comp Physiol A. ,  1977, vol.  116 (pg.  247- 255) Google Scholar CrossRef Search ADS   Ludwig J,  Oliver D,  Frank G,  Klocker N,  Gummer AW,  Fakler B.  Reciprocal electromechanical properties of rat prestin: the motor molecule from rat outer hair cells,  Proc Natl Acad Sci U S A. ,  2001, vol.  98 (pg.  4178- 4183) Google Scholar CrossRef Search ADS PubMed  Manley GA. ,  Peripheral hearing mechanisms in reptiles and birds ,  1990 New York Springer-Verlag Manley GA.  Cochlear mechanisms from a phylogenetic viewpoint,  Proc Natl Acad Sci U S A. ,  2000, vol.  97 (pg.  11736- 11743) Google Scholar CrossRef Search ADS PubMed  Mann DA,  Higgs DM,  Tavolga WN,  Souza MJ,  Popper AN.  Ultrasound detection by clupeiform fishes,  J Acoust Soc Am. ,  2001, vol.  109 (pg.  3048- 3054) Google Scholar CrossRef Search ADS PubMed  Masterton B,  Heffner H,  Ravizza R.  The evolution of human hearing,  J Acoust Soc Am. ,  1969, vol.  45 (pg.  966- 985) Google Scholar CrossRef Search ADS PubMed  Matsuda K,  Zheng J,  Du GG,  Klocker N,  Madison LD,  Dallos P.  N-linked glycosylation sites of the motor protein prestin: effects on membrane targeting and electrophysiological function,  J Neurochem. ,  2004, vol.  89 (pg.  928- 938) Google Scholar CrossRef Search ADS PubMed  McGuire RM,  Liu H,  Pereira FA,  Raphael RM.  Cysteine mutagenesis reveals transmembrane residues associated with charge translocation in prestin,  J Biol Chem. ,  2010, vol.  285 (pg.  3103- 3113) Google Scholar CrossRef Search ADS PubMed  Mills DM,  Shepherd RK.  Distortion product otoacoustic emission and auditory brainstem responses in the echidna (Tachyglossus aculeatus),  J Assoc Res Otolaryngol. ,  2001, vol.  2 (pg.  130- 146) Google Scholar CrossRef Search ADS PubMed  Mitchell C,  Vernon J,  Herman P.  What does the lemur really hear?,  J Acoust Soc Am. ,  1971, vol.  50 (pg.  710- 711) Google Scholar CrossRef Search ADS PubMed  Mohl B.  Auditory sensitivity of the common seal in air and water,  J Aud Res. ,  1968, vol.  8 (pg.  27- 38) Murphy WJ,  Pevzner PA,  O'Brien SJ.  Mammalian phylogenomics comes of age,  Trends Genet. ,  2004, vol.  20 (pg.  631- 639) Google Scholar CrossRef Search ADS PubMed  Navaratnam D,  Bai JP,  Samaranayake H,  Santos-Sacchi J.  N-terminal-mediated homomultimerization of prestin, the outer hair cell motor protein,  Biophys J. ,  2005, vol.  89 (pg.  3345- 3352) Google Scholar CrossRef Search ADS PubMed  Nei M,  Kumar S. ,  Molecular evolution and phylogenetics ,  2000 New York Oxford University Press Nienhuys TG,  Clark GM.  Critical bands following the selective destruction of cochlear inner and outer hair cells,  Acta Oto-laryngol. ,  1979, vol.  88 (pg.  350- 358) Google Scholar CrossRef Search ADS   Okoruwa OE,  Weston MD,  Sanjeevi DC,  Millemon AR,  Fritzsch B,  Hallworth R,  Beisel KW.  Evolutionary insights into the unique electromotility motor of mammalian outer hair cells,  Evol Dev. ,  2008, vol.  10 (pg.  300- 315) Google Scholar CrossRef Search ADS PubMed  Oliver D,  He DZ,  Klocker N,  Ludwig J,  Schulte U,  Waldegger S,  Ruppersberg JP,  Dallos P,  Fakler B.  Intracellular anions as the voltage sensor of prestin, the outer hair cell motor protein,  Science ,  2001, vol.  292 (pg.  2340- 2343) Google Scholar CrossRef Search ADS PubMed  Pasqualetto E,  Aiello R,  Gesiot L,  Bonetto G,  Bellanda M,  Battistutta R.  Structure of the cytosolic portion of the motor protein prestin and functional role of the STAS domain in SLC26/SulP anion transporters,  J Mol Biol. ,  2010, vol.  400 (pg.  448- 462) Google Scholar CrossRef Search ADS PubMed  Popper AN.  Hair cell heterogeneity and ultrasonic hearing: recent advances in understanding fish hearing,  Philos Trans R Soc Lond B Biol Sci. ,  2000, vol.  355 (pg.  1277- 1280) Google Scholar CrossRef Search ADS PubMed  Popper AN,  Fay RR.  Evolution of the ear and hearing: issues and questions,  Brain Behav Evol. ,  1997, vol.  50 (pg.  213- 221) Google Scholar CrossRef Search ADS PubMed  Rajagopalan L,  Patel N,  Madabushi S, et al.  (12 co-authors) Essential helix interactions in the anion transporter domain of prestin revealed by evolutionary trace analysis,  J Neurosci. ,  2006, vol.  26 (pg.  12727- 12734) Google Scholar CrossRef Search ADS PubMed  Reimer K.  Ontogeny of hearing in the marsupial, Monodelphis domestica, as revealed by brainstem auditory evoked potentials,  Hear Res. ,  1995, vol.  92 (pg.  143- 150) Google Scholar CrossRef Search ADS PubMed  Retzius G. ,  Das gehororgan der wirbeltiere: I Das gehororgan der fische und amphibien ,  1881 Stockholm (Sweden) Samson and Wallin Retzius G. ,  Das gehororgan der wirbeltiere: II. Das gehororgan der amphibien ,  1884 Stockholm (Sweden) Samson and Wallin Ryan A.  Hearing sensitivity of the Mongolian gerbil, Meriones unguiculatis,  J Acoust Soc Am. ,  1976, vol.  59 (pg.  1222- 1226) Google Scholar CrossRef Search ADS PubMed  Ryan MJ,  Tuttle MD,  Barclay MR.  Behavioral responses of the frog-eating bat, Trachops cirrhosus, to sonic frequencies,  J Comp Physiol A. ,  1983, vol.  150 (pg.  413- 418) Google Scholar CrossRef Search ADS   Saitou N,  Nei M.  The neighbor-joining method: a new method for reconstructing phylogenetic trees,  Mol Biol Evol. ,  1987, vol.  4 (pg.  406- 425) Google Scholar PubMed  Santos-Sacchi J.  Reversible inhibition of voltage-dependent outer hair cell motility and capacitance,  J Neurosci. ,  1991, vol.  11 (pg.  3096- 3110) Google Scholar PubMed  Schaechinger TJ,  Oliver D.  Nonmammalian orthologs of prestin (SLC26A5) are electrogenic divalent/chloride anion exchangers,  Proc Natl Acad Sci U S A. ,  2007, vol.  104 (pg.  7693- 7698) Google Scholar CrossRef Search ADS PubMed  Shi P,  Zhang J,  Yang H,  Zhang YP.  Adaptive diversification of bitter taste receptor genes in mammalian evolution,  Mol Biol Evol. ,  2003, vol.  20 (pg.  805- 814) Google Scholar CrossRef Search ADS PubMed  Shindyalov IN,  Bourne PE.  Protein structure alignment by incremental combinatorial extension (CE) of the optimal path,  Protein Eng. ,  1998, vol.  11 (pg.  739- 747) Google Scholar CrossRef Search ADS PubMed  Sindic A,  Chang MH,  Mount DB,  Romero MF.  Renal physiology of SLC26 anion exchangers,  Curr Opin Nephrol Hypertens. ,  2007, vol.  16 (pg.  484- 490) Google Scholar CrossRef Search ADS PubMed  Tamura K,  Peterson D,  Peterson N,  Stecher G,  Nei M,  Kumar S.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,  Mol Biol Evol. ,  2011, vol.  28 (pg.  2731- 2739) Google Scholar CrossRef Search ADS PubMed  Tan X,  Pecka JL,  Tang J,  Okoruwa OE,  Zhang Q,  Beisel KW,  He DZ.  From zebrafish to mammal: functional evolution of prestin, the motor protein of cochlear outer hair cells,  J Neurophysiol. ,  2011, vol.  105 (pg.  36- 44) Google Scholar CrossRef Search ADS PubMed  Tanaka T,  Nei M.  Positive Darwinian selection observed at the variable-region genes of immunoglobulins,  Mol Biol Evol. ,  1989, vol.  6 (pg.  447- 459) Google Scholar PubMed  Teeling EC,  Springer MS,  Madsen O,  Bates P,  O'Brien SJ,  Murphy WJ.  A molecular phylogeny for bats illuminates biogeography and the fossil record,  Science ,  2005, vol.  307 (pg.  580- 584) Google Scholar CrossRef Search ADS PubMed  Thomson KS.  The evolution of the tetrapod middle ear in the rhipidistian-amphibian transition,  Am Zool. ,  1966, vol.  6 (pg.  379- 397) Google Scholar CrossRef Search ADS PubMed  Tung CH,  Yang JM.  fastSCOP: a fast web server for recognizing protein structural domains and SCOP superfamilies,  Nucleic Acids Res. ,  2007, vol.  35 (pg.  W438- W443) Google Scholar CrossRef Search ADS PubMed  Warren WC,  Hillier LW,  Marshall Graves JA, et al.  (106 co-authors) Genome analysis of the platypus reveals unique signatures of evolution,  Nature ,  2008, vol.  453 (pg.  175- 183) Google Scholar CrossRef Search ADS PubMed  Werner G. ,  Das labyrinth der wirbeltiere ,  1960 Jena (Germany) Fischer Verlag  p.309. Wier CC,  Jesteadt W,  Green DM.  Frequency discrimination as a function of frequency and sensation level,  J Acoust Soc Am. ,  1977, vol.  61 (pg.  178- 184) Google Scholar CrossRef Search ADS PubMed  Wu G,  Robertson DH,  Brooks CL3rd,  Vieth M.  Detailed analysis of grid-based molecular docking: a case study of CDOCKER-A CHARMm-based MD docking algorithm,  J Comput Chem. ,  2003, vol.  24 (pg.  1549- 1562) Google Scholar CrossRef Search ADS PubMed  Wu X,  Gao J,  Guo Y,  Zuo J.  Hearing threshold elevation precedes hair-cell loss in prestin knockout mice,  Brain Res Mol Brain Res. ,  2004, vol.  126 (pg.  30- 37) Google Scholar CrossRef Search ADS PubMed  Yang H,  Shi P,  Zhang YP,  Zhang J.  Composition and evolution of the V2r vomeronasal receptor gene repertoire in mice and rats,  Genomics ,  2005, vol.  86 (pg.  306- 315) Google Scholar CrossRef Search ADS PubMed  Yang Z,  Nielsen R.  Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages,  Mol Biol Evol. ,  2002, vol.  19 (pg.  908- 917) Google Scholar CrossRef Search ADS PubMed  Yang Z,  Wong WS,  Nielsen R.  Bayes empirical Bayes inference of amino acid sites under positive selection,  Mol Biol Evol. ,  2005, vol.  22 (pg.  1107- 1118) Google Scholar CrossRef Search ADS PubMed  Zhang J.  Frequent false detection of positive selection by the likelihood method with branch-site models,  Mol Biol Evol. ,  2004, vol.  21 (pg.  1332- 1339) Google Scholar CrossRef Search ADS PubMed  Zhang J,  Nielsen R,  Yang Z.  Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level,  Mol Biol Evol. ,  2005, vol.  22 (pg.  2472- 2479) Google Scholar CrossRef Search ADS PubMed  Zhang J,  Rosenberg HF,  Nei M.  Positive Darwinian selection after gene duplication in primate ribonuclease genes,  Proc Natl Acad Sci U S A. ,  1998, vol.  95 (pg.  3708- 3713) Google Scholar CrossRef Search ADS PubMed  Zheng J,  Du GG,  Matsuda K,  Orem A,  Aguinaga S,  Deak L,  Navarrete E,  Madison LD,  Dallos P.  The C-terminus of prestin influences nonlinear capacitance and plasma membrane targeting,  J Cell Sci. ,  2005, vol.  118 (pg.  2987- 2996) Google Scholar CrossRef Search ADS PubMed  Zheng J,  Shen W,  He DZ,  Long KB,  Madison LD,  Dallos P.  Prestin is the motor protein of cochlear outer hair cells,  Nature ,  2000, vol.  405 (pg.  149- 155) Google Scholar CrossRef Search ADS PubMed  © The Author 2012. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com TI - Hearing Aid for Vertebrates via Multiple Episodic Adaptive Events on Prestin Genes JF - Molecular Biology and Evolution DO - 10.1093/molbev/mss087 DA - 2012-03-13 UR - https://www.deepdyve.com/lp/oxford-university-press/hearing-aid-for-vertebrates-via-multiple-episodic-adaptive-events-on-AYHZGTlp0F SP - 2187 EP - 2198 VL - 29 IS - 9 DP - DeepDyve ER -