Search

The ITAM-bearing transmembrane signaling subunits (TSS) are indispensable components of activating leukocyte receptor complexes. The TSS-encoding genes map to paralogous chromosomal regions, which are thought to arise from ancient genome tetraploidization(s). To assess a possible role of tetraploidization in the TSS evolution, we studied TSS and other functionally linked genes in the amphibian species Xenopus laevis whose genome was duplicated about 40 MYR ago. We found that X. laevis has retained a duplicated set of sixteen TSS genes, all except one being transcribed. Furthermore, duplicated TCRα loci and genes encoding TSS-coupling protein kinases have also been retained. No clear evidence for functional divergence of the TSS paralogs was obtained from gene expression and sequence analyses. We suggest that the main factor of maintenance of duplicated TSS genes in X. laevis was a protein dosage effect and that this effect might have facilitated the TSS set expansion in early vertebrates.

Feline calicivirus (FCV) is a common cause of mild to severe upper respiratory tract disease (URTD) in cats. FCV strain 21223 was isolated from a kitten with severe pneumonia in a disease outbreak with unusually high mortality (35 %) that occurred in a Missouri feline colony in 1995-1996. Phylogenetic analysis of the genome sequence of strain 21223 indicated the emergence of a new FCV strain. Analysis of the full-length genome sequence of a closely related (99.5 % nucleotide identity) strain, 3786, obtained from an asymptomatic animal in the same colony four months later, showed the presence of seven amino acid substitutions, with six of them located in the VP1 capsid sequence encoded by ORF2. Comparative analysis of the E-region sequences (426-521 aa ORF2) presumably involved in virus-host cell receptor interactions did not identify amino acid substitutions unique to the virulent strain. We determined the complete genome sequences of four virus isolates that were collected in regional catteries in the months following the outbreak that were associated with different manifestations of the disease (URTD, chronic stomatitis, and gingivitis). We show that genetically distinct FCV strains were cocirculating in the area, and no apparent correlation could be made between overall sequence and observed disease.

Enzymatic deamination of bases in DNA or RNA leads to an alteration of flow of genetic information. Adenosine deaminases edit RNA (ADARs, TADs). Specialized cytidine deaminases are involved in RNA/DNA editing in lipid metabolism (APOBEC1) and in innate (APOBEC3 family) and humoral (AID) immunity. APOBEC2 is required for proper muscle development and, along with AID, was implicated in demethylation of DNA. The functions of APOBEC4, APOBEC5, and other deaminases recently discovered by bioinformatics approaches are unknown. What is the basis for the diverse biological functions of enzymes with similar enzyme structure and the same principal enzymatic reaction? AID, APOBEC1, lamprey CDA1, and APOBEC3G enzymes cause uracil DNA glycosylase-dependent induction of mutations when overproduced ectopically in bacteria or yeast. APOBEC2, on the contrary, is nonmutagenic. We studied the effects of the expression of various deaminases in yeast and bacteria. The mutagenic specificities of four deaminases, hAID, rAPOBEC1, hAPOBEC3G, and lamprey CDA1, are strikingly different. This suggests the existence of an intrinsic component of deaminase targeting. The expression of yeast CDD1 and TAD2/TAD3, human APOBEC4, Xanthomonas oryzae APOBEC5, and deaminase encoded by Micromonas sp. gene MICPUN_56782 was nonmutagenic. A lack of a mutagenic effect for Cdd1 is expected because the enzyme functions in the salvage of pyrimidine nucleotides, and it is evolutionarily distant from RNA/DNA editing enzymes. The reason for inactivity of deaminases grouped with APOBEC2 is not obvious from their structures. This can not be explained by protein insolubility and peculiarities of cellular distribution and requires further investigation.

The variable lymphocyte receptors (VLRs) of jawless vertebrates such as lamprey and hagfish are composed of highly diverse modular leucine-rich repeats. Each lymphocyte assembles a unique VLR by rearrangement of the germline gene. In the lamprey genome, we identify here about 850 distinct cassettes encoding leucine-rich repeat modules that serve as sequence templates for the hypervariable VLR repertoires. The data indicate a gene conversion-like process in VLR diversification. Genomic analysis suggested a link between the VLR and platelet glycoprotein receptors. Lamprey lymphocytes express two putative deaminases of the AID-APOBEC family that may be involved in VLR diversification, as indicated by in vitro mutagenesis and recombination assays. Vertebrate acquired immunity could have therefore originated from lymphocyte receptor diversification by an ancestral AID-like DNA cytosine deaminase.

Kindlins are essential for integrin-mediated cell adhesion. This study focuses on the evolutionary origin and subsequent functional specialization of kindlins and their role in evolutionary adaptation of cell adhesiveness in multicellular organisms., Kindlins are integrin-interacting proteins essential for integrin-mediated cell adhesiveness. In this study, we focused on the evolutionary origin and functional specialization of kindlins as a part of the evolutionary adaptation of cell adhesive machinery. Database searches revealed that many members of the integrin machinery (including talin and integrins) existed before kindlin emergence in evolution. Among the analyzed species, all metazoan lineages—but none of the premetazoans—had at least one kindlin-encoding gene, whereas talin was present in several premetazoan lineages. Kindlin appears to originate from a duplication of the sequence encoding the N-terminal fragment of talin (the talin head domain) with a subsequent insertion of the PH domain of separate origin. Sequence analysis identified a member of the actin filament–associated protein 1 (AFAP1) superfamily as the most likely origin of the kindlin PH domain. The functional divergence between kindlin paralogues was assessed using the sequence swap (chimera) approach. Comparison of kindlin 2 (K2)/kindlin 3 (K3) chimeras revealed that the F2 subdomain, in particular its C-terminal part, is crucial for the differential functional properties of K2 and K3. The presence of this segment enables K2 but not K3 to localize to focal adhesions. Sequence analysis of the C-terminal part of the F2 subdomain of K3 suggests that insertion of a variable glycine-rich sequence in vertebrates contributed to the loss of constitutive K3 targeting to focal adhesions. Thus emergence and subsequent functional specialization of kindlins allowed multicellular organisms to develop additional tissue-specific adaptations of cell adhesiveness.

Using iterative database searches, we identified a new subfamily of the AID/APOBEC family of RNA/DNA editing cytidine deaminases. The new subfamily, which is represented by readily identifiable orthologs in mammals, chicken, and frog, but not fishes, was designated APOBEC4. The zinc-coordinating motifs involved in catalysis and the secondary structure of the APOBEC4 deaminase domain are evolutionarily conserved, suggesting that APOBEC4 proteins are active polynucleotide (deoxy)cytidine deaminases. In reconstructed maximum likelihood phylogenetic trees, APOBEC4 forms a distinct clade with a high statistical support. APOBEC4 and APOBEC1 are joined in a moderately supported cluster clearly separated from AID, APOBEC2 and APOBEC3 subfamilies. In mammals, APOBEC4 is expressed primarily in testis which suggests the possibility that it is an editing enzyme for mRNAs involved in spermatogenesis.

Most of the eukaryotic protein-coding genes are interrupted by multiple introns. A substantial fraction of introns occupy the same position in orthologous genes from distant eukaryotes, such as plants and animals, and consequently are inferred to have been inherited from the common ancestor of these organisms. In contrast to these conserved introns, many other introns appear to have been gained during evolution of each major eukaryotic lineage. The mechanism(s) of insertion of new introns into genes remains unknown. Because the nucleotides that flank splice junctions are nonrandom, it has been proposed that introns are preferentially inserted into specific target sequences termed protosplice sites. However, it remains unclear whether the consensus nucleotides flanking the splice junctions are remnants of the original protosplice sites or if they evolved convergently after intron insertion. Here, we directly address the existence of protosplice sites by examining the context of introns inserted within codons that encode amino acids conserved in all eukaryotes and accordingly are not subject to selection for splicing efficiency. We show that introns are either predominantly inserted into specific protosplice sites, which have the consensus sequence (A/C)AG|Gt, or that they are inserted randomly but are preferentially fixed at such sites.

We present a classification analysis of the mutation spectra of the p53 gene and construct maps of hotspots for the germline (Li-Fraumein syndrome), different types of tumors and their derived cell lines. While spectra from solid tumors share common hotspots with the germline spectrum, they also contain unique sets of somatic hotspots that are not observed in the germline. All these hotspots correspond to amino acid replacements in the DNA-binding interface of p53. The mutation spectra of lymphomas and cell lines derived from lymphomas and lung cancers contained few hotspots compared to solid tumors. Thus, the distribution of hotspots in the p53 gene appears to depend on the tumor type and cell growth conditions; this specificity is missed by the bulk hotspot analysis. A negative correlation was detected between the amino acid replacement propensity in tumors and evolutionary variability: the hotspots are located in the positions that are highly conserved in p53 and its paralogs, p63 and p73. In all the mutation spectra, substitutions leading to amino acid replacements strongly dominate over silent substitutions, indicating that functional sites evolving under strong purifying selection are subject to intensive positive selection in p53-dependent tumors. These results are compatible with the gain-of-function concept of the role of p53 in tumorigenesis.

The mechanisms and evolutionary dynamics of intron insertion and loss in eukaryotic genes remain poorly understood. Reconstruction of parsimonious scenarios of gene structure evolution in paralogous gene families in animals and plants revealed numerous gains and losses of introns. In all analyzed lineages, the number of acquired new introns was substantially greater than the number of lost ancestral introns. This trend held even for lineages in which vertical evolution of genes involved more intron losses than gains, suggesting that gene duplication boosts intron insertion. However, dating gene duplications and the associated intron gains and losses based on the molecular clock assumption showed that very few, if any, introns were gained during the last approximately 100 million years of animal and plant evolution, in agreement with previous conclusions reached through analysis of orthologous gene sets. These results are generally compatible with the emerging notion of intensive insertion and loss of introns during transitional epochs in contrast to the relative quiet of the intervening evolutionary spans.

We have purified wild type and exonuclease-deficient four-subunit DNA polymerase epsilon (Pol epsilon) complex from Saccharomyces cerevisiae and analyzed the fidelity of DNA synthesis by the two enzymes. Wild type Pol epsilon synthesizes DNA accurately, generating single-base substitutions and deletions at average error rates of/=2 x 10-5 and/=5 x 10-7, respectively. Pol epsilon lacking 3' --5' exonuclease activity is less accurate to a degree suggesting that wild type Pol epsilon proofreads at least 92% of base substitution errors and at least 99% of frameshift errors made by the polymerase. Surprisingly the base substitution fidelity of exonuclease-deficient Pol epsilon is severalfold lower than that of proofreading-deficient forms of other replicative polymerases. Moreover the spectrum of errors shows a feature not seen with other A, B, C, or X family polymerases: a high proportion of transversions resulting from T.dTTP, T.dCTP, and C.dTTP mispairs. This unique error specificity and amino acid sequence alignments suggest that the structure of the polymerase active site of Pol epsilon differs from those of other B family members. We observed both similarities and differences between the spectrum of substitutions generated by proofreading-deficient Pol epsilon in vitro and substitutions occurring in vivo in a yeast strain defective in Pol epsilon proofreading and DNA mismatch repair. We discuss the implications of these findings for the role of Pol epsilon polymerase activity in DNA replication.

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.

Prochlorococcus marinus , the dominant photosynthetic organism in the ocean, is found in two main ecological forms: high-light-adapted genotypes in the upper part of the water column and low-light-adapted genotypes at the bottom of the illuminated layer. P. marinus SS120, the complete genome sequence reported here, is an extremely low-light-adapted form. The genome of P. marinus SS120 is composed of a single circular chromosome of 1,751,080 bp with an average G+C content of 36.4%. It contains 1,884 predicted protein-coding genes with an average size of 825 bp, a single rRNA operon, and 40 tRNA genes. Together with the 1.66-Mbp genome of P. marinus MED4, the genome of P. marinus SS120 is one of the two smallest genomes of a photosynthetic organism known to date. It lacks many genes that are involved in photosynthesis, DNA repair, solute uptake, intermediary metabolism, motility, phototaxis, and other functions that are conserved among other cyanobacteria. Systems of signal transduction and environmental stress response show a particularly drastic reduction in the number of components, even taking into account the small size of the SS120 genome. In contrast, housekeeping genes, which encode enzymes of amino acid, nucleotide, cofactor, and cell wall biosynthesis, are all present. Because of its remarkable compactness, the genome of P. marinus SS120 might approximate the minimal gene complement of a photosynthetic organism.

Noncoding DNA in the human–mouse orthologous intergenic regions contains ‘islands' of conserved sequences, the functions of which remain largely unknown. We hypothesized that some of these regions might be matrix–scaffold attachment regions, MARs (or S/MARs). MARs comprise one of the few classes of eukaryotic noncoding DNA with an experimentally characterized function, being involved in the attachment of chromatin to the nuclear matrix, chromatin remodeling and transcription regulation. To test our hypothesis, we analyzed the co-occurrence of predicted MARs with highly conserved noncoding DNA regions in human–mouse genomic alignments. We found that 11% of the conserved noncoding DNA consists of predicted MARs. Conversely, more than half of the predicted MARs co-occur with one or more independently identified conserved sequence blocks. An excess of conserved predicted MARs is seen in intergenic regions preceding 5′ ends of genes, suggesting that these MARs are primarily involved in transcriptional control.

2'-5' Oligoadenylate synthetases (OAS) are a family of enzymes, which are best known for their important role in interferon-dependent antiviral mechanisms, but are also involved in the regulation of apoptosis, cell growth and differentiation in vertebrates. These enzymes bind double-stranded RNA and catalyze the synthesis of 2'-5' oligoadenylates from ATP. Several 2'-5' oligoadenylate synthetase-like proteins, which lack the ability to synthesize 2'-5' A, have been recently identified in humans and mice; the functions of these inactivated OAS derivatives remain unknown. Examination of phylogenetic trees shows that OAS inactivation in mammals occurred on several independent occasions. Comparative sequence analysis of OAS, poly(A)-polymerases, TRF4/sigma-family polymerases, archaeal CCA-adding enzymes and uridilyltransferases from trypanosomes resulted in the identification of a C-terminal domain, which is conserved in all these enzymes and is distinct from the nucleotidyltransferase domain. Secondary structure prediction shows that this domain has a four-helix core, which is most closely related to the ATP-cone domain, a regulatory nucleotide-binding domain present in ribonucleotide reductases and several other enzymes and transcription regulators. These observations, taken together with the experimental evidence of nuclease activity in the TRF4/sigma-family of polymerases, suggest that the C-terminal domain of OAS and their homologs might have nuclease activity. The putative nuclease domain is preferentially conserved in OAS derivatives that lack an active nucleotidyltransferase domain and, as indicated by the analysis of the ratio of synonymous to non-synonymous substitutions, appears to be subject to purifying selection in these proteins. In contrast, phylogenetic analysis provided evidence of episodic positive selection in the mouse OAS-like proteins with inactivated nucleotidyltransferase domains, which suggests that some of these proteins might have distinct antiviral functions.

Human centromeres remain poorly characterized regions of the human genome despite their importance for the maintenance of chromosomes. In part this is due to the difficulty of cloning of highly repetitive DNA fragments and distinguishing chromosome-specific clones in a genomic library. In this work we report the highly selective isolation of human centromeric DNA using transformation-associated recombination (TAR) cloning. A TAR vector with alphoid DNA monomers as targeting sequences was used to isolate large centromeric regions of human chromosomes 2, 5, 8, 11, 15, 19, 21 and 22 from human cells as well as monochromosomal hybrid cells. The alphoid DNA array was also isolated from the 12 Mb human mini-chromosome DeltaYq74 that contained the minimum amount of alphoid DNA required for proper chromosome segregation. Preliminary results of the structural analyses of different centromeres are reported in this paper. The ability of the cloned human centromeric regions to support human artificial chromosome (HAC) formation was assessed by transfection into human HT1080 cells. Centromeric clones from DeltaYq74 did not support the formation of HACs, indicating that the requirements for the existence of a functional centromere on an endogenous chromosome and those for forming a de novo centromere may be distinct. A construct with an alphoid DNA array from chromosome 22 with no detectable CENP-B motifs formed mitotically stable HACs in the absence of drug selection without detectable acquisition of host DNAs. In summary, our results demonstrated that TAR cloning is a useful tool for investigating human centromere organization and the structural requirements for formation of HAC vectors that might have a potential for therapeutic applications.

Mutation frequencies vary along a nucleotide sequence, and nucleotide positions with an exceptionally high mutation frequency are called hotspots. Mutation hotspots in DNA often reflect intrinsic properties of the mutation process, such as the specificity with which mutagens interact with nucleic acids and the sequence-specificity of DNA repair/replication enzymes. They might also reflect structural and functional features of target protein or RNA sequences in which they occur. The determinants of mutation frequency and specificity are complex and there are many analytical methods for their study. This paper discusses computational approaches to analysing mutation spectra (distribution of mutations along the target genes) that include many detectable (mutable) positions. The following methods are reviewed: mutation hotspot prediction; pairwise and multiple comparisons of mutation spectra; derivation of a consensus sequence; and analysis of correlation between nucleotide sequence features and mutation spectra. Spectra of spontaneous and induced mutations are used for illustration of the complexities and pitfalls of such analyses. In general, the DNA sequence context of mutation hotspots is a fingerprint of interactions between DNA and DNA repair/ replication/modification enzymes, and the analysis of hotspot context provides evidence of such interactions.

In overlapping genes, the same DNA sequence codes for two proteins using different reading frames. Analysis of overlapping genes can help in understanding the mode of evolution of a coding region from noncoding DNA. We identified 71 pairs of convergent genes, with overlapping 3' ends longer than 15 nucleotides, that are conserved in at least two prokaryotic genomes. Among the overlap regions, we observed a statistically significant bias towards the 123:132 phase (i.e. the second codon base in one gene facing the degenerate third position in the second gene). This phase ensures the least mutual constraint on nonconservative amino acid replacements in both overlapping coding sequences. The excess of this phase is compatible with directional (positive) selection acting on the overlapping coding regions. This could be a general evolutionary mode for genes emerging from noncoding sequences, in which the protein sequence has not been subject to selection.

SEL1L is a recently cloned and organ-specific expressing human gene whose function is still at an embryonic stage but displays several interesting characteristics, among which a remarkable cross-species conservation. During evolution, the gene structural complexity increased, suggesting a diversification of its function; however, several amino acid motifs remain perfectly conserved from the bacteria to the human protein. SEL1L is the human ortholog of the C. elegans gene sel-1; the latter is implicated in the negative regulation of LIN-12/GLP-1/Notch receptor proteins. These receptor proteins play fundamental roles in signal transduction pathways and are key players in cell fate determination during the development of various organs. Studies in model organisms, such as C. elegans, helped to illuminate fundamental mechanisms involved in normal cellular functions and human diseases. This paper describes the conserved nature of SEL1L across a wide range of species suggesting, that the encoded protein most likely exerts a very important biological function; it may belong to a subclass of genes considered to be "essential."

During a systematic analysis of conserved gene context in prokaryotic genomes, a previously undetected, complex, partially conserved neighborhood consisting of more than 20 genes was discovered in most Archaea (with the exception of Thermoplasma acidophilum and Halobacterium NRC-1) and some bacteria, including the hyperthermophiles Thermotoga maritima and Aquifex aeolicus. The gene composition and gene order in this neighborhood vary greatly between species, but all versions have a stable, conserved core that consists of five genes. One of the core genes encodes a predicted DNA helicase, often fused to a predicted HD-superfamily hydrolase, and another encodes a RecB family exonuclease; three core genes remain uncharacterized, but one of these might encode a nuclease of a new family. Two more genes that belong to this neighborhood and are present in most of the genomes in which the neighborhood was detected encode, respectively, a predicted HD-superfamily hydrolase (possibly a nuclease) of a distinct family and a predicted, novel DNA polymerase. Another characteristic feature of this neighborhood is the expansion of a superfamily of paralogous, uncharacterized proteins, which are encoded by at least 20–30% of the genes in the neighborhood. The functional features of the proteins encoded in this neighborhood suggest that they comprise a previously undetected DNA repair system, which, to our knowledge, is the first repair system largely specific for thermophiles to be identified. This hypothetical repair system might be functionally analogous to the bacterial–eukaryotic system of translesion, mutagenic repair whose central components are DNA polymerases of the UmuC-DinB-Rad30-Rev1 superfamily, which typically are missing in thermophiles.

The recent sequencing of the complete genome of the fruit fly Drosophila melanogaster has yielded about 30% of the predicted genes with no obvious counterparts in other organisms. These rapidly evolving genes remain largely unexplored. Here, we present evidence for a striking variability in an important Drosophila cell cycle regulator encoded by the gene roughex (rux) in closely related fly species. The unusual level of Rux protein variability indicates that there are very low overall constraints on amino acid substitutions. Despite the lack of sequence similarity, certain common features, including the presence of a C-terminal nuclear localization signal and a functionally important N-terminal RXL cyclin-binding motif, exist between Rux and cyclin-dependent kinase inhibitors of the Cip/Kip family. These results indicate that even some genes involved in key regulatory processes in eukaryotes evolve at extremely high rates.

We describe here the error specificity of mammalian DNA polymerase eta (pol eta), an enzyme that performs translesion DNA synthesis and may participate in somatic hypermutation of immunoglobulin genes. Both mouse and human pol eta lack intrinsic proofreading exonuclease activity and both copy undamaged DNA inaccurately. Analysis of more than 1500 single-base substitutions by human pol eta indicates that error rates for all 12 mismatches are high and variable depending on the composition and symmetry of the mismatch and its location. pol eta also generates tandem base substitutions at an unprecedented rate, and kinetic analysis indicates that it extends a tandem double mismatch about as efficiently as other replicative enzymes extend single-base mismatches. This ability to use an aberrant primer terminus and the high rate of single and double-base substitutions support the idea that pol eta may forego strict shape complementarity in order to facilitate highly efficient lesion bypass. Relaxed discrimination is further indicated by pol eta infidelity for a wide variety of nucleotide deletion and addition errors. The nature and location of these errors suggest that some may be initiated by strand slippage, while others result from additional mechanisms.

The presence of repetitive elements can create serious problems for sequence analysis, especially in the case of homology searches in nucleotide sequence databases. Repetitive elements should be treated carefully by using special programs and databases. In this paper, various aspects of SINE (short interspersed repetitive element) identification, analysis and evolution are discussed.

Several LINE-1s have been isolated and characterized from genomic DNA of the vole, Microtus kirgisorum. Blot hybridization revealed specific restriction patterns of L1 elements in vole genomes. Rehybridization of the genomic blot with a cloned 5'-end fragment revealed two major bands indicating the presence of two different L1 subfamilies. The copy numbers are estimated for different parts of M. kirgisorum L1 elements. Data also demonstrate that most vole L1 elements are truncated at the 5'-end; however, in contrast to mouse, the ORF1 copy number is higher in vole. A difference between the substitution rates of the ORF1 5'-region (approximately 330 nucleotides) and the rest of the L1 coding regions is revealed.

Base analog 6- N -hydroxylaminopurine is a potent mutagen in variety of prokaryotic and eukaroytic organisms. In the review, we discuss recent results of the studies of HAP mutagenic activity, genetic control and specificity in bacteria and yeast with the emphasis to the mechanisms protecting living cells from mutagenic and toxic effects of this base analog.

A repetitive DNA sequence, MS2, was isolated from EcoRI-digested genomic DNA of the vole, Microtus subarvalis. The fragment was cloned and sequenced. Sequence analysis of this 1194-bp fragment revealed a 156-bp region demonstrating a 55% homology with the mouse B1 repeat. The remaining MS2 sequence shows no significant homology with other known Gen-Bank sequences. The results of in situ hybridization of MS2 on vole metaphase chromosomes indicate the fragment is confined to heterochromatin blocks of the sex chromosomes in all but one species (M. arvalis). Distribution of MS2 sequences provides evidence for heterogeneity of the giant heterochromatin blocks of the XY Chromosomes (Chrs) in voles, for the unique cluster-like localization of MS2 within these blocks.

We report a new approach based on the Monte Carlo method, to analyze gene conversion. With this, we have examined 235 somatic mutations in the chicken Vlambda gene and found that about 75% of somatic mutations significantly correlated with donor sequences in 25 pseudo Vlambda genes (set C) versus about 25% that did not (set N). The RGYW and TAA consensus sequences of mutational hotspots were earlier revealed in mammalian V genes. Analysis for correlation between somatic mutations in the Vlambda gene and the consensus sequences showed that the somatic mutations of set N were correlating with the consensus sequences (P(WWrandom)0.01) and the somatic mutations of set C were not. Based on further statistical analysis, we suggest that there must be at least two mechanisms responsible for somatic hypermutagenesis in the Vlambda gene: gene conversion and another, which accounts for the elevated frequencies of somatic mutations at the RGYW and TAA consensus sequences.

Among thousands of long non-coding RNAs (lncRNAs) only a small subset is functionally characterized and the functional annotation of lncRNAs on the genomic scale remains inadequate. In this study we computationally characterized two functionally different parts of human lncRNAs transcriptome based on their ability to bind the polycomb repressive complex, PRC2. This classification is enabled by the fact that while all lncRNAs constitute a diverse set of sequences, the classes of PRC2-binding and PRC2 non-binding lncRNAs possess characteristic combinations of sequence-structure patterns and, therefore, can be separated within the feature space. Based on the specific combination of features, we built several machine-learning classifiers and identified the SVM-based classifier as the best performing. We further showed that the SVM-based classifier is able to generalize on the independent data sets. We observed that this classifier, trained on the human lncRNAs, can predict up to 59.4% of PRC2-binding lncRNAs in mice. This suggests that, despite the low degree of sequence conservation, many lncRNAs play functionally conserved biological roles.

The deaminase-like fold includes, in addition to nucleic acid/nucleotide deaminases, several catalytic domains such as the JAB domain, and others involved in nucleotide and ADP-ribose metabolism. Using sensitive sequence and structural comparison methods, we develop a comprehensive natural classification of the deaminase-like fold and show that its ancestral version was likely to operate on nucleotides or nucleic acids. Consequently, we present evidence that a specific group of JAB domains are likely to possess a DNA repair function, distinct from the previously known deubiquitinating peptidase activity. We also identified numerous previously unknown clades of nucleic acid deaminases. Using inference based on contextual information, we suggest that most of these clades are toxin domains of two distinct classes of bacterial toxin systems, namely polymorphic toxins implicated in bacterial interstrain competition and those that target distantly related cells. Genome context information suggests that these toxins might be delivered via diverse secretory systems, such as Type V, Type VI, PVC and a novel PrsW-like intramembrane peptidase-dependent mechanism. We propose that certain deaminase toxins might be deployed by diverse extracellular and intracellular pathogens as also endosymbionts as effectors targeting nucleic acids of host cells. Our analysis suggests that these toxin deaminases have been acquired by eukaryotes on several independent occasions and recruited as organellar or nucleo-cytoplasmic RNA modifiers, operating on tRNAs, mRNAs and short non-coding RNAs, and also as mutators of hyper-variable genes, viruses and selfish elements. This scenario potentially explains the origin of mutagenic AID/APOBEC-like deaminases, including novel versions from Caenorhabditis, Nematostella and diverse algae and a large class of fast-evolving fungal deaminases. These observations greatly expand the distribution of possible unidentified mutagenic processes catalyzed by nucleic acid deaminases.

Background Accurate estimation of the divergence time of the extant eukaryotes is a fundamentally important but extremely difficult problem owing primarily to gross violations of the molecular clock at long evolutionary distances and the lack of appropriate calibration points close to the date of interest. These difficulties are intrinsic to the dating of ancient divergence events and are reflected in the large discrepancies between estimates obtained with different approaches. Estimates of the age of Last Eukaryotic Common Ancestor (LECA) vary approximately twofold, from ~1,100 million years ago (Mya) to ~2,300 Mya. Results We applied the genome-wide analysis of rare genomic changes associated with conserved amino acids (RGC_CAs) and used several independent techniques to obtain date estimates for the divergence of the major lineages of eukaryotes with calibration intervals for insects, land plants and vertebrates. The results suggest an early divergence of monocot and dicot plants, approximately 340 Mya, raising the possibility of plant-insect coevolution. The divergence of bilaterian animal phyla is estimated at ~400-700 Mya, a range of dates that is consistent with cladogenesis immediately preceding the Cambrian explosion. The origin of opisthokonts (the supergroup of eukaryotes that includes metazoa and fungi) is estimated at ~700-1,000 Mya, and the age of LECA at ~1,000-1,300 Mya. We separately analyzed the red algal calibration interval which is based on single fossil. This analysis produced time estimates that were systematically older compared to the other estimates. Nevertheless, the majority of the estimates for the age of the LECA using the red algal data fell within the 1,200-1,400 Mya interval. Conclusion The inference of a "young LECA" is compatible with the latest of previously estimated dates and has substantial biological implications. If these estimates are valid, the approximately 1 to 1.4 billion years of evolution of eukaryotes that is open to comparative-genomic study probably was preceded by hundreds of millions years of evolution that might have included extinct diversity inaccessible to comparative approaches. Reviewers This article was reviewed by William Martin, Herve Philippe (nominated by I. King Jordan), and Romain Derelle.

Retroposition, a leading mechanism for gene duplication, is an important process shaping the evolution of genomes. Retrogenes are also involved in the gene structure evolution as a major player in the process of intron deletion. Here, we demonstrate the role of retrogenes in intron gain in mammals. We identified one case of “intronization,” the transformation of exonic sequences into an intron, in the primate specific retrogene RNF113B and two independent “intronization” events in the retrogene DCAF12L2, one in the common ancestor of primates and rodents and another one in the rodent lineage. Intron gain resulted from the origin of new splice variants, and both genes have two transcript forms, one with retained intron and one with the intron spliced out. Evolution of these genes, especially RNF113B, has been very dynamic and has been accompanied by several additional events including parental gene loss, secondary retroposition, and exaptation of transposable elements.

–derived SINE, specific forthree rodent families: Muridae, Cricetidae, and Spalacidae. Theseelements are composed of three regions: a 58 RNA-related region(RNA-polymerase III promoter), a tRNA-unrelated region, and a38 AT-rich region (polyadenylation signals); the total length of aB2 element is approximately 200 bp (Deininger 1989; Weiner etal. 1986; Rogers 1985; Krayev et al. 1982).Previously, Kass et al. (1997) and Serdobova and Kramerov(1998) concluded B2 elements do not exist in the human genome.Our computer analysis of sequences from nucleotide databases(Altschul et al. 1990; Benson et al. 1999) revealed several humansequences with significant homology to the rodent B2 consensussequence (Bains and Temple-Smith 1989; Jurka et al. 1992; Okadaand Oshima 1995). However, some of these sequences were ob-tained from human–rodent hybrid cells and may represent con-tamination or transposition within the cells of the hybrids. Toverify the presence of B2 SINEs in the human genome, we de-signed oligonucleotide primers from the mouse B2 element con-sensus sequence. The 58 primer identified the tRNA-related region,and the 38 primer identified the 38 region immediately upstreamfrom the AT-rich sequences of the B2 consensus sequence (Fig. 1).Human genomic DNA from different sources (placenta, blood,and commercial placental DNA) was amplified, cloned, and se-quenced. The average divergence of these sequences from therodent B2 consensus sequences (Kass et al. 1997) is about 6%.However, only one specific subclass of human B2 elements, am-plified with mouse B2 primer sequences, was analyzed, and con-clusions about the divergence of all human B2 elements cannot beattempted. Examples of potential human B2 elements from sub-telomeric DNA sequences (Ning and Flint 1996) with significantly

Rats of the OXYS strain are sensitive to oxidative stress and serve as a biological model of premature aging. We have compared spectra of somatic mutations in a control region of mtDNA from the liver of the OXYS rat strain and of Wistar rats as a control. The majority of nucleotide substitutions in the mutation spectra were represented by transitions: 94 and 97% in the OXYS and Wistar rats, respectively. It was shown that 40% of somatic mutations in the control region of mtDNA from Wistar rats were significantly consistent with the model of dislocation mutagenesis. No statistical support for this model was found for mutations in the control region of mtDNA from OXYS rats. The mutation frequency in the ETAS section was higher in the OXYS strain rats than in Wistar rats. These results suggest different mechanisms of mutagenesis in the two rat strains under study.

Background Evolution of DNA polymerases, the key enzymes of DNA replication and repair, is central to any reconstruction of the history of cellular life. However, the details of the evolutionary relationships between DNA polymerases of archaea and eukaryotes remain unresolved. Results We performed a comparative analysis of archaeal, eukaryotic, and bacterial B-family DNA polymerases, which are the main replicative polymerases in archaea and eukaryotes, combined with an analysis of domain architectures. Surprisingly, we found that eukaryotic Polymerase ε consists of two tandem exonuclease-polymerase modules, the active N-terminal module and a C-terminal module in which both enzymatic domains are inactivated. The two modules are only distantly related to each other, an observation that suggests the possibility that Pol ε evolved as a result of insertion and subsequent inactivation of a distinct polymerase, possibly, of bacterial descent, upstream of the C-terminal Zn-fingers, rather than by tandem duplication. The presence of an inactivated exonuclease-polymerase module in Pol ε parallels a similar inactivation of both enzymatic domains in a distinct family of archaeal B-family polymerases. The results of phylogenetic analysis indicate that eukaryotic B-family polymerases, most likely, originate from two distantly related archaeal B-family polymerases, one form giving rise to Pol ε, and the other one to the common ancestor of Pol α, Pol δ, and Pol ζ. The C-terminal Zn-fingers that are present in all eukaryotic B-family polymerases, unexpectedly, are homologous to the Zn-finger of archaeal D-family DNA polymerases that are otherwise unrelated to the B family. The Zn-finger of Polε shows a markedly greater similarity to the counterpart in archaeal PolD than the Zn-fingers of other eukaryotic B-family polymerases. Conclusion Evolution of eukaryotic DNA polymerases seems to have involved previously unnoticed complex events. We hypothesize that the archaeal ancestor of eukaryotes encoded three DNA polymerases, namely, two distinct B-family polymerases and a D-family polymerase all of which contributed to the evolution of the eukaryotic replication machinery. The Zn-finger might have been acquired from PolD by the B-family form that gave rise to Pol ε prior to or in the course of eukaryogenesis, and subsequently, was captured by the ancestor of the other B-family eukaryotic polymerases. The inactivated polymerase-exonuclease module of Pol ε might have evolved by fusion with a distinct polymerase, rather than by duplication of the active module of Pol ε, and is likely to play an important role in the assembly of eukaryotic replication and repair complexes. Reviewers This article was reviewed by Patrick Forterre, Arcady Mushegian, and Chris Ponting. For the full reviews, please go to the Reviewers' Reports section.

Based on the analysis of some immunoglobulin V-gene sequences, somatic mutations are assumed to occur by correction of complementary violations in heteroduplexes formed by direct or inverted repeats. Correlation between somatic mutations and repeats is investigated by a statistical weights method in 323 somatic substitutions in 14 V-genes. Assuming absence of correlation, the probability of observing data in the sample would be very low (0.00004). This result supports the idea that somatic mutations may arise from heteroduplex repair. The high frequency of these mutations in complementarity-determining regions (CDRs) of V-genes may be due to a high concentration of repeats in these regions. Analysis of somatic substitutions has revealed that stabilizing selection seems to provide conservation of framework regions (FRs) (which leads to preservation of the protein's three-dimensional structure). Positive selection may be provided by B-lymphocyte proliferation with large changes in CDRs.

The eukaryotic DNA polymerase delta (Pol delta) participates in genome replication, homologous recombination, DNA repair and damage tolerance. Regulation of the plethora of Pol delta functions depends on the interaction between the second (p50) and third (p66) non-catalytic subunits. We report the crystal structure of p50*p66(N) complex featuring oligonucleotide binding and phosphodiesterase domains in p50 and winged helix-turn-helix N-terminal domain in p66. Disruption of the interaction between the yeast orthologs of p50 and p66 by strategic amino acid changes leads to cold-sensitivity, sensitivity to hydroxyurea and to reduced UV mutagenesis, mimicking the phenotypes of strains where the third subunit of Pol delta is absent. The second subunits of all B family replicative DNA polymerases in archaea and eukaryotes, except Pol delta, share a three-domain structure similar to p50*p66(N), raising the possibility that a portion of the gene encoding p66 was derived from the second subunit gene relatively late in evolution.

The two types of eukaryotic spliceosomal introns, U2 and U12, possess different splice signals and are excised by distinct spliceosomes. The nature of the primordial introns remains uncertain. A comparison of the amino acid distributions at insertion sites of introns that retained their positions throughout eukaryotic evolution with the distributions for human and Arabidopsis thaliana U2 and U12 introns reveals close similarity with U2 but not U12. Thus, the primordial spliceosomal introns were, most likely, U2-type.

Activation-induced deaminase (AID) is the master regulator of class switch recombination (CSR) and somatic hypermutation (SHM), but the mechanisms regulating AID function are obscure. The differential pattern of switch plasmid activity in three IgM(+)/AID(+) and two IgG(+)/AID(+) B cell lines prompted an analysis of global gene expression to discover the origin of these cells. Gene profiling suggested that the IgG(+)/AID(+) B cell lines derived from germinal center B cells. Analysis of SHM potential demonstrates that the IgVkappa domains are inducibly diversified at high rate during in vitro culture. The mutation spectra focused to A:T base pairs, revealing a component of the hypermutation program that occurs preferentially during phase 2 of SHM. The A:T error spectra were analyzed and were not characteristic of polymerase eta activity. A differential pattern of three consensus motifs used for A:T base substitutions was observed in WT and Poleta-, Msh2- and Msh6-deficient B cells. Strikingly, mutations in our B cell lines recapitulated the mutable motif profile for Poleta and Msh2 deficiency, respectively, and suggest that an additional pathway for the generation of A:T mutations in SHM is conserved in mouse and human.

In the positive-sense RNA genome of Beet yellows Closterovirus (BYV), the 3′-terminal open reading frames (ORFs) 2–8 are expressed as a nested set of subgenomic (sg) RNAs. ORFs 2–6, coding for the structural and movement proteins, form a ‘five-gene block’ conserved in closteroviruses. We mapped the 5′-end of the ORF 4 sgRNA, which encodes the p64 protein, at adenosine-11169 in the BYV genome. This completes the mapping of the transcription start sites for the five-gene block sgRNAs of BYV. Computer-assisted analysis of the sequences upstream of BYV ORFs 2, 3, 4, 5, and 6 revealed two conserved motifs, which might constitute the subgenomic promoter elements. These motifs are conserved in the equivalent positions upstream of three orthologous genes of Citrus tristeza Closterovirus and two orthologous genes of Beet yellow stunt Closterovirus .

As the number of sequenced genomes from diverse walks of life rapidly increases, phylogenetic analysis is entering a new era: reconstruction of the evolutionary history of organisms on the basis of full-scale comparison of their genomes. In addition to brute force, genome-wide analysis of alignments, rare genomic changes (RGCs) that are thought to comprise derived shared characters of individual clades are increasingly used in genome-wide phylogenetic studies. We propose a new type of RGCs designated RGC_CAMs (after Conserved Amino acids-Multiple substitutions), which are inferred using a genome-scale analysis of protein and underlying nucleotide sequence alignments. The RGC_CAM approach utilizes amino acid residues conserved in major eukaryotic lineages, with the exception of a few species comprising a putative clade, and selects for phylogenetic inference only those amino acid replacements that require 2 or 3 nucleotide substitutions, in order to reduce homoplasy. The RGC_CAM analysis was combined with a procedure for rigorous statistical testing of competing phylogenetic hypotheses. The RGC_CAM method is shown to be robust to branch length differences and taxon sampling. When applied to animal phylogeny, the RGC_CAM approach strongly supports the coelomate clade that unites chordates with arthropods as opposed to the ecdysozoan (molting animals) clade. This conclusion runs against the view of animal evolution that is currently prevailing in the evo-devo community. The final solution to the coelomate-ecdysozoa controversy will require a much larger set of complete genome sequences representing diverse animal taxa. It is expected that RGC_CAM and other RGC-based methods will be crucial for these future, definitive phylogenetic studies.

Many fundamental questions concerning the emergence and subsequent evolution of eukaryotic exon–intron organization are still unsettled. Genome-scale comparative studies, which can shed light on crucial aspects of eukaryotic evolution, require adequate computational tools. We describe novel computational methods for studying spliceosomal intron evolution. Our goal is to give a reliable characterization of the dynamics of intron evolution. Our algorithmic innovations address the identification of orthologous introns, and the likelihood-based analysis of intron data. We discuss a compression method for the evaluation of the likelihood function, which is noteworthy for phylogenetic likelihood problems in general. We prove that after O(nl) preprocessing time, subsequent evaluations take O(nl/log l) time almost surely in the Yule–Harding random model of n-taxon phylogenies, where l is the input sequence length. We illustrate the practicality of our methods by compiling and analyzing a data set involving 18 eukaryotes, which is more than in any other study to date. The study yields the surprising result that ancestral eukaryotes were fairly intron-rich. For example, the bilaterian ancestor is estimated to have had more than 90% as many introns as vertebrates do now. Availability: The Java implementations of the algorithms are publicly available from the corresponding author's site http://www.iro.umontreal.ca/~csuros/introns/. Data are available on request. Contact: csuros@iro.umontreal.ca

The functioning of the eukaryotic genome depends on efficient and accurate DNA replication and repair. The process of replication is complicated by the ongoing decomposition of DNA and damage of the genome by endogenous and exogenous factors. DNA damage can alter base coding potential resulting in mutations, or block DNA replication, which can lead to double-strand breaks (DSB) and to subsequent chromosome loss. Replication is coordinated with DNA repair systems that operate in cells to remove or tolerate DNA lesions. DNA polymerases can serve as sensors in the cell cycle checkpoint pathways that delay cell division until damaged DNA is repaired and replication is completed. Eukaryotic DNA template-dependent DNA polymerases have different properties adapted to perform an amazingly wide spectrum of DNA transactions. In this review, we discuss the structure, the mechanism, and the evolutionary relationships of DNA polymerases and their possible functions in the replication of intact and damaged chromosomes, DNA damage repair, and recombination.

Background Carcinogenesis typically involves multiple somatic mutations in caretaker (DNA repair) and gatekeeper (tumor suppressors and oncogenes) genes. Analysis of mutation spectra of the tumor suppressor that is most commonly mutated in human cancers, p53, unexpectedly suggested that somatic evolution of the p53 gene during tumorigenesis is dominated by positive selection for gain of function. This conclusion is supported by accumulating experimental evidence of evolution of new functions of p53 in tumors. These findings prompted a genome-wide analysis of possible positive selection during tumor evolution. Methods A comprehensive analysis of probable somatic mutations in the sequences of Expressed Sequence Tags (ESTs) from malignant tumors and normal tissues was performed in order to access the prevalence of positive selection in cancer evolution. For each EST, the numbers of synonymous and non-synonymous substitutions were calculated. In order to identify genes with a signature of positive selection in cancers, these numbers were compared to: i) expected numbers and ii) the numbers for the respective genes in the ESTs from normal tissues. Results We identified 112 genes with a signature of positive selection in cancers, i.e., a significantly elevated ratio of non-synonymous to synonymous substitutions, in tumors as compared to 37 such genes in an approximately equal-sized EST collection from normal tissues. A substantial fraction of the tumor-specific positive-selection candidates have experimentally demonstrated or strongly predicted links to cancer. Conclusion The results of EST analysis should be interpreted with extreme caution given the noise introduced by sequencing errors and undetected polymorphisms. Furthermore, an inherent limitation of EST analysis is that multiple mutations amenable to statistical analysis can be detected only in relatively highly expressed genes. Nevertheless, the present results suggest that positive selection might affect a substantial number of genes during tumorigenic somatic evolution.

Instead of the immunoglobulin-type antigen receptors of jawed vertebrates, jawless fish have variable lymphocyte receptors (VLRs), which consist of leucine-rich repeat (LRR) modules. Somatic diversification of the VLR gene is shown here to occur through a multistep assembly of LRR modules randomly selected from a large bank of flanking cassettes. The predicted concave surface of the VLR is lined with hypervariable positively selected residues, and computational analysis suggests a repertoire of about 10 14 unique receptors. Lamprey immunized with anthrax spores responded with the production of soluble antigen-specific VLRs. These findings reveal that two strikingly different modes of antigen recognition through rearranged lymphocyte receptors have evolved in the jawless and jawed vertebrates.

By comparing sequences of human, mouse and rat orthologous genes, we show that in 5′-untranslated regions (5′-UTRs) of mammalian cDNAs but not in 3′-UTRs or coding sequences, AUG is conserved to a significantly greater extent than any of the other 63 nt triplets. This effect is likely to reflect, primarily, bona fide evolutionary conservation, rather than cDNA annotation artifacts, because the excess of conserved upstream AUGs (uAUGs) is seen in 5′-UTRs containing stop codons in-frame with the start AUG and many of the conserved AUGs are found in different frames, consistent with the location in authentic non-coding sequences. Altogether, conserved uAUGs are present in at least 20–30% of mammalian genes. Qualitatively similar results were obtained by comparison of orthologous genes from different species of the yeast genus Saccharomyces. Together with the observation that mammalian and yeast 5′-UTRs are significantly depleted in overall AUG content, these findings suggest that AUG triplets in 5′-UTRs are subject to the pressure of purifying selection in two opposite directions: the uAUGs that have no specific function tend to be deleterious and get eliminated during evolution, whereas those uAUGs that do serve a function are conserved. Most probably, the principal role of the conserved uAUGs is attenuation of translation at the initiation stage, which is often additionally regulated by alternative splicing in the mammalian 5′-UTRs. Consistent with this hypothesis, we found that open reading frames starting from conserved uAUGs are significantly shorter than those starting from non-conserved uAUGs, possibly, owing to selection for optimization of the level of attenuation.

DNA polymerase (pol) η participates in hypermutation of A:T bases in Ig genes because humans deficient for the polymerase have fewer substitutions of these bases. To determine whether polymerase η is also responsible for the well-known preference for mutations of A vs T on the nontranscribed strand, we sequenced variable regions from three patients with xeroderma pigmentosum variant (XP-V) disease, who lack polymerase η. The frequency of mutations in the intronic region downstream of rearranged JH4 gene segments was similar between XP-V and control clones; however, there were fewer mutations of A:T bases and correspondingly more substitutions of C:G bases in the XP-V clones (p < 10−7). There was significantly less of a bias for mutations of A compared with T nucleotides in the XP-V clones compared with control clones, whereas the frequencies for mutations of C and G were identical in both groups. An analysis of mutations in the WA sequence motif suggests that polymerase η generates more mutations of A than T on the nontranscribed strand. This in vivo data from polymerase η-deficient B cells correlates well with the in vitro specificity of the enzyme. Because polymerase η inserts more mutations opposite template T than template A, it would generate more substitutions of A on the newly synthesized strand.

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and repair/replication/modification enzymes, and the analysis of hotspot context provides evidence of such interactions. The hotspots might also reflect structural and functional features of the respective DNA sequences and provide information about natural selection. We discuss analysis of 8-oxoguanine-induced mutations in pro- and eukaryotic genes, polymorphic positions in the human mitochondrial DNA and mutations in the HIV-1 retrovirus. Comparative analysis of 8-oxoguanine-induced mutations and spontaneous mutation spectra suggested that a substantial fraction of spontaneous A x T--C x T mutations is caused by 8-oxoGTP in nucleotide pools. In the case of human mitochondrial DNA, significant differences between molecular mechanisms of mutations in hypervariable segments and coding part of DNA were detected. Analysis of mutations in the HIV-1 retrovirus suggested a complex interplay between molecular mechanisms of mutagenesis and natural selection.

Summary To study spontaneous base substitutions in human mitochondrial DNA (mtDNA), we reconstructed the mutation spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. Classification analysis revealed numerous mutation hotspots in HVS I and II mutation spectra. Statistical analysis suggested that strand dislocation mutagenesis, operating in monotonous runs of nucleotides, plays an important role in generating base substitutions in the mtDNA control region. The frequency of mutations compatible with the primer strand dislocation in the HVS I region was almost twice as high as that for template strand dislocation. Frequencies of mutations compatible with the primer and template strand dislocation models are almost equal in the HVS II region. Further analysis of strand dislocation models suggested that an excess of pyrimidine transitions in mutation spectra, reconstructed on the basis of the L-strand sequence, is caused by an excess of both L-strand pyrimidine transitions and H-strand purine transitions. In general, no significant bias toward parent H-strand-specific dislocation mutagenesis was found in the HVS I and II regions.

Human sperm protein associated with the nucleus on the X chromosome ( SPANX ) genes comprise a gene family with five known members ( SPANX - A1 , - A2 , - B , - C , and - D ), encoding cancer/testis-specific antigens that are potential targets for cancer immunotherapy. These highly similar paralogous genes cluster on the X chromosome at Xq27. We isolated and sequenced primate genomic clones homologous to human SPANX . Analysis of these clones and search of the human genome sequence revealed an uncharacterized group of genes, SPANX - N , which are present in all primates as well as in mouse and rat. In humans, four SPANX - N genes comprise a series of tandem duplicates at Xq27; a fifth member of this subfamily is located at Xp11. Similarly to SPANX - A / D , human SPANX - N genes are expressed in normal testis and some melanoma cell lines; testis-specific expression of SPANX is also conserved in mouse. Analysis of the taxonomic distribution of the long and short forms of the intron indicates that SPANX - N is the ancestral form, from which the SPANX - A / D subfamily evolved in the common ancestor of the hominoid lineage. Strikingly, the coding sequences of the SPANX genes evolved much faster than the intron and the 5′ untranslated region. There is a strong correlation between the rates of evolution of synonymous and nonsynonymous codon positions, both of which are accelerated 2-fold or more compared to the noncoding sequences. Thus, evolution of the SPANX family appears to have involved positive selection that affected not only the protein sequence but also the synonymous sites in the coding sequence.

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.