Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Extreme morphological and ecological homoplasy in tropical salamanders

Fossorial salamanders typically have elongate and attenuated heads and bodies, diminutive limbs, hands and feet, and extremely elongate tails. Batrachoseps from California, Lineatriton from eastern México, and Oedipina from southern México to Ecuador, all members of the family Plethodontidae, tribe Bolitoglossini, resemble one another in external morphology, which has evolved independently. Whereas Oedipina and Batrachoseps are elongate because there are more trunk vertebrae, a widespread homoplasy (parallelism) in salamanders, the genus Lineatriton is unique in having evolved convergently by an alternate “giraffe-neck” developmental program. Lineatriton has the same number of trunk vertebrae as related, nonelongated taxa, but individual trunk vertebrae are elongated. A robust phylogenetic hypothesis, based on sequences of three mtDNA genes, finds Lineatriton to be deeply nested within a clade characterized by generalized ecology and morphology. Lineatriton lineolus, the only currently recognized taxon in the genus, shows unanticipated genetic diversity. Surprisingly, geographically separated populations of L. lineolus are not monophyletic, but are sister taxa of different species of the morphologically generalized genus Pseudoeurycea. Lineatriton, long thought to be a unique monospecific lineage, is polyphyletic. Accordingly, the specialized morphology of Lineatriton displays homoplasy at two hierarchical levels: (i) with respect to other elongate lineages in the family (convergence), and (ii) within what is currently recognized as a single taxon (parallelism). These evolutionary events are of adaptive significance because to invade the lowland tropics salamanders must be either arboreal or fossorial; the repeated evolution of elongation and attenuation has led to multiple lowland invasions.

Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

How the owl resolves auditory coding ambiguity

The barn owl (Tyto alba) uses interaural time difference (ITD) cues to localize sounds in the horizontal plane. Low-order binaural auditory neurons with sharp frequency tuning act as narrow-band coincidence detectors; such neurons respond equally well to sounds with a particular ITD and its phase equivalents and are said to be phase ambiguous. Higher-order neurons with broad frequency tuning are unambiguously selective for single ITDs in response to broad-band sounds and show little or no response to phase equivalents. Selectivity for single ITDs is thought to arise from the convergence of parallel, narrow-band frequency channels that originate in the cochlea. ITD tuning to variable bandwidth stimuli was measured in higher-order neurons of the owl’s inferior colliculus to examine the rules that govern the relationship between frequency channel convergence and the resolution of phase ambiguity. Ambiguity decreased as stimulus bandwidth increased, reaching a minimum at 2–3 kHz. Two independent mechanisms appear to contribute to the elimination of ambiguity: one suppressive and one facilitative. The integration of information carried by parallel, distributed processing channels is a common theme of sensory processing that spans both modality and species boundaries. The principles underlying the resolution of phase ambiguity and frequency channel convergence in the owl may have implications for other sensory systems, such as electrolocation in electric fish and the computation of binocular disparity in the avian and mammalian visual systems.

Functional organization of the hippocampal memory system

In humans declarative or explicit memory is supported by the hippocampus and related structures of the medial temporal lobe working in concert with the cerebral cortex. This paper reviews our progress in developing an animal model for studies of cortical–hippocampal interactions in memory processing. Our findings support the view that the cortex maintains various forms of memory representation and that hippocampal structures extend the persistence and mediate the organization of these codings. Specifically, the parahippocampal region, through direct and reciprocal interconnections with the cortex, is sufficient to support the convergence and extended persistence of cortical codings. The hippocampus itself is critical to the organization cortical representations in terms of relationships among items in memory and in the flexible memory expression that is the hallmark of declarative memory.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.

The chromosomal distribution of mouse odorant receptor genes.

Odorant receptors (ORs) on nasal olfactory sensory neurons are encoded by a large multigene family. Each member of the family is expressed in a small percentage of neurons that are confined to one of several spatial zones in the nose but are randomly distributed throughout that zone. This pattern of expression suggests that when the sensory neuron selects which OR gene to express it may be confined to a particular zonal gene set of several hundred OR genes but select from among the members of that set via a stochastic mechanism. Both locus-dependent and locus-independent models of OR gene choice have been proposed. To investigate the feasibility of these models, we determined the chromosomal locations of 21 OR genes expressed in four different spatial zones. We found that OR genes are clustered within multiple loci that are broadly distributed in the genome. These loci lie within paralogous chromosomal regions that appear to have arisen by duplications of large chromosomal domains followed by extensive gene duplication and divergence. Our studies show that OR genes expressed in the same zone map to numerous loci; moreover, a single locus can contain genes expressed in different zones. These findings raise the possibility that OR gene choice may be locus-independent or involve consecutive stochastic choices.

Determining divergence times with a protein clock: Update and reevaluation

A recent study of the divergence times of the major groups of organisms as gauged by amino acid sequence comparison has been expanded and the data have been reanalyzed with a distance measure that corrects for both constraints on amino acid interchange and variation in substitution rate at different sites. Beyond that, the availability of complete genome sequences for several eubacteria and an archaebacterium has had a great impact on the interpretation of certain aspects of the data. Thus, the majority of the archaebacterial sequences are not consistent with currently accepted views of the Tree of Life which cluster the archaebacteria with eukaryotes. Instead, they are either outliers or mixed in with eubacterial orthologs. The simplest resolution of the problem is to postulate that many of these sequences were carried into eukaryotes by early eubacterial endosymbionts about 2 billion years ago, only very shortly after or even coincident with the divergence of eukaryotes and archaebacteria. The strong resemblances of these same enzymes among the major eubacterial groups suggest that the cyanobacteria and Gram-positive and Gram-negative eubacteria also diverged at about this same time, whereas the much greater differences between archaebacterial and eubacterial sequences indicate these two groups may have diverged between 3 and 4 billion years ago.

Convergence of two repressors through heterodimer formation of androgen receptor and testicular orphan receptor-4: A unique signaling pathway in the steroid receptor superfamily

The androgen receptor (AR) binds to androgen response elements and regulates target genes via a mechanism involving coregulators. Here we demonstrate that the AR can interact with the testicular orphan receptor-4 (TR4) and function as a repressor to down-regulate the TR4 target genes by preventing the TR4 binding to its target DNA. Interestingly, the heterodimerization of AR and TR4 also allows TR4 to repress AR target gene expression. Simultaneous exposure to both receptors therefore could result in bidirectional suppression of their target genes. Together, these data demonstrate that the coupling of two different receptors, through the heterodimerization of AR and TR4, is a unique signaling pathway in the steroid receptor superfamily, which may facilitate further understanding of the complicated androgen action in prostate cancer or libido.

An ancestral MADS-box gene duplication occurred before the divergence of plants and animals

Changes in genes encoding transcriptional regulators can alter development and are important components of the molecular mechanisms of morphological evolution. MADS-box genes encode transcriptional regulators of diverse and important biological functions. In plants, MADS-box genes regulate flower, fruit, leaf, and root development. Recent sequencing efforts in Arabidopsis have allowed a nearly complete sampling of the MADS-box gene family from a single plant, something that was lacking in previous phylogenetic studies. To test the long-suspected parallel between the evolution of the MADS-box gene family and the evolution of plant form, a polarized gene phylogeny is necessary. Here we suggest that a gene duplication ancestral to the divergence of plants and animals gave rise to two main lineages of MADS-box genes: TypeI and TypeII. We locate the root of the eukaryotic MADS-box gene family between these two lineages. A novel monophyletic group of plant MADS domains (AGL34 like) seems to be more closely related to previously identified animal SRF-like MADS domains to form TypeI lineage. Most other plant sequences form a clear monophyletic group with animal MEF2-like domains to form TypeII lineage. Only plant TypeII members have a K domain that is downstream of the MADS domain in most plant members previously identified. This suggests that the K domain evolved after the duplication that gave rise to the two lineages. Finally, a group of intermediate plant sequences could be the result of recombination events. These analyses may guide the search for MADS-box sequences in basal eukaryotes and the phylogenetic placement of new genes from other plant species.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Activity-dependent CREB phosphorylation: Convergence of a fast, sensitive calmodulin kinase pathway and a slow, less sensitive mitogen-activated protein kinase pathway

The cAMP-responsive element binding protein (CREB), a key regulator of gene expression, is activated by phosphorylation on Ser-133. Several different protein kinases possess the capability of driving this phosphorylation, making it a point of potential convergence for multiple intracellular signaling cascades. Previous work in neurons has indicated that physiologic synaptic stimulation recruits a fast calmodulin kinase IV (CaMKIV)-dependent pathway that dominates early signaling to CREB. Here we show in hippocampal neurons that the fast, CaMK-dependent pathway can be followed by a slower pathway that depends on Ras/mitogen-activated protein kinase (MAPK), along with CaMK. This pathway was blocked by dominant-negative Ras and was specifically recruited by depolarizations that produced strong intracellular Ca2+ transients. When both pathways were recruited, phosphorylated CREB (pCREB) formation was overwhelmingly dominated by the CaMK pathway between 0 and 10 min, and by the MAPK pathway at 60 min, whereas the two pathways acted in concert at 30 min. The Ca2+ signals that produced only rapid CaMK signaling to pCREB or both rapid CaMK and slow MAPK signaling deviated significantly for only ≈1 min, yet their differential impact on pCREB extended over a much longer period, between 20 and 60 min and beyond, which is of likely significance for gene expression. The CaMK-dependent MAPK pathway may inform the nucleus about stimulus amplitude. In contrast, the CaMKIV pathway may be well suited to conveying information on the precise timing of localized synaptic stimuli, befitting its greater speed and sensitivity, whereas the previously described calcineurin pathway may carry information about stimulus duration.

Duplication and functional divergence in the chalcone synthase gene family of Asteraceae: evolution with substrate change and catalytic simplification.

Plant-specific polyketide synthase genes constitute a gene superfamily, including universal chalcone synthase [CHS; malonyl-CoA:4-coumaroyl-CoA malonyltransferase (cyclizing) (EC 2.3.1.74)] genes, sporadically distributed stilbene synthase (SS) genes, and atypical, as-yet-uncharacterized CHS-like genes. We have recently isolated from Gerbera hybrida (Asteraceae) an unusual CHS-like gene, GCHS2, which codes for an enzyme with structural and enzymatic properties as well as ontogenetic distribution distinct from both CHS and SS. Here, we show that the GCHS2-like function is encoded in the Gerbera genome by a family of at least three transcriptionally active genes. Conservation within the GCHS2 family was exploited with selective PCR to study the occurrence of GCHS2-like genes in other Asteraceae. Parsimony analysis of the amplified sequences together with CHS-like genes isolated from other taxa of angiosperm subclass Asteridae suggests that GCHS2 has evolved from CHS via a gene duplication event that occurred before the diversification of the Asteraceae. Enzyme activity analysis of proteins produced in vitro indicates that the GCHS2 reaction is a non-SS variant of the CHS reaction, with both different substrate specificity (to benzoyl-CoA) and a truncated catalytic profile. Together with the recent results of Durbin et al. [Durbin, M. L., Learn, G. H., Jr., Huttley, G. A. & Clegg, M. T. (1995) Proc. Natl. Acad. Sci. USA 92, 3338-3342], our study confirms a gene duplication-based model that explains how various related functions have arisen from CHS during plant evolution.

Extraordinary divergence and positive Darwinian selection in a fusagenic protein coating the acrosomal process of abalone spermatozoa.

During fertilization in marine invertebrates, fusion between sperm and egg cell membranes occurs at the tip of the sperm acrosomal process. In abalone sperm the acrosomal process is coated with an 18-kDa protein. In situ, this protein has no effect on the egg vitelline envelope, but in vitro it is a potent fusagen of liposomes. Thus, the 18-kDa protein may mediate membrane fusion between the gametes, a step in gamete recognition known to restrict heterospecific fertilization in other species. The cDNA and deduced amino acid sequences of the 18-kDa protein were determined for five species of California abalone. The deduced amino acid sequences exhibit extraordinary divergence; the percent identity varies from 27% to 87%. Analysis of nucleotide substitution shows extremely high frequencies of amino acid-altering substitution compared to silent substitution, demonstrating that positive Darwinian selection promotes the divergence of this protein. However, amino acid replacement is conservative with respect to size and polarity of residue. The data support the developing idea that in free-spawning marine invertebrates, the proteins mediating fertilization may be subjected to intense, and as yet unknown, selective forces. The extraordinary divergence of fertilization proteins may be related to the establishment of barriers to heterospecific fertilization.