Positive selection during the diversification of class I vomeronasal receptor-like (V1RL) genes,putative pheromone receptor genes,in human and primate evolution

Vomeronasal receptors are the major receptors for pheromones in vertebrates, and five putative type 1 vomeronasal receptors (V1RL) have been identified in humans. The evolution of the V1RL1 gene in non-human primates, and patterns of selection on V1RL genes, were investigated. The presumed ortholog of V1RL1 was sequenced from 13 species of nonhuman primate, and in eight of these species V1RL1 was a pseudogene. Phylogenetic reconstructions reveal that V1RL1 pseudogene formation occurred independently in multiple primate lineages. Using maximum likelihood estimates of dN/dS ratios in PAML, we show that V1RL genes have evolved under neutral evolution in lineages in which they became a pseudogene. In contrast, among lineages in which V1RL genes contain an open reading frame, the majority of sites are under purifying selection and a minority are under significant positive selection. These results provide an interesting case where all three categories of selection can be teased apart in the same data set using maximum likelihood methods. The finding of positive selection on V1RL genes during primate evolution provides indirect support for the hypothesis that V1RL genes have a function in species-specific pheromone detection in primates.     

Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Balancing selection acting over the evolutionary history of a lineage can result in the retention of alleles among species for longer than expected under neutral evolution. The associated pattern of trans-species polymorphism, in which similar or even identical alleles are shared among species, is often used to infer that balancing selection has occurred. The genes of the major histocompatibility complex (MHC) are thought to be subject to balancing selection that maintains alleles associated with response to specific pathogens. To explore the role of balancing selection in shaping MHC diversity in ctenomyid rodents, we examined allelic variability at the class II DRB and DQA loci in 18 species in the genus Ctenomys. Previous studies of four of these species had revealed significant within-population evidence of positive selection on MHC loci. The current study expands upon these analyses to (1) evaluate among-species evidence of positive selection and (2) explore the potential for balancing selection on MHC genes. Interspecific nucleotide sequence variation revealed significant evidence of positive selection on the DRB and DQA loci. At the same time, comparisons of phylogenetic trees for these MHC loci with a putative species tree based on mitochondrial sequence data revealed multiple examples of trans-specific polymorphism, including sharing of identical DRB and DQA alleles among distantly related species of Ctenomys. These findings suggest that MHC genes in these animals have historically been subject to balancing selection and yield new insights into the complex suite of forces shaping MHC diversity in free-living vertebrates.     

Isolation and chromosomal localization of the human En-2 gene

By low stringency hybridization we have isolated from a human cosmid genomic library sequences homologous with a probe from the Drosophila engrailed gene. Partial nucleotide sequence analysis shows a consensus splice acceptor site followed by an open reading frame (ORF) that can encode 104 amino acids; the first 94 amino acids have 71% identity with the Drosophila engrailed protein. The shared region contains a homeo domain and is within the region of engrailed shared with the Drosophila invected gene and the mouse En-1 and En-2 genes. At the amino acid level, the human sequence is 85% identical with the mouse En-1 gene and 100% identical with the mouse En-2 gene. Hybridization against a panel of human-hamster somatic cell hybrids maps this human En-2 gene to chromosome 7, and regional mapping by in situ hybridization to human chromosomes localizes it to region 7q36 at the end of the long arm.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A Novel Human Gene Encoding a G-Protein-Coupled Receptor (GPR15) Is Located on Chromosome 3

We used sequence similarities among G-protein-coupled receptor genes to discover a novel receptor gene. Using primers based on conserved regions of the opioid-related receptors, we isolated a PCR product that was used to locate the full-length coding region of a novel human receptor gene, which we have namedGPR15.A comparison of the amino acid sequence of the receptor encoded byGPR15with other receptors revealed that it shared sequence identity with the angiotensin II AT1 and AT2 receptors, the interleukin 8b receptor, and the orphan receptors GPR1 and AGTL1.GPR15was mapped to human chromosome 3q11.2–q13.1.     

甜瓜抗霜霉病基因同源序列克隆与分析

采用RT—PCR扩增的方法,从高抗霜霉病甜瓜品种‘日本安农二号’中克隆到约3kb的cDNA片段（命名为MRGH-D,该基因是一个连续的通读编码框,编码1007个氨基酸。推测的蛋白质分子量为113．7kDa,等电点为7．88,蛋白质预测无跨膜区。根据推测的氨基酸序列,该基因属于TIR—NBS—LRR类抗病基因,具有TIR-NBS—LRR类抗病基因所有的保守结构域。核苷酸序列和氨基酸序列同源性分析结果显示,MRGH-J与甜瓜抗病基因的同源序列MRGHl2及抗霜霉病相关基因mp-19均具有高达99％的同源性,推测该基因可能在甜瓜抗霜霉病中起作用。     

“两广二号”家蚕抗病毒基因Bmlipase-1和BmSP-2克隆与原核表达

本研究以优良杂交品种"两广二号"家蚕为试材,克隆了该杂交品种家蚕两个抗家蚕核型多角体病毒(BmNPV)基因:脂肪酶基因Bmlipase-1和丝氨酸蛋白酶基因BmSP-2,测序并分别与不同品种蚕的同源基因序列进行比较。结果显示,"两广二号"家蚕Bmlipase-1基因ORF长度为885bp,编码294个氨基酸,BmSP-2扩增长度为855bp,编码284个氨基酸;它们的核苷酸和推导氨基酸序列同源性皆达92%以上,Bmlipase-1更保守,同源性大于99%";两广二号"家蚕的Bmlipase-1基因脂肪酶活化部位和BmSP-2基因酶催化三联体位点的氨基酸残基与不同品种蚕的完全相同。以上结果说明这两个抗病毒基因在蚕的遗传进化过程中高度保守,提示其可能在机体消化或者免疫防御方面起着重要生理作用。将这两个抗病毒基因在大肠杆菌BL21中进行融合表达,获得的融合Bmlipase-1和BmSP-2蛋白分子量分别为47kD和42kD左右。     

Adaptive evolution of fertilization proteins within a genus: variation in ZP2 and ZP3 in deer mice (Peromyscus)

Rapid evolution of reproductive proteins has been documented in a wide variety of taxa. In internally fertilized species, knowledge about the evolutionary dynamics of these proteins between closely related taxa is primarily limited to accessory gland proteins in the semen of Drosophila. Investigation of additional taxa and functional classes of proteins is necessary in order to determine if there is a general pattern of adaptive evolution of reproductive proteins between recently diverged species. We performed an evolutionary analysis of 2 egg coat proteins, ZP2 and ZP3, in 15 species of deer mice (genus Peromyscus). Both of these proteins are involved in egg-sperm binding, a critical step in maintaining species-specific fertilization. Here, we show that Zp2 and Zp3 gene trees are not consistent with trees based on nonreproductive genes, Mc1r and Lcat, where species formed monophyletic clades. In fact, for both of the reproductive genes, intraspecific amino acid variation was extensive and alleles were sometimes shared across species. We document positive selection acting on ZP2 and ZP3 and identify specific amino acid sites that are likely targets of selection using both maximum likelihood approaches and patterns of parallel amino acid change. In ZP3, positively selected sites are clustered in and around the region implicated in sperm binding in Mus, suggesting changes may impact egg-sperm binding and fertilization potential. Finally, we identify lineages with significantly elevated rates of amino acid substitution using a Bayesian mapping approach. These findings demonstrate that the pattern of adaptive reproductive protein evolution found at higher taxonomic levels can be documented between closely related mammalian species, where reproductive isolation has evolved recently.     

Genomic organization and characterization of two vomeronasal 1 receptor-like genes (ora1 and ora2) in Atlantic salmon Salmo salar

Olfactory receptors are encoded by three large multigene superfamilies (OR, V1R and V2R) in mammals. Fish do not possess a vomeronasal system; therefore, it has been proposed that their V1R-like genes be classified as olfactory receptors related to class A G protein-coupled receptors (ora). Unlike mammalian genomes, which contain more than a hundred V1R genes, the five species of teleost fish that have been investigated to date appear to have six ora genes (ora1-6) except for pufferfish that have lost ora1. The common ancestor of salmonid fishes is purported to have undergone a whole genome duplication. As salmonids have a life history that requires the use of olfactory cues to navigate back to their natal habitats to spawn, we set out to determine if ora1 or ora2 is duplicated in a representative species, Atlantic salmon (Salmo salar). We used an oligonucleotide probe designed from a conserved sequence of several teleost ora2 genes to screen an Atlantic salmon BAC library (CHORI-214). Hybridization-positive BACs belonged to a single fingerprint contig of the Atlantic salmon physical map. All were also positive for ora2 by PCR. One of these BACs was chosen for further study, and shotgun sequencing of this BAC identified two V1R-like genes, ora1 and ora2, that are in a head-to-head conformation as is seen in some other teleosts. The gene products, ora1 and ora2, are highly conserved among teleosts. We only found evidence for a single ora1-2 locus in the Atlantic salmon genome, which was mapped to linkage group 6. Fluorescent in situ hybridization (FISH) analysis placed ora1-2 on chromosome 12. Conserved synteny was found surrounding the ora1 and ora2 genes in Atlantic salmon, medaka and three-spined stickleback, but not zebrafish.     

Protamine 3 shows evidence of weak, positive selection in mouse species (genus Mus)--but it is not a protamine

Protamines are short and highly basic sperm-specific nuclear proteins that replace somatic histones during spermiogenesis in a process that is crucial for sperm formation and function. Many mammals have two protamine genes (PRM1 and PRM2) located in a gene cluster, which appears to evolve fast. Another gene in this cluster (designated protamine 3 [PRM3]) encodes a protein that is conserved among mammals but that does not seem to be involved in chromatin condensation. We have compared protein sequences and amino acid compositions of protamines in this gene cluster, searched for evidence of positive selection of PRM3, and examined whether sexual selection (sperm competition) may drive the evolution of the PRM3 gene. Nucleotide and amino acid analyses of mouse sequences revealed that PRM3 was very different from PRM1 and from both the precursor and the mature sequences of PRM2. Among 10 mouse species, PRM3 showed weak evidence of positive selection in two species, but there was no clear association with levels of sperm competition. In analyses from among mammalian species, no evidence of positive selection was found in PRM3. We conclude that PRM3 exhibits several clear differences from other protamines and, furthermore, that it cannot be regarded as a true protamine.     

Evaluation of candidate F-box genes for the pollen S of gametophytic self-incompatibility in the Pyrinae (Rosaceae) on the basis of their phylogenomic context

The recent analysis of the S-locus region of apple and Japanese pear, two species of Pyrinae (Rosaceae), suggested multiple and different F-box genes (called SFBBs) as candidates for the male determinant (pollen S) of RNase-based gametophytic self-incompatibility in these two species. Here, we followed a phylogenetic approach to take advantage of the pattern of molecular evolution of the S-locus of Pyrinae in characterizing SFBB homologs belonging to S-haplotypes of apple and three species of Pyrus (European, Japanese, and Chinese pears). Our results suggested that the S-locus region of Pyrinae contains no less than six SFBB members and that its structure seems to be rather conserved between apple and pear species. In accordance with the prevailing theory on S-haplotype evolution, the pollen S is expected to have coevolved with the S-RNase and to show some common features derived from the long-term evolution under frequency-dependent balancing selection, i.e., high sequence diversity, evidence of positive selection, and shared ancestral polymorphisms. Using this conceptual framework, we present evidence that some SFBB genes may be better candidates for pollen S in Pyrinae than others. Overall, the SFBB genes analyzed exhibited much lower sequence diversity than their associated S-RNases; likewise, they showed little or no evidence of positive selection. However, evidence of coevolution with the S-RNase clearly emerged for two of them. Altogether our results suggested different evolutionary histories for different SFBBs putatively derived from their distinct involvement in self-incompatibility.     

Identification by sequence analysis of a second rat brain cDNA encoding the dopamine (D2) receptor

A rat brain cDNA library constructed in lambda ZAP II was screened with three oligonucleotide probes based on the reported coding region of the D2 receptor gene, RGB-2. A complete cDNA clone, D2(8)-1, showing positive signals with the three probes was subsequently identified by restriction analysis and dideoxy sequence analysis to be a variant of the RGB-2 gene. Comparison of the two genes revealed almost complete homology except that D2(8)-1 contains an 87 bp insert within the protein coding region and 265 additional nucleotides 5' upstream from the 5' end reported for RGB-2. It is suggested that at least two mRNA species encoding for D2 receptors exist in rat brain, possibly resulting from alternative splicing of RNA.     

Duplicate genes for Fe-containing superoxide dismutase in Streptomyces coelicolor A3(2)

Streptomyces coelicolor Müller contains two types of superoxide dismutase (SOD) containing Ni (encoded by sodN) or Fe (encoded by sodF). Unlike a single species of Fe-containing SOD in Müller strain, multiple forms of FeSODs were detected in S. coelicolor A3(2) strain by activity staining and Western blot analysis. Genomic Southern hybridization suggested the presence of at least two copies of the sodF-like gene in A3(2). Two different genes for FeSOD (sodF1 and sodF2) were isolated from the phage library of A3(2) genome. The nucleotide sequence of the sodF1 coding region was identical with the unique sodF gene from Müller while that of sodF2 shared 88% identity. The gene products of sodF1 and sodF2 were identified by activity staining and immunoblot analysis. Expression from the sodF1 gene was repressed by nickel as sensitively as Müller sodF, suggesting the presence of Ni-responsive regulatory site within the region shared by the two genes. Among 12 other Streptomyces species examined, only S. fradiae contained two FeSOD-like polypeptides. We postulate that the additional copy of the sodF gene (sodF2) was provided by the horizontal transfer from remotely related bacteria.     

Positive selection in phytotoxic protein-encoding genes of Botrytis species

Evolutionary patterns of sequence divergence were analyzed in genes from the fungal genus Botrytis (Ascomycota), encoding phytotoxic proteins homologous to a necrosis and ethylene-inducing protein from Fusarium oxysporum. Fragments of two paralogous genes (designated NEP1 and NEP2) were amplified from all known Botrytis species and sequenced. NEP1 sequences of two Botrytis species contain premature stop codons, indicating that they may be non-functional. Both paralogs of all species encode proteins with a remarkably similar predicted secondary structure, however, they contain different types of post-translational modification motifs, which are conserved across the genus. While both NEP genes are, overall, under purifying selection, we identified a number of amino acids under positive selection based on inference using maximum likelihood models. Positively selected amino acids in NEP1 were not under selection in corresponding positions in NEP2. The biological significance of positively selected residues and the role of NEP proteins in pathogenesis remain to be resolved.     

Salmonid Opsin Sequences Undergo Positive Selection and Indicate an Alternate Evolutionary Relationship in Oncorhynchus

Positive selection can be demonstrated by statistical analysis when non-synonymous nucleotide substitutions occur more frequently than synonymous substitutions (dN>dS). This pattern of sequence evolution has been observed in the rhodopsin gene of cichlids. Mutations in opsin genes resulting in amino acid (AA) replacement appear to be associated with the evolution of specific color patterns and the evolution of courtship behaviors. Within fish, AA replacements in opsin proteins have improved vision at great depths and have occurred in deep-sea species. Salmonids experience diverse photic environments during their life history. Furthermore, sexual selection has resulted in species-specific male and female coloration during spawning. To look for evidence of positive selection in salmonid opsins, we sequenced the RH1, RH2, LWS, SWS1, and SWS2 genes from six Pacific salmon species as well as the Atlantic salmon. These salmonids include landlocked and migratory species and species that vary in their coloration during spawning. In each opsin gene comparison from all species sampled, traditional dN:dS analysis did not indicate positive selection. However, the more sensitive Creevey–McInerney statistical analysis indicates that RH1 and RH2 experienced positive selection early in the evolution and speciation of salmonids.     

Adaptive evolution in the glucose transporter 4 gene Slc2a4 in Old World fruit bats (family: Pteropodidae)

Frugivorous and nectarivorous bats are able to ingest large quantities of sugar in a short time span while avoiding the potentially adverse side-effects of elevated blood glucose. The glucose transporter 4 protein (GLUT4) encoded by the Slc2a4 gene plays a critical role in transmembrane skeletal muscle glucose uptake and thus glucose homeostasis. To test whether the Slc2a4 gene has undergone adaptive evolution in bats with carbohydrate-rich diets in relation to their insect-eating sister taxa, we sequenced the coding region of the Slc2a4 gene in a number of bat species, including four Old World fruit bats (Pteropodidae) and three New World fruit bats (Phyllostomidae). Our molecular evolutionary analyses revealed evidence that Slc2a4 has undergone a change in selection pressure in Old World fruit bats with 11 amino acid substitutions detected on the ancestral branch, whereas, no positive selection was detected in the New World fruit bats. We noted that in the former group, amino acid replacements were biased towards either Serine or Isoleucine, and, of the 11 changes, six were specific to Old World fruit bats (A133S, A164S, V377F, V386I, V441I and G459S). Our study presents preliminary evidence that the Slc2a4 gene has undergone adaptive changes in Old World fruit bats in relation to their ability to meet the demands of a high sugar diet.     

稻属中一个高度保守和组成型表达酪氨酸／丝氨酸－酸酸双特性蛋白激酶基因的克隆与分析

从水稻中克隆了一个在稻属植物中高度保守和组成型表达的丝氨酸／苏氨酸蛋白激酶基因（OsSTK）。该基因包含两个外显子和一个114bp的小内含子序列,预测编码一个419个氨基酸的蛋白质。该基因推导的氨基酸序列与其它已知序列的一致性均低于52％。利用从不同种和类型的野生稻克隆的部分该基因序列构建的系统树与野生稻的分类和进化关系相一致。OSPKN-端拥有一段富含丝氨酸、碱性氨基酸和带电荷氨基酸的特异性导肽序列,其中包含“GDGDGDGDG”短重复序列。由于该基因蛋白激酶结构域中的VIb,VIII和XI亚结构域中同时具有酪氨酸蛋白激酶和丝氨酸／苏氨酸蛋白激酶的特性,推测该基因可能同时具有催化酪氨酸和丝氨酸、苏氨酸磷酸化的双重功能。     

Two different genes encoding acetylcholinesterase existing in cotton aphid (Aphis gossypii).

Two acetylcholinesterase (AChE) genes, Ace1 and Ace2, have been cloned from cotton aphid, Aphis gossypii Glover, using the rapid amplification of cDNA ends (RACE) technique. To the best of our knowledge, this should be the first direct molecular evidence that multiple AChE genes exist in insects. The Ace1 gene was successfully amplified along its full length of 2371 bp. The open reading frame is 2031 bp long and encodes 676 amino acids (GenBank accession No. AF502082). The Ace2 gene was amplified as a mega-fragment of 2130 bp lacking part of 5'-end untranslated region (UTR). The open reading frame is 1992 bp long and ecodes a protein of 664 amino acids (GenBank accession No. AF502081). Both genes have the conserved amino acids and features shared by the AChE family, but share only 35% identity in amino acid sequence. The Ace1 gene is highly homologous to the AChE gene of Schizaphis graminum (AF321574) with 95% identity, and Ace2 to that of Myzus persicae (AF287291) with 92% identity. Phylogenetic analysis showed that the two cloned AChEs of A. gossypii are different in evolution. The phylogenetic tree generated by the PHYLIP program package inferred that AChE2 of A. gossypii is a more ancestral form of AChE. Homology modeling of structures using Torpedo californica (2ACE_) and Drosophila melanogaster (1Q09:A) native acetylcholinesterase structure as main template indicated that the two AChEs of Aphis gossypii might have different three-dimensional structures. Alternative splicing of Ace1 near the 5'-end resulting in two proteins differing by the presence or absence of a fragment of four amino acids is also reported.