The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana

Next-generation sequencing as a tool to study microbial evolution

Thanks to their short generation times and large population sizes, microbes evolve rapidly. Evolutionary biologists have exploited this to observe evolution in real time. The falling costs of whole-genome sequencing using next-generation technologies now mean that it is realistic to use this as a tool to study this rapid microbial evolution both in the laboratory and in the wild. Such experiments are being used to accurately estimate the rates of mutation, reveal the genetic targets and dynamics of natural selection, uncover the correlation (or lack thereof) between genetic and phenotypic change, and provide data to test long-standing evolutionary hypotheses. These advances have important implications for our understanding of the within- and between-host evolution of microbial pathogens.

Genome studies at the PAG 2011 conference

The contents of the plenary lectures presented at the Plant and Animal Genome (PAG) meeting in January 2011 are summarized in order to provide some insights into the advances in plant, animal and microbe genome studies as they impact on our understanding of complex biological systems. The areas of biology covered include the dynamics of genome change, biological recognition processes and the new processes that underpin investment in science. This overview does not attempt to summarize the diversity of activities that are covered during the PAG through workshops, posters and the suppliers of cutting-edge technologies, but reviews major advances in specific research areas.

Reversion-reporter transgenes to analyze all six base-substitution pathways in Arabidopsis

To expand the repertoire of Arabidopsis (Arabidopsis thaliana) mutation-reporter transgenes, we constructed six mutant alleles in the same codon of the β-glucuronidase-encoding GUS transgene. Each allele reverts to GUS+ only via a particular one of the six transition/transversion pathways. AcV5 epitope tags, fused carboxyl terminal to the inactive GUS- proteins, enabled semiquantitative immunoassays in plant protein extracts. Spontaneous G:C→T:A transversions, previously not measured using reporter transgenes, were quite frequent. This may reflect mispairing of adenine with 8-oxoguanine in DNA attacked by endogenous oxyradicals. Spontaneous G:C→A:T was modest and other reversions were relatively low, as reported previously. Frequencies of ultraviolet C-induced TT→TC and TC→TT reversions were both high. With increased transgene copy number, spontaneous G:C→T:A reversions increased but ultraviolet C-induced reversions decreased. Frequencies of some reversion events were reduced among T4 versus T3 generation plants. Based on these and other analyses of sources of experimental variation, we propose guidelines for the employment of these lines to study genotoxic stress in planta.

Evolution and functional diversification of MIRNA genes

MicroRNAs (miRNAs) are small regulatory RNAs found in diverse eukaryotic lineages. In plants, a minority of annotated MIRNA gene families are conserved between plant families, while the majority are family- or species-specific, suggesting that most known MIRNA genes arose relatively recently in evolutionary time. Given the high proportion of young MIRNA genes in plant species, new MIRNA families are likely spawned and then lost frequently. Unlike highly conserved, ancient miRNAs, young miRNAs are often weakly expressed, processed imprecisely, lack targets, and display patterns of neutral variation, suggesting that young MIRNA loci tend to evolve neutrally. Genome-wide analyses from several plant species have revealed that variation in miRNA foldback expression, structure, processing efficiency, and miRNA size have resulted in the unique functionality of MIRNA loci and resulting miRNAs. Additionally, some miRNAs have evolved specific properties and functions that regulate other transcriptional or posttranscriptional silencing pathways. The evolution of miRNA processing and functional diversity underscores the dynamic nature of miRNA-based regulation in complex regulatory networks.

The composition and origins of genomic variation among individuals of the soybean reference cultivar Williams 82

Soybean (Glycine max) is a self-pollinating species that has relatively low nucleotide polymorphism rates compared with other crop species. Despite the low rate of nucleotide polymorphisms, a wide range of heritable phenotypic variation exists. There is even evidence for heritable phenotypic variation among individuals within some cultivars. Williams 82, the soybean cultivar used to produce the reference genome sequence, was derived from backcrossing a Phytophthora root rot resistance locus from the donor parent Kingwa into the recurrent parent Williams. To explore the genetic basis of intracultivar variation, we investigated the nucleotide, structural, and gene content variation of different Williams 82 individuals. Williams 82 individuals exhibited variation in the number and size of introgressed Kingwa loci. In these regions of genomic heterogeneity, the reference Williams 82 genome sequence consists of a mosaic of Williams and Kingwa haplotypes. Genomic structural variation between Williams and Kingwa was maintained between the Williams 82 individuals within the regions of heterogeneity. Additionally, the regions of heterogeneity exhibited gene content differences between Williams 82 individuals. These findings show that genetic heterogeneity in Williams 82 primarily originated from the differential segregation of polymorphic chromosomal regions following the backcross and single-seed descent generations of the breeding process. We conclude that soybean haplotypes can possess a high rate of structural and gene content variation, and the impact of intracultivar genetic heterogeneity may be significant. This detailed characterization will be useful for interpreting soybean genomic data sets and highlights important considerations for research communities that are developing or utilizing a reference genome sequence.

Transposable elements and small RNAs contribute to gene expression divergence between Arabidopsis thaliana and Arabidopsis lyrata

Transposable elements (TEs) are often the primary determinant of genome size differences among eukaryotes. In plants, the proliferation of TEs is countered through epigenetic silencing mechanisms that prevent mobility. Recent studies using the model plant Arabidopsis thaliana have revealed that methylated TE insertions are often associated with reduced expression of nearby genes, and these insertions may be subject to purifying selection due to this effect. Less is known about the genome-wide patterns of epigenetic silencing of TEs in other plant species. Here, we compare the 24-nt siRNA complement from A. thaliana and a closely related congener with a two- to threefold higher TE copy number, Arabidopsis lyrata. We show that TEs--particularly siRNA-targeted TEs--are associated with reduced gene expression within both species and also with gene expression differences between orthologs. In addition, A. lyrata TEs are targeted by a lower fraction of uniquely matching siRNAs, which are associated with more effective silencing of TE expression. Our results suggest that the efficacy of RNA-directed DNA methylation silencing is lower in A. lyrata, a finding that may shed light on the causes of differential TE proliferation among species.

Population genetics of genomics-based crop improvement methods

Many genome-wide association studies (GWAS) in humans are concluding that, even with very large sample sizes and high marker densities, most of the genetic basis of complex traits may remain unexplained. At the same time, recent research in plant GWAS is showing much greater success with fewer resources. Both GWAS and genomic selection (GS), a method for predicting phenotypes by the use of genome-wide marker data, are receiving considerable attention among plant breeders. In this review we explore how differences in population genetic histories, as well as past selection for traits of interest, have produced trait architectures and patterns of linkage disequilibrium (LD) that frequently differ dramatically between domesticated plants and humans, making detection of quantitative trait loci (QTL) effects in crops more rewarding and less costly than in humans.

The case for resequencing studies of Arabidopsis thaliana accessions: mining the dark matter of natural genetic variation

Ultra-high-throughput sequencing (UHTS) techniques are evolving rapidly and may soon become an affordable and routine tool for sequencing plant DNA, even in smaller plant biology labs. Here we review recent insights into intraspecific genome variation gained from UHTS, which offers a glimpse of the rather unexpected levels of structural variability among Arabidopsis thaliana accessions. The challenges that will need to be addressed to efficiently assemble and exploit this information are also discussed.

Functional conservation of nucleosome formation selectively biases presumably neutral molecular variation in yeast genomes

One prominent pattern of mutational frequency, long appreciated in comparative genomics, is the bias of purine/pyrimidine conserving substitutions (transitions) over purine/pyrimidine altering substitutions (transversions). Traditionally, this transitional bias has been thought to be driven by the underlying rates of DNA mutation and/or repair. However, recent sequencing studies of mutation accumulation lines in model organisms demonstrate that substitutions generally do not accumulate at rates that would indicate a transitional bias. These observations have called into question a very basic assumption of molecular evolution; that naturally occurring patterns of molecular variation in noncoding regions accurately reflect the underlying processes of randomly accumulating neutral mutation in nuclear genomes. Here, in Saccharomyces yeasts, we report a very strong inverse association (r = −0.951, P < 0.004) between the genome-wide frequency of substitutions and their average energetic effect on nucleosome formation, as predicted by a structurally based energy model of DNA deformation around the nucleosome core. We find that transitions occurring at sites positioned nearest the nucleosome surface, which are believed to function most importantly in nucleosome formation, alter the deformation energy of DNA to the nucleosome core by only a fraction of the energy changes typical of most transversions. When we examined the same substitutions set against random background sequences as well as an existing study reporting substitutions arising in mutation accumulation lines of Saccharomyces cerevisiae, we failed to find a similar relationship. These results support the idea that natural selection acting to functionally conserve chromatin organization may contribute significantly to genome-wide transitional bias, even in noncoding regions. Because nucleosome core structure is highly conserved across eukaryotes, our observations may also help to further explain locally elevated transition bias at CpG islands, which are known to destabilize nucleosomes at vertebrate promoters.

Phylogenomic analysis of the uracil-DNA glycosylase superfamily

The spontaneous deamination of cytosine produces uracil mispaired with guanine in DNA, which will produce a mutation, unless repaired. In all domains of life, uracil-DNA glycosylases (UDGs) are responsible for the elimination of uracil from DNA. Thus, UDGs contribute to the integrity of the genetic information and their loss results in mutator phenotypes. We are interested in understanding the role of UDG genes in the evolutionary variation of the rate and the spectrum of spontaneous mutations. To this end, we determined the presence or absence of the five main UDG families in more than 1,000 completely sequenced genomes and analyzed their patterns of gene loss and gain in eubacterial lineages. We observe nonindependent patterns of gene loss and gain between UDG families in Eubacteria, suggesting extensive functional overlap in an evolutionary timescale. Given that UDGs prevent transitions at G:C sites, we expected the loss of UDG genes to bias the mutational spectrum toward a lower equilibrium G + C content. To test this hypothesis, we used phylogenetically independent contrasts to compare the G + C content at intergenic and 4-fold redundant sites between lineages where UDG genes have been lost and their sister clades. None of the main UDG families present in Eubacteria was associated with a higher G + C content at intergenic or 4-fold redundant sites. We discuss the reasons of this negative result and report several features of the evolution of the UDG superfamily with implications for their functional study. uracil-DNA glycosylase, mutation rate evolution, mutational bias, GC content, DNA repair, mutator gene.

Atypical E2F activity coordinates PHR1 photolyase gene transcription with endoreduplication onset

Because of their sessile life style, plants have evolved the ability to adjust to environmentally harsh conditions. An important aspect of stress adaptation involves the reprogramming of the cell cycle to ensure optimal growth. The atypical E2F transcription factor DP-E2F-like 1 (E2Fe/DEL1) had been found previously to be an important regulator of the endocycle onset. Here, a novel role for E2Fe/DEL1 was identified as a transcriptional repressor of the type-II cyclobutane pyrimidine dimer-photolyase DNA repair gene PHR1. Upon ultraviolet-B (UV-B) treatment, plants knocked out for E2Fe/DEL1 had improved DNA repair abilities when compared with control plants, whereas those overexpressing it performed less well. Better DNA repair allowed E2Fe/DEL1 knockout plants to resume endoreduplication faster than control plants, contributing in this manner to UV-B radiation resistance by compensating the stress-induced reduction in cell number by ploidy-dependent cell growth. As E2Fe/DEL1 levels decreased upon UV-B treatment, we hypothesize that the coordinated transcriptional induction of PHR1 with the endoreduplication onset contributes to the adaptation of plants exposed to UV-B stress.

Type II protein arginine methyltransferase 5 (PRMT5) is required for circadian period determination in Arabidopsis thaliana

Posttranslational modification is an important element in circadian clock function from cyanobacteria through plants and mammals. For example, a number of key clock components are phosphorylated and thereby marked for subsequent ubiquitination and degradation. Through forward genetic analysis we demonstrate that protein arginine methyltransferase 5 (PRMT5; At4g31120) is a critical determinant of circadian period in Arabidopsis. PRMT5 is coregulated with a set of 1,253 genes that shows alterations in phase of expression in response to entrainment to thermocycles versus photocycles in constant temperature. PRMT5 encodes a type II protein arginine methyltransferase that catalyzes the symmetric dimethylation of arginine residues (Rsme2). Rsme2 modification has been observed in many taxa, and targets include histones, components of the transcription complex, and components of the spliceosome. Neither arginine methylation nor PRMT5 has been implicated previously in circadian clock function, but the period lengthening associated with mutational disruption of prmt5 indicates that Rsme2 is a decoration important for the Arabidopsis clock and possibly for clocks in general.

High throughput DNA sequencing: The new sequencing revolution

Improvements in technology have rapidly changed the field of DNA sequencing. These improvements are boosted by bio-medical research. Plant science has benefited from this breakthrough, and a number of plant genomes are now available, new biological questions can be approached and new breeding strategies can be designed. The first part of this review aims to briefly describe the principles of the new sequencing methods, many of which are already used in plant laboratories. The second part summarizes the state of plant genome sequencing and illustrates the achievements in the last few years. Although already impressive, these results represent only the beginning of a new genomic era in plant science. Finally we describe some of the exciting discoveries in the structure and evolution of plant genomes made possible by genome sequencing in terms of biodiversity, genome expression and epigenetic regulations. All of these findings have already influenced plant breeding and biodiversity protection. Finally we discuss current trends, challenges and perspectives.

The population genomics of plant adaptation

There has been an enormous increase in the amount of data on DNA sequence polymorphism available for many organisms in the last decade. New sequencing technologies provide great potential for investigating natural selection in plants using population genomic approaches. However, plant populations frequently show significant departures from the assumptions of standard models used to detect selection and many forms of directional selection do not fit with classical population genetics theory. Here, we explore the extent to which plant populations show departures from standard model assumptions, and the implications this has for detecting selection on molecular variation. A growing number of multilocus studies of nucleotide variation suggest that changes in population size, particularly bottlenecks, and strong subdivision may be common in plants. This demographic variation presents important challenges for models used to infer selection. In addition, selection from standing genetic variation and multiple independent adaptive substitutions can further complicate efforts to understand the nature of selection. We discuss emerging patterns from plant studies and propose that, rather than treating population history as a nuisance variable when testing for selection, the interaction between demography and selection is of fundamental importance for evolutionary studies of plant populations using molecular data.

Germ cell mutagens: risk assessment challenges in the 21st century

Heritable mutations may result in a wide variety of detrimental outcomes, from embryonic lethality to genetic disease in the offspring. Despite this, today's commonly used test batteries do not include assays for germ cell mutation. Current challenges include a lack of practical assays and concrete evidence for human germline mutagens, and large data gaps that often impede risk assessment. Moreover, most regulatory assessments are based on the assumption that somatic cell mutation assays also protect the germline by default, which has not been adequately confirmed. The field is also faced with new challenges aimed at dramatically reducing animal testing, and attempts to rapidly classify thousands of chemicals using high throughput in vitro assays. These approaches may not adequately capture effects that may be particular to gametes, since many aspects of the germline are unique. In light of these challenges, an urgent need exists to develop new approaches to evaluate the potential of toxicants to cause germline mutation. The application of new technologies will greatly enhance our understanding of mutation in humans exposed to environmental mutagens. However, we must be poised to collect and interpret these data, and facilitate risk translation to regulators and the public. Genetic toxicologists must also become actively involved in the development of high-throughput tools to study germline mutation. Appropriate attention to these areas will result in the development of policies that prioritize the protection of the germline and future generations from DNA sequence mutations.

Evidence that mutation is universally biased towards AT in bacteria

Mutation is the engine that drives evolution and adaptation forward in that it generates the variation on which natural selection acts. Mutation is a random process that nevertheless occurs according to certain biases. Elucidating mutational biases and the way they vary across species and within genomes is crucial to understanding evolution and adaptation. Here we demonstrate that clonal pathogens that evolve under severely relaxed selection are uniquely suitable for studying mutational biases in bacteria. We estimate mutational patterns using sequence datasets from five such clonal pathogens belonging to four diverse bacterial clades that span most of the range of genomic nucleotide content. We demonstrate that across different types of sites and in all four clades mutation is consistently biased towards AT. This is true even in clades that have high genomic GC content. In all studied cases the mutational bias towards AT is primarily due to the high rate of C/G to T/A transitions. These results suggest that bacterial mutational biases are far less variable than previously thought. They further demonstrate that variation in nucleotide content cannot stem entirely from variation in mutational biases and that natural selection and/or a natural selection-like process such as biased gene conversion strongly affect nucleotide content.

The origin of the non-recombining region of sex chromosomes in Carica and Vasconcellea

Carica and Vasconcellea are two closely related sister genera in the family Caricaceae, and were once classified as two sections under Carica. Sex chromosomes have been found in papaya and originated approximately 2-3 million years ago. The objectives of this study were to determine whether sex chromosomes have evolved in Vasconcellea. Six X/Y gene pairs were cloned, sequenced and analyzed from three dioecious, one trioecious and one monoecious species of Vasconcellea. The isolation of distinctive X and Y alleles in dioecious and trioecious species of Vasconcellea demonstrated that sex chromosomes have evolved in this genus. Phylogenetic analyses indicated a monophyletic relationship between the X/Y alleles of Carica and those of Vasconcellea. Distinctive clusters of X/Y alleles were documented in V. parviflora and V. pulchra for all available gene sequences, and in V. goudatinana and V. cardinamarcensis for some X/Y alleles. The X and Y alleles within each species shared most single nucleotide polymorphism haplotypes that differed from other species. Limited evidence of gene conversion was documented among the X/Y alleles of some species, but was not sufficient to cause the evolutionary patterns reported herein. The Carica and Vasconcellea sex chromosomes may have originated from the same autosomes bearing the X allelic form that still exist in the monoecious species V. monoica, and have evolved independently after the speciation event that separated Carica from Vasconcellea. Within Vasconcellea, sex chromosomes have evolved at the species level, at least for some species.

Are we Genomic Mosaics? Variations of the Genome of Somatic Cells can Contribute to Diversify our Phenotypes

Theoretical and experimental evidences support the hypothesis that the genomes and the epigenomes may be different in the somatic cells of complex organisms. In the genome, the differences range from single base substitutions to chromosome number; in the epigenome, they entail multiple postsynthetic modifications of the chromatin. Somatic genome variations (SGV) may accumulate during development in response both to genetic programs, which may differ from tissue to tissue, and to environmental stimuli, which are often undetected and generally irreproducible. SGV may jeopardize physiological cellular functions, but also create novel coding and regulatory sequences, to be exposed to intraorganismal Darwinian selection. Genomes acknowledged as comparatively poor in genes, such as humans', could thus increase their pristine informational endowment. A better understanding of SGV will contribute to basic issues such as the "nature vs nurture" dualism and the inheritance of acquired characters. On the applied side, they may explain the low yield of cloning via somatic cell nuclear transfer, provide clues to some of the problems associated with transdifferentiation, and interfere with individual DNA analysis. SGV may be unique in the different cells types and in the different developmental stages, and thus explain the several hundred gaps persisting in the human genomes "completed" so far. They may compound the variations associated to our epigenomes and make of each of us an "(epi)genomic" mosaic. An ensuing paradigm is the possibility that a single genome (the ephemeral one assembled at fertilization) has the capacity to generate several different brains in response to different environments.

Evolution of the mutation rate

Understanding the mechanisms of evolution requires information on the rate of appearance of new mutations and their effects at the molecular and phenotypic levels. Although procuring such data has been technically challenging, high-throughput genome sequencing is rapidly expanding knowledge in this area. With information on spontaneous mutations now available in a variety of organisms, general patterns have emerged for the scaling of mutation rate with genome size and for the likely mechanisms that drive this pattern. Support is presented for the hypothesis that natural selection pushes mutation rates down to a lower limit set by the power of random genetic drift rather than by intrinsic physiological limitations, and that this has resulted in reduced levels of replication, transcription, and translation fidelity in eukaryotes relative to prokaryotes.

MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana

MicroRNAs (miRNAs) are short regulatory RNAs processed from partially self-complementary foldbacks within longer MIRNA primary transcripts. Several MIRNA families are conserved deeply through land plants, but many are present only in closely related species or are species specific. The finding of numerous evolutionarily young MIRNA, many with low expression and few if any targets, supports a rapid birth-death model for MIRNA evolution. A systematic analysis of MIRNA genes and families in the close relatives, Arabidopsis thaliana and Arabidopsis lyrata, was conducted using both whole-genome comparisons and high-throughput sequencing of small RNAs. Orthologs of 143 A. thaliana MIRNA genes were identified in A. lyrata, with nine having significant sequence or processing changes that likely alter function. In addition, at least 13% of MIRNA genes in each species are unique, despite their relatively recent speciation (approximately 10 million years ago). Alignment of MIRNA foldbacks to the Arabidopsis genomes revealed evidence for recent origins of 32 families by inverted or direct duplication of mostly protein-coding gene sequences, but less than half of these yield miRNA that are predicted to target transcripts from the originating gene family. miRNA nucleotide divergence between A. lyrata and A. thaliana orthologs was higher for young MIRNA genes, consistent with reduced purifying selection compared with deeply conserved MIRNA genes. Additionally, target sites of younger miRNA were lost more frequently than for deeply conserved families. In summary, our systematic analyses emphasize the dynamic nature of the MIRNA complement of plant genomes.

Local-scale patterns of genetic variability, outcrossing, and spatial structure in natural stands of Arabidopsis thaliana

As Arabidopsis thaliana is increasingly employed in evolutionary and ecological studies, it is essential to understand patterns of natural genetic variation and the forces that shape them. Previous work focusing mostly on global and regional scales has demonstrated the importance of historical events such as long-distance migration and colonization. Far less is known about the role of contemporary factors or environmental heterogeneity in generating diversity patterns at local scales. We sampled 1,005 individuals from 77 closely spaced stands in diverse settings around Tübingen, Germany. A set of 436 SNP markers was used to characterize genome-wide patterns of relatedness and recombination. Neighboring genotypes often shared mosaic blocks of alternating marker identity and divergence. We detected recent outcrossing as well as stretches of residual heterozygosity in largely homozygous recombinants. As has been observed for several other selfing species, there was considerable heterogeneity among sites in diversity and outcrossing, with rural stands exhibiting greater diversity and heterozygosity than urban stands. Fine-scale spatial structure was evident as well. Within stands, spatial structure correlated negatively with observed heterozygosity, suggesting that the high homozygosity of natural A. thaliana may be partially attributable to nearest-neighbor mating of related individuals. The large number of markers and extensive local sampling employed here afforded unusual power to characterize local genetic patterns. Contemporary processes such as ongoing outcrossing play an important role in determining distribution of genetic diversity at this scale. Local "outcrossing hotspots" appear to reshuffle genetic information at surprising rates, while other stands contribute comparatively little. Our findings have important implications for sampling and interpreting diversity among A. thaliana accessions.

Rate, molecular spectrum, and consequences of human mutation

Although mutation provides the fuel for phenotypic evolution, it also imposes a substantial burden on fitness through the production of predominantly deleterious alleles, a matter of concern from a human-health perspective. Here, recently established databases on de novo mutations for monogenic disorders are used to estimate the rate and molecular spectrum of spontaneously arising mutations and to derive a number of inferences with respect to eukaryotic genome evolution. Although the human per-generation mutation rate is exceptionally high, on a per-cell division basis, the human germline mutation rate is lower than that recorded for any other species. Comparison with data from other species demonstrates a universal mutational bias toward A/T composition, and leads to the hypothesis that genome-wide nucleotide composition generally evolves to the point at which the power of selection in favor of G/C is approximately balanced by the power of random genetic drift, such that variation in equilibrium genome-wide nucleotide composition is largely defined by variation in mutation biases. Quantification of the hazards associated with introns reveals that mutations at key splice-site residues are a major source of human mortality. Finally, a consideration of the long-term consequences of current human behavior for deleterious-mutation accumulation leads to the conclusion that a substantial reduction in human fitness can be expected over the next few centuries in industrialized societies unless novel means of genetic intervention are developed.