Increased positive selection pressure within the complementarity determining regions of the T‐cell receptor β gene in New World monkeys

Because of the long‐term co‐evolution of TCR and MHC molecules, numerous nucleotide substitutions have accumulated within the domains of TCRβ genes. We previously found that nonsynonymous nucleotide substitutions occurred more frequently in complementarity determining region (CDR)β than in CDRα, even though only a limited number of common marmoset (Callithrix jacchus) and human T‐cell receptor β variable (TRBV) sequences were compared. This interesting finding raised the question of whether the increased selective pressure within CDRβ was species‐specific. In this study, we identified 21 TRBV region sequences from the common marmoset and performed comparative sequence analyses of the T‐cell receptor α variable (TRAV) and TRBV regions from human, chimpanzee, rhesus monkey, cotton‐top tamarin, Ma's night monkey, and common marmoset. The ratios of the number of nonsynonymous nucleotide substitutions per site (d_N) to the d_S values (d_N/d_S) were less than 1 within the framework regions (FRs) of TRAV and TRBV region sequences, suggesting that purifying selection is largely dominant within the FRs. In contrast, the d_N values were statistically significantly greater for CDRβ than for CDRα only in New World monkeys. Also, increased d_N/d_S ratios (d_N/d_S>1) were observed within CDRβ between humans and New World monkeys and, interestingly, between New World monkeys, which share a relatively recent common ancestor. Moreover, phylogenetic analysis by maximum likelihood analysis provided firm evidence to support that positive selection occurred within CDRβ along New World monkey lineages. These results suggest that increased positive selection pressure within CDRβ is common in New World monkeys rather than being species‐specific. This study provides an intriguing insight into the co‐evolution of TCR and MHC molecules within primates. Am. J. Primatol. 73:1082–1092, 2011. © 2011 Wiley‐Liss, Inc.     

Ancestral Population Sizes and Species Divergence Times in the Primate Lineage on the Basis of Intron and BAC End Sequences

The effective sizes of ancestral populations and species divergence times of six primate species (humans, chimpanzees, gorillas, orangutans, and representatives of Old World monkeys and New World monkeys) are estimated by applying the two-species maximum likelihood (ML) method to intron sequences of 20 different loci. Examination of rate heterogeneity of nucleotide substitutions and intragenic recombination identifies five outrageous loci (ODC1, GHR, HBE, INS, and HBG). The estimated ancestral polymorphism ranges from 0.21 to 0.96% at major divergences in primate evolution. One exceptionally low polymorphism occurs when African and Asian apes diverged. However, taking into consideration the possible short generation times in primate ancestors, it is concluded that the ancestral population size in the primate lineage was no smaller than that of extant humans. Furthermore, under the assumption of 6 million years (myr) divergence between humans and chimpanzees, the divergence time of humans from gorillas, orangutans, Old World monkeys, and New World monkeys is estimated as 7.2, 18, 34, and 65 myr ago, respectively, which are generally older than traditional estimates. Beside the intron sequences, three other data sets of orthologous sequences are used between the human and the chimpanzee comparison. The ML application to these data sets including 58,156 random BAC end sequences (BES) shows that the nucleotide substitution rate is as low as 0.6–0.8 × 10^–9 per site per year and the extent of ancestral polymorphism is 0.33–0.51%. With such a low substitution rate and short generation time, the relatively high extent of polymorphism suggests a fairly large effective population size in the ancestral lineage common to humans and chimpanzees.[Reviewing Editor: Dr. Magnus Nordborg]     

Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon

Polymorphism of the nucleotide sequences encoding 149 amino acids of linked major histocompatibility complex (Mhc) class II 131 and 132 peptides, and of the intervening intron (548–773 base pairs), was examined within and among seven Pacific salmon (Oncorhynchus) species. Levels of nucleotide diversity were higher for theB1 sequence than forB2 or the intron in comparisons both within and between species. For the codons of the peptide binding region of the BI sequence, the level of nonsynonymous nucleotide substitution (d_N) exceeded the level of synonymous substitution (d_S) by a factor of ten for within-species comparisons, and by a factor of four for between-species comparisons. The excess of d_N indicates that balancing selection maintains diversity at this salmonidMhc class II locus, as is common forMhc loci in other vertebrates. Levels of nucleotide diversity for both the exon and intron sequences were greater among than within species, and there were numerous species-specific nucleotides present in both the coding and noncoding regions. Thus, neighbor-joining analysis of both the intron and exon regions provided phylogenies in which the sequences clustered strongly by species. There was little evidence of shared ancestral (trans-species) polymorphism in the exon phylogeny, and the intron phylogeny depicted standard relationships among the Pacific salmon species. The lack of shared allelicB1 lineages in these closely related species may result from severe bottlenecks that occurred during speciation or during the ice ages that glaciated the rim of the north Pacific Ocean approximately every 100 000 years in the Pleistocene.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers U34692-U34720     

Germline repertoire of the immunoglobulin V H 3 family in rhesus monkeys

 To facilitate molecular studies of antibody responses in rhesus monkeys (Macaca mulatta), we cloned and sequenced germline segments from its largest and most diverse immunoglobulin heavy-chain gene family, V _H 3. Using a PCR-based approach, we characterized 29 sequences, 20 with open reading frames (ORFs) and 9 pseudogenes. The leader sequences, introns, exons, and recombination signal sequences of M. mulatta V _H 3 gene segments are not strictly identical to those of humans, but the mature coding regions demonstrate, on average, greater than 90% sequence similarity. Although the framework regions are more highly conserved, the complementarity-determining regions (CDRs) also show strong similarities, and their predicted three-dimensional structures resemble those of their human homologues. In one instance, homologous macaque and human CDR1 sequences were 100% identical at the nucleotide level, and some CDR2s shared nucleotide identity as high as 96.5%. However, some rhesus V _H 3 ORFs have unusual structural features, including atypical CDR lengths and uncommon amino acids at structurally crucial positions. The similarity of rhesus and human V _H 3 homologues reinforces the notion that humoral immunity in this nonhuman primate species is an appropriate system for modeling human antibody responses. Received: 10 August 1999 / Revised: 30 December 1999     

Expression Patterns of Duplicate Genes in the Developing Root in Arabidopsis thaliana

Data on gene expression in the development of the root in Arabidopsis thaliana were used to test for expression profile differences among multi-gene families and to examine the extent to which expression differences accompanied coding sequences divergence within families. Significant differences among families were observed on two principal axes, accounting for over 80% of the variance in the expression data. The number of synonymous nucleotide substitutions per synonymous site (d_S) and the number of nonsynonymous nucleotide substitutions per nonsynonymous site (d_N) were estimated between the members of two-member families (N=428) and between phylogenetically independent sister pairs (N=190) of sequences within larger families. Ribosomal proteins and a few other proteins were exceptional in showing highly divergent expression patterns in spite of very low levels of amino acid sequence divergence, as indicated by the low d_N relative to d_S. However, the majority of gene duplicates showed relatively high levels of amino acid sequence divergence without appreciable change in expression pattern in the cell types analyzed. Reviewing Editor:Dr. Manyuan Long     

Comprehensive analysis and characterization of the TCR α chain sequences in the common marmoset

The common marmoset (Callithrix jacchus) is useful as a nonhuman primate model of human diseases. Although the marmoset model has great potential for studying autoimmune diseases and immune responses against pathogens, little information is available regarding the genes involved in adaptive immunity. Here, we identified one TCR α constant (TRAC), 46 TRAJ (joining), and 35 TRAV (variable) segments from marmoset cDNA. Marmoset TRAC, TRAJ, and TRAV shared 80%, 68–100%, and 79–98% identity with their human counterparts at the amino acid level, respectively. The amino acid sequences were less conserved in TRAC than in TCRβ chain constant (TRBC). Comparative analysis of TRAV between marmosets and humans showed that the rates of synonymous substitutions per site (d _S ) were not significantly different between the framework regions (FRs) and complementarity determining regions (CDRs), whereas the rates of nonsynonymous substitutions per site (d _N ) were significantly lower in the FRs than in CDRs. Interestingly, the d _N values of the CDRs were greater for TRBV than TRAV. These results suggested that after the divergence of Catarrhini from Platyrrhini, amino acid substitutions were decreased in the FRs by purifying selection and occurred more frequently in CDRβ than in CDRα by positive selection, probably depending on structural and functional constraints. This study provides not only useful information facilitating the investigation of adaptive immunity using the marmoset model but also new insight into the molecular evolution of the TCR heterodimer in primate species.     

Nucleotide diversity and molecular evolution of the PSY1 gene in Zea mays compared to some other grass species

Phytoene synthase (PSY), which is encoded by the phytoene synthase 1 (PSY1) gene, is the first rate-limiting enzyme in the plant carotenoid biosynthetic pathway. In order to examine the genetic diversity and evolution pattern of PSY1 within the Andropogoneae, sequences of 76 accessions from 5 species (maize, teosinte, tripsacum, coix, and sorghum) of the Andropogoneae were tested, along with 4 accessions of rice (Oryza sativa L.) included as outliers. Both the number and the order of exons and introns were relatively conserved across the species tested. Three domains were identified in the coding sequence, including signal peptide (SP), PSY, and highly conserved squalene synthase (SQS) domain. Although no positive selection signal was detected at an overall coding level among all species tested, the SP domain and the region upstream of the SQS–PSY domain appear to have undergone rapid evolution, as evidenced by a high d _N/d _S ratio (>1.0). At the nucleotide level, positive selection and balancing selection were detected only among the yellow maize germplasm and the white maize germplasm, respectively. The phylogenetic tree based on full-length sequences of PSY1-like regions supported the monophyletic theory of the Andropogoneae and the closest relationship between Zea and Tripsacum among the Andropogoneae. Coix, which was theorized to have a closer relationship with maize due to similarities in morphology and chromosome number, has been shown in this study to have diverged relatively early from the other Andropogoneae, including maize.     

Primate evolution at the DNA level and a classification of hominoids

Summary The genetic distances among primate lineages estimated from orthologous noncoding nucleotide sequences of -type globin loci and their flanking and intergenic DNA agree closely with the distances (delta T₅₀H values) estimated by cross hybridization of total genomic single-copy DNAs. These DNA distances and the maximum parsimony tree constructed for the nucleotide sequence orthologues depict a branching pattern of primate lineages that is essentially congruent with the picture from phylogenetic analyses of morphological characters. The molecular evidence, however, resolves ambiguities in the morphological picture and provides an objective view of the cladistic position of humans among the primates. The molecular data group humans with chimpanzees in subtribe Hominina, with gorillas in tribe Hominini, orangutans in subfamily Homininae, gibbons in family Hominidae, Old World monkeys in infraorder Catarrhini, New World monkeys in semisuborder Anthropoidea, tarsiers in suborder Haplorhini, and strepsirhines (lemuriforms and lorisiforms) in order Primates. A seeming incongruency between organismal and molecular levels of evolution, namely that morphological evolution appears to have speeded up in higher primates, especially in the lineage to humans, while molecular evolution has slowed down, may have the trivial explanation that relatively small genetic changes may sometimes result in marked phenotypic changes.     

Structural and evolutionary comparisons of four alleles of the mouse immunoglobulin kappa chain gene,Igk-VSer

The mouseIgK-VSer gene encodes an immunoglobulin light chain variable region which gives rise to two phenotypic polymorphisms of mouse chains. The nucleotide sequences of coding and flanking regions of theIgk-VSer ^c andIgk-VSer ^d alleles found in recently inbred strains of wild mice are compared with those of theIgk-VSer ^a andIgk-VSer ^b alleles described previously. Results suggest that the gene is evolving randomly and that framework 2 and complentarity determining region 2 are preserved, presumably for overall light chain structure. Results indicate that all four allels have an octamer motif upstream of the gene which should be functional and allow prediction of whether or not the product of the germ line gene will be detectable as either the I_B-peptide or Ef1^a phenotypic polymorphism. Southern hybridization of genomic DNA using as probe a 1-kbXba I-Xba I fragment located approximately 4 kb upstream of the BALB/cIgk-VSer ^b coding region demonstrated the presence of homologous DNA in mice bearing theIgk-VSer ^a allele and absence from mice bearing theIgk-VSer ^c andIgk-VSer ^d alleles. Nucleotide sequence comparison of BALB/c and SK/CamRk (Igk-VSer ^d ) DNA in this region demonstrated that BALB/c contained an insertion 2.4 kb in length which was absent from SK/CamRk. Both strains contain DNA homologous to the reverse complement of the mouse Bam5 repetitive element at the point of the insertion, with BALB/c containing approximately 70 nucleotides more of the element than SK/CamRk. Surprisingly, the strains containing DNA related to theXba I-Xba I probe are not those determined to be the most similar by nucleotide sequence comparisons and by the Phylogenetic Analysis Using Parsimony program. The evolutionary relationship of the alleles and a possible basis for the inconsistency presented by theXba I-Xba I fragment-related DNA are discussed.     

ASPM and the Evolution of Cerebral Cortical Size in a Community of New World Monkeys

The ASPM (abnormal spindle-like microcephaly associated) gene has been proposed as a major determinant of cerebral cortical size among primates, including humans. Yet the specific functions of ASPM and its connection to human intelligence remain controversial. This debate is limited in part by a taxonomic focus on Old World monkeys and apes. Here we expand the comparative context of ASPM sequence analyses with a study of New World monkeys, a radiation of primates in which enlarged brain size has evolved in parallel in spider monkeys (genus Ateles) and capuchins (genus Cebus). The primate community of Costa Rica is perhaps a model system because it allows for independent pairwise comparisons of smaller- and larger-brained species within two taxonomic families. Accordingly, we analyzed the complete sequence of exon 18 of ASPM in Ateles geoffroyi, Alouatta palliata, Cebus capucinus, and Saimiri oerstedii. As the analysis of multiple species in a genus improves phylogenetic reconstruction, we also analyzed eleven published sequences from other New World monkeys. Our exon-wide, lineage-specific analysis of eleven genera and the ratio of rates of nonsynonymous to synonymous substitutions (d_N/d_S) on ASPM revealed no detectable evidence for positive selection in the lineages leading to Ateles or Cebus, as indicated by d_N/d_S ratios of <1.0 (0.6502 and 0.4268, respectively). Our results suggest that a multitude of interacting genes have driven the evolution of larger brains among primates, with different genes involved in this process in different encephalized lineages, or at least with evidence for positive selection not readily apparent for the same genes in all lineages. The primate community of Costa Rica may serve as a model system for future studies that aim to elucidate the molecular mechanisms underlying cognitive capacity and cortical size.     

Patterns of nucleotide substitutions inferred from the phylogenies of the class I major histocompatibility complex genes

Summary Patterns of nucleotide substitutions in human major histocompatibility complex (MHC) class I genes were estimated by using phylogenetic trees of DNA sequences. The pattern is defined as a set of 12 parameters, each of which represents the relative frequency of substitutions from a particular nucleotide to another. The pattern at the antigen recognition sites (ARS) in functional MHC genes was remarkably different from that at the remaining coding region (non-ARS). In particular, the proportion of transitions among all the nucleotide substitutions (P _s) was extremely low at the third codon positions of ARS. In the HLA-A genes, P _s at the third codon positions was only 6% in ARS, whereas it was 69% in non-ARS. In HLA-B, the corresponding values were 30% in ARS and 80% in non-ARS, respectively. On the other hand, P _s in a class I pseudogene (HLA-H) was 57%, which was in good agreement with P _s in other pseudogenes. Because pseudogenes are selectively neutral, the pattern in pseudogenes is regarded as the pattern of spontaneous substitution mutations. In general, the pattern in functional genes that are subject to selective forces deviates from the pattern in pseudogenes. At the third codon positions in coding regions, transitions scarcely cause amino acid replacements, whereas about half of transversions do cause replacements. Accordingly, P _s at the third codon positions decreases if amino acid replacements are accelerated by natural selection but increases if amino acids are conserved by functional constraint. Our observations imply that the ARS region is subject to natural selection favoring amino acid replacements, whereas the non-ARS region is subject to functional constraint. Offprint requests to: T. Gojobori     

Detecting positive selection within genomes: the problem of biased gene conversion

The identification of loci influenced by positive selection is a major goal of evolutionary genetics. A popular approach is to perform scans of alignments on a genome-wide scale in order to find regions evolving at accelerated rates on a particular branch of a phylogenetic tree. However, positive selection is not the only process that can lead to accelerated evolution. Notably, GC-biased gene conversion (gBGC) is a recombination-associated process that results in the biased fixation of G and C nucleotides. This process can potentially generate bursts of nucleotide substitutions within hotspots of meiotic recombination. Here, we analyse the results of a scan for positive selection on genes on branches across the primate phylogeny. We show that genes identified as targets of positive selection have a significant tendency to exhibit the genomic signature of gBGC. Using a maximum-likelihood framework, we estimate that more than 20 per cent of cases of significantly elevated non-synonymous to synonymous substitution rates ratio (d_N/d_S), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high d_N/d_S (more than 2). Our results indicate that gBGC significantly affects the evolution of coding sequences in primates, often leading to patterns of evolution that can be mistaken for positive selection.     

Distribution of nucleotide differences between two randomly chosen cistrons in a subdivided population: The finite island model

Using the island model of finite size, the distributions as well as the means and variances of d_w and d_b are obtained, where d_w is the number of nucleotide differences between two cistrons randomly chosen from the same colony and d_b is the corresponding number between two cistrons randomly chosen from different colonies. The rate for the means to approach equilibrium is independent of mutation while that for the variances is somewhat retarded by mutation. At the steady state, the mean of d_w is independent of population subdivision and migration rate, as long as there is migration. It has been shown that the actual genic variation in a colony may be much larger than that revealed by the heterozygosity in the colony.     

C4 gene polymorphism in primates: evolution,generation, and Chido and Rodgers antigenicity

Eleven new C4d genomic primate sequences of the fourth complement factor (C4) have been obtained. Seven of them belong to five species not yet explored for this gene: Pan paniscus (pygmy chimpanzee), Cercopithecus aethiops (green monkey), Macaca mulatta (rhesus monkey), Macaca fascicularis (cynomolgus), and Saguinus oedipus (cotton top tamarin). The New World monkeys (tamarins, four individuals) sequenced for C4 have a single C4d sequence only, which shows a B isotypic specificity and a Rodgers 3 (Rg3), Chido 1 (Ch1) antigenicity. Rg3 and Ch1 could thus be the oldest Rg/Ch specificity (at least 50 million years old) and Rg1, Rg2, Ch3, and Ch6 could be more recent human-specific antigens. Mechanisms of C4d polymorphism generation were analyzed by compiling all the presently available sequences. Examples of both point mutations and crossing-over events among C4d primate sequences could be detected. The problem of a possible trans-species inheritance of C4d polymorphism was addressed and two apparently contradicting dendrograms were obtained. One of them, constructed by using both exon and intron sequences, does not support trans-species evolution, but supports the proposed theory of extensive homogenization of the C4 genes occurring within each species, because alleles from each primate species cluster together. Another completely different dendrogram, obtained by using exon sequences only, suggests the existence of trans-species evolution for C4d polymorphism, because alleles belonging to different species cluster together in a way similar to that found for HLA class I or II alleles. However, orangutan sequences group together in both kinds of C4d sequence dendrograms and seem to have arisen from an ancestor different from that of chimpanzee, gorilla and man C4d sequences. Finally, further data have been obtained that support trans-species conservation of A-ness and B-ness and the existence of trans-specifically conserved allelic motifs, both in intronic and exonic sequences.     

Rapid Evolution of the Ribonuclease A Superfamily: Adaptive Expansion of Independent Gene Clusters in Rats and Mice

The two eosinophil ribonucleases, eosinophil-derived neurotoxin (EDN/RNase 2) and eosinophil cationic protein (ECP/RNase 3), are among the most rapidly evolving coding sequences known among primates. The eight mouse genes identified as orthologs of EDN and ECP form a highly divergent, species-limited cluster. We present here the rat ribonuclease cluster, a group of eight distinct ribonuclease A superfamily genes that are more closely related to one another than they are to their murine counterparts. The existence of independent gene clusters suggests that numerous duplications and diversification events have occurred at these loci recently, sometime after the divergence of these two rodent species (∼10–15 million years ago). Nonsynonymous substitutions per site (d _N) calculated for the 64 mouse/rat gene pairs indicate that these ribonucleases are incorporating nonsilent mutations at accelerated rates, and comparisons of nonsynonymous to synonymous substitution (d _N / d _S) suggest that diversity in the mouse ribonuclease cluster is promoted by positive (Darwinian) selection. Although the pressures promoting similar but clearly independent styles of rapid diversification among these primate and rodent genes remain uncertain, our recent findings regarding the function of human EDN suggest a role for these ribonucleases in antiviral host defense. Received: 8 April 1999 / Accepted: 22 June 1999     

Molecular organization of the Escherichia coli gab cluster: nucleotide sequence of the structural genes gabD and gabP and expression of the GABA permease gene

We have determined the nucleotide sequences of two structural genes of the Escherichia coli gab cluster, which encodes the enzymes of the 4-aminobutyrate degradation pathway: gabD, coding for succinic semialdehyde dehydrogenase (SSDH, EC 1.2.1.16) and gabP, coding for the 4-aminobutyrate (GABA) transport carrier (GABA permease). We have previously reported the nucleotide sequence of the third structural gene of the cluster, gabT, coding for glutamate: succinic semialdehyde transaminase (EC 2.6.1.19). All three gab genes are transribed unidirectionally and their orientation within the cluster is 5-gabD-gabT-gabP-3. gabT and gabP are separated by an intergenic region of 234-bp, which contains three repetetive extragenic palindromic (REP) sequences. The gabD gene consists of 1,449 nucleotides specifying a protein of 482 amino acids with a molecular mass of 51.7 kDa. The protein shows significant homologies to the NAD⁺-dependent aldehyde dehydrogenase (EC 1.2.1.3) from Aspergillus nidulans and several mammals, and to the tumor associated NADP⁺-dependent aldehyde dehydrogenase (EC 1.2.1.4) from rat. The permease gene gabP comprises 1,401 nucleotides coding a highly hydrophobic protein of 466 amino acids with a molecular mass of 51.1 kDa. The GABA permease shows features typical for an integral membrane protein and is highly homologous to the aromatic acid carrier from E. coli, the proline, arginine and histidine permeases from Saccharomyces cerevisiae and the proline transport protein from A. nidulans. Uptake of GABA was increased ca. 5-fold in transformants of E. coli containing gabP plasmids. Strong overexpression of the gabP gene under control of the isopropyl-2-d-thiogalactoside (IPTG) inducible tac promoter, however, resulted in a severe growth inhibition of the transformed strains. The GABA carrier was characterized using moderately overexpressing transformants. The K _m of GABA uptake was found to be 11.8 M and the V_max 0.33 nmol/min · mg cells. Uptake of GABA was stimulated by ammonium sulfate and abolished by 2,4-dinitrophenol. Aspartate competed with GABA for uptake.     

Amino Acid Repeats Cause Extraordinary Coding Sequence Variation in the Social Amoeba Dictyostelium discoideum

Protein sequences are normally the most conserved elements of genomes owing to purifying selection to maintain their functions. We document an extraordinary amount of within-species protein sequence variation in the model eukaryote Dictyostelium discoideum stemming from triplet DNA repeats coding for long strings of single amino acids. D. discoideum has a very large number of such strings, many of which are polyglutamine repeats, the same sequence that causes various human neurological disorders in humans, like Huntington’s disease. We show here that D. discoideum coding repeat loci are highly variable among individuals, making D. discoideum a candidate for the most variable proteome. The coding repeat loci are not significantly less variable than similar non-coding triplet repeats. This pattern is consistent with these amino-acid repeats being largely non-functional sequences evolving primarily by mutation and drift.     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Nature of allelic sequence polymorphism at the KIR3DL3 locus

KIR3DL3 is a framework gene of the Leukocyte Receptor Complex, present in all individuals and haplotypes analysed to date. We describe 17 novel KIR3DL3 alleles, including seven single nucleotide polymorphic (SNP) positions within the coding region. Sequence variation within introns included a VNTR within intron 1. As KIR3DL3 mRNA is known to be expressed in decidual NK cells, we investigated the impact of KIR3DL3 allelic variation on pre-eclampsia. No statistical difference in allele frequency or polymorphism was observed between pre-eclampsia patient and control cohorts. Linkage disequilibrium (LD) analysis of exonic SNPs suggested that recombination may be a mechanism of generating sequence diversity within KIR3DL3. A potential recombination hotspot was located within intron 5. A strong LD was detected between polymorphism in exon 6 of KIR3DL3 and the KIR gene −2DL3 or -2DS2 loci, which define the centromeric end of two main haplotypes (A and B) of the KIR cluster. Comparison of primate KIR sequences indicated that the Ig domains of KIR3DL3 are highly conserved between chimpanzee, gorilla and humans. Investigation of KIR3DL3 dN/dS ratios indicated a greater level of synonymous mutations consistent with purifying selection, although positive selection was detected acting on two sites within the stem region.Electronic supplementary material Supplementary material is available for this article at and is accessible for authorized users.