Gene duplications and losses among vertebrate deoxyribonucleoside kinases of the non-TK1 Family

ABSTRACT

Deoxyribonucleoside kinases (dNKs) salvage deoxyribonucleosides (dNs) and catalyze the rate limiting step of this salvage pathway by converting dNs into corresponding monophosphate forms. These enzymes serve as an excellent model to study duplicated genes and their evolutionary history. So far, among vertebrates only four mammalian dNKs have been studied for their substrate specificity and kinetic properties. However, some vertebrates, such as fish, frogs, and birds, apparently possess a duplicated homolog of deoxycytidine kinase (dCK). In this study, we characterized a family of dCK/deoxyguanosine kinase (dGK)-like enzymes from a frog Xenopus laevis and a bird Gallus gallus. We showed that X. laevis has a duplicated dCK gene and a dGK gene, whereas G. gallus has a duplicated dCK gene but has lost the dGK gene. We cloned, expressed, purified, and subsequently determined the kinetic parameters of the dCK/dGK enzymes encoded by these genes. The two dCK enzymes in G. gallus have broader substrate specificity than their human or X. laevis counterparts. Additionally, the duplicated dCK enzyme in G. gallus might have become mitochondria. Based on our study we postulate that changing and adapting substrate specificities and subcellular localization are likely the drivers behind the evolution of vertebrate dNKs.     

Phylogeny reconstruction using duplicate genes

In this paper, we propose a new method (uninode coding) for coding duplicate (paralogous) genes to infer species trees. Uninode coding incorporates data from duplicated and unduplicated gene copies in phylogenetic analyses of taxa. Uninode coding utilizes global parsimony through the inclusion of both duplicated and unduplicated gene copies, allows one to code all data sources from a taxon into a single terminal, and overcomes problems of character dependence among duplicated and unduplicated gene copies. We present an example of uninode coding using the phytochrome A and phytochrome C data from a study by Donoghue and Mathews.     

Population-genetic models of the fates of duplicate genes

The ultimate fate of a duplicated gene is that it either silenced through inactivating mutations or both copies are maintained by selection. This later fate can occur via neofunctionalization wherein one copy acquires a new function or by subfunctionalization wherein the original function of the gene is partitioned across both copies. The relative probabilities of these three different fates involve often very subtle iterations between of population size, mutation rate, and selection. All three of these fates are critical to the expansion and diversification of gene families.     

A three-step model for the rearrangement of the chloroplast trnK-psbA region of the gymnosperm Pinus contorta.

A region of the Pinus contorta chloroplast genome which contains a duplication of the psbA gene was characterized. From previous experiments it was known that the two copies of the psbA gene were located approximately 3.3 kilobase pairs (kbp) apart, that they had the same orientation and that one endpoint of the duplication was 19 base pairs (bp) downstream of the psbA stop codon. In order to determine the size and additional genetic content of the duplicated segment, both copies as well as the intervening DNA were sequenced completely. It was found that the duplicated segment was 1969 bp long, that the two copies were completely identical and were separated by 2431 bp. The duplicated segment carried, in addition to psbA, the 3' exon of the trnK gene, which was partially included in a 124 bp direct repeat. The translocated copy of the duplicated segment was found to be inserted upstream of the trnK(UUU) gene and was immediately followed by a repeated 41 bp stretch from the psbA coding region. The trnK gene was split by a 2509 bp intron which contained an open reading frame of 515 codons. Sequence comparisons of the duplicated segment and its flanking DNA to the corresponding regions of P. sylvestris, a species which lacks the rearrangements found in P. contorta, made it possible to identify 3-9 bp homologies within which recombinations had occurred. A model was derived which would accommodate the conversion of a trnK-psbA locus of the ancestral P. sylvestris-like organization into the rearranged structure found in P. contorta.     

Molecular evolution of a duplication: the sex-peptide (Acp70A) gene region of Drosophila subobscura and Drosophila madeirensis

In Drosophila melanogaster, the Acp70A gene, which is involved in the postmating reactions of the female, is a single-copy gene. However, in Drosophila subobscura, the gene is duplicated and both copies are transcribed. To study the molecular evolution of the duplication, a 2.1- kb fragment encompassing both copies of the duplication was sequenced for 10 lines of D. subobscura and one line of Drosophila madeirensis. Estimates of the divergence between the two copies of the duplicated region and between the two species studied, D. subobscura and D. madeirensis, revealed that both copies of the Acp70a gene had evolved independently since their duplication. The ratio of nonsynonymous to silent divergence between copies was generally higher than one. The McDonald and Kreitman test revealed an excess of nonsynonymous changes fixed since the duplication and before the split of the D. subobscura and D. madeirensis lineages. These results point to natural selection driving protein evolution after the duplication. Specifically, adaptive evolution appears to have caused the initial differentiation between copies of the N-terminal parts of the proteins, while purifying selection could be responsible for the high conservation of the C- terminal parts.      

A mammalian protein kinase with potential for serine/threonine and tyrosine phosphorylation is related to cell cycle regulators.

In a screen of mouse erythroleukemia cDNA expression libraries with anti-phosphotyrosine antibodies, designed to isolate tyrosine kinase coding sequences, we identified several cDNAs encoding proteins identical or very similar to known protein-tyrosine kinases. However, two frequently isolated cDNAs, clk and nek, encode proteins which are most closely related to protein kinases involved in regulating progression through the cell cycle, and contain motifs generally considered diagnostic of protein-serine/threonine kinases. The clk gene product contains a C-terminal cdc2-like kinase domain, most similar to the FUS3 catalytic domain. The Clk protein, expressed in bacteria, becomes efficiently phosphorylated in vitro on tyrosine as well as serine/threonine, and phosphorylates the exogenous substrate poly(glu, tyr) on tyrosine. Direct biochemical evidence indicates that both protein-tyrosine and protein-serine/threonine kinase activities are intrinsic to the Clk catalytic domain. These results suggest the existence of a novel class of protein kinases, with an unusual substrate specificity, which may be involved in cell cycle control.     

Complex arrangement of genes within a 220-kb region of double-duplicated DNA on human 2q37.1

Gene duplication events are followed by divergence of initially identical gene copies, due to the subsequent accumulation of mutations. These mutations tend to be degenerative and may lead to either nonfunctionalization or subfunctionalization of the gene copies. Here we report the molecular characterization of a 220-kb genomic DNA fragment from human 2q37.1, in which a double duplication and a partial triplication event has taken place. As a result, this region contains four copies of alkaline phosphatase (P), four copies of the ECEL1 gene (X), two copies of a newly identified gene (N), and two copies of a cholinergic receptor subunit (R), in the order N-P-X-P-X-P-X-N-P-X-R-R. While three of the four ECEL1 copies, one copy of the phosphatase gene and one copy of the newly identified gene have lost their function, three phosphatase gene copies and the two receptor subunits are still functionally active and thus may provide an example for subfunctionalization of duplicated genes.     

Differential selection after duplication in mammalian developmental genes

Gene duplication provides the opportunity for subsequent refinement of distinct functions of the duplicated copies. Either through changes in coding sequence or changes in regulatory regions, duplicate copies appear to obtain new or tissue-specific functions. If this divergence were driven by natural selection, we would expect duplicated copies to have differentiated patterns of substitutions. We tested this hypothesis using genes that duplicated before the human/mouse split and whose orthologous relations were clear. The null hypothesis is that the number of amino acid changes between humans and mice was distributed similarly across different paralogs. We used a method modified from Tang and Lewontin to detect heterogeneity in the amino acid substitution pattern between those different paralogs. Our results show that many of the paralogous gene pairs appear to be under differential selection in the human/mouse comparison. The properties that led to diversification appear to have arisen before the split of the human and mouse lineages. Further study of the diverged genes revealed insights regarding the patterns of amino acid substitution that resulted in differences in function and/or expression of these genes. This approach has utility in the study of newly identified members of gene families in genomewide data mining and for contrasting the merits of alternative hypotheses for the evolutionary divergence of function of duplicated genes.     

The coalescent with selection on copy number variants

We develop a coalescent-based simulation tool to generate patterns of single nucleotide polymorphisms (SNPs) in a wide region encompassing both the original and duplicated genes. Selection on the new duplicated copy and interlocus gene conversion between the two copies are incorporated. This simulation enables us to explore how selection on duplicated copies affects the pattern of SNPs. The fixation of an advantageous duplicated copy causes a strong reduction in polymorphism not only in the duplicated copy but also in its flanking regions, which is a typical signature of a selective sweep by positive selection. After fixation, polymorphism gradually increases by accumulating neutral mutations and eventually reaches the equilibrium value if there is no gene conversion. When gene conversion is active, the number of SNPs in the duplicated copy quickly increases by transferring SNPs from the original copy; therefore, the time when we can recognize the signature of selection is decreased. Because this effect of gene conversion is restricted only to the duplicated region, more power to detect selection is expected if a flanking region to the duplicated copy is used.     

Alternate use of divergent forms of an ancient exon in the fructose-1,6-bisphosphate aldolase gene of Drosophila melanogaster.

The fructose-1,6-bisphosphate aldolase gene of Drosophila melanogaster contains three divergent copies of an evolutionarily conserved 3' exon. Two mRNAs encoding aldolase contain three exons and differ only in the poly(A) site. The first exon is small and noncoding. The second encodes the first 332 amino acids, which form the catalytic domain, and is homologous to exons 2 through 8 of vertebrates. The third exon encodes the last 29 amino acids, thought to control substrate specificity, and is homologous to vertebrate exon 9. A third mRNA substitutes a different 3' exon (4a) for exon 3 and encodes a protein very similar to aldolase. A fourth mRNA begins at a different promoter and shares the second exon with the aldolase messages. However, two exons, 3a and 4a, together substitute for exon 3. Like exon 4a, exon 3a is homologous to terminal aldolase exons. The exon 3a-4a junction is such that exon 4a would be translated in a frame different from that which would produce a protein with similarity to aldolase. The putative proteins encoded by the third and fourth mRNAs are likely to be aldolases with altered substrate specificities, illustrating alternate use of duplicated and diverged exons as an evolutionary mechanism for adaptation of enzymatic activities.     

Docking domains and substrate-specificity determination for MAP kinases

Signalling specificity in eukaryotic cells is maintained by several mechanisms. One mechanism by which mitogen-activated protein (MAP) kinases ensure their specificity of action is by interacting with their substrates through docking domains. These docking domains recruit the kinases to the correct substrates and enhance their fidelity and efficiency of action. Additional specificity determinants in the substrates serve to enhance the specificity of substrate phosphorylation by MAP kinases further.     

How Does the Cell Count the Number of Ectopic Copies of a Gene in the Premeiotic Inactivation Process Acting in Ascobolus Immersus?

Repeated genes, artificially introduced in Ascobolus immersus by integrative transformation, are frequently inactivated during the sexual phase. Inactivation is observed in about 50% of meioses if duplicated genes are at ectopic chromosomal locations, and in 90% of meioses if genes are tandemly repeated. Inactivation is associated with extensive methylation of the cytosine residues of the duplicated sequences and is induced in the still haploid nuclei of the dikaryotic cell which will undergo karyogamy and subsequent meiosis. Only repeated sequences become methylated. This raises the intriguing question of how the premeiotic inactivation machinery is informed that a nucleus contains multiple copies of a gene. By using in crosses recombinant strains of A. immersus in which either one, two or three genetically independent copies of the exogenous amdS gene from Aspergillus nidulans had been introduced, we could follow the premeiotic inactivation of each one of the ectopic amdS copies. This led us to propose that a prerequisite for inactivation is a premeiotic pairing of repeated sequences and that each copy can undergo successive cycles of pairing. In fact, once methylated, a copy can pair with a still unmethylated copy, so that an uneven number of copies can be subject to inactivation.     

The Evolutionary Dynamics of Genetic Incompatibilities Introduced by Duplicated Genes in Arabidopsis thaliana

Although gene duplications provide genetic backup and allow genomic changes under relaxed selection, they may potentially limit gene flow. When different copies of a duplicated gene are pseudofunctionalized in different genotypes, genetic incompatibilities can arise in their hybrid offspring. Although such cases have been reported after manual crosses, it remains unclear whether they occur in nature and how they affect natural populations. Here, we identified four duplicated-gene based incompatibilities including one previously not reported within an artificial Arabidopsis intercross population. Unexpectedly, however, for each of the genetic incompatibilities we also identified the incompatible alleles in natural populations based on the genomes of 1,135 Arabidopsis accessions published by the 1001 Genomes Project. Using the presence of incompatible allele combinations as phenotypes for GWAS, we mapped genomic regions that included additional gene copies which likely rescue the genetic incompatibility. Reconstructing the geographic origins and evolutionary trajectories of the individual alleles suggested that incompatible alleles frequently coexist, even in geographically closed regions, and that their effects can be overcome by additional gene copies collectively shaping the evolutionary dynamics of duplicated genes during population history.     

Molecular conversion of NAD kinase to NADH kinase through single amino acid residue substitution

NAD kinase phosphorylates NAD+ to form NADP+ and is strictly specific to NAD+, whereas NADH kinase phosphorylates both NAD+ and NADH, thereby showing relaxed substrate specificity. Based on their primary and tertiary structures, the difference in the substrate specificities between NAD and NADH kinases was proposed to be caused by one aligned residue: Gly or polar amino acid (Gln or Thr) in five NADH kinases and a charged amino acid (Arg) in two NAD kinases. The substitution of Arg with Gly in the two NAD kinases relaxed the substrate specificity (i.e. converted the NAD kinases to NADH kinases). The substitution of Arg in one NAD kinase with polar amino acids also relaxed the substrate specificity, whereas substitution with charged and hydrophobic amino acids did not show a similar result. In contrast, the substitution of Gly with Arg in one NADH kinase failed to convert it to NAD kinase. These results suggest that a charged or hydrophobic amino acid residue in the position of interest is crucial for strict specificity of NAD kinases to NAD+, whereas Gly or polar amino acid residue is not the sole determinant for the relaxed substrate specificity of NADH kinases. The significance of the conservation of the residue at the position in 207 NAD kinase homologues is also discussed.     

A few amino acid substitutions can convert deoxyribonucleoside kinase specificity from pyrimidines to purines

In mammals, the four native deoxyribonucleosides are phosphorylated to the corresponding monophosphates by four deoxyribonucleoside kinases, which have specialized substrate specificities. These four enzymes are likely to originate from a common progenitor kinase. Insects appear to have only one multisubstrate deoxyribonucleoside kinase (dNK, EC 2.7.1.145), which prefers pyrimidine nucleosides, but can also phosphorylate purine substrates. When the structures of the human deoxyguanosine kinase (dGK, EC 2.7.1.113) and the dNK from Drosophila melanogaster were compared, a limited number of amino acid residues were identified and proposed to be responsible for the substrate specificity. Three of these key residues in Drosophila dNK were then mutagenized and the mutant enzymes were characterized regarding their ability to phosphorylate native deoxyribonucleosides and nucleoside analogs. The mutations converted the dNK substrate specificity from predominantly pyrimidine specific into purine specific. A similar scenario could have been followed during the evolution of kinases. Upon gene duplication of the progenitor kinase, only a limited number of single amino acid changes has taken place in each copy and resulted in substrate-specialized enzymes.     

A diffusion model for the fate of tandem gene duplicates in diploids

Suppose one chromosome in one member of a population somehow acquires a duplicate copy of the gene, fully linked to the original gene's locus. Preservation is the event that eventually every chromosome in the population is a descendant of the one which initially carried the duplicate. For a haploid population in which the absence of all copies of the gene is lethal, the probability of preservation has recently been estimated via a diffusion approximation. That approximation is shown to carry over to the case of diploids and arbitrary strong selection against the absence of the gene. The techniques used lead to some new results. In the large population limit, it is shown that the relative probability that descendants of a small number of individuals carrying multiple copies of the gene fix in the population is proportional to the number of copies carried. The probability of preservation is approximated when chromosomes carrying two copies of the gene are subject to additional, fully non-functionalizing mutations, thereby modelling either an additional cost of replicating a longer genome, or a partial duplication of the gene. In the latter case the preservation probability depends only on the mutation rate to null for the duplicated portion of the gene.     

Gene function, gene networks and the fate of duplicated genes

For both copies of a duplicated gene to become fixed in a population and subsequently maintained, selection must favour individuals with both genes over individuals with one. Here I review and assess some of the proposed ways that gene structure and function might affect the likelihood of both copies acquiring distinct functions and therefore positive selection. In particular I focus on the interacting pathways of genes that make up gene networks, and how these may affect genes duplicated both singly and en masse. Using the Wnt and hedgehog pathways as examples and data from developmental and genome analyses, I show that, while some of these theories may genuinely reflect what has occurred in animal evolution, there are still insufficient data to rigorously assess their relative importance. This, however, is likely to change in the near future.     

Analysis of two genes encoding heat-labile toxins and located on a single Ent plasmid from Escherichia coli.

A clinical isolate of Escherichia coli harbored two copies of the heat-labile toxin (LT)-encoding gene (elt) on a 157-kb plasmid. The arrangement of the gene copies is different from the cholera toxin-encoding gene duplication described for some strains of Vibrio cholerae. The nucleotide sequences of the elt alleles are not identical (differing by 2 bp) and the duplicated region flanking the alleles extends 314 bp on one side and 1122 bp on the other side of each copy. Different partial copies of IS600 were identified 280 bp 3' to the stop codon of each allele. Partial and complete copies of other IS were also identified near the elt alleles. The data suggest that the regions surrounding the genes are hot spots for IS which would account for the observed heterogeneity in DNA flanking elt genes.     

Duplication and Gene Conversion in the Drosophila melanogaster Genome

Using the genomic sequences of Drosophila melanogaster subgroup, the pattern of gene duplications was investigated with special attention to interlocus gene conversion. Our fine-scale analysis with careful visual inspections enabled accurate identification of a number of duplicated blocks (genomic regions). The orthologous parts of those duplicated blocks were also identified in the D. simulans and D. sechellia genomes, by which we were able to clearly classify the duplicated blocks into post- and pre-speciation blocks. We found 31 post-speciation duplicated genes, from which the rate of gene duplication (from one copy to two copies) is estimated to be 1.0×10⁻⁹ per single-copy gene per year. The role of interlocus gene conversion was observed in several respects in our data: (1) synonymous divergence between a duplicated pair is overall very low. Consequently, the gene duplication rate would be seriously overestimated by counting duplicated genes with low divergence; (2) the sizes of young duplicated blocks are generally large. We postulate that the degeneration of gene conversion around the edges could explain the shrinkage of “identifiable” duplicated regions; and (3) elevated paralogous divergence is observed around the edges in many duplicated blocks, supporting our gene conversion–degeneration model. Our analysis demonstrated that gene conversion between duplicated regions is a common and genome-wide phenomenon in the Drosophila genomes, and that its role should be especially significant in the early stages of duplicated genes. Based on a population genetic prediction, we applied a new genome-scan method to test for signatures of selection for neofunctionalization and found a strong signature in a pair of transporter genes.