Evolutionary History of Cathepsin L (L-like) Family Genes in Vertebrates

Cathepsin L family, an important cysteine protease found in lysosomes, is categorized into cathepsins B, F, H, K, L, S, and W in vertebrates. This categorization is based on their sequence alignment and traditional functional classification, but the evolutionary relationship of family members is unclear. This study determined the evolutionary relationship of cathepsin L family genes in vertebrates through phylogenetic construction. Results showed that cathepsins F, H, S and K, and L and V were chronologically diverged. Tandem-repeat duplication was found to occur in the evolutionary history of cathepsin L family. Cathepsin L in zebrafish, cathepsins S and K in xenopus, and cathepsin L in mice and rats underwent evident tandem-repeat events. Positive selection was detected in cathepsin L-like members in mice and rats, and amino acid sites under positive selection pressure were calculated. Most of these sites appeared at the connection of secondary structures, suggesting that the sites may slightly change spatial structure. Severe positive selection was also observed in cathepsin V (L2) of primates, indicating that this enzyme had some special functions. Our work provided a brief evolutionary history of cathepsin L family and differentiated cathepsins S and K from cathepsin L based on vertebrate appearance. Positive selection was the specific cause of differentiation of cathepsin L family genes, confirming that gene function variation after expansion events was related to interactions with the environment and adaptability.     

Evolutionary Patterns of Gene Families Generated in the Early Stage of Vertebrates

In this paper we have analyzed 49 vertebrate gene families that were generated in the early stage of vertebrates and/or shortly before the origin of vertebrates, each of which consists of three or four member genes. We have dated the first (T₁) and second (T₂) gene duplications of 26 gene families with 3 member genes. The means of T₁ (594 mya) and T₂ (488 mya) are largely consistent to a well-cited version of two-round (2R) genome duplication theory. Moreover, in most cases, the time interval between two successive gene duplications is large enough that the fate of duplicate genes generated by the first gene duplication was likely to be determined before the second one took place. However, the phylogenetic pattern of 23 gene families with 4 members is complicated; only 5 of them are predicted by 2R model, but 11 families require an additional gene (or genome) duplication. For the rest (7 families), at least one gene duplication event had occurred before the divergence between vertebrate and Drosophila, indicating a possible misleading of the 4:1 rule (member gene ratio between vertebrates and invertebrates). Our results show that Ohno's 2R conjecture is valid as a working hypothesis for providing a most parsimonious explanation. Although for some gene families, additional gene duplication is needed, the credibility of the third genome duplication (3R) remains to be investigated. Received: 13 December 1999 / Accepted: 7 April 2000     

Axon Guidance Events in the Wiring of the Mammalian Olfactory System

The detection of odorant signals from the environment and the generation of appropriate behavioral outputs in response to these signals rely on the olfactory system. Olfactory sensory neurons (OSNs) of the olfactory epithelium are located in the nasal cavity and project axons that synapse onto dendrites of second-order neurons in the olfactory bulb (OB) that in turn relay the information gathered to higher order regions of the brain. The connections formed are remarkably accurate such that axons of OSNs expressing the same olfactory receptor innervate specific glomeruli within the complex three-dimensional structure that represents the OB. The molecular determinants that control this complex process are beginning to be identified. In this review, we discuss the role of various families of axon guidance cues and of recently characterized families of adhesion molecules in the formation of stereotypic connections in the olfactory system of mice. Cho and Prince contributed equally.     

Elongation Factor 1A Family Regulates the Recycling of the M4 Muscarinic Acetylcholine Receptor

In this study, we tested the hypothesis that the elongation 1A (eEF1A) family regulates the cell surface density of the M4 subtype of the muscarinic acetylcholine receptors (mAChR) following agonist-induced internalization. Here, we show that mouse brains lacking eEF1A2 have no detectable changes in M4 expression or localization. We, however, did discover that eEF1A1, the other eEF1A isoform, is expressed in adult neurons contrary to previous reports. This novel finding suggested that the lack of change in M4 expression and distribution in brains lacking eEF1A2 might be due to compensatory effects of eEF1A1. Supporting this theory, we demonstrate that the overexpression of either eEF1A1 or eEF1A2 inhibits M4 recovery to the cell surface after agonist-induced internalization in PC12 cells. Furthermore, eEF1A1 or eEF1A2 had no effect on the recovery of the M1 subtype in PC12 cells. These results demonstrate the novel ability of the eEF1A family to specifically regulate the M4 mAChR.     

Stable Expression of a Synthetic Gene for the Human Motilin Receptor: Use in an Aequorin-Based Receptor Activation Assay

A synthetic gene for the human motilin receptor containing 33 unique restriction sites was designed and stably coexpressed in HEK293 cells with the bioluminescent Ca(2+) indicator protein aequorin. The dose-dependent response of the receptor to motilin was demonstrated using transient transfections, and a stable cell line was selected. [(125)I]Motilin binding was used to estimate receptor expression level for the stable cell line, and titration of a membrane preparation indicated a K(d) value of 0.8 nM. The same cell line was used to evaluate a panel of erythromycin-derived agonists and provided EC(50) values for receptor activation that agree closely with data obtained in contractility assays. The peptide antagonist ANQ11125 (Phe3Leu13 motilin 1-14) inhibited motilin induced response with a K(i) value of 10 nM. The system is well-suited for the screening of compound libraries and receptor mutagenesis studies.     

Evolutionary Genomics of Weedy Rice in the USA

Red rice is an interfertiie, weedy form of cultivated rice (Oryza sativa L.) that competes aggressively with the cropin the southern US, reducing yields and contaminating harvests. No wild Oryza species occur In North America andthe weed has been proposed to have evolved through multiple mechanisms, including "de-domestication" of UScrop cultivars, accidental introduction of Asian weeds, and hybridization between US crops and Asian wild/weedyOryza strains. The phenotype of US red rice ranges from "crop mimics", which share some domestication traitswith the crop, to strains closely resembling Asian wild Oryza species. Assessments of genetic diversity haveindicated that many weed strains are closely related to Asian taxa (including indica and aus rice varieties, whichhave never been cultivated in the US, and the Asian crop progenitor O. rufipogon), whereas others show geneticsimilarity to the tropical japonica varieties cultivated in the southern US. Herein, we review what is known aboutthe evolutionary origins and genetic diversity of US red rice and describe an ongoing research project to furthercharacterize the evolutionary genomics of this aggressive weed.     

High Incidence of Interchromosomal Transpositions in the Evolutionary History of a Subset of Or Genes in Drosophila

In insects, the odorant receptor (Or) multigene family is an intermediate-sized family with genes present in all chromosomes, indicating that duplication followed by interchromosomal transposition played an important role in the early stages of the family evolution. Here, we have explored the occurrence of interchromosomal transpositions in more recent stages through the comparative analysis of a subset of Or genes in Drosophila, where the gene content of chromosomal arms is highly conserved. The studied subset consisted of 11 Or genes located on the left arm of chromosome 3 (Muller’s D element) in D. melanogaster. Our study focused on the number and chromosomal arm location of these members of the family across the 12 Drosophila species with complete genome sequences. In contrast to previous results from in situ hybridization comparative mapping that were mainly based on single-copy genes, our study, based on members of a multigene family of moderate size, revealed repeated interchromosomal transposition events and a complex history of some of the studied genes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phylogeny and Evolution of the SERA Multigene Family in the Genus <Emphasis Type="Italic">Plasmodium</Emphasis>

The serine repeat antigen gene family of Plasmodium falciparum (Pf-SERA) consists of nine gene members. By sequence similarity search, 45 genes were identified to be homologous to the Pf-SERA genes in the ongoing seven Plasmodium genome sequencing project databases for the species: P. reichenowi, P. vivax, P. knowlesi, P. yoelii, P. berghei, P. chabaudi, and P. gallinaceum. In combination with additional PCR-based sequencing, we found that almost all SERA genes in each species were aligned in a tandem cluster and sandwiched between two conserved hypothetical protein genes, except for P. reichenowi, which could not be confirmed. The minimum and maximum numbers of clustered genes were 2 and 12 for P. gallinaceum and P. vivax, respectively. The best tree of the maximum likelihood analysis demonstrated that all Plasmodium SERA homologues, except for SERA1 of P. gallinaceum (Pg-SERA1), can be classified into four groups, represented by Pf-SERA5, Pf-SERA6, Pf-SERA7, and Pf-SERA8. Genes in the Pf-SERA8 group, although highly divergent and distantly related to the sequences of other groups, were not pseudogenes. P. berghei SERA5, the counterpart of Pf-SERA8, was expressed in the mosquito stage. P. gallinaceum lacks the orthologues to Pf-SERA5, Pf-SERA6, and Pf-SERA7, suggesting that P. gallinaceum diverged from a common ancestor of all eight Plasmodium species examined before gene duplication(s) occurred to generate these paralogous groups. Here, we reveal an evolutionary trail of SERA gene cluster in the genus Plasmodium and discuss a phylogeny of Plasmodium species from the viewpoint of the evolution of a multigene family.     

Isolation and Characterization of Five Actin cDNAs from the Cestode Diphyllobothrium dendriticum: A Phylogenetic Study of the Multigene Family

Five cDNAs (pDidact2–pDidact6), representing different actin genes, were isolated from a Diphyllobothrium dendriticum cDNA library, and the DNA as well as the putative amino acid sequences were determined. The corresponding Didact2 and Didact4 genes code for peptides 376 amino acids long, with molecular weights 41,772 and 41,744 Da, respectively, while the deduced Didact3 protein is 377 amino acids long and weighs 41,912 Da. The pDidact5 and -6 cDNAs lack nucleotides corresponding to three to six amino acids at the amino-terminus. Two of the five cDNAs contain the conventional AATAAA as the putative polyadenylation signal, one has the common variant ATTAAA, whereas the hexanucleotide AATAGA is found 15 and 18 nucleotides, respectively, upstream of the poly(A) site in two of the cDNAs. Phylogenetic studies including 102 actin protein sequences revealed that there are at least four different types of cestode actins. In this study three of these types were found to be expressed in the adult D. dendriticum tapeworm. Structurally the cestode actin groupings differ from each other to an extent seen only among the metazoan actins between the vertebrate muscle and cytoplasmic isoforms. In the phylogenetic trees constructed, cestode actins were seen to map to two different regions, one on the border of the metazoan actins and the other within this group. It is, however, difficult to say whether the cestode actins branched off early in the metazoan evolution or if this position in the phylogenetic tree only reflects upon differences in evolutionary rate. Received: 19 June 1996 / Accepted: 20 August 1996     

Molecular evolution of the nicotinic acetylcholine receptor: An example of multigene family in excitable cells

An extensive phylogenetic analysis of the nicotinic-acetylcholine-receptor subunit gene family has been performed by cladistic and phenetic methods. The conserved parts of amino acid sequences have been analyzed by CLUSTAL V and PHYLIP software. The structure of the genes was also taken in consideration. The results show that a first gene duplication may have occurred before the appearance of Bilateria. Three subfamilies then appeared: I-the neuronal -bungarotoxin binding-site subunits (7, 8); III-the neuronal nicotinic subunits (2–6, 2–4), which also contain the muscle acetylcholine-binding subunit (1); and IV—the muscle non- subunits (1, , ). The Insecta subunits (subfamily II) could be orthologous to family III and IV. Several tissular switches of expression from neuron to muscle and the converse can be inferred from the extant expression of subunits and the reconstructed trees. The diversification of the neuronal nicotinic subfamily begins in the stem lineage of chordates, the last duplications occurring shortly before the onset of the mammalian lineage. Such evolution parallels the increase in complexity of the cholinergic systems.Abbreviations -Bgt -bungarotoxin - ACh acetylcholine - MP maximum of parsimony - MYA million years ago - NJ neighbor-joining - nAChR nicotinic acetylcholine receptor Correspondence to: N. Le Novère     

Evolutionary Aspects of Functional and Pseudogene Members of the Phytochrome Gene Family in Scots Pine

According to the neutral theory of evolution, mutation and genetic drift are the only forces that shape unconstrained, neutral, gene evolution. Thus, pseudogenes (which often evolve neutrally) provide opportunities to obtain direct estimates of mutation rates that are not biased by selection, and gene families comprising functional and pseudogene members provide useful material for both estimating neutral mutation rates and identifying sites that appear to be under positive or negative selection pressures. Conifers could be very useful for such analyses since they have large and complex genomes. There is evidence that pseudogenes make significant contributions to the size and complexity of gene families in pines, although few studies have examined the composition and evolution of gene families in conifers. In this work, I examine the complexity and rates of mutation of the phytochrome gene family in Pinus sylvestris and show that it includes not only functional genes but also pseudogenes. As expected, the functional PHYO does not appear to have evolved neutrally, while phytochrome pseudogenes show signs of unconstrained evolution.     

Evolution of the bovine lysozyme gene family: Changes in gene expression and reversion of function

Recruitment of lysozyme to a digestive function in ruminant artiodactyls is associated with amplification of the gene. At least four of the approximately ten genes are expressed in the stomach, and several are expressed in nonstomach tissues. Characterization of additional lysozymelike sequences in the bovine genome has identified most, if not all, of the members of this gene family. There are at least six stomachlike lysozyme genes, two of which are pseudogenes. The stomach lysozyme pseudogenes show a pattern of concerted evolution similar to that of the functional stomach genes. At least four nonstomach lysozyme genes exist. The nonstomach lysozyme genes are not monophyletic. A gene encoding a tracheal lysozyme was isolated, and the stomach lysozyme of advanced ruminants was found to be more closely related to the tracheal lysozyme than to the stomach lysozyme of the camel or other nonstomach lysozyme genes of ruminants. The tracheal lysozyme shares with stomach lysozymes of advanced ruminants the deletion of amino acid 103, and several other adaptive sequence characteristics of stomach lysozymes. I suggest here that tracheal lysozyme has reverted from a functional stomach lysozyme. Tracheal lysozyme then represents a second instance of a change in lysozyme gene expression and function within ruminants. Correspondence to: D.M. Irwin     

Evolutionary analysis of rhodopsin and cone pigments: connecting the three-dimensional structure with spectral tuning and signal transfer

Extensive sequence data and structural sampling of expressed proteins from different species lead to the idea that entire molecules or specific domain folds belong to large superfamilies of proteins. A subset of G protein-coupled receptors, one of the largest families involved in cellular signaling, rod and cone opsins are involved in phototransduction in photoreceptor cells. Here, the evolutionary analysis of opsin sequences and structures predicts key residues involved in the transmission of the signal from the binding site of the chromophore to the cytoplasmic surface and residues that are involved in the spectral tuning of opsins to short wavelengths of light.     

Developmental genetic profiles of glutamate receptor system,neuromodulator system,protector of normal tissue and mitochondria,and reelin in marmoset cortex: Potential molecular mechanisms of pruning phase of spines in primate synaptic formation process during the end of infancy and prepuberty (II)

This is the second report of a series paper, which reports molecular mechanisms underlying the occurrence of pruning spine phase after rapid spinogenesis phase in neonates and young infant in the primate brain. We performed microarray analysis between the peak of spine numbers [postnatal 3 months (M)] and spine pruning (postnatal 6 M) in prefrontal, inferior temporal, and primary visual cortices of the common marmoset (Callithrix jacchus). The pruning phase is not clearly defined in rodents but is in primates including the marmoset. The differentially expressed genes between 3 M and 6 M in all three cortical areas were selected by two-way analysis of variance. The list of selected genes was analyzed by canonical pathway analysis using “Ingenuity Pathway Analysis of complex omics data” (IPA; Ingenuity Systems, Qiagen, Hilden, Germany). In this report, we discuss these lists of genes for the glutamate receptor system, G-protein-coupled neuromodulator system, protector of normal tissue and mitochondria, and reelin. (1) Glutamate is a common neurotransmitter. Its receptors AMPA1, GRIK1, and their scaffold protein DLG4 decreased as spine numbers decreased. Instead, GRIN3 (NMDA receptor) increased, suggesting that strong NMDA excitatory currents may be required for a single neuron to receive sufficient net synaptic activity in order to compensate for the decrease in synapse. (2) Most of the G protein-coupled receptor genes (e.g., ADRA1D, HTR2A, HTR4, and DRD1) in the selected list were upregulated at 6 M. The downstream gene ROCK2 in these receptor systems plays a role of decreasing synapses, and ROCK2 decreased at 6 M. (3) Synaptic phagosytosis by microglia with complement and other cytokines could cause damage to normal tissue and mitochondria. SOD1, XIAP, CD46, and CD55, which play protective roles in normal tissue and mitochondria, showed higher expression at 6 M than at 3 M, suggesting that normal brain tissue is more protected at 6 M. (4) Reelin has an important role in cortical layer formation. In addition, RELN and three different pathways of reelin were expressed at 6 M, suggesting that new synapse formation decreased at that age. Moreover, if new synapses were formed, their positions were free and probably dependent on activity.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

Characterization of Calcium-activated Chloride Channels in Patches Excised from the Dendritic Knob of Mammalian Olfactory Receptor Neurons

We investigated the properties of calcium-activated chloride channels in inside-out membrane patches from the dendritic knobs of acutely dissociated rat olfactory receptor neurons. Patches typically contained large calcium-activated currents, with total conductances in the range 30–75 nS. The dose response curve for calcium exhibited an EC₅₀ of about 26 μm. In symmetrical NaCl solutions, the current-voltage relationship reversed at 0 mV and was linear between −80 and +70 mV. When the intracellular NaCl concentration was progressively reduced from 150 to 25 mM, the reversal potential changed in a manner consistent with a chloride-selective conductance. Indeed, modeling these data with the Goldman-Hodgkin-Katz equation revealed a P_Na/P_Cl of 0.034. The halide permeability sequence was P_Cl > P_F > P_I > P_Br indicating that permeation through the channel was dominated by ion binding sites with a high field strength. The channels were also permeable to the large organic anions, SCN⁻, acetate⁻, and gluconate⁻, with the permeability sequence P_Cl > P_SCN > P_acetate > P_gluconate. Significant permeation to gluconate ions suggested that the channel pore had a minimum diameter of at least 5.8 A. Received: 16 April 1997/Revised: 3 October 1997     

Fifteen Million Years of Evolution in the Oryza Genus Shows Extensive Gene Family Expansion

In analyzing gene families in the whole-genome sequences available for O. sativa （AA）, O. glaberrima （AA）, and O. brachyantha （FF）, we observed large size expansions in the AA genomes compared to FF genomes for the superfamilies F-box and NB-ARC, and five additional families： the Aspartic proteases, BTB/POZ proteins （BTB）, Glutaredoxins, Trypsin a-amylase inhibitor proteins, and Zf-Dof proteins. Their evolutionary dynamic was investigated to understand how and why such important size variations are observed between these closely related species. We show that expansions resulted from both amplification, largely by tandem duplications, and contraction by gene losses. For the F-box and NB-ARC gene families, the genes conserved in all species were under strong purifying selection while expanded orthologous genes were under more relaxed purifying selection. In F-box, NB-ARC, and BTB, the expanded groups were enriched in genes with little evidence of expression, in comparison with conserved groups. We also detected 87 loci under positive selection in the expanded groups. These results show that most of the duplicated copies in the expanded groups evolve neutrally after duplication because of functional redundancy but a fraction of these genes were preserved following neofunctionalization. Hence, the lineage-specific expansions observed between Oryza species were partly driven by directional selection.     

Concentration Dependence of Sodium Permeation and Sodium Ion Interactions in the Cyclic AMP-gated Channels of Mammalian Olfactory Receptor Neurons

The dependence of currents through the cyclic nucleotide-gated (CNG) channels of mammalian olfactory receptor neurons (ORNs) on the concentration of NaCl was studied in excised inside-out patches from their dendritic knobs using the patch-clamp technique. With a saturating concentration (100 μm) of adenosine 3′, 5′-cyclic monophosphate (cAMP), the changes in the reversal potential of macroscopic currents were studied at NaCl concentrations from 25 to 300 mm. In symmetrical NaCl solutions without the addition of divalent cations, the current-voltage relations were almost linear, reversing close to 0 mV. When the external NaCl concentration was maintained at 150 mm and the internal concentrations were varied, the reversal potentials of the cAMP-activated currents closely followed the Na⁺ equilibrium potential indicating that P _Cl/P _Na≈ 0. However, at low external NaCl concentrations (≤100 mm) there was some significant chloride permeability. Our results further indicated that Na⁺ currents through these channels: (i) did not obey the independence principle; (ii) showed saturation kinetics with K _ms in the range of 100–150 mm and (iii) displayed a lack of voltage dependence of conductance in asymmetric solutions that suggested that ion-binding sites were situated midway along the channel. Together, these characteristics indicate that the permeation properties of the olfactory CNG channels are significantly different from those of photoreceptor CNG channels. Received: 7 November 1996/Revised: 24 March 1997