Power analysis of tests for loss of selective constraint in cave crayfish and nonphotosynthetic plant lineages

Loss of selective constraint on a gene may be expected following changes in the environment or life history that render its function unnecessary. The long-term persistence of protein-coding genes after the loss of known functional necessity can occur by chance or because of selective maintenance of an unknown gene function. The selective maintenance of an alternative gene function is not demonstrated by the failure of statistical tests to reject the hypothesis that there has been no change in the degree of constraint on the evolution of coding genes. Maintenance may be inferred, however, when power analyses of such tests demonstrate that there has been a sufficient number of nucleotide substitutions to detect the loss of selective constraint. Here, we describe a power analysis for tests of loss of constraint on protein-coding genes. The power analysis was applied to loss-of-constraint tests for opsin gene evolution in cave-dwelling crayfish and rbcL evolution in nonphotosynthetic parasitic plants. The power of previously applied tests for loss of constraint on cave crayfish opsin genes was insufficient to distinguish between chance retention and selective maintenance of opsin genes. However, the power of codon-based likelihood ratio tests for change in d(N)/d(S) (=omega) (nonsynonymous to synonymous change) did have sufficient power to detect a loss of constraint on rbcL associated with a loss of photosynthesis in most examples but failed to detect such a change in three independent lineages. We conclude that rbcL has been selectively maintained in these holoparasitic plant lineages. This conclusion suggests that either these taxa are photosynthetic for at least a part of their life or rbcL may have an unknown function in these plants unrelated to photosynthesis.     

Dynamics of neuroreceptor function in a chemical synapse

A model for the dynamics of functioning of neuroreceptors in chemical synapses is proposed. It is shown that the functioning of a neuroreceptor by the action of mediator is determined by correlation of the mechanical elasticity of the receptor with its electrical polarizability and the presence of spiralized fragments. It is shown in terms of the quantum-mechanical approach that, as the charged groups of the mediator come close to the C-end of the neuroreceptor, a shift of potential surface occurs, which corresponds to the deformation of the receptor. The microscopic approach enables one to describe more exactly the dynamics of functioning of neuroreceptors than the macroscopic approach; however, its implementation requires the use of powerful computers. If only qualitative assessment is needed, the macroscopic approach would suffice. It is assumed that in terms of the model, the saltatory (jump-like) conduction of the spike along the axon can be explained.     

Evolution and phylogenetic utility of the period gene in Lepidoptera

Evolution and phylogenetic utility of the period gene are explored through sequence analysis of a relatively conserved 909-bp fragment in 26 lepidopteran species. Taxa range from tribes to superfamilies, primarily within the putative clade Macrolepidotera plus near outgroups, and include both strongly established and problematic groupings. Their divergence dates probably range from the late Cretaceous through much of the Tertiary. Comparisons within the same set of closely related species show that amino acid substitutions in period occur 4.9 and 44 times as frequently as they do in two other nuclear genes--dopa decarboxylase and elongation factor-1 alpha, respectively. In contrast, rates of observed synonymous substitution are within 60% of each other for these three genes. Synonymous changes in period approach saturation by the family level, whereas nonsynonymous and amino acid divergences across the Macrolepidoptera are less than half the maximal values reported for this gene. Phylogenetic analyses of period strongly supported groupings at the family level and below. In contrast to previous analyses at this level with other nuclear genes, much of the information lies in nonsynonymous change. Relationships up to the superfamily level were recovered with decreasing effectiveness, and little, if any, signal was apparent regarding relationships among superfamilies. This could reflect rapid radiation of the superfamilies, however, rather than saturation in the period locus; thus, period, in combination with other genes, remains a plausible candidate for approaching the difficult problems of lepidopteran family and superfamily relationships.      

Charge scan reveals an extended region at the intracellular end of the GABA receptor pore that can influence ion selectivity

Selective permeability is a fundamental property of ion channels. The Cys-loop receptor superfamily is composed of both excitatory (ACh, 5-HT) and inhibitory (GABA, glycine) neurotransmitter-operated ion channels. In the GABA receptor, it has been previously shown that the charge selectivity of the integral pore can be altered by a single mutation near the intracellular end of the second transmembrane-spanning domain (TM2). We have extended these findings and now show that charge selectivity of the anionic rho1 GABA receptor can be influenced by the introduction of glutamates, one at a time, over an 8-amino acid stretch (-2' to 5') in the proposed intracellular end of TM2 and the TM1-TM2 intracellular linker. Depending on the position, glutamate substitutions in this region produced sodium to chloride permeability ratios (P(Na)+(/Cl)-) varying from 0.64 to 3.4 (wild type P(Na)+(/Cl)- = 0). In addition to providing insight into the mechanism of ion selectivity, this functional evidence supports a model proposed for the homologous nicotinic acetylcholine receptor in which regions of the protein, in addition to TM2, form the ion pathway.     

A conserved gene cluster as a putative functional unit in insect innate immunity

The Nimrod gene superfamily is an important component of the innate immune response. The majority of its member genes are located in close proximity within the Drosophila melanogaster genome and they lie in a larger conserved cluster (“Nimrod cluster”), made up of non-related groups (families, superfamilies) of genes. This cluster has been a part of the Arthropod genomes for about 300-350 million years. The available data suggest that the Nimrod cluster is a functional module of the insect innate immune response.     

Mining of assembled expressed sequence tag (EST) data for protein families: application to the G protein-coupled receptor superfamily

The availability of large expressed sequence tag (EST) databases has led to a revolution in the way new genes are identified. Mining of these databases using known protein sequences as queries is a powerful technique for discovering orthologous and paralogous genes. The scientist is often confronted, however, by an enormous amount of search output owing to the inherent redundancy of EST data. In addition, high search sensitivity often cannot be achieved using only a single member of a protein superfamily as a query. In this paper a technique for addressing both of these issues is described. Assembled EST databases are queried with every member of a protein superfamily, the results are integrated and false positives are pruned from the set. The result is a set of assemblies enriched in members of the protein superfamily under consideration. The technique is applied to the G protein-coupled receptor (GPCR) superfamily in the construction of a GPCR Resource. A novel full-length human GPCR identified from the GPCR Resource is presented, illustrating the utility of the method.     

The origin and evolution of protein superfamilies.

The organization of proteins into superfamilies based primarily on their sequences is introduced: examples are given of the methods used to cluster the related sequences and to elucidate the evolutionary history of the corresponding genes within each superfamily. Within the framework of this organization, the amount of sequence information currently and potentially available in all living forms can be discussed. The 116 superfamilies already sampled reflect possibly 10% of the total number. There are related proteins from many species in all of these superfamilies, suggesting that the origin of a new superfamily is rare indeed. The proteins so far sequenced are so rigorously conserved by the evolutionary process that we would expect to recognize as related descendants of any protein found in the ancestral vertebrate. The evolutionary history of the thyrotropin-gonadotropin beta chain superfamily is discussed in detail as an example. Some proteins are so constrained in structure that related forms can be recognized in prokaryotes and eukaryotes. Evolution in these superfamilies can be traced back close to the origin of life itself. From the evolutionary tree of the c-type cytochromes the identity of the prokaryote types involved in the symbiotic origin of mitochondria and chloroplasts begins to emerge.     

Maximum-likelihood analysis of molecular adaptation in abalone sperm lysin reveals variable selective pressures among lineages and sites

Maximum-likelihood models of codon substitution were used to analyze sperm lysin genes of 25 abalone (HALIOTIS:) species to identify lineages and amino acid sites under diversifying selection. The models used the nonsynonymous/synonymous rate ratio (omega = d(N)/d(S)) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio tests suggested significant variation in selective pressure among lineages. The variable selective pressure provided an explanation for the previous observation that the omega ratio is >1 in comparisons of closely related species and <1 in comparisons of distantly related species. Computer simulations demonstrated that saturation of nonsynonymous substitutions and constraint on lysin structure were unlikely to account for the observed pattern. Lineages linking closely related sympatric species appeared to be under diversifying selection, while lineages separating distantly related species from different geographic locations were associated with low evolutionary rates. The selective pressure indicated by the omega ratio was found to vary greatly among amino acid sites in lysin. Sites under potential diversifying selection were identified. Ancestral lysins were inferred to trace the route of evolution at individual sites and to provide lysin sequences for future laboratory studies.     

A network of conserved intramolecular contacts defines the off-state of the transmembrane switch mechanism in a seven-transmembrane receptor

Activation of the rhodopsin-like 7-transmembrane (7-TM) receptors requires switching interhelical constraints that stabilize the inactive state to a new set of contacts in the activated state, which binds the cognate G-protein. The free energy to drive this is provided by agonist binding, which has higher affinity to the active than to the inactive conformation. We have sought specific interhelical constraint contacts, using the M(1) muscarinic acetylcholine receptor as a model. Histidine substitutions of particular groups of amino acids, in transmembrane domains 3, 6, and 7, created high-affinity Zn(2+) binding sites, demonstrating the close proximity of their side chains in the inactive state. Alanine point substitutions have shown the effect of weakening the individual intramolecular contacts. In each case, the acetylcholine affinity was increased, implying promotion of the activated state. These amino acids are highly conserved throughout the 7-TM receptor superfamily. We propose that they form an important part of a network of conserved interhelical contacts that defines the off-state of a general transmembrane switch mechanism.     

An aminopeptidase,ARTS-1, is required for interleukin-6 receptor shedding

Aminopeptidase regulator of TNFR1 shedding (ARTS-1) binds to the type I tumor necrosis factor receptor (TNFR1) and promotes receptor shedding. Because hydroxamic acid-based metalloprotease inhibitors prevent shedding of both TNFR1 and the interleukin-6 receptor (IL-6Ralpha), we hypothesized that ARTS-1 might also regulate shedding of IL-6Ralpha, a member of the type I cytokine receptor superfamily that is structurally different from TNFR1. Reciprocal co-immunoprecipitation experiments identified that membrane-associated ARTS-1 directly binds to a 55-kDa IL-6Ralpha, a size consistent with soluble IL-6Ralpha generated by ectodomain cleavage of the membrane-bound receptor. Furthermore, ARTS-1 promoted IL-6Ralpha shedding, as demonstrated by a direct correlation between increased membrane-associated ARTS-1 protein, increased IL-6Ralpha shedding, and decreased membrane-associated IL-6Ralpha in cell lines overexpressing ARTS-1. The absence of basal IL-6Ralpha shedding from arts-1 knock-out cells identified that ARTS-1 was required for constitutive IL-6Ralpha shedding. Furthermore, the mechanism of constitutive IL-6Ralpha shedding requires ARTS-1 catalytic activity. Thus, ARTS-1 promotes the shedding of two cytokine receptor superfamilies, the type I cytokine receptor superfamily (IL-6Ralpha) and the TNF receptor superfamily (TNFR1). We propose that ARTS-1 is a multifunctional aminopeptidase that may modulate inflammatory events by promoting IL-6Ralpha and TNFR1 shedding.     

The vestigial olfactory receptor subgenome of odontocete whales: phylogenetic congruence between gene-tree reconciliation and supermatrix methods

The macroevolutionary transition of whales (cetaceans) from a terrestrial quadruped to an obligate aquatic form involved major changes in sensory abilities. Compared to terrestrial mammals, the olfactory system of baleen whales is dramatically reduced, and in toothed whales is completely absent. We sampled the olfactory receptor (OR) subgenomes of eight cetacean species from four families. A multigene tree of 115 newly characterized OR sequences from these eight species and published data for Bos taurus revealed a diverse array of class II OR paralogues in Cetacea. Evolution of the OR gene superfamily in toothed whales (Odontoceti) featured a multitude of independent pseudogenization events, supporting anatomical evidence that odontocetes have lost their olfactory sense. We explored the phylogenetic utility of OR pseudogenes in Cetacea, concentrating on delphinids (oceanic dolphins), the product of a rapid evolutionary radiation that has been difficult to resolve in previous studies of mitochondrial DNA sequences. Phylogenetic analyses of OR pseudogenes using both gene-tree reconciliation and supermatrix methods yielded fully resolved, consistently supported relationships among members of four delphinid subfamilies. Alternative minimizations of gene duplications, gene duplications plus gene losses, deep coalescence events, and nucleotide substitutions plus indels returned highly congruent phylogenetic hypotheses. Novel DNA sequence data for six single-copy nuclear loci and three mitochondrial genes (> 5000 aligned nucleotides) provided an independent test of the OR trees. Nucleotide substitutions and indels in OR pseudogenes showed a very low degree of homoplasy in comparison to mitochondrial DNA and, on average, provided more variation than single-copy nuclear DNA. Our results suggest that phylogenetic analysis of the large OR superfamily will be effective for resolving relationships within Cetacea whether supermatrix or gene-tree reconciliation procedures are used.     

Superfamily active site templates

We show that three-dimensional signatures consisting of only a few functionally important residues can be diagnostic of membership in superfamilies of enzymes. Using the enolase superfamily as a model system, we demonstrate that such a signature, or template, can identify superfamily members in structural databases with high sensitivity and specificity. This is remarkable because superfamilies can be highly diverse, with members catalyzing many different overall reactions; the unifying principle can be a conserved partial reaction or chemical capability. Our definition of a superfamily thus hinges on the disposition of residues involved in a conserved function, rather than on fold similarity alone. A clear advantage of basing structure searches on such active site templates rather than on fold similarity is the specificity with which superfamilies with distinct functional characteristics can be identified within a large set of proteins with the same fold, such as the (beta/alpha)8 barrels. Preliminary results are presented for an additional group of enzymes with a different fold, the haloacid dehalogenase superfamily, suggesting that this approach may be generally useful for assigning reading frames of unknown function to specific superfamilies and thereby allowing inference of some of their functional properties.     

Neurochemical and Molecular Pharmacological Aspects of the GABAB Receptor

Metabotropic gamma-aminobutyric acid (GABA)B receptors are known to modulate the synaptic release of various neurotransmitters in the nervous system. Activation of GABA(B) receptor induces the inhibition of adenylyl cyclase activity, while it does not stimulate the formation of inositol phosphates. Activation of a potassium conductance and suppression of a calcium conductance are also recognized, similarly to some of G protein-coupled receptors. Recent molecular cloning has revealed that GABA(B) receptor possesses a large extracellular domain including the binding site for GABA and seven transmembrane domains. Their molecular structures in the brain are unique and interesting because of heterodimerization consisting of two distinct genes: GABABR1 and GABABR2. Such assembled receptors can be classified as a novel type of the metabotropic receptor superfamily.     

Evolution of protein superfamilies and bacterial genome size

We present the structural annotation of 56 different bacterial species based on the assignment of genes to 816 evolutionary superfamilies in the CATH domain structure database. These assignments have enabled us to analyse the recurrence of specific superfamilies within and across the genomes. We have selected the superfamilies that have a very broad representation and therefore appear to be universally distributed in a significant number of bacterial lineages. Occurrence profiles of these universally distributed superfamilies are compared with genome size in order to estimate the correlation between superfamily duplication and the increase in proteome size. This distinguishes between those size-dependent superfamilies where frequency of occurrence is highly correlated with increase in genome size, and size-independent superfamilies where no correlation is observed. Consideration of the size correlation and the ratio between the mean and the standard deviations for all the superfamily profiles allows more detailed subdivisions and classification of superfamilies. For example, within the size-independent superfamilies, we distinguished a group that are distributed evenly amongst all the genomes. Within the size-dependent superfamilies we differentiated two groups: linearly distributed and non-linearly distributed. Functional annotation using the COG database was performed for all superfamilies in each of these groups, and this revealed significant differences amongst the three sets of superfamilies. Evenly distributed, size-independent domains are shown to be involved primarily in protein translation and biosynthesis. For the size-dependent superfamilies, linearly distributed superfamilies are involved mainly in metabolism, and non-linearly distributed superfamily domains are involved principally in gene regulation.     

Conversion of the ion selectivity of the 5-HT(3a) receptor from cationic to anionic reveals a conserved feature of the ligand-gated ion channel superfamily.

The 5-hydroxytryptamine(3) (5-HT(3)) receptor is a member of a superfamily of ligand-gated ion channels, which includes nicotinic acetylcholine, gamma-aminobutyric acid, and glycine receptors. The receptors are either cation or anion selective, leading to their distinctive involvement in either excitatory or inhibitory neurotransmission. Using a combination of site-directed mutagenesis and electrophysiological characterization of homomeric 5-HT(3A) receptors expressed in HEK293 cells, we have identified a set of mutations that convert the ion selectivity of the 5-HT(3A) receptor from cationic to anionic; these were substitution of V13'T in M2 together with neutralization of glutamate residues (E-1'A) and the adjacent insertion of a proline residue (P-1') in the M1-M2 loop. Mutant receptors showed significant chloride permeability (P(Cl)/P(Na) = 12.3, P(Na)/P(Cl) = 0.08), whereas WT receptors are predominantly permeable to sodium (P(Na)/P(Cl) > 20, P(Cl)/P(Na) < 0.05). Since the equivalent mutations have previously been shown to convert alpha7 nicotinic acetylcholine receptors from cationic to anionic (Galzi J.-L., Devillers-Thiery, A, Hussy, N., Bertrand, S. Changeux, J. P., and Bertrand, D. (1992) Nature 359, 500-505) and, recently, the converse mutations have allowed the construction of a cation selective glycine receptor (Keramidas, A., Moorhouse, A. J., French, C. R., Schofield, P. R., and Barry, P. H. (2000) Biophys. J. 78, 247-259), it appears that the determinants of ion selectivity represent a conserved feature of the ligand-gated ion channel superfamily.     

Conversion of the ion selectivity of the 5-HT(3a) receptor from cationic to anionic reveals a conserved feature of the ligand-gated ion channel superfamily

The 5-hydroxytryptamine(3) (5-HT(3)) receptor is a member of a superfamily of ligand-gated ion channels, which includes nicotinic acetylcholine, gamma-aminobutyric acid, and glycine receptors. The receptors are either cation or anion selective, leading to their distinctive involvement in either excitatory or inhibitory neurotransmission. Using a combination of site-directed mutagenesis and electrophysiological characterization of homomeric 5-HT(3A) receptors expressed in HEK293 cells, we have identified a set of mutations that convert the ion selectivity of the 5-HT(3A) receptor from cationic to anionic; these were substitution of V13'T in M2 together with neutralization of glutamate residues (E-1'A) and the adjacent insertion of a proline residue (P-1') in the M1-M2 loop. Mutant receptors showed significant chloride permeability (P(Cl)/P(Na) = 12.3, P(Na)/P(Cl) = 0.08), whereas WT receptors are predominantly permeable to sodium (P(Na)/P(Cl) > 20, P(Cl)/P(Na) < 0.05). Since the equivalent mutations have previously been shown to convert alpha7 nicotinic acetylcholine receptors from cationic to anionic (Galzi J.-L., Devillers-Thiery, A, Hussy, N., Bertrand, S. Changeux, J. P., and Bertrand, D. (1992) Nature 359, 500-505) and, recently, the converse mutations have allowed the construction of a cation selective glycine receptor (Keramidas, A., Moorhouse, A. J., French, C. R., Schofield, P. R., and Barry, P. H. (2000) Biophys. J. 78, 247-259), it appears that the determinants of ion selectivity represent a conserved feature of the ligand-gated ion channel superfamily.     

Photoactivated derivatives of peptide and protein ligands in the studies of neuroreceptors

The application of photoactivatable derivatives for studies of different types of neuroreceptors belonging to two superfamilies, G-protein dependent receptors and ligand-gated ionic channels, is discussed. Studies of the structure of voltage-gated and ion-gated channels with the use of specific photoactivatable derivatives of neurotoxins are also described. Possibilities and prospects for the application of photoactivatable ligands in studies of the spatial structure of neuroreceptors are reviewed.     

Systematic Analysis of γ-Aminobutyric Acid (GABA) Metabolism and Function in the Social Amoeba Dictyostelium discoideum

While GABA has been suggested to regulate spore encapsulation in the social amoeba Dictyostelium discoideum, the metabolic profile and other potential functions of GABA during development remain unclear. In this study, we investigated the homeostasis of GABA metabolism by disrupting genes related to GABA metabolism and signaling. Extracellular levels of GABA are tightly regulated during early development, and GABA is generated by the glutamate decarboxylase, GadB, during growth and in early development. However, overexpression of the prespore-specific homologue, GadA, in the presence of GadB reduces production of extracellular GABA. Perturbation of extracellular GABA levels delays the process of aggregation. Cytosolic GABA is degraded by the GABA transaminase, GabT, in the mitochondria. Disruption of a putative vesicular GABA transporter (vGAT) homologue DdvGAT reduces secreted GABA. We identified the GABA_B receptor-like family member GrlB as the major GABA receptor during early development, and either disruption or overexpression of GrlB delays aggregation. This delay is likely the result of an abolished pre-starvation response and late expression of several “early” developmental genes. Distinct genes are employed for GABA generation during sporulation. During sporulation, GadA alone is required for generating GABA and DdvGAT is likely responsible for GABA secretion. GrlE but not GrlB is the GABA receptor during late development.