Indirect mechanisms of genomic instability and the biological significance of mutations at tandem repeat loci

Radiation induction of genomic instability has two features: induction of untargeted mutation and delayed mutation. These phenomena have been studied mostly in tissue culture cells, but analyses have also been conducted in whole body systems. The study of response in whole body systems frequently applies repeat sequences as markers to detect mutations. These studies have generated conflicting findings. In addition, lack of knowledge of the mechanisms involved in repeat mutation confounds the interpretation of the biological significance of increased rates of repeat mutation. In this review, some of the existing controversies of genomic instability are discussed in relation to the mechanism of repeat mutation. Analyses of published and unpublished studies indicate a mechanistic similarity between radiation-induced genomic instability at repeat loci and dynamic mutations of triplet repeats. Because of their repetitive nature, repeat sequences frequently block progression of replication forks and are consequently resolved by slippage and/or recombination. Irradiation of cells induces S checkpoints and promotes slippage/recombination mediated repeat mutations. Thus, genomic instability at repeat loci might be viewed as a consequence of cellular attempts to restore the stability of replication in the face of the stalled replication fork; this process can occur both spontaneously as well as after exposure to radiation.     

Automated validation of genetic variants from large databases: ensuring that variant references refer to the same genomic locations

SUMMARY: Accurate annotations of genomic variants are necessary to achieve full-genome clinical interpretations that are scientifically sound and medically relevant. Many disease associations, especially those reported before the completion of the HGP, are limited in applicability because of potential inconsistencies with our current standards for genomic coordinates, nomenclature and gene structure. In an effort to validate and link variants from the medical genetics literature to an unambiguous reference for each variant, we developed a software pipeline and reviewed 68 641 single amino acid mutations from Online Mendelian Inheritance in Man (OMIM), Human Gene Mutation Database (HGMD) and dbSNP. The frequency of unresolved mutation annotations varied widely among the databases, ranging from 4 to 23%. A taxonomy of primary causes for unresolved mutations was produced. AVAILABILITY: This program is freely available from the web site (http://safegene.hms.harvard.edu/aa2nt/).     

MutDB: annotating human variation with functionally relevant data

SUMMARY: We have developed a resource, MutDB (http://mutdb.org/), to aid in determining which single nucleotide polymorphisms (SNPs) are likely to alter the function of their associated protein product. MutDB contains protein structure annotations and comparative genomic annotations for 8000 disease-associated mutations and SNPs found in the UCSC Annotated Genome and the human RefSeq gene set. MutDB provides interactive mutation maps at the gene and protein levels, and allows for ranking of their predicted functional consequences based on conservation in multiple sequence alignments. AVAILABILITY: http://mutdb.org/ Supplementary information: http://mutdb.org/about/about.html     

The evolution of mutation rates: separating causes from consequences

Natural selection can adjust the rate of mutation in a population by acting on allelic variation affecting processes of DNA replication and repair. Because mutation is the ultimate source of the genetic variation required for adaptation, it can be appealing to suppose that the genomic mutation rate is adjusted to a level that best promotes adaptation. Most mutations with phenotypic effects are harmful, however, and thus there is relentless selection within populations for lower genomic mutation rates. Selection on beneficial mutations can counter this effect by favoring alleles that raise the mutation rate, but the effect of beneficial mutations on the genomic mutation rate is extremely sensitive to recombination and is unlikely to be important in sexual populations. In contrast, high genomic mutation rates can evolve in asexual populations under the influence of beneficial mutations, but this phenomenon is probably of limited adaptive significance and represents, at best, a temporary reprieve from the continual selection pressure to reduce mutation. The physiological cost of reducing mutation below the low level observed in most populations may be the most important factor in setting the genomic mutation rate in sexual and asexual systems, regardless of the benefits of mutation in producing new adaptive variation. Maintenance of mutation rates higher than the minimum set by this "cost of fidelity" is likely only under special circumstances.     

OmicBrowse: a browser of multidimensional omics annotations

OmicBrowse is a browser to explore multiple datasets coordinated in the multidimensional omic space integrating omics knowledge ranging from genomes to phenomes and connecting evolutional correspondences among multiple species. OmicBrowse integrates multiple data servers into a single omic space through secure peer-to-peer server communications, so that a user can easily obtain an integrated view of distributed data servers, e.g. an integrated view of numerous whole-genome tiling-array data retrieved from a user's in-house private-data server, along with various genomic annotations from public internet servers. OmicBrowse is especially appropriate for positional-cloning purposes. It displays both genetic maps and genomic annotations within wide chromosomal intervals and assists a user to select candidate genes by filtering their annotations or associated documents against user-specified keywords or ontology terms. We also show that an omic-space chart effectively represents schemes for integrating multiple datasets of multiple species. Availability: OmicBrowse is developed by the Genome-Phenome Superbrain Project and is released as free open-source software under the GNU General Public License at http://omicspace.riken.jp.     

RNA mutations of prox1 detected in human esophageal cancer cells by the shifted termination assay

We have recently reported a novel finding that a candidate tumor suppressor gene prox1 suffered adenosine-to-inosine (A-to-I) RNA mutation without genomic mutation in a subset of human cancer cells and lost its function. Hence, screening of mutations in both cDNA and genomic DNA could be important in the analysis of causes for cancers. Here, we applied a sensitive, accurate, and simple method, called shifted termination assay (STA) for detection of an A-to-I RNA mutation (R334G) in prox1. We prepared PCR-amplified samples containing the target base of RNA mutation from cDNAs and genomic DNAs of various cell lines and clinical samples, to demonstrate that the STA method can be used to identify not only genomic mutations but also RNA mutations more effectively compared to sequencing. By means of STA, we found prox1 R334G RNA mutations but not genomic DNA mutations in 4 of 8 cases of esophageal cancers. This method can help us to detect RNA mutation effectively and progress research of a potential oncogenic principle.     

A workflow for mutation extraction and structure annotation

Rich information on point mutation studies is scattered across heterogeneous data sources. This paper presents an automated workflow for mining mutation annotations from full-text biomedical literature using natural language processing (NLP) techniques as well as for their subsequent reuse in protein structure annotation and visualization. This system, called mSTRAP (Mutation extraction and STRucture Annotation Pipeline), is designed for both information aggregation and subsequent brokerage of the mutation annotations. It facilitates the coordination of semantically related information from a series of text mining and sequence analysis steps into a formal OWL-DL ontology. The ontology is designed to support application-specific data management of sequence, structure, and literature annotations that are populated as instances of object and data type properties. mSTRAPviz is a subsystem that facilitates the brokerage of structure information and the associated mutations for visualization. For mutated sequences without any corresponding structure available in the Protein Data Bank (PDB), an automated pipeline for homology modeling is developed to generate the theoretical model. With mSTRAP, we demonstrate a workable system that can facilitate automation of the workflow for the retrieval, extraction, processing, and visualization of mutation annotations -- tasks which are well known to be tedious, time-consuming, complex, and error-prone. The ontology and visualization tool are available at (http://datam.i2r.a-star.edu.sg/mstrap).     

GeCo++: a C++ library for genomic features computation and annotation in the presence of variants

We propose a C++ class library developed to the purpose of making the implementation of sequence analysis algorithms easier and faster when genomic annotations and variations need to be considered. The library provides a class hierarchy to seamlessly bind together annotations of genomic elements to sequences and to algorithm results; it allows to evaluate the effect of mutations/variations in terms of both element position shifts and of algorithm results, limiting recalculation to the minimum. Particular care has been posed to keep memory and time overhead into acceptable limits. AVAILABILITY AND IMPLEMENTATION: A complete tutorial as well as a detailed doxygen generated documentation and source code is freely available at http://bioinformatics.emedea.it/geco, under the GPL license. The library was written in standard ISO C++, and does not depend on external libraries.     

Reconstructing the ancestral germ line methylation state of young repeats

One of the key objectives of comparative genomics is the characterization of the forces that shape genomes over the course of evolution. In the last decades, evidence has been accumulated that for vertebrate genomes also epigenetic modifications have to be considered in this context. Especially, the elevated mutation frequency of 5-methylcytosine (5mC) is assumed to facilitate the depletion of CpG dinucleotides in species that exhibit global DNA methylation. For instance, the underrepresentation of CpG dinucleotides in many mammalian genomes is attributed to this effect, which is only neutralized in so-called CpG islands (CGIs) that are preferentially unmethylated and thus partially protected from rapid CpG decay. For primate-specific CpG-rich transposable elements from the ALU family, it is unclear whether their elevated CpG frequency is caused by their small age or by the absence of DNA methylation. In consequence, these elements are often misclassified in CGI annotations. We present a method for the estimation of germ line methylation from pairwise ancestral-descendant alignments. The approach is validated in a simulation study and tested on DNA repeats from the AluSx family. We conclude that a predicted unmethylated state in the germ line is highly correlated with epigenetic activity of the respective genomic region. Thus, CpG-rich repeats can be facilitated as in silico probes for the epigenetic potential of their genomic neighborhood.     

DNannotator: Annotation software tool kit for regional genomic sequences

Sequence annotation is essential for genomics-based research. Investigators of a specific genomic region who have developed abundant local discoveries such as genes and genetic markers, or have collected annotations from multiple resources, can be overwhelmed by the difficulty in creating local annotation and the complexity of integrating all the annotations. Presenting such integrated data in a form suitable for data mining and high-throughput experimental design is even more daunting. DNannotator, a web application, was designed to perform batch annotation on a sizeable genomic region. It takes annotation source data, such as SNPs, genes, primers, and so on, prepared by the end-user and/or a specified target of genomic DNA, and performs de novo annotation. DNannotator can also robustly migrate existing annotations in GenBank format from one sequence to another. Annotation results are provided in GenBank format and in tab-delimited text, which can be imported and managed in a database or spreadsheet and combined with existing annotation as desired. Graphic viewers, such as Genome Browser or Artemis, can display the annotation results. Reference data (reports on the process) facilitating the user's evaluation of annotation quality are optionally provided. DNannotator can be accessed at http://sky.bsd.uchicago.edu/DNannotator.htm.     

The Chlamydomonas Genome Project,version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

A major update to the Chlamydomonas reference genome and structural annotations is presented, revealing large genomic changes that have occurred in laboratory culture.

IN A NUTSHELL Background: Chlamydomonas reinhardtii (Chlamydomonas) is an important reference organism. The first draft genome for the species was sequenced 20 years ago. The current assembly, version 5, contains many gaps and some misassemblies. Although the structural annotations are generally of high quality, some gene models are missing, and many genes have misleading names. Finally, the reference genome has always been based on CC-503, a mating-type plus strain that was mutagenized to achieve a cell wall-less phenotype. Question: We aimed to update the Chlamydomonas reference genome and annotations using long-read sequencing. We also sequenced a second strain, the mating-type minus CC-4532. Via comparison, we tested whether the mutagenesis of CC-503 resulted in any gene-disrupting structural mutations. Findings: We produced highly contiguous genome assemblies for CC-503 and CC-4532 that were validated by linkage analyses. Most of the filled gaps were in genic regions, leading to substantial improvements in the gene models. Gene symbols were manually overhauled, providing a reliable nomenclature that can serve as a template for green algal genome projects. We discovered that the CC-503 genome harbors many large mutations and is unstable, making it an unsuitable reference. Genomic instability may stem from the deletion of the helicase RECQ3. We also fnd that the CC-4532 genome is experiencing an ongoing proliferation of transposable elements. We expect that all strains carry some large laboratory mutations. Next steps: The CC-4532 genome is presented as the new Chlamydomonas reference. The assembly and annotation bring great advancements and will serve as a foundation for future updates. However, no single strain can provide a perfect reference, and we anticipate a truly representative Chlamydomonas pan-genome.     

The UCSC Genome Browser Database

The University of California Santa Cruz (UCSC) Genome Browser Database is an up to date source for genome sequence data integrated with a large collection of related annotations. The database is optimized to support fast interactive performance with the web-based UCSC Genome Browser, a tool built on top of the database for rapid visualization and querying of the data at many levels. The annotations for a given genome are displayed in the browser as a series of tracks aligned with the genomic sequence. Sequence data and annotations may also be viewed in a text-based tabular format or downloaded as tab-delimited flat files. The Genome Browser Database, browsing tools and downloadable data files can all be found on the UCSC Genome Bioinformatics website (http://genome.ucsc.edu), which also contains links to documentation and related technical information.     

Genomic instability induced by mutations in Saccharomyces cerevisiae POL1

Mutations of chromosome replication genes can be one of the early events that promote genomic instability. Among genes that are involved in chromosomal replication, DNA polymerase alpha is essential for initiation of replication and lagging-strand synthesis. Here we examined the effect of two mutations in S. cerevisiae POL1, pol1-1 and pol1-17, on a microsatellite (GT)(16) tract. The pol1-17 mutation elevated the mutation rate 13-fold by altering sequences both inside and downstream of the (GT)(16) tract, whereas the pol1-1 mutation increased the mutation rate 54-fold by predominantly altering sequences downstream of the (GT)(16) tract in a RAD52-dependent manner. In a rad52 null mutant background pol1-1 and pol1-17 also exhibited different plasmid and chromosome loss phenotypes. Deletions of mismatch repair (MMR) genes induce a differential synergistic increase in the mutation rates of pol1-1 and pol1-17. These findings suggest that perturbations of DNA replication in these two pol1 mutants are caused by different mechanisms, resulting in various types of mutations. Thus, mutations of POL1 can induce a variety of mutator phenotypes and can be a source of genomic instability in cells.     

Avoiding Dangerous Missense: Thermophiles Display Especially Low Mutation Rates

Rates of spontaneous mutation have been estimated under optimal growth conditions for a variety of DNA-based microbes, including viruses, bacteria, and eukaryotes. When expressed as genomic mutation rates, most of the values were in the vicinity of 0.003–0.004 with a range of less than two-fold. Because the genome sizes varied by roughly 10⁴-fold, the mutation rates per average base pair varied inversely by a similar factor. Even though the commonality of the observed genomic rates remains unexplained, it implies that mutation rates in unstressed microbes reach values that can be finely tuned by evolution. An insight originating in the 1920s and maturing in the 1960s proposed that the genomic mutation rate would reflect a balance between the deleterious effect of the average mutation and the cost of further reducing the mutation rate. If this view is correct, then increasing the deleterious impact of the average mutation should be countered by reducing the genomic mutation rate. It is a common observation that many neutral or nearly neutral mutations become strongly deleterious at higher temperatures, in which case they are called temperature-sensitive mutations. Recently, the kinds and rates of spontaneous mutations were described for two microbial thermophiles, a bacterium and an archaeon. Using an updated method to extrapolate from mutation-reporter genes to whole genomes reveals that the rate of base substitutions is substantially lower in these two thermophiles than in mesophiles. This result provides the first experimental support for the concept of an evolved balance between the total genomic impact of mutations and the cost of further reducing the basal mutation rate.     

Detection of eight BRCA1 mutations in 10 breast/ovarian cancer families, including 1 family with male breast cancer.

Genetic epidemiological evidence suggests that mutations in BRCA1 may be responsible for approximately one half of early onset familial breast cancer and the majority of familial breast/ovarian cancer. The recent cloning of BRCA1 allows for the direct detection of mutations, but the feasibility of presymptomatic screening for cancer susceptibility is unknown. We analyzed genomic DNA from one affected individual from each of 24 families with at least three cases of ovarian or breast cancer, using SSCP assays. Variant SSCP bands were subcloned and sequenced. Allele-specific oligonucleotide hybridization was used to verify sequence changes and to screen DNA from control individuals. Six frameshift and two missense mutations were detected in 10 different families. A frameshift mutation was detected in a male proband affected with both breast and prostate cancer. A 40-bp deletion was detected in a patient who developed intra-abdominal carcinomatosis 1 year after prophylactic oophorectomy. Mutations were detected throughout the gene, and only one was detected in more than a single family. These results provide further evidence that inherited breast and ovarian cancer can occur as a consequence of a wide array of BRCA1 mutations. These results suggests that development of a screening test for BRCA1 mutations will be technically challenging. The finding of a mutation in a family with male breast cancer, not previously thought to be related to BRCA1, also illustrates the potential difficulties of genetic counseling for individuals known to carry mutations.     

Tissue Culture-Induced Heritable Genomic Variation in Rice,and Their Phenotypic Implications

Background

Somaclonal variation generally occurs in plants regenerated from tissue culture. However, fundamental issues regarding molecular characteristics, mutation rates and mutation spectra of plant somatic variation as well as their phenotypic relevance have been addressed only recently. Moreover, these studies have reported highly discrepant results in different plant species and even in the same plant genotype.

Methodology/principal findings

We investigated heritable genomic variation induced by tissue culture in rice by whole genome re-sequencing of an extensively selfed somaclonal line (TC-reg-2008) and its wild type (WT) donor (cv. Hitomebore). We computed the overall mutation rate, single nucleotide polymorphisms (SNPs), small scale insertions/deletions (Indels) and mobilization of transposable elements (TEs). We assessed chromosomal distribution of the various types of genomic variations, tested correlations between SNPs and Indels, and examined concomitancy between TE activity and its cytosine methylation states. We also performed gene ontology (GO) analysis of genes containing nonsynonymous mutations and large-effect mutations, and assayed effects of the genomic variations on phenotypes under both normal growing condition and several abiotic stresses. We found that heritable somaclonal genomic variation occurred extensively in rice. The genomic variations distributed non-randomly across each of the 12 rice chromosomes, and affected a large number of functional genes. The phenotypic penetrance of the genomic variations was condition-dependent.

Conclusions/significance

Tissue culture is a potent means to generate heritable genetic variations in rice, which bear distinct difference at least in space (chromosomal distribution) from those occurred under natural settings. Our findings have provided new information regarding the mutation rate and spectrum as well as chromosomal distribution pattern of somaclonal variation in rice. Our data also suggest that rice possesses a strong capacity to canalize genetic variations under normal growing conditions to maintain phenotypic robustness, which however can be released by certain abiotic stresses to generate variable phenotypes.     

Mutation accumulation in space and the maintenance of sexual reproduction

The maintenance of sexual reproduction remains one of the major puzzles of evolutionary biology, since, all else being equal, an asexual mutant should have a twofold fitness advantage over the sexual wildtype. Most theories suggest that sex helps either to purge deleterious mutations, or to adapt to changing environments. Both mechanisms have their limitations if they act in isolation because they require either high genomic mutation rates or very virulent pathogens, and it is therefore often thought that they must act together to maintain sex. Typically, however, these theories have in common that they are not based on spatial processes. Here, we show that local dispersal and local competition can explain the maintenance of sexual reproduction as a means of purging deleterious mutations. Using a spatially explicit individual-based model, we find that even with reasonably low genomic mutation rates and large total population sizes, asexual clones cannot invade a sexual population. Our results demonstrate how spatial processes affect mutation accumulation such that it can fully erode the twofold benefit of asexuality faster than an asexual clone can take over a sexual population. Thus, the cost of sex is generally overestimated in models that ignore the effects of space on mutation accumulation.     

Is the Mutation Rate Lower in Genomic Regions of Stronger Selective Constraints?

A study of the plant Arabidopsis thaliana detected lower mutation rates in genomic regions where mutations are more likely to be deleterious, challenging the principle that mutagenesis is blind to its consequence. To examine the generality of this finding, we analyze large mutational data from baker''s yeast and humans. The yeast data do not exhibit this trend, whereas the human data show an opposite trend that disappears upon the control of potential confounders. We find that the Arabidopsis study identified substantially more mutations than reported in the original data-generating studies and expected from Arabidopsis'' mutation rate. These extra mutations are enriched in polynucleotide tracts and have relatively low sequencing qualities so are likely sequencing errors. Furthermore, the polynucleotide “mutations” can produce the purported mutational trend in Arabidopsis. Together, our results do not support lower mutagenesis of genomic regions of stronger selective constraints in the plant, fungal, and animal models examined.