A comprehensive examination of protein sequences for evidence of internal gene duplication

Summary We have implemented a routine procedure for screening protein sequences for evidence of intragenic duplications. We tested 163 protein sequences representing 116 superfamilies of unrelated proteins. Twenty superfamilies contain proteins with internal gene duplications. The intragenic duplications detected can be divided into two major types. (1) One or more duplications of all or part of a gene produce a protein with two or several detectable regions of sequence homology. Sequences from 18 superfamilies contained this type of duplication. (2) Repeated reduplication of a small DNA segment can produce a protein that is repetitive over most of its length. Three superfamilies contain such repetitive sequences. We also investigated the limits of detection of ancient duplications using sequences derived by random mutation of a model sequence consisting of ten 10-residue repeats. The original repetitive nature of the sequence was usually detected after 250 point mutations even though the ancestral segment could not be accurately reconstructed.     

Amazonian Amerindians exhibit high variability of KIR profiles

Natural killer cell immunoglobulin-like receptors (KIRs) mediate cell lysis through the recognition of human leukocyte antigen class I complexes in target cells, playing an important role in innate immune response. In this context, disease-based selective pressures could be relevant, leaving signatures detected by population studies. However, most population studies on KIR variability have focused on Europe and Asia, while Americas, Oceania, and Africa remain poorly studied. The aim of this study was to analyze the variability of KIR genes in Amerindian tribes from the Amazon region to infer about their evolutionary history. KIR profiles were estimated in 40 individuals from six Amazonian Amerindian tribes using single specific primer polymerase chain reaction. Twenty-five different profiles were identified, and surprisingly, the haplogroup A frequency was the lowest observed in human populations (16%). Results showed also that KIR variability was higher in this group in contrast to Venezuelan Amerindians. Principal components analysis evidenced that Amerindians formed a separated group from other worldwide populations and showed a higher intraethnic differentiation in comparison to other ethnic groups. Such pattern may reflect small effective size and intense genetic drift. However, because of the role of KIR in immune response, selective pressures cannot be entirely ruled out.     

Phylogenies from amino acid sequences aligned with gaps: The problem of gap weighting

Summary The common but generally overlooked problem of how best to construct phylogenies from orthologous amino acid sequences, when their alignment requires the placement therein of gaps denoting insertions/deletions in the evolutionary history of their genes since their common ancestor, has been studied. Three diverse methods were examined: 1. each missing residue in a gap is weighted as equivalent to the average number of minimum nucleotide replacements in known conjugate amino acid pairs of those same two sequences, which weight necessarily differs for each pair of sequences; 2. each missing residue in a gap is weighted as equivalent to a fixed number of nucleotide replacements; and 3. each gap, regardless of length, is weighted as equivalent to a fixed number of nucleotide replacements. For the flavodoxins, each method yielded a different best tree and suggests that the choice of method may be crucial. For the plant ferredoxins, all methods give results inconsistent with botanical classification and suggests the sequences may not all be orthologous. For the bacterial ferredoxins, the method was less germane than the actual weight used, five different best trees being obtained depending upon the weight. The best tree for all ferredoxins (prokaryotic plus eukaryotic) combined proved to be greatly dependent upon the gap locations with several reasonable alignments yielding different best trees. They also suggest that functional equivalence may well prove to be a poor guide to which residues have a common ancestral codon. The rubredoxin sequences show that a partial internal gene duplication occurred in thePseudomonas line, probably very soon after its divergence from the other genera. Together, the results clearly indicate that the phylogenetic answer one gets may greatly depend upon how one treats the gaps but they fail to indicate what treatment may be best. This results partly from the fact that the phylogenies of the taxa represented are not known with sufficient confidence to be sure when the procedures are performing best.     

Impact of KIR/HLA ligand combinations on immune responses in malignant melanoma

Tumor growth and dissemination depend partly on the reactivity of natural killer (NK) cells and T cells expressing NK-associated receptors. Their effector functions are regulated by an array of activating and inhibitory cell surface receptors with MHC class I ligand specificity, such as the killer immunoglobulin-like receptors (KIRs). Given the extensive genomic diversity of KIRs and their HLA ligands, it is reasonable to speculate that HLA, KIR gene variations and specific KIR-ligand combinations will have an impact on disease susceptibility and/or progression. Here, we discuss how KIR genotypes and KIR/HLA immunogenetic profiles may be involved in tumorigenesis, especially in malignant melanoma (MM). A hypothetical model of the impact of KIR/ligand combinations on immune responses in MM is proposed.     

An approach to delineate primers for a group of poorly conserved sequences incorporating the common motif region

Glutathione synthetase (gshB) has previously been reported to confer tolerance to acidic soil condition in Rhizobium species. Cloning the gene coding for this enzyme necessitates the designing of proper primer sets which in turn depends on the identification of high quality sequence similarity in multiple global alignments. In this experiment, a group of homologous gene sequences related to gshB gene (accession no: gi-86355669:327589-328536) of Rhizobium etli CFN 42, were extracted from NCBI nucleotide sequence databases using BLASTN and were analyzed for designing degenerate primers. However, the T-coffee multiple global alignment results did not show any block of conserved region for the above sequence set to design the primers. Therefore, we attempted to identify the location of common motif region based on multiple local alignments employing the MEME algorithm supported with MAST and Primer3. The results revealed some common motif regions that enabled us to design the primer sets for related gshB gene sequences. The result will be validated in wet lab.     

Methods for assessing gene content diversity of KIR with examples from a global set of populations

A number of statistical methods are widely used to describe allelic variation at specific genetic loci and its implication on the evolutionary history of these loci. Although the methods were developed primarily to study allelic variation at loci that are virtually always present in the genome, they are often applied to data of gene content variation (i.e., presence/absence of multiple homologous genes) at the killer cell immunoglobulin-like receptor (KIR) gene cluster. In this paper, we discuss methodological issues involved in the analysis of gene content variation data in the KIR region and also its covariation with polymorphism at the human leukocyte antigen class I loci, which encode ligands for KIR. A comparison of several statistical methods and measures (gene frequency, haplotype frequency, and linkage disequilibrium estimation) using the Centre d’Etude du Polymorphisme Humain data will be provided using KIR haplotypes that have been determined by segregation analysis, noting the strengths and weaknesses of the methods when only the presence/absence data is considered. Finally, application of these methods to a set of globally distributed populations is described (see Single et al., Nat Genet 39:1114–1119, 2007) in order to illustrate the challenges faced when inferring the joint effects of natural selection and demographic history on these immune-related genes.     

Evidence for balancing selection acting on KIR2DL4 genotypes in rhesus macaques of Indian origin

The interaction of killer-cell immunoglobulin-like receptors (KIR) and their respective major histocompatibility complex (MHC) ligands can alter the activation state of the natural killer (NK) cell. In both humans and rhesus macaques, particular types of non-classical MHC class I molecules are predominantly expressed on the trophoblast. In humans, human leukocyte antigen G has been demonstrated to act as a ligand for KIR2DL4, present on all NK cells, whereas Mamu-AG may execute a similar function in rhesus macaques. During primate evolution, orthologues of KIR2DL4 appear to have been highly conserved, suggesting strong purifying selection. A cohort of 112 related and unrelated rhesus macaques of mostly Indian origin were selected to study their KIR2DL4 genes for the occurrence of polymorphism. Comparison of the proximal region provided evidence for strong conservative selection acting on the exons encoding the Ig domains. As is found in humans, in the Indian rhesus macaque population, two different KIR2DL4 entities are encountered, which differ for their intra-cellular signalling motifs. One genotype contains a complex mutation in the distal region of exon 9, which negates a serine/threonine kinase site. Furthermore, both allelic entities are present in a distribution, which suggests that balancing selection is operating on these two distinct forms of KIR2DL4. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Regulation of Human Cytidine Triphosphate Synthetase 2 by Phosphorylation

Cytidine triphosphate synthetase (CTPS) is the rate-limiting enzyme in de novo CTP synthesis and is required for the formation of RNA, DNA, and phospholipids. This study determined the kinetic properties of the individual human CTPS isozymes (hCTPS1 and hCTPS2) and regulation through substrate concentration, oligomerization, and phosphorylation. Kinetic analysis demonstrated that both hCTPS1 and hCTPS2 were maximally active at physiological concentrations of ATP, GTP, and glutamine, whereas the K_m and IC₅₀ values for the substrate UTP and the product CTP, respectively, were close to their physiological concentrations, indicating that the intracellular concentrations of UTP and CTP may precisely regulate hCTPS activity. Low serum treatment increased hCTPS2 phosphorylation, and five probable phosphorylation sites were identified in the hCTPS2 C-terminal domain. Metabolic labeling of hCTPS2 with [³²P]H₃PO₄ demonstrated that Ser⁵⁶⁸ and Ser⁵⁷¹ were two major phosphorylation sites, and additional studies demonstrated that Ser⁵⁶⁸ was phosphorylated by casein kinase 1 both in vitro and in vivo. Interestingly, mutation of Ser⁵⁶⁸ (S568A) but not Ser⁵⁷¹ significantly increased hCTPS2 activity, demonstrating that Ser⁵⁶⁸ is a major inhibitory phosphorylation site. The S568A mutation had a greater effect on the glutamine than ammonia-dependent activity, indicating that phosphorylation of this site may influence the glutaminase domain of hCTPS2. Deletion of the C-terminal regulatory domain of hCTPS1 also greatly increased the V_max of this enzyme. In summary, this is the first study to characterize the kinetic properties of hCTPS1 and hCTPS2 and to identify Ser⁵⁶⁸ as a major site of CTPS2 regulation by phosphorylation.     

A mutation in <Emphasis Type="Italic">KIR3DS1</Emphasis> that results in truncation and lack of cell surface expression

The KIR gene cluster exhibits a high degree of polymorphism in terms of gene content as well as allelic polymorphism, and data suggest that it is evolving rapidly. The KIR3DL1 locus is one of the most polymorphic loci within this cluster and is unique in that it encodes an activating receptor KIR3DS1, as well as multiple inhibitory KIR3DL1 allotypes. Because KIR3DS1 has been implicated in a number of diseases, we tested for the presence of KIR3DS1 variants that might affect its expression and activating capacity. Preliminary FACS analysis indicated that indeed some individuals with the KIR3DS1 allele showed no cell surface expression of the molecule. Sequencing analysis identified a variant with a complex deletion/substitution mutation in exon 4 (which encodes the D1 extracellular domain), resulting in a premature stop codon. We subsequently genotyped 3,960 unrelated individuals and determined the frequencies of this allele across geographically distinct world populations. The data indicate that the null KIR3DS1 allele is uncommon, arose on a single haplotype, and spread across geographically distinct populations. An erratum to this article can be found at     

The genomic organization and evolution of the natural killer immunoglobulin-like receptor (KIR) gene cluster

Natural killer (NK) immunoglobulin-like receptors (KIRs) are a family of polymorphic receptors which interact with specific motifs on HLA class I molecules and modulate NK cytolytic activity. In this study, we analyzed a recently sequenced subgenomic region on chromosome 19q13.4 containing eight members of the KIR receptor repertoire. Six members are clustered within a 100-kb continuous sequence. These genes include a previously unpublished member of the KIR gene family 2DS6, as well as 2DL1, 2DL4, 3DL1, 2DS4, 3DL2, from centromere to telomere. Two additional KIR genes, KIRCI and 2DL3, which may be located centromeric of this cluster were also analyzed. We show that the KIR genes have undergone repeated gene duplications. Diversification between the genes has occurred postduplication primarily as a result of retroelement indels and gene truncation. Using pre- and postduplication Alu sequences identified within these genes as evolutionary molecular clocks, the evolution and duplication of this gene cluster is estimated to have occurred 30–45 million years ago, during primate evolution. A proposed model of the duplication history of the KIR gene family leading to their present organization is presented. Received: 25 November 1999 / Revised: 10 January 2000     

Association of supervillin with KIR2DL1 regulates the inhibitory signaling of natural killer cells

Inhibitory signaling is crucial in the regulation of the cytotoxicity of natural killer (NK) cells. Here, we show that KIR2DL1, an inhibitory receptor of NK cells, associates with supervillin, an F-actin binding protein. Interaction of supervillin with KIR2DL1 is dependent on the KIR2DL1 receptor stimulation and requires the phosphorylation of tyrosines in both ITIM motifs. “Knockdown” of expression of supervillin by RNA interference (RNAi) restores the KIR2DL1-suppressed cytotoxicity of NK cells. Inhibition of supervillin by RNAi also enhances the polarization of cytolytic granules (both granzyme B and perforin) to the synapse formed between YTS-GFP-KIR2DL1 NK cells and 721.221-HLA-Cw4 target cells. Further study reveals that supervillin is required for KIR2DL1-mediated inhibition of Vav1 and ERK phoshorylation. Moreover, we have found that binding of supervillin with KIR2DL1 facilitates the recruitment of SHPs especially SHP-2 to KIR2DL1 receptor. Thus, our findings demonstrate that supervillin is a novel molecule that associates with KIR2DL1 receptor and regulates the inhibitory signaling in NK cells.     

Dipeptide frequencies in proteins and the CpG deficiency in vertebrate DNA

Summary Analysis of vertebrate protein sequences totalling 4040 residues shows that amino acids with a high proportion of codons ending in C occur with significantly reduced frequency before amino acids whose codons start with G. This effect is not shown by control bacterial protein sequences. The consequent implication of shortage of XXC. GXX codon pairs in vertebrate messenger RNA is discussed in relation to the extreme rarity of the base doublet CpG in vertebrate DNA.     

Chain-terminating natural mutations affect the function of activating KIR receptors 3DS1 and 2DS3

To determine the nucleotide polymorphism of activating killer-cell immunoglobulin-like receptors (aKIR) 3DS1 and 2DS3, we developed a novel direct-sequencing method and analyzed DNA samples of 175 KIR3DS1(+) individuals and 72 KIR2DS3(+) individuals from the white population. The putative ligand-binding extracellular immunoglobulin (Ig)-like domains of these aKIR receptors are highly conserved, a scenario contrary to inhibitory KIRs that recognize polymorphic human leukocyte antigen (HLA) class I molecules. Null alleles 3DS1*049N and 2DS3*003N that do not express cell-surface receptors were discovered, and they occur commonly in whites (3DS1*049N = 2%; 2DS3*003N = 0.8%). Sequence-specific polymerase chain reaction (PCR) detecting these null alleles is negative with DNA from nonwhite subjects, suggesting that these null alleles are specific to whites and probably originated after the colonization of modern humans in Europe.     

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.     

Coding DNA sequences: statistical distributions

The base distributions in coding DNA sequences (CDS) are investigated. We explore the scaling properties of the 4-dimensional directed random walk and compare them with that for the DNA sequences. Inference from these observation are, however, contradicted by alternate analysis using factorial moments. To resolve this conflict we look directly at the nucleotide base distributions. In all the cases the base distributions change from gaussian to non-gaussian as the scale size is increased. The CDS, therefore, have nucleotide distributions different from the random.     

Genotypic diversity of the Killer Cell Immunoglobulin-like Receptors (KIR) and their HLA class I Ligands in a Saudi population

Killer Cell Immunoglobulin-like Receptors (KIR) have been used as good markers for the study of genetic predisposition in many diseases and in human genetic population dynamics. In this context, we have investigated the genetic diversity of KIR genes and their main HLA class I ligands in Saudi population and compared the data with other studies of neighboring populations. One hundred and fourteen randomly selected healthy Saudi subjects were genotyped for the presence or absence of 16 KIR genes and their HLA-C1, -C2, -Bw4^Thr80 and Bw4^Ile80 groups, using a PCR-SSP technique. The results show the occurrence of the framework genes (3DL2, 3DL3 and 2DL4) and the pseudogenes (2DP1 and 3DP1) at highest frequencies. All inhibitory KIR (iKIR) genes appeared at higher frequencies than activating genes (aKIR), except for 2DS4 with a frequency of 90.35%. A total of 55 different genotypes were observed appearing at different frequencies, where 12 are considered novel. Two haplotypes were characterized, AA and Bx (BB and AB), which were observed in 24.5% and 75.5% respectively of the studied group. The frequencies of iKIR + HLA associations were found to be much higher than aKIR + HLA. KIR genes frequencies in the Saudi population are comparable with other Middle Eastern and North African populations.