Nonstructural p26 proteins encoded by the 3'-proximal genes of velariviruses and criniviruses are orthologs.

The 3'-most genes in RNA-2 of the Crinivirus genus members (family Closteroviridae) code for non-structural p26 proteins that share amino acid sequence similarity [Stewart LR, Hwang MS, Falk BW (2009) Virus Res 145:293-299]. In this study, sensitive bioinformatic tools have been used to identify the homologous p26 proteins encoded by the 3' genes in monopartite genomes of the members of Velarivirus, another Closteroviridae genus, and mint vein banding-associated virus, an unassigned member of the family. The p26 proteins showed similarity in their predicted secondary structures, but an amino acid sequence alignment showed no strictly conserved positions, thus indicating a high plasticity of these non-structural proteins. The implications of the sequence analysis for possible functions of the crinivirus and velarivirus p26 proteins are discussed.

Retention of duplicated ITAM-containing transmembrane signaling subunits in the tetraploid amphibian species Xenopus laevis.

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.

Novel aspects of Kindlin-3 function in humans based on a new case of leukocyte adhesion deficiency III.

BACKGROUND: Kindlin-3 is a novel integrin activator in hematopoietic cells, and its deficiency leads to immune problems and severe bleeding, known as leukocyte adhesion deficiency III (LAD-III). Our current understanding of Kindlin-3 function primarily relies on analysis of animal models or cell lines. OBJECTIVES: To understand the functions of Kindlin-3 in human primary blood cells. PATIENTS/METHODS: We analyzed primary and immortalized hematopoietic cells obtained from a new LAD-III patient with immune problems, bleeding, a history of anemia, and abnormally shaped red blood cells. RESULTS: The patient's white blood cells (WBCs) and platelets showed defects in agonist-induced integrin activation and botrocetin-induced platelet agglutination. Primary leukocytes from this patient exhibited abnormal activation of ß(1) integrin. Integrin activation defects were responsible for the observed deficiency in the botrocetin-induced platelet response. Analysis of patient genomic DNA revealed a novel mutation in the Kindlin3 gene. The mutation abolished Kindlin-3 expression in primary WBCs and platelets, owing to abnormal splicing. Kindlin-3 is expressed in red blood cells (RBCs), and its deficiency is proposed to lead to abnormally shaped RBCs. Immortalized patient WBCs expressed a truncated form of Kindlin-3 that was not sufficient to support integrin activation. Expression of Kindlin-3 cDNA in immortalized patient WBCs rescued integrin activation defects, whereas overexpression of the truncated form did not. CONCLUSIONS: Kindlin-3 deficiency impairs integrin function, including activation of ß(1) integrin. Abnormalities in glycoprotein Ib-IX function in Kindlin-3-deficient platelets are secondary to integrin defects. The region of Kindlin-3 encoded by exon 11 is crucial for its ability to activate integrins in humans.

Mutator effects and mutation signatures of editing deaminases produced in bacteria and yeast.

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.

Analysis of mitochondrial DNA somatic mutations in OXYS and Wistar strain rats.

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.

[Vertebrate immunity: mutator proteins and their evolution].

M.E. Lobashev has brilliantly postulated in 1947 that error-prone repair contribute to mutations in cells. This was shown to be true once the mechanisms of UV mutagenesis in Escherichia coli were deciphered. Induced mutations are generated during error-prone SOS DNA repair with the involvement of inaccurate DNA polymerases belonging to the Y family. Currently, several distinct mutator enzymes participating in spontaneous and induced mutagenesis have been identified. Upon induction of these proteins, mutation rates increase by several orders of magnitude. These proteins regulate the mutation rates in evolution and in ontogeny during immune response. In jawed vertebrates, somatic hypermutagenesis occurs in the variable regions of immunoglobulin genes, leading to affinity maturation of antibodies. The process is initiated by cytidine deamination in DNA to uracil by AID (Activation-Induced Deaminase). Further repair of uracil-containing DNA through proteins that include the Y family DNA polymerases causes mutations, induce gene conversion, and class switch recombination. In jawless vertebrates, the variable lymphocyte receptors (VLR) serve as the primary molecules for adaptive immunity. Generation of mature VLRs most likely depends on agnathan AID-like deaminases. AID and its orthologs in lamprey (PmCDA1 and PMCDA2) belong to the AID/APOBEC family of RNA/DNA editing cytidine deaminases. This family includes enzymes with different functions: APOBEC1 edits RNA, APOBEC3 restricts retroviruses. The functions of APOBEC2 and APOBEC4 have not been yet determined. Here, we report a new member of the AID/APOBEC family, APOBEC5, in the bacterium Xanthomonas oryzae. The widespread presence of RNA/DNA editing deaminases suggests that they are an ancient means of generating genetic diversity.

Mutagenesis by transient misalignment in the human mitochondrial DNA control region.

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.

Cloning of human centromeres by transformation-associated recombination in yeast and generation of functional human artificial chromosomes.

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.

[Context organization of mRNA 5'-untranslated regions of higher plants]. / Kontekstnaia organizatsiia 5'-netransliruemykh raionov mRNK genov vysshikh rastenii.

Computer analysis of nucleotide sequences of 5'-untranslated regions (5'-UTR) of higher plants mRNA adopted from the EMBL nucleotide sequence databank was carried out. It was demonstrated that the average nucleotide frequencies of the leader sequences and adjacent regions of basal promoters are similar, whereas introns and 3'-UTR have a higher content of T and a lower content of C. A particular 5'-UTR contextual feature is a misbalance in the content of complementary nucleotides; probably a stable secondary structure adversely affects the translation properties of the leader sequence. About 20% of 5'-UTR contain AUG triplets, which is twice the earlier estimate. Considered are the properties of leader open reading frames (uORF), the possible causes of their high content in 5'-UTRs of eukaryotic mTNAs, and correlations between the features of uORFs and of the protein-encoding sequence of the gene. It is demonstrated that in effectively translated mRNAs the leader AUG triplets are more frequently located in a nonoptimal context, whereas terminating codons of uORFs more frequently exist in the optimal one. A hypothesis is put forward that the efficiency of termination at the uORF stop codon might substantially interfere with the mRNA translation activity.

Genome trees constructed using five different approaches suggest new major bacterial clades.

BACKGROUND: The availability of multiple complete genome sequences from diverse taxa prompts the development of new phylogenetic approaches, which attempt to incorporate information derived from comparative analysis of complete gene sets or large subsets thereof. Such attempts are particularly relevant because of the major role of horizontal gene transfer and lineage-specific gene loss, at least in the evolution of prokaryotes. RESULTS: Five largely independent approaches were employed to construct trees for completely sequenced bacterial and archaeal genomes: i) presence-absence of genomes in clusters of orthologous genes; ii) conservation of local gene order (gene pairs) among prokaryotic genomes; iii) parameters of identity distribution for probable orthologs; iv) analysis of concatenated alignments of ribosomal proteins; v) comparison of trees constructed for multiple protein families. All constructed trees support the separation of the two primary prokaryotic domains, bacteria and archaea, as well as some terminal bifurcations within the bacterial and archaeal domains. Beyond these obvious groupings, the trees made with different methods appeared to differ substantially in terms of the relative contributions of phylogenetic relationships and similarities in gene repertoires caused by similar life styles and horizontal gene transfer to the tree topology. The trees based on presence-absence of genomes in orthologous clusters and the trees based on conserved gene pairs appear to be strongly affected by gene loss and horizontal gene transfer. The trees based on identity distributions for orthologs and particularly the tree made of concatenated ribosomal protein sequences seemed to carry a stronger phylogenetic signal. The latter tree supported three potential high-level bacterial clades,: i) Chlamydia-Spirochetes, ii) Thermotogales-Aquificales (bacterial hyperthermophiles), and ii) Actinomycetes-Deinococcales-Cyanobacteria. The latter group also appeared to join the low-GC Gram-positive bacteria at a deeper tree node. These new groupings of bacteria were supported by the analysis of alternative topologies in the concatenated ribosomal protein tree using the Kishino-Hasegawa test and by a census of the topologies of 132 individual groups of orthologous proteins. Additionally, the results of this analysis put into question the sister-group relationship between the two major archaeal groups, Euryarchaeota and Crenarchaeota, and suggest instead that Euryarchaeota might be a paraphyletic group with respect to Crenarchaeota. CONCLUSIONS: We conclude that, the extensive horizontal gene flow and lineage-specific gene loss notwithstanding, extension of phylogenetic analysis to the genome scale has the potential of uncovering deep evolutionary relationships between prokaryotic lineages.

Rapid evolution of a cyclin A inhibitor gene, roughex, in Drosophila.

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.

Presence of ATG triplets in 5' untranslated regions of eukaryotic cDNAs correlates with a 'weak' context of the start codon.

MOTIVATION: The context of the start codon (typically, AUG) and the features of the 5' Untranslated Regions (5' UTRs) are important for understanding translation regulation in eukaryotic mRNAs and for accurate prediction of the coding region in genomic and cDNA sequences. The presence of AUG triplets in 5' UTRs (upstream AUGs) might effect the initiation rate and, in the context of gene prediction, could reduce the accuracy of the identification of the authentic start. To reveal potential connections between the presence of upstream AUGs and other features of 5' UTRs, such as their length and the start codon context, we undertook a systematic analysis of the available eukaryotic 5' UTR sequences. RESULTS: We show that a large fraction of 5' UTRs in the available cDNA sequences, 15-53% depending on the organism, contain upstream ATGs. A negative correlation was observed between the information content of the translation start signal and the length of the 5' UTR. Similarly, a negative correlation exists between the 'strength' of the start context and the number of upstream ATGs. Typically, cDNAs containing long 5' UTRs with multiple upstream ATGs have a 'weak' start context, and in contrast, cDNAs containing short 5' UTRs without ATGs have 'strong' starts. These counter-intuitive results may be interpreted in terms of upstream AUGs having an important role in the regulation of translation efficiency by ensuring low basal translation level via double negative control and creating the potential for additional regulatory mechanisms. One of such mechanisms, supported by experimental studies of some mRNAs, includes removal of the AUG-containing portion of the 5' UTR by alternative splicing. AVAILABILITY: An ATG_ EVALUATOR program is available upon request or at www.itba.mi.cnr.it/webgene. CONTACT: rogozin@ncbi.nlm.nih.gov, milanesi@itba.mi.cnr.it.

Error rate and specificity of human and murine DNA polymerase eta.

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.

[Study of the DNA primary structure effect on induction of mutations by alkylating agents]. / Issledovanie vliianiia osobennostei pervichnoi struktury DNK na vozniknovenie mutatsii, indutsirovannykh alkiliruuiushchimi agentami.

Analysis and comparison of mutation spectra is one of the major tasks of molecular biology, since mutation spectra often reveal important properties of various mutagens and proteins involved in the repair/replication systems. Mutability is known to vary significantly along the nucleotide sequence. Mutations are abundant at certain positions (mutation hotspots). In this work, we applied regression analysis based on the basic logic patterns to understand the role of the nucleotide sequence context in mutation induction. The spectra of mutations induced by various alkylating agents were studied. The nucleotide bases at positions -2, -1, +1 and +2 were shown to have the most significant effect in G:C-->A:T replacements.

Characterization of the genomic Xist locus in rodents reveals conservation of overall gene structure and tandem repeats but rapid evolution of unique sequence.

The Xist locus plays a central role in the regulation of X chromosome inactivation in mammals, although its exact mode of action remains to be elucidated. Evolutionary studies are important in identifying conserved genomic regions and defining their possible function. Here we report cloning, sequence analysis, and detailed characterization of the Xist gene from four closely related species of common vole (field mouse), Microtus arvalis. Our analysis reveals that there is overall conservation of Xist gene structure both between different vole species and relative to mouse and human Xist/XIST. Within transcribed sequence, there is significant conservation over five short regions of unique sequence and also over Xist-specific tandem repeats. The majority of unique sequences, however, are evolving at an unexpectedly high rate. This is also evident from analysis of flanking sequences, which reveals a very high rate of rearrangement and invasion of dispersed repeats. We discuss these results in the context of Xist gene function and evolution.

Comparative study and prediction of DNA fragments associated with various elements of the nuclear matrix.

Scaffold/matrix-associated region (S/MAR) sequences are DNA regions that are attached to the nuclear matrix, and participate in many cellular processes. The nuclear matrix is a complex structure consisting of various elements. In this paper we compared frequencies of simple nucleotide motifs in S/MAR sequences and in sequences extracted directly from various nuclear matrix elements, such as nuclear lamina, cores of rosette-like structures, synaptonemal complex. Multivariate linear discriminant analysis revealed significant differences between these sequences. Based on this result we have developed a program, ChrClass (Win/NT version, ftp.bionet.nsc.ru/pub/biology/chrclass/chrclass.zip), for the prediction of the regions associated with various elements of the nuclear matrix in a query sequence. Subsequently, several test samples were analyzed by using two S/MAR prediction programs (a ChrClass and MAR-Finder) and a simple MRS criterion (S/MAR recognition signature) indicating the presence of S/MARs. Some overlap between the predictions of all MAR prediction tools has been found. Simultaneous use of the ChrClass, MRS criterion and MAR-Finder programs may help to obtain a more clearcut picture of S/MAR distribution in a query sequence. In general, our results suggest that the proportion of missed S/MARs is lower for ChrClass, whereas the proportion of wrong S/MARs is lower for MAR-Finder and MRS.

Somatic mutation hotspots correlate with DNA polymerase eta error spectrum.

Mutational spectra analysis of 15 immunoglobulin genes suggested that consensus motifs RGYW and WA were universal descriptors of somatic hypermutation. Highly mutable sites, "hotspots", that matched WA were preferentially found in one DNA strand and RGYW hotspots were found in both strands. Analysis of base-substitution hotspots in DNA polymerase error spectra showed that 33 of 36 hotspots in the human polymerase eta spectrum conformed to the WA consensus. This and four other characteristics of polymerase eta substitution specificity suggest that errors introduced by this enzyme during synthesis of the nontranscribed DNA strand in variable regions may contribute to strand-specific somatic hypermutagenesis of immunoglobulin genes at A-T base pairs.

Genome alignment, evolution of prokaryotic genome organization, and prediction of gene function using genomic context.

Gene order in prokaryotes is conserved to a much lesser extent than protein sequences. Only several operons, primarily those that code for physically interacting proteins, are conserved in all or most of the bacterial and archaeal genomes. Nevertheless, even the limited conservation of operon organization that is observed can provide valuable evolutionary and functional clues through multiple genome comparisons. A program for constructing gapped local alignments of conserved gene strings in two genomes was developed. The statistical significance of the local alignments was assessed using Monte Carlo simulations. Sets of local alignments were generated for all pairs of completely sequenced bacterial and archaeal genomes, and for each genome a template-anchored multiple alignment was constructed. In most pairwise genome comparisons, <10% of the genes in each genome belonged to conserved gene strings. When closely related pairs of species (i.e., two mycoplasmas) are excluded, the total coverage of genomes by conserved gene strings ranged from <5% for the cyanobacterium Synechocystis sp to 24% for the minimal genome of Mycoplasma genitalium, and 23% in Thermotoga maritima. The coverage of the archaeal genomes was only slightly lower than that of bacterial genomes. The majority of the conserved gene strings are known operons, with the ribosomal superoperon being the top-scoring string in most genome comparisons. However, in some of the bacterial-archaeal pairs, the superoperon is rearranged to the extent that other operons, primarily those subject to horizontal transfer, show the greatest level of conservation, such as the archaeal-type H+-ATPase operon or ABC-type transport cassettes. The level of gene order conservation among prokaryotic genomes was compared to the cooccurrence of genomes in clusters of orthologous genes (COGs) and to the conservation of protein sequences themselves. Only limited correlation was observed between these evolutionary variables. Gene order conservation shows a much lower variance than the cooccurrence of genomes in COGs, which indicates that intragenome homogenization via recombination occurs in evolution much faster than intergenome homogenization via horizontal gene transfer and lineage-specific gene loss. The potential of using template-anchored multiple-genome alignments for predicting functions of uncharacterized genes was quantitatively assessed. Functions were predicted or significantly clarified for approximately 90 COGs (approximately 4% of the total of 2414 analyzed COGs). The most significant predictions were obtained for the poorly characterized archaeal genomes; these include a previously uncharacterized restriction-modification system, a nuclease-helicase combination implicated in DNA repair, and the probable archaeal counterpart of the eukaryotic exosome. Multiple genome alignments are a resource for studies on operon rearrangement and disruption, which is central to our understanding of the evolution of prokaryotic genomes. Because of the rapid evolution of the gene order, the potential of genome alignment for prediction of gene functions is limited, but nevertheless, such predictions information significantly complements the results obtained through protein sequence and structure analysis.