Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2.

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

The 29-nucleotide deletion in SARS-CoV: truncated versions of ORF8 are under purifying selection.

BACKGROUND: Accessory proteins have diverse roles in coronavirus pathobiology. One of them in SARS-CoV (the causative agent of the severe acute respiratory syndrome outbreak in 2002-2003) is encoded by the open reading frame 8 (ORF8). Among the most dramatic genomic changes observed in SARS-CoV isolated from patients during the peak of the pandemic in 2003 was the acquisition of a characteristic 29-nucleotide deletion in ORF8. This deletion cause splitting of ORF8 into two smaller ORFs, namely ORF8a and ORF8b. Functional consequences of this event are not entirely clear. RESULTS: Here, we performed evolutionary analyses of ORF8a and ORF8b genes and documented that in both cases the frequency of synonymous mutations was greater than that of nonsynonymous ones. These results suggest that ORF8a and ORF8b are under purifying selection, thus proteins translated from these ORFs are likely to be functionally important. Comparisons with several other SARS-CoV genes revealed that another accessory gene, ORF7a, has a similar ratio of nonsynonymous to synonymous mutations suggesting that ORF8a, ORF8b, and ORF7a are under similar selection pressure. CONCLUSIONS: Our results for SARS-CoV echo the known excess of deletions in the ORF7a-ORF7b-ORF8 complex of accessory genes in SARS-CoV-2. A high frequency of deletions in this gene complex might reflect recurrent searches in "functional space" of various accessory protein combinations that may eventually produce more advantageous configurations of accessory proteins similar to the fixed deletion in the SARS-CoV ORF8 gene.

Low-Density Lipoprotein Receptor (LDLR) Is Involved in Internalization of Lentiviral Particles Pseudotyped with SARS-CoV-2 Spike Protein in Ocular Cells.

Here, we present evidence that caveolae-mediated endocytosis using LDLR is the pathway for SARS-CoV-2 virus internalization in the ocular cell line ARPE-19. Firstly, we found that, while Angiotensin-converting enzyme 2 (ACE2) is expressed in these cells, blocking ACE2 by antibody treatment did not prevent infection by SARS-CoV-2 spike pseudovirions, nor did antibody blockade of extracellular vimentin and other cholesterol-rich lipid raft proteins. Next, we implicated the role of cholesterol homeostasis in infection by showing that incubating cells with different cyclodextrins and oxysterol 25-hydroxycholesterol (25-HC) inhibits pseudovirion infection of ARPE-19. However, the effect of 25-HC is likely not via cholesterol biosynthesis, as incubation with lovastatin did not appreciably affect infection. Additionally, is it not likely to be an agonistic effect of 25-HC on LXR receptors, as the LXR agonist GW3965 had no significant effect on infection of ARPE-19 cells at up to 5 µM GW3965. We probed the role of endocytic pathways but determined that clathrin-dependent and flotillin-dependent rafts were not involved. Furthermore, 20 µM chlorpromazine, an inhibitor of clathrin-mediated endocytosis (CME), also had little effect. In contrast, anti-dynamin I/II antibodies blocked the entry of SARS-CoV-2 spike pseudovirions, as did dynasore, a noncompetitive inhibitor of dynamin GTPase activity. Additionally, anti-caveolin-1 antibodies significantly blocked spike pseudotyped lentiviral infection of ARPE-19. However, nystatin, a classic inhibitor of caveolae-dependent endocytosis, did not affect infection while indomethacin inhibited only at 10 µM at the 48 h time point. Finally, we found that anti-LDLR antibodies block pseudovirion infection to a similar degree as anti-caveolin-1 and anti-dynamin I/II antibodies, while transfection with LDLR-specific siRNA led to a decrease in spike pseudotyped lentiviral infection, compared to scrambled control siRNAs. Thus, we conclude that SARS-CoV-2 spike pseudovirion infection in ARPE-19 cells is a dynamin-dependent process that is primarily mediated by LDLR.

Viral discovery and diversity in trypanosomatid protozoa with a focus on relatives of the human parasite Leishmania.

Knowledge of viral diversity is expanding greatly, but many lineages remain underexplored. We surveyed RNA viruses in 52 cultured monoxenous relatives of the human parasite Leishmania (Crithidia and Leptomonas), as well as plant-infecting PhytomonasLeptomonas pyrrhocoris was a hotbed for viral discovery, carrying a virus (Leptomonas pyrrhocoris ostravirus 1) with a highly divergent RNA-dependent RNA polymerase missed by conventional BLAST searches, an emergent clade of tombus-like viruses, and an example of viral endogenization. A deep-branching clade of trypanosomatid narnaviruses was found, notable as Leptomonas seymouri bearing Narna-like virus 1 (LepseyNLV1) have been reported in cultures recovered from patients with visceral leishmaniasis. A deep-branching trypanosomatid viral lineage showing strong affinities to bunyaviruses was termed "Leishbunyavirus" (LBV) and judged sufficiently distinct to warrant assignment within a proposed family termed "Leishbunyaviridae" Numerous relatives of trypanosomatid viruses were found in insect metatranscriptomic surveys, which likely arise from trypanosomatid microbiota. Despite extensive sampling we found no relatives of the totivirus Leishmaniavirus (LRV1/2), implying that it was acquired at about the same time the Leishmania became able to parasitize vertebrates. As viruses were found in over a quarter of isolates tested, many more are likely to be found in the >600 unsurveyed trypanosomatid species. Viral loss was occasionally observed in culture, providing potentially isogenic virus-free lines enabling studies probing the biological role of trypanosomatid viruses. These data shed important insights on the emergence of viruses within an important trypanosomatid clade relevant to human disease.

Analysis of Stop Codons within Prokaryotic Protein-Coding Genes Suggests Frequent Readthrough Events.

Nonsense mutations turn a coding (sense) codon into an in-frame stop codon that is assumed to result in a truncated protein product. Thus, nonsense substitutions are the hallmark of pseudogenes and are used to identify them. Here we show that in-frame stop codons within bacterial protein-coding genes are widespread. Their evolutionary conservation suggests that many of them are not pseudogenes, since they maintain dN/dS values (ratios of substitution rates at non-synonymous and synonymous sites) significantly lower than 1 (this is a signature of purifying selection in protein-coding regions). We also found that double substitutions in codons-where an intermediate step is a nonsense substitution-show a higher rate of evolution compared to null models, indicating that a stop codon was introduced and then changed back to sense via positive selection. This further supports the notion that nonsense substitutions in bacteria are relatively common and do not necessarily cause pseudogenization. In-frame stop codons may be an important mechanism of regulation: Such codons are likely to cause a substantial decrease of protein expression levels.

The Functional Consequences of the Novel Ribosomal Pausing Site in SARS-CoV-2 Spike Glycoprotein RNA.

The SARS-CoV-2 Spike glycoprotein (S protein) acquired a unique new 4 amino acid -PRRA- insertion sequence at amino acid residues (aa) 681-684 that forms a new furin cleavage site in S protein as well as several new glycosylation sites. We studied various statistical properties of the -PRRA- insertion at the RNA level (CCUCGGCGGGCA). The nucleotide composition and codon usage of this sequence are different from the rest of the SARS-CoV-2 genome. One of such features is two tandem CGG codons, although the CGG codon is the rarest codon in the SARS-CoV-2 genome. This suggests that the insertion sequence could cause ribosome pausing as the result of these rare codons. Due to population variants, the Nextstrain divergence measure of the CCU codon is extremely large. We cannot exclude that this divergence might affect host immune responses/effectiveness of SARS-CoV-2 vaccines, possibilities awaiting further investigation. Our experimental studies show that the expression level of original RNA sequence "wildtype" spike protein is much lower than for codon-optimized spike protein in all studied cell lines. Interestingly, the original spike sequence produces a higher titer of pseudoviral particles and a higher level of infection. Further mutagenesis experiments suggest that this dual-effect insert, comprised of a combination of overlapping translation pausing and furin sites, has allowed SARS-CoV-2 to infect its new host (human) more readily. This underlines the importance of ribosome pausing to allow efficient regulation of protein expression and also of cotranslational subdomain folding.

Mutational signatures and mutable motifs in cancer genomes.

Cancer is a genetic disorder, meaning that a plethora of different mutations, whether somatic or germ line, underlie the etiology of the 'Emperor of Maladies'. Point mutations, chromosomal rearrangements and copy number changes, whether they have occurred spontaneously in predisposed individuals or have been induced by intrinsic or extrinsic (environmental) mutagens, lead to the activation of oncogenes and inactivation of tumor suppressor genes, thereby promoting malignancy. This scenario has now been recognized and experimentally confirmed in a wide range of different contexts. Over the past decade, a surge in available sequencing technologies has allowed the sequencing of whole genomes from liquid malignancies and solid tumors belonging to different types and stages of cancer, giving birth to the new field of cancer genomics. One of the most striking discoveries has been that cancer genomes are highly enriched with mutations of specific kinds. It has been suggested that these mutations can be classified into 'families' based on their mutational signatures. A mutational signature may be regarded as a type of base substitution (e.g. C:G to T:A) within a particular context of neighboring nucleotide sequence (the bases upstream and/or downstream of the mutation). These mutational signatures, supplemented by mutable motifs (a wider mutational context), promise to help us to understand the nature of the mutational processes that operate during tumor evolution because they represent the footprints of interactions between DNA, mutagens and the enzymes of the repair/replication/modification pathways.

Dopamine response gene pathways in dorsal striatum MSNs from a gene expression viewpoint: cAMP-mediated gene networks.

BACKGROUND: Medium spiny neurons (MSNs) comprise the main body (95% in mouse) of the dorsal striatum neurons and represent dopaminoceptive GABAergic neurons. The cAMP (cyclic Adenosine MonoPhosphate)-mediated cascade of excitation and inhibition responses observed in MSN intracellular signal transduction is crucial for neuroscience research due to its involvement in the motor and behavioral functions. In particular, all types of addictions are related to MSNs. Shedding the light on the mechanics of the above-mentioned cascade is of primary importance for this research domain. RESULTS: A mouse model of chronic social conflicts in daily agonistic interactions was used to analyze dorsal striatum neurons genes implicated in cAMP-mediated phosphorylation activation pathways specific for MSNs. Based on expression correlation analysis, we succeeded in dissecting Drd1- and Drd2-dopaminoceptive neurons (D1 and D2, correspondingly) gene pathways. We also found that D1 neurons genes clustering are split into two oppositely correlated states, passive and active ones, the latter apparently corresponding to D1 firing stage upon protein kinase A (PKA) activation. We observed that under defeat stress in chronic social conflicts the loser mice manifest overall depression of dopamine-mediated MSNs activity resulting in previously reported reduced motor activity, while the aggressive mice with positive fighting experience (aggressive mice) feature an increase in both D1-active phase and D2 MSNs genes expression leading to hyperactive behavior pattern corresponded by us before. Based on the alternative transcript isoforms expression analysis, it was assumed that many genes (Drd1, Adora1, Pde10, Ppp1r1b, Gnal), specifically those in D1 neurons, apparently remain transcriptionally repressed via the reversible mechanism of promoter CpG island silencing, resulting in alternative promoter usage following profound reduction in their expression rate. CONCLUSION: Based on the animal stress model dorsal striatum pooled tissue RNA-Seq data restricted to cAMP related genes subset we elucidated MSNs steady states exhaustive projection for the first time. We correspond the existence of D1 active state not explicitly outlined before, and connected with dynamic dopamine neurotransmission cycles. Consequently, we were also able to indicate an oscillated postsynaptic dopamine vs glutamate action pattern in the course of the neurotransmission cycles.

Involvement of DNA mismatch repair in the maintenance of heterochromatic DNA stability in Saccharomyces cerevisiae.

Heterochromatin contains a significant part of nuclear DNA. Little is known about the mechanisms that govern heterochromatic DNA stability. We show here that in the yeast Saccharomyces cerevisiae (i) DNA mismatch repair (MMR) is required for the maintenance of heterochromatic DNA stability, (ii) MutLα (Mlh1-Pms1 heterodimer), MutSα (Msh2-Msh6 heterodimer), MutSß (Msh2-Msh3 heterodimer), and Exo1 are involved in MMR at heterochromatin, (iii) Exo1-independent MMR at heterochromatin frequently leads to the formation of Pol ζ-dependent mutations, (iv) MMR cooperates with the proofreading activity of Pol ε and the histone acetyltransferase Rtt109 in the maintenance of heterochromatic DNA stability, (v) repair of base-base mismatches at heterochromatin is less efficient than repair of base-base mismatches at euchromatin, and (vi) the efficiency of repair of 1-nt insertion/deletion loops at heterochromatin is similar to the efficiency of repair of 1-nt insertion/deletion loops at euchromatin.

Crossing fitness valleys via double substitutions within codons.

BACKGROUND: Single nucleotide substitutions in protein-coding genes can be divided into synonymous (S), with little fitness effect, and non-synonymous (N) ones that alter amino acids and thus generally have a greater effect. Most of the N substitutions are affected by purifying selection that eliminates them from evolving populations. However, additional mutations of nearby bases potentially could alleviate the deleterious effect of single substitutions, making them subject to positive selection. To elucidate the effects of selection on double substitutions in all codons, it is critical to differentiate selection from mutational biases. RESULTS: We addressed the evolutionary regimes of within-codon double substitutions in 37 groups of closely related prokaryotic genomes from diverse phyla by comparing the fractions of double substitutions within codons to those of the equivalent double S substitutions in adjacent codons. Under the assumption that substitutions occur one at a time, all within-codon double substitutions can be represented as "ancestral-intermediate-final" sequences (where "intermediate" refers to the first single substitution and "final" refers to the second substitution) and can be partitioned into four classes: (1) SS, S intermediate-S final; (2) SN, S intermediate-N final; (3) NS, N intermediate-S final; and (4) NN, N intermediate-N final. We found that the selective pressure on the second substitution markedly differs among these classes of double substitutions. Analogous to single S (synonymous) substitutions, SS double substitutions evolve neutrally, whereas analogous to single N (non-synonymous) substitutions, SN double substitutions are subject to purifying selection. In contrast, NS show positive selection on the second step because the original amino acid is recovered. The NN double substitutions are heterogeneous and can be subject to either purifying or positive selection, or evolve neutrally, depending on the amino acid similarity between the final or intermediate and the ancestral states. CONCLUSIONS: The results of the present, comprehensive analysis of the evolutionary landscape of within-codon double substitutions reaffirm the largely conservative regime of protein evolution. However, the second step of a double substitution can be subject to positive selection when the first step is deleterious. Such positive selection can result in frequent crossing of valleys on the fitness landscape.

Palmitoylation of Metazoan Carotenoid Oxygenases.

Abundant in nature, carotenoids are a class of fat-soluble pigments with a polyene tetraterpenoid structure. They possess antioxidant properties and their consumption leads to certain health benefits in humans. Carotenoid cleavage oxygenases (CCOs) are a superfamily of enzymes which oxidatively cleave carotenoids and they are present in all kingdoms of life. Complexity of CCO evolution is high. For example, in this study we serendipitously found a new family of eukaryotic CCOs, the apocarotenoid oxygenase-like (ACOL) family. This family has several members in animal genomes and lacks the animal-specific amino acid motif PDPCK. This motif is likely to be associated with palmitoylation of some animal CCOs. We recently demonstrated that two mammalian members of the carotenoid oxygenase family retinal pigment epithelial-specific 65 kDa protein (RPE65) and beta-carotene oxygenase 2 (BCO2) are palmitoylated proteins. Here we used the acyl-resin-assisted capture (acyl-RAC) method to demonstrate protein palmitoylation and immunochemistry to localize mouse BCO2 (mBCO2) in COS7 cell line in the absence and presence of its substrate ß-carotene. We demonstrate that mBCO2 palmitoylation depends on the evolutionarily conserved motif PDPCK and that metazoan family members lacking the motif (Lancelet beta-carotene oxygenase-like protein (BCOL) and Acropora ACOL) are not palmitoylated. Additionally, we observed that the palmitoylation status of mBCO2 and its membrane association depend on the presence of its substrate ß-carotene. Based on our results we conclude that most metazoan carotenoid oxygenases retain the evolutionarily conserved palmitoylation PDPCK motif to target proteins to internal membranes depending on substrate status. Exceptions are in the secreted BCOL subfamily and the strictly cytosolic ancient ACOL subfamily of carotenoid oxygenases.

Aberrant RNA splicing is the major pathogenic effect in a knock-in mouse model of the dominantly inherited c.1430A>G human RPE65 mutation.

Human RPE65 mutations cause a spectrum of retinal dystrophies that result in blindness. While RPE65 mutations have been almost invariably recessively inherited, a c.1430A>G (p.(D477G)) mutation has been reported to cause autosomal dominant retinitis pigmentosa (adRP). To study the pathogenesis of this human mutation, we have replicated the mutation in a knock-in (KI) mouse model using CRISPR/Cas9-mediated genome editing. Significantly, in contrast to human patients, heterozygous KI mice do not exhibit any phenotypes in visual function tests. When raised in regular vivarium conditions, homozygous KI mice display relatively undisturbed visual functions with minimal retinal structural changes. However, KI/KI mouse retinae are more sensitive to light exposure and exhibit signs of degenerative features when subjected to light stress. We find that instead of merely producing a missense mutant protein, the A>G nucleotide substitution greatly affects appropriate splicing of Rpe65 mRNA by generating an ectopic splice site in comparable context to the canonical one, thereby disrupting RPE65 protein expression. Similar splicing defects were also confirmed for the human RPE65 c.1430G mutant in an in vitro Exontrap assay. Our data demonstrate that a splicing defect is associated with c.1430G pathogenesis, and therefore provide insights in the therapeutic strategy for human patients.

Exploring background mutational processes to decipher cancer genetic heterogeneity.

Much remains unknown about the progression and heterogeneity of mutational processes in different cancers and their diagnostic and clinical potential. A growing body of evidence supports mutation rate dependence on the local DNA sequence context for various types of mutations. We propose several tools for the analysis of cancer context-dependent mutations, which are implemented in an online computational framework MutaGene. The framework explores DNA context-dependent mutational patterns and underlying somatic cancer mutagenesis, analyzes mutational profiles of cancer samples, identifies the combinations of underlying mutagenic processes including those related to infidelity of DNA replication and repair machinery, and various other endogenous and exogenous mutagenic factors. As a result, the combination of mutagenic processes can be identified in any query sample with subsequent comparison to mutational profiles derived from malignant and benign samples. In addition, mutagen or cancer-specific mutational background models are applied to calculate expected DNA and protein site mutability to decouple relative contributions of mutagenesis and selection in carcinogenesis, thus elucidating the site-specific driving events in cancer. MutaGene is freely available at https://www.ncbi.nlm.nih.gov/projects/mutagene/.

Evolutionary switches between two serine codon sets are driven by selection.

Serine is the only amino acid that is encoded by two disjoint codon sets so that a tandem substitution of two nucleotides is required to switch between the two sets. Previously published evidence suggests that, for the most evolutionarily conserved serines, the codon set switch occurs by simultaneous substitution of two nucleotides. Here we report a genome-wide reconstruction of the evolution of serine codons in triplets of closely related species from diverse prokaryotes and eukaryotes. The results indicate that the great majority of codon set switches proceed by two consecutive nucleotide substitutions, via a threonine or cysteine intermediate, and are driven by selection. These findings imply a strong pressure of purifying selection in protein evolution, which in the case of serine codon set switches occurs via an initial deleterious substitution quickly followed by a second, compensatory substitution. The result is frequent reversal of amino acid replacements and, at short evolutionary distances, pervasive homoplasy.

Disruption of Transcriptional Coactivator Sub1 Leads to Genome-Wide Re-distribution of Clustered Mutations Induced by APOBEC in Active Yeast Genes.

Mutations in genomes of species are frequently distributed non-randomly, resulting in mutation clusters, including recently discovered kataegis in tumors. DNA editing deaminases play the prominent role in the etiology of these mutations. To gain insight into the enigmatic mechanisms of localized hypermutagenesis that lead to cluster formation, we analyzed the mutational single nucleotide variations (SNV) data obtained by whole-genome sequencing of drug-resistant mutants induced in yeast diploids by AID/APOBEC deaminase and base analog 6-HAP. Deaminase from sea lamprey, PmCDA1, induced robust clusters, while 6-HAP induced a few weak ones. We found that PmCDA1, AID, and APOBEC1 deaminases preferentially mutate the beginning of the actively transcribed genes. Inactivation of transcription initiation factor Sub1 strongly reduced deaminase-induced can1 mutation frequency, but, surprisingly, did not decrease the total SNV load in genomes. However, the SNVs in the genomes of the sub1 clones were re-distributed, and the effect of mutation clustering in the regions of transcription initiation was even more pronounced. At the same time, the mutation density in the protein-coding regions was reduced, resulting in the decrease of phenotypically detected mutants. We propose that the induction of clustered mutations by deaminases involves: a) the exposure of ssDNA strands during transcription and loss of protection of ssDNA due to the depletion of ssDNA-binding proteins, such as Sub1, and b) attainment of conditions favorable for APOBEC action in subpopulation of cells, leading to enzymatic deamination within the currently expressed genes. This model is applicable to both the initial and the later stages of oncogenic transformation and explains variations in the distribution of mutations and kataegis events in different tumor cells.

Comparative mutational analyses of influenza A viruses.

The error-prone RNA-dependent RNA polymerase (RdRP) and external selective pressures are the driving forces for RNA viral diversity. When confounded by selective pressures, it is difficult to assess if influenza A viruses (IAV) that have a wide host range possess comparable or distinct spontaneous mutational frequency in their RdRPs. We used in-depth bioinformatics analyses to assess the spontaneous mutational frequencies of two RdRPs derived from human seasonal (A/Wuhan/359/95; Wuhan) and H5N1 (A/Vietnam/1203/04; VN1203) viruses using the mini-genome system with a common firefly luciferase reporter serving as the template. High-fidelity reverse transcriptase was applied to generate high-quality mutational spectra which allowed us to assess and compare the mutational frequencies and mutable motifs along a target sequence of the two RdRPs of two different subtypes. We observed correlated mutational spectra (τ correlation P < 0.0001), comparable mutational frequencies (H3N2:5.8 ± 0.9; H5N1:6.0 ± 0.5), and discovered a highly mutable motif "(A)AAG" for both Wuhan and VN1203 RdRPs. Results were then confirmed with two recombinant A/Puerto Rico/8/34 (PR8) viruses that possess RdRP derived from Wuhan or VN1203 (RG-PR8×Wuhan(PB2, PB1, PA, NP) and RG-PR8×VN1203(PB2, PB1, PA, NP)). Applying novel bioinformatics analysis on influenza mutational spectra, we provide a platform for a comprehensive analysis of the spontaneous mutation spectra for an RNA virus.

Genome-wide mutation avalanches induced in diploid yeast cells by a base analog or an APOBEC deaminase.

Genetic information should be accurately transmitted from cell to cell; conversely, the adaptation in evolution and disease is fueled by mutations. In the case of cancer development, multiple genetic changes happen in somatic diploid cells. Most classic studies of the molecular mechanisms of mutagenesis have been performed in haploids. We demonstrate that the parameters of the mutation process are different in diploid cell populations. The genomes of drug-resistant mutants induced in yeast diploids by base analog 6-hydroxylaminopurine (HAP) or AID/APOBEC cytosine deaminase PmCDA1 from lamprey carried a stunning load of thousands of unselected mutations. Haploid mutants contained almost an order of magnitude fewer mutations. To explain this, we propose that the distribution of induced mutation rates in the cell population is uneven. The mutants in diploids with coincidental mutations in the two copies of the reporter gene arise from a fraction of cells that are transiently hypersensitive to the mutagenic action of a given mutagen. The progeny of such cells were never recovered in haploids due to the lethality caused by the inactivation of single-copy essential genes in cells with too many induced mutations. In diploid cells, the progeny of hypersensitive cells survived, but their genomes were saturated by heterozygous mutations. The reason for the hypermutability of cells could be transient faults of the mutation prevention pathways, like sanitization of nucleotide pools for HAP or an elevated expression of the PmCDA1 gene or the temporary inability of the destruction of the deaminase. The hypothesis on spikes of mutability may explain the sudden acquisition of multiple mutational changes during evolution and carcinogenesis.

MitoCOGs: clusters of orthologous genes from mitochondria and implications for the evolution of eukaryotes.

BACKGROUND: Mitochondria are ubiquitous membranous organelles of eukaryotic cells that evolved from an alpha-proteobacterial endosymbiont and possess a small genome that encompasses from 3 to 106 genes. Accumulation of thousands of mitochondrial genomes from diverse groups of eukaryotes provides an opportunity for a comprehensive reconstruction of the evolution of the mitochondrial gene repertoire. RESULTS: Clusters of orthologous mitochondrial protein-coding genes (MitoCOGs) were constructed from all available mitochondrial genomes and complemented with nuclear orthologs of mitochondrial genes. With minimal exceptions, the mitochondrial gene complements of eukaryotes are subsets of the superset of 66 genes found in jakobids. Reconstruction of the evolution of mitochondrial genomes indicates that the mitochondrial gene set of the last common ancestor of the extant eukaryotes was slightly larger than that of jakobids. This superset of mitochondrial genes likely represents an intermediate stage following the loss and transfer to the nucleus of most of the endosymbiont genes early in eukaryote evolution. Subsequent evolution in different lineages involved largely parallel transfer of ancestral endosymbiont genes to the nuclear genome. The intron density in nuclear orthologs of mitochondrial genes typically is nearly the same as in the rest of the genes in the respective genomes. However, in land plants, the intron density in nuclear orthologs of mitochondrial genes is almost 1.5-fold lower than the genomic mean, suggestive of ongoing transfer of functional genes from mitochondria to the nucleus. CONCLUSIONS: The MitoCOGs are expected to become an important resource for the study of mitochondrial evolution. The nearly complete superset of mitochondrial genes in jakobids likely represents an intermediate stage in the evolution of eukaryotes after the initial, extensive loss and transfer of the endosymbiont genes. In addition, the bacterial multi-subunit RNA polymerase that is encoded in the jakobid mitochondrial genomes was replaced by a single-subunit phage-type RNA polymerase in the rest of the eukaryotes. These results are best compatible with the rooting of the eukaryotic tree between jakobids and the rest of the eukaryotes. The land plants are the only eukaryotic branch in which the gene transfer from the mitochondrial to the nuclear genome appears to be an active, ongoing process.

The vast, conserved mammalian lincRNome.

We compare the sets of experimentally validated long intergenic non-coding (linc)RNAs from human and mouse and apply a maximum likelihood approach to estimate the total number of lincRNA genes as well as the size of the conserved part of the lincRNome. Under the assumption that the sets of experimentally validated lincRNAs are random samples of the lincRNomes of the corresponding species, we estimate the total lincRNome size at approximately 40,000 to 50,000 species, at least twice the number of protein-coding genes. We further estimate that the fraction of the human and mouse euchromatic genomes encoding lincRNAs is more than twofold greater than the fraction of protein-coding sequences. Although the sequences of most lincRNAs are much less strongly conserved than protein sequences, the extent of orthology between the lincRNomes is unexpectedly high, with 60 to 70% of the lincRNA genes shared between human and mouse. The orthologous mammalian lincRNAs can be predicted to perform equivalent functions; accordingly, it appears likely that thousands of evolutionarily conserved functional roles of lincRNAs remain to be characterized.

Genetic characterization of feline calicivirus strains associated with varying disease manifestations during an outbreak season in Missouri (1995-1996).

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.