Mutation spectra in Salmonella of analogues of MX: implications of chemical structure for mutational mechanisms

We determined the mutation spectra in Salmonella of four chlorinated butenoic acid analogues (BA-1 through BA-4) of the drinking water mutagen 3-chloro-4-(dichloromethyl)-5-hydroxy-2(5H)-furanone (MX) and compared the results with those generated previously by us for MX and a related compound, MCF. We then considered relationships between the properties of mutagenic potency and mutational specificity for these six chlorinated butenoic acid analogues. In TA98, the three most potent mutagens, BA-3, BA-4, MX, and the organic extract, all induced large percentages of complex frameshifts (33-67%), which distinguish these agents from any other class of compound studied previously. In TA100, which has only GC sites for mutation recovery, >71% of the mutations induced by all of the agents were GC-->TA transversions. The availability of both GC and TA sites for mutation in TA104 resulted in greater distinctions in mutational specificity than in TA100. MX targeted GC sites almost exclusively (98%); the structurally similar BA-4 and BA-2 produced mutations at similar frequencies at both GC and AT sites; and the structurally similar BA-3 and BA-1 induced most mutations at AT sites (69%). Thus, large variations in structural properties influencing relative mutagenic potency appeared to be distinct from the more localized similar structural features influencing mutagenic specificity in TA104. Among a set of physicochemical properties examined for the six butenoic acids, a significant correlation was found between pK(a) and mutagenic potency in TA100, even when the unionized fraction of the activity dose was considered. In addition, a correlation in CLOGP for BA-1 to BA-4 suggested a role for bioavailability in determining mutagenic potency. These results illustrate the potential value of structural analyses for exploring the relationship between chemical structure and mutational mechanisms. To our knowledge, this is the first study in which such analyses have been applied to structural analogues for which both mutagenic potency and mutation spectra date were available.     

Comparison of chaetognath mitochondrial genomes and phylogenetical implications

Two complete mitochondrial genomes (mtDNAs) of chaetognaths, Spadella cephaloptera and Paraspadella gotoi, have been recently published. These genomes are highly unusual. They are the smallest metazoan mtDNAs so far known; atp6 and atp8 genes are missing; lastly, our reanalysis has evidenced that, contrarily to what has been previously published for one sequence, both contain a unique transfer RNA (tRNA(Met)) evidencing that both have the same gene content. Indeed, even if the gene order seems very different, two gene blocks are conserved. In addition, comparison of gene arrangement suggests phylogenetical relationships between chaetognaths and some lophotrochozoa like annelids and molluscs.     

Analysis of phylogenetically reconstructed mutational spectra in human mitochondrial DNA control region

Analysis of mutations in mitochondrial DNA is an important issue in population and evolutionary genetics. To study spontaneous base substitutions in human mitochondrial DNA we reconstructed the mutational spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. An excess of pyrimidine transitions was found both in HVS I and II regions. By means of classification analysis numerous mutational hotspots were revealed in these spectra. Context analysis of hotspots revealed a complex influence of neighboring bases on mutagenesis in the HVS I region. Further statistical analysis suggested that a transient misalignment dislocation mutagenesis operating in monotonous runs of nucleotides play an important role for generating base substitutions in mitochondrial DNA and define context properties of mtDNA. Our results suggest that dislocation mutagenesis in HVS I and II is a fingerprint of errors produced by DNA polymerase gamma in the course of human mitochondrial DNA replication     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.     

Comparative analysis of the hypervariable segment 1 mutational spectra in human mitochondrial DNA phylogeographic groups

Variability of the mtDNA hypervariable segment 1 (HVS 1) nucleotide sequences belonging to 88 phylogeographic clusters characteristic for human populations of Africa, West and East Eurasia was analyzed. Statistically significant differences between distribution of mutations in mitochondrial gene pools of the human continental groups were revealed. The list of the HVS 1 nucleotide positions characterizing by instability explained by the model of mtDNA strands dislocation during the replication process is suggested. It was shown that DNA strands dislocation during mtDNA replication is one of the key mechanisms of the context-dependent mtDNA mutagenesis during the regional differentiation of human populations.     

Bayesian analysis of mutational spectra

Studies that examine both the frequency of gene mutation and the pattern or spectrum of mutational changes can be used to identify chemical mutagens and to explore the molecular mechanisms of mutagenesis. In this article, we propose a Bayesian hierarchical modeling approach for the analysis of mutational spectra. We assume that the total number of independent mutations and the numbers of mutations falling into different response categories, defined by location within a gene and/or type of alteration, follow binomial and multinomial sampling distributions, respectively. We use prior distributions to summarize past information about the overall mutation frequency and the probabilities corresponding to the different mutational categories. These priors can be chosen on the basis of data from previous studies using an approach that accounts for heterogeneity among studies. Inferences about the overall mutation frequency, the proportions of mutations in each response category, and the category-specific mutation frequencies can be based on posterior distributions, which incorporate past and current data on the mutant frequency and on DNA sequence alterations. Methods are described for comparing groups and for assessing dose-related trends. We illustrate our approach using data from the literature.     

Altering the laccase functionality by in vivo assembly of mutant libraries with different mutational spectra

The generation of diversity for directed protein evolution experiments shows an important bottleneck in the in vitro random mutagenesis protocols. Most of them are biased towards specific changes that eventually confer a predicted and conservative mutational spectrum, limiting the exploration of the vast protein space. The current work describes a simple methodology to in vivo recombine mutant libraries with different nucleotide bias created by in vitro methods. This in vivo assembly was based on the accurate physiology of Saccharomyces cerevisiae, which as host, provided its high homologous recombination frequency to shuffle the libraries in a nonmutagenic way. The fungal thermophilic laccase from Myceliophthora thermophila expressed in S. cerevisiae was submitted to this protocol under the selective pressure of high concentrations of organic solvents. Mutant 2E9 with approximately 3-fold better kinetics than parent type showed two consecutive amino acid changes (G614D -GGC/GAC- and E615K -GAG/AAG-) because of the in vivo shuffling of the mutant libraries. Both mutations are located in the C-terminal tail that is specifically processed at the Golgi during the maturation of the protein by the Kex2 protease. Notoriously, the oxygen consumption at the T2/T3 trinuclear copper cluster was altered and the catalytic copper at the T1 site was perturbed showing differences in its redox potential and geometry. The change in the isoelectric point of C-terminal extension upon mutations seems to affect the folding of the protein at the posttranslational processing steps providing new insights in the significance of the C-terminal tail for the functionality of the ascomycete laccases.     

Direct measurement of mutational spectra in humans

By combining high fidelity in vitro DNA amplification and mutant DNA sequence separation by denaturing gradient gel electrophoresis, we are able to directly observe mutational hotspots in human genomic DNA. Our technological development has progressed through the stage of identifying mutant sequences in independently derived, 6-thioguanine-resistant human B cells. We are now analyzing uncloned, complex populations derived from several thousand 6-thioguanine-resistant cells and report preliminary data concerning the mutational spectra of benzo[a]pyrene diol epoxide and ultraviolet light in exon 3 of the hypoxanthine-guanine phosphoribosyltransferase gene. In addition, the approach appears to be general for any gene sequence for which a means to select mutants exists. The more global need to eliminate phenotypic selection is, however, our primary impetus. Our analysis leads us to conclude that no known in vitro DNA polymerase has sufficient fidelity to permit direct observation of unselected mutants. Therefore, an additional change in technology will be necessary to observe nonselected mutant DNA sequences at the low frequencies found in human tissues.     

Massive gene rearrangements of mitochondrial genomes and implications for the phylogeny of Trichoptera (Insecta)

Mitochondrial genomes have been widely used for phylogenetic reconstruction and evolutionary analysis in various groups of Insecta. Gene rearrangements in the mitogenome can be informative characters for phylogenetic reconstruction and adaptive evolution. Trichoptera is one of the most important groups of aquatic insects. Prior to this study, complete mitogenomes from Trichoptera were restricted to eight families, resulting in a biased view of their mitogenome structure and evolution. Here, we assemble new mitogenomes for 66 species by high-throughput sequencing. The mitogenomes of 19 families and 47 genera are documented for the first time. Combined with 16 previously published mitogenomes of Trichoptera, we find 14 kinds of gene rearrangement patterns novel for Trichoptera, including rearrangement of protein-coding genes, tRNAs and control regions. Simultaneously, we provide evidence for the occurrence of tandem duplication and non-random loss events in the mitogenomes of three families. Phylogenetic analyses show that Hydroptilidae was recovered as a sister group to Annulipalpia. The increased nucleotide substitution rate and adaptive evolution may have affected the mitochondrial gene rearrangements in Trichoptera. Our study offers new insights into the mechanisms and patterns of mitogenome rearrangements in Insecta at large and into the usefulness of mitogenomic gene order as a phylogenetic marker within Trichoptera.     

Evolutionary implications of analyses of complete mitochondrial genomes across order Zoantharia (Cnidaria: Hexacorallia)

Cnidarians are early-diverging metazoans, but evolutionary aspects of some taxa are still poorly understood, as in the order Zoantharia (Anthozoa: Hexacorallia). Zoantharians have been divided into two suborders based on the arrangement of the fifth septae as complete (Macrocnemina) or incomplete (Brachycnemina). Previous molecular phylogenetic analyses have indicated the need for re-evaluation as Macrocnemina has been found to be paraphyletic. Despite many phylogenetic studies, the recovery of complete mitochondrial genomes (mt-genomes) for systematic and evolutionary studies of zoantharians has been limited. The present study represents the first to sequence the complete mt-genomes of members of eight of nine zoantharian families. Although all examined mt-genomes had the same gene order arrangement, there were variations among mt-genomes' sizes, nucleotide substitution rates, and introns. Only two species did not have the cox1 intron, which harbors a gene coding a homing endonuclease of the LAGLIDADG type. Our mitogenomic analyses also showed relatively high nucleotide diversity in mt-DNA regions other than the standard regions traditionally considered for DNA barcoding of this group. Phylogenetic analyses using 13 mt-genome protein-coding genes recovered a fully resolved tree with clear separation between macrocnemic representatives. Ancestral state reconstruction analyses revealed three main transitions in arrangement of the marginal musculature through the evolutionary history of the order. An “early” transition from reticulate mesogleal to a cteniform endodermal arrangement was followed by transitions that occurred in the common ancestor of the Brachycnemina and family Hydrozoanthidae. Our results indicate the need for clarification of higher-level phylogeny and taxonomy of Zoantharia.     

5-fluorouracil forward mutation assay in Salmonella: determination of mutational target and spontaneous mutational spectra

A forward mutation assay in Salmonella typhimurium that selects for 5-fluoruracil (FU) resistance has been developed. The two genes possibly involved in FU resistance, the uracil phosphoribosyl transferase gene (upp) and the uracil transport protein (uraA), have been cloned from S. typhimurium and sequenced. One hundred percent of FU-resistant clones display sequence changes in the upp gene, indicating that its loss is the major mechanism involved in FU resistance. The spontaneous mutational spectra at the upp locus were then determined in two S. typhimurium strains, FU100 and FU1535, that differ only in the presence of pKM101 plasmid. The pKM101 plasmid provides error-prone replicative bypass of DNA lesions and renders FU100 more susceptible to induced mutagenesis. Fluctuation analysis of FU-resistant clones demonstrated a 10-fold higher spontaneous mutation rate at the upp locus in FU100 relative to FU1535. Over 300 independent FU-resistant clones were then used to generate the spectra at the upp locus in both the strains. Approximately 40% of all the mutations were base substitutions, present at the same relative percentage in both the strains. Frameshift mutations also accounted for approximately 40% of the total; however, their incidence was slightly elevated in FU100. The remaining mutations were larger insertions and deletions, which were both slightly elevated in FU1535. pKM101 significantly elevated the rate of all classes of mutations at the upp locus, with profound effects on A:T to T:A transversions and -2-base frameshift mutations. These initial mutational spectra at the upp locus reveal 147 mutable sites, or 23% of the total 627-base coding sequence and suggest that the target can detect a diverse spectrum of mutagenic events.     

Complete mitochondrial genomes of three Oxya grasshoppers (Orthoptera) and their implications for phylogenetic reconstruction

Oxya is a genus of grasshoppers (Orthoptera: Acridoidea) attacking rice and other gramineous plants in Africa and Asia. In the present study, we characterized complete mitochondrial genomes (mitogenomes) of three species, Oxya japonica japonica (15,427 bp), Oxya hainanensis (15,443 bp) and Oxya agavisa robusta (15,552 bp) collected from China. The three mitogenomes contained a typical gene set of metazoan mitogenomes and shared the same gene order with other Acridid grasshoppers, including the rearrangement of tRNA^Asp and tRNA^Lys. Analyses of pairwise genetic distances showed that ATP8 was the least conserved gene, while COI the most conserved. To determine the position of Oxya grasshoppers in the phylogeny of Acrididae, we reconstructed phylogenetic trees among 64 species from across 11 subfamilies using nucleotide sequences of mitogenomes. While the tree confirms traditional classifications of Acrididae at major higher-levels, it suggests a few modifications for classifications at lower-levels.     

In vitro and in vivo mutational analysis of the 3'-terminal regions of hepatitis e virus genomes and replicons

Hepatitis E virus (HEV) replication is not well understood, mainly because the virus does not infect cultured cells efficiently. However, Huh-7 cells transfected with full-length genomes produce open reading frame 2 protein, indicative of genome replication (6). To investigate the role of 3'-terminal sequences in RNA replication, we constructed chimeric full-length genomes with divergent 3'-terminal sequences of genotypes 2 and 3 replacing that of genotype 1 and transfected them into Huh-7 cells. The production of viral proteins by these full-length chimeras was indistinguishable from that of the wild type, suggesting that replication was not impaired. In order to better quantify HEV replication in cell culture, we constructed an HEV replicon with a reporter (luciferase). Luciferase production was cap dependent and RNA-dependent RNA polymerase dependent and increased following transfection of Huh-7 cells. Replicons harboring the 3'-terminal intergenotypic chimera sequences were also assayed for luciferase production. In spite of the large sequence differences among the 3' termini of the viruses, replication of the chimeric replicons was surprisingly similar to that of the parental replicon. However, a single unique nucleotide change within a predicted stem structure at the 3' terminus substantially reduced the efficiency of replication: RNA replication was partially restored by a covariant mutation. Similar patterns of replication were obtained when full-length genomes were inoculated into rhesus macaques, suggesting that the in vitro system could be used to predict the effect of 3'-terminal mutations in vivo. Incorporation of the 3'-terminal sequences of the swine strain of HEV into the genotype 1 human strain did not enable the human strain to infect swine.     

The mitochondrial genomes of Cambaroides similis and Procambarus clarkii (Decapoda: Astacidea: Cambaridae): the phylogenetic implications for Reptantia

Kim, S., Park, M.‐H., Jung, J.‐H., Ahn, D.‐H., Sultana, T., Kim, S., Park, J.‐K., Choi, H.‐G. & Min, G.‐S. (2012). The mitochondrial genomes of Cambaroides similis and Procambarus clarkii (Decapoda: Astacidea: Cambaridae): the phylogenetic implications for Reptantia. —Zoologica Scripta, 41, 281–292. We determined the complete mitochondrial (mt) genome sequences of two northern hemisphere freshwater crayfish species, Cambaroides similis and Procambarus clarkii (Decapoda: Astacidea: Cambaridae). These species have an identical gene order with typical metazoan mt genome compositions. However, their gene arrangement was very distinctive compared with the pan‐crustacean ground pattern because of the presence of a long inverted block, which included 19 coding genes and a control region (CR). Because the CR was inverted, their nucleotide frequencies showed a reversed strand‐specific bias compared with the other decapods. Based on a comparative analysis of mt genome arrangements between southern and northern hemisphere crayfish and their putative close marine relative (Homarus americanus, a true clawed lobster), we postulated that the ancestor of freshwater crayfish had a typical pan‐crustacean mtDNA gene order, similar to its marine relatives. Based on this assumption, we traced the most parsimonious gene rearrangement scenario of the northern hemisphere crayfish. In a phylogenetic study on the infraordinal relationships in reptan decapods, the lineage Lineata [Thalassinidea (Brachyura, Anomura)] was well supported, while the infraorder positions of Achelata and Astacidea remained unidentified.     

Regression trees for analysis of mutational spectra in nucleotide sequences.

MOTIVATION: The study and comparison of mutational spectra is an important problem in molecular biology, because these spectra often reveal important features of the action of various mutagens and the functioning of repair/replication enzymes. As is known, mutability varies significantly along nucleotide sequences: mutations often concentrate at certain positions in a sequence, otherwise termed 'hotspots'. RESULTS: Herein, we propose a regression analysis method based on the use of regression trees in order to analyse the influence of nucleotide context on the occurrence of such hotspots. The REGRT program developed has been tested on simulated and real mutational spectra. For the G:C-->T:A mutational spectra induced by Sn1 alkylating agents (nine spectra), the prediction accuracy was 0. 99. AVAILABILITY: The REGRT program is available upon request from V.Berikov.     

Composition strand asymmetries in prokaryotic genomes: mutational bias and biased gene orientation

Most prokaryotic genomes display strand compositional asymmetries, but the reasons for these biases remain unclear. When the distribution of gene orientation is biased, as it often is, this may induce a bias in composition, as codon frequencies are not identical. We show here that this effect can be estimated and removed, and that the residual base skews are the highest at third base codon positions and lower at first and second positions. This strongly suggests that compositional asymmetries result from 1) a replication-related mutational bias that is filtered through selective pressure and/or from 2) an uneven distribution of gene orientation. In most cases, the mutational bias alters the codon usage and amino acid frequencies of the leading and the lagging strand. However, these features are not ubiquitous amongst prokaryotes, and the biological reasons for them remain to be found.