The Role of Small RNA-Based Epigenetic Silencing for Purifying Selection on Transposable Elements in Capsella grandiflora

To avoid negative effects of transposable element (TE) proliferation, plants epigenetically silence TEs using a number of mechanisms, including RNA-directed DNA methylation. These epigenetic modifications can extend outside the boundaries of TE insertions and lead to silencing of nearby genes, resulting in a trade-off between TE silencing and interference with nearby gene regulation. Therefore, purifying selection is expected to remove silenced TE insertions near genes more efficiently and prevent their accumulation within a population. To explore how effects of TE silencing on gene regulation shapes purifying selection on TEs, we analyzed whole genome sequencing data from 166 individuals of a large population of the outcrossing species Capsella grandiflora. We found that most TEs are rare, and in chromosome arms, silenced TEs are exposed to stronger purifying selection than those that are not silenced by 24-nucleotide small RNAs, especially with increasing proximity to genes. An age-of-allele test of neutrality on a subset of TEs supports our inference of purifying selection on silenced TEs, suggesting that our results are robust to varying transposition rates. Our results provide new insights into the processes affecting the accumulation of TEs in an outcrossing species and support the view that epigenetic silencing of TEs results in a trade-off between preventing TE proliferation and interference with nearby gene regulation. We also suggest that in the centromeric and pericentromeric regions, the negative aspects of epigenetic TE silencing are missing.

Related papers:

title: Barbara McClintock and the discovery of jumping genes
abstract: Barbara McClintock’s life shows us how important it is to nurture original and unconventional thinking in science if we are to get out of the rut of ordinariness. After a long period of relative neglect, she was awarded the Nobel Prize in 1983 for her work on genetic instability (transposition).

title: On the role of unequal exchange in the containment of transposable element copy number.
abstract: A population genetics model of the role of asymmetric pairing and unequal exchange in the stabilization of transposable element copy number in natural populations is proposed and analysed. Monte Carlo simulations indicate that the approximations incorporated into the analysis are robust in the relevant parameter ranges. Given several simple assumptions concerning transposition and excision, equal and unequal exchange, and chromosome structure, predictions of the relative numbers of transposable elements in various regions of the Drosophila melanogaster genome are compared to the observed distribution of roo / B104 elements across chromosomal regions with differing rates of exchange, and between X chromosomes and autosomes. There is no indication of an accumulation of elements in the distal regions of chromosomes, which is expected if unequal exchange is reduced concomitantly with normal crossing over in the distal regions. There is, however, an indication of an excess of elements relative to physical length in the proximal regions of the chromosomes, which also have restricted crossing over. This observation is qualitatively consistent with the model's predictions. The observed distribution of elements between the mid-sections of the X chromosomes and autosomes is consistent with the predictions of one of two models of unequal exchange.

title: The role of host factors in the population dynamics of selfish transposable elements.
abstract: Previous models of the evolution of selfish transposable genetic elements have failed to include the possibility that transposition may be limited by shortage of a host-encoded factor. The titration of host factors may be important in limiting the rate of transpositional increase in these elements. This will be exacerbated if multiple copies of the host factor protein must bind simultaneously to the target element. In the case of the Drosophila melanogaster P transposable element, which can exist as autonomous and as non-autonomous copies, there is evidence that a host-encoded protein, IRBP, is required for the transposition process. We have produced a specific model of the invasion of a host population by the P element, in which we have incorporated the requirement for the multiple binding of a host factor. We find that, for the P family, in which it is apparently transposition itself that creates selective harm to the host, the effect of selection in the context of host factor limitation is to drive up copy number. This can result in a novel high copy number-low transposition state. We also find that host factor limitation reinforces the tendency for transposable elements that create sterility to be replaced by their deletion derivatives.

title: Epigenetic silencing of transposable elements: a trade-off between reduced transposition and deleterious effects on neighboring gene expression.
abstract: Transposable elements (TEs) are ubiquitous genomic parasites. The deleterious consequences of the presence and activity of TEs have fueled debate about the evolutionary forces countering their expansion. Purifying selection is thought to purge TE insertions from the genome, and TE sequences are targeted by hosts for epigenetic silencing. However, the interplay between epigenetic and evolutionary forces countering TE expansion remains unexplored. Here we analyze genomic, epigenetic, and population genetic data from Arabidopsis thaliana to yield three observations. First, gene expression is negatively correlated with the density of methylated TEs. Second, the signature of purifying selection is detectable for methylated TEs near genes but not for unmethylated TEs or for TEs far from genes. Third, TE insertions are distributed by age and methylation status, such that older, methylated TEs are farther from genes. Based on these observations, we present a model in which host silencing of TEs near genes has deleterious effects on neighboring gene expression, resulting in the preferential loss of methylated TEs from gene-rich chromosomal regions. This mechanism implies an evolutionary tradeoff in which the benefit of TE silencing imposes a fitness cost via deleterious effects on the expression of nearby genes.

title: How important are transposons for plant evolution?
abstract: For decades, transposable elements have been known to produce a wide variety of changes in plant gene expression and function. This has led to the idea that transposable element activity has played a key part in adaptive plant evolution. This Review describes the kinds of changes that transposable elements can cause, discusses evidence that those changes have contributed to plant evolution and suggests future strategies for determining the extent to which these changes have in fact contributed to plant adaptation and evolution. Recent advances in genomics and phenomics for a range of plant species, particularly crops, have begun to allow the systematic assessment of these questions.

title: Transposable element influences on gene expression in plants.
abstract: Transposable elements (TEs) comprise a major portion of many plant genomes and bursts of TE movements cause novel genomic variation within species. In order to maintain proper gene function, plant genomes have evolved a variety of mechanisms to tolerate the presence of TEs within or near genes. Here, we review our understanding of the interactions between TEs and gene expression in plants by assessing three ways that transposons can influence gene expression. First, there is growing evidence that TE insertions within introns or untranslated regions of genes are often tolerated and have minimal impact on expression level or splicing. However, there are examples in which TE insertions within genes can result in aberrant or novel transcripts. Second, TEs can provide novel alternative promoters, which can lead to new expression patterns or original coding potential of an alternate transcript. Third, TE insertions near genes can influence regulation of gene expression through a variety of mechanisms. For example, TEs may provide novel cis-acting regulatory sites behaving as enhancers or insert within existing enhancers to influence transcript production. Alternatively, TEs may change chromatin modifications in regions near genes, which in turn can influence gene expression levels. Together, the interactions of genes and TEs provide abundant evidence for the role of TEs in changing basic functions within plant genomes beyond acting as latent genomic elements or as simple insertional mutagens. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

title: Reconstructing de novo silencing of an active plant retrotransposon
abstract: Transposable elements (TEs) contribute to genome size, organization and evolution. In plants, their activity is primarily controlled by transcriptional gene silencing (TGS), usually investigated at steady states, reflecting how long-established silent conditions are maintained, faithfully reiterated or temporarily modified. How active, invasive TEs are detected and silenced de novo in plants remains largely unknown. Using inbred lineages of hybrid Arabidopsis thaliana epigenomes combining wild-type and mutant chromosomes, we have deciphered the sequence of physiological and molecular events underlying the de novo invasion, proliferation and eventual demise of the single-copy endogenous retrotransposon Evadé (EVD). We show how this reconstructed TE burst causes widespread genome diversification and de novo epiallelism that could serve as sources for selectable and potentially adaptive traits.

title: Silencing of active transposable elements in plants.
abstract: In plant genomes the vast majority of transposable elements (TEs) are found in a transcriptionally silenced state that is epigenetically propagated from generation to generation. Although the mechanism of this maintenance of silencing has been well studied, it is now clear that the pathways responsible for maintaining TEs in a silenced state differ from the pathways responsible for initially targeting the TE for silencing. Recently, attention in this field has focused on investigating the molecular mechanisms that initiate and establish TE silencing. Here we review the current models of how TEs are triggered for silencing, the data supporting each model, and the key future questions in this fast moving field.

title: RNA-Directed DNA Methylation: The Evolution of a Complex Epigenetic Pathway in Flowering Plants.
abstract: RNA-directed DNA methylation (RdDM) is an epigenetic process in plants that involves both short and long noncoding RNAs. The generation of these RNAs and the induction of RdDM rely on complex transcriptional machineries comprising two plant-specific, RNA polymerase II (Pol II)-related RNA polymerases known as Pol IV and Pol V, as well as a host of auxiliary factors that include both novel and refashioned proteins. We present current views on the mechanism of RdDM with a focus on evolutionary innovations that occurred during the transition from a Pol II transcriptional pathway, which produces mRNA precursors and numerous noncoding RNAs, to the Pol IV and Pol V pathways, which are specialized for RdDM and gene silencing. We describe recently recognized deviations from the canonical RdDM pathway, discuss unresolved issues, and speculate on the biological significance of RdDM for flowering plants, which have a highly developed Pol V pathway.

title: Heritable transposon silencing initiated by a naturally occurring transposon inverted duplication
abstract: It has been suggested that gene silencing evolved as a defense against genomic parasites such as transposons. This idea is based on analysis of mutations that reactivate transposons that are stably silenced: they affect maintenance rather than initiation of silencing. Here we describe the cloning and characterization of a naturally occurring locus able to heritably silence the otherwise highly active MuDR transposon in maize. This locus, Mu killer (Muk), results from the inverted duplication of a partially deleted autonomous MuDR element located at the breakpoint of a genomic deletion. Muk produces a hybrid hairpin transcript that is processed into small RNAs, which are amplified when the target MuDR transcript is present. Muk provides the first example of a naturally occurring transposon derivative capable of initiating the heritable silencing of an active transposon family. Further, transposon-generated inverted duplications may be important for the generation of double-stranded RNAs used in gene silencing.

title: Effects of recombination rate and gene density on transposable element distributions in Arabidopsis thaliana.
abstract: Transposable elements (TEs) comprise a major component of eukaryotic genomes, and exhibit striking deviations from random distribution across the genomes studied, including humans, flies, nematodes, and plants. Although considerable progress has been made in documenting these patterns, the causes are subject to debate. Here, we use the genome sequence of Arabidopsis thaliana to test for the importance of competing models of natural selection against TE insertions. We show that, despite TE accumulation near the centromeres, recombination does not generally correlate with TE abundance, suggesting that selection against ectopic recombination does not influence TE distribution in A. thaliana. In contrast, a consistent negative correlation between gene density and TE abundance, and a strong under-representation of TE insertions in introns suggest that selection against TE disruption of gene expression is playing a more important role in A. thaliana. High rates of self-fertilization may reduce the importance of recombination rate in genome structuring in inbreeding organisms such as A. thaliana and Caenorhabditis elegans.

title: Demography and weak selection drive patterns of transposable element diversity in natural populations of Arabidopsis lyrata
abstract: Transposable elements (TEs) are the major component of most plant genomes, and characterizing their population dynamics is key to understanding plant genome complexity. Yet there have been few studies of TE population genetics in plant systems. To study the roles of selection, transposition, and demography in shaping TE population diversity, we generated a polymorphism dataset for six TE families in four populations of the flowering plant Arabidopsis lyrata. The TE data indicated significant differentiation among populations, and maximum likelihood procedures suggested weak selection. For strongly bottlenecked populations, the observed TE band-frequency spectra fit data simulated under neutral demographic models constructed from nucleotide polymorphism data. Overall, we propose that TEs are subjected to weak selection, the efficacy of which varies as a function of demographic factors. Thus, demographic effects could be a major factor driving distributions of TEs among plant lineages.

title: The evolution of transposable elements in natural populations of self-fertilizing Arabidopsis thaliana and its outcrossing relative Arabidopsis lyrata
abstract: BackgroundTransposable Elements (TEs) make up the majority of plant genomes, and thus understanding TE evolutionary dynamics is key to understanding plant genome evolution. Plant reproductive systems are diverse and mating type variation is one factor among many hypothesized to influence TE evolutionary dynamics. Here, we collected a large TE-display data set in self-fertilizing Arabidopsis thaliana, and compared it to data gathered in outcrossing Arabidopsis lyrata. We analyzed seven TE families in four natural populations of each species to tease apart the effects of mating system, demography, transposition, and selection in determining patterns of TE diversity.ResultsMeasures of TE band differentiation were largely consistent across TE families. However, patterns of diversity in A. thaliana Ac elements differed significantly from that other TEs, perhaps signaling a lack of recent transposition. Across TE families, we estimated higher allele frequencies and lower selection coefficients on A. thaliana TE insertions relative to A. lyrata TE insertions.ConclusionsThe differences in TE distributions between the two Arabidopsis species represents a synthesis of evolutionary forces that include the transposition dynamics of individual TE families and the demographic histories of populations. There are also species-specific differences that could be attributed to the effects of mating system, including higher overall allele frequencies in the selfing lineage and a greater proportion of among population TE diversity in the outcrossing lineage.

title: An Age-of-Allele Test of Neutrality for Transposable Element Insertions
abstract: How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, previous studies have used models of transposition–selection equilibrium that assume a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence. By conditioning on the age of an individual TE insertion allele (inferred by the number of unique substitutions that have occurred within the particular TE sequence since insertion), we determine the probability distribution of the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this nonequilibrium neutral model, we are able to explain ∼80% of the variance in TE insertion allele frequencies based on age alone. Controlling for both nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications, or other copy number variants.

title: Ancestral repeats have shaped epigenome and genome composition for millions of years in Arabidopsis thaliana
abstract: Little is known about the evolution of repeated sequences over long periods of time. Using two independent approaches, we show that the majority of the repeats found in the Arabidopsis thaliana genome are ancient and likely to derive from the retention of fragments deposited during ancestral bursts that occurred early in the Brassicaceae evolution. We determine that the majority of young repeats are found in pericentromeric domains, while older copies are frequent in the gene-rich regions. Our results further suggest that the DNA methylation of repeats through small RNA-mediated pathways can last over prolonged periods of time. We also illustrate the way repeated sequences are composted by mutations towards genomic dark matter over time, probably driven by the deamination of methylcytosines, which also have an impact on epigenomic landscapes. Overall, we show that the ancient proliferation of repeat families has long-term consequences on A. thaliana biology and genome composition.

title: Recent speciation associated with the evolution of selfing in Capsella
abstract: The evolution from outcrossing to predominant self-fertilization represents one of the most common transitions in flowering plant evolution. This shift in mating system is almost universally associated with the “selfing syndrome,” characterized by marked reduction in flower size and a breakdown of the morphological and genetic mechanisms that prevent self-fertilization. In general, the timescale in which these transitions occur, and the evolutionary dynamics associated with the evolution of the selfing syndrome are poorly known. We investigated the origin and evolution of selfing in the annual plant Capsella rubella from its self-incompatible, outcrossing progenitor Capsella grandiflora by characterizing multilocus patterns of DNA sequence variation at nuclear genes. We estimate that the transition to selfing and subsequent geographic expansion have taken place during the past 20,000 years. This transition was probably associated with a shift from stable equilibrium toward a near-complete population bottleneck causing a major reduction in effective population size. The timing and severe founder event support the hypothesis that selfing was favored during colonization as new habitats emerged after the last glaciation and the expansion of agriculture. These results suggest that natural selection for reproductive assurance can lead to major morphological evolution and speciation on relatively short evolutionary timescales.

title: Genome-wide evidence for efficient positive and purifying selection in Capsella grandiflora, a plant species with a large effective population size.
abstract: Recent studies comparing genome-wide polymorphism and divergence in Drosophila have found evidence for a surprisingly high proportion of adaptive amino acid fixations, but results for other taxa are mixed. In particular, few studies have found convincing evidence for adaptive amino acid substitution in plants. To assess the generality of this finding, we have sequenced 257 loci in the outcrossing crucifer Capsella grandiflora, which has a large effective population size and low population structure. Using a new method that jointly infers selective and demographic effects, we estimate that 40% of amino acid substitutions were fixed by positive selection in this species, and we also infer a low proportion of slightly deleterious amino acid mutations. We contrast these estimates with those for a similar data set from the closely related Arabidopsis thaliana and find significantly higher rates of adaptive evolution and fewer nearly neutral mutations in C. grandiflora. In agreement with results for other taxa, genes involved in reproduction show the strongest evidence for positive selection in C. grandiflora. Taken together, these results imply that both positive and purifying selection are more effective in C. grandiflora than in A. thaliana, consistent with the contrasting demographic history and effective population sizes of these species.

title: Hybrid origins and the earliest stages of diploidization in the highly successful recent polyploid Capsella bursa-pastoris
abstract: Whole genome duplication events have occurred repeatedly during flowering plant evolution, and there is growing evidence for predictable patterns of gene retention and loss following polyploidization. Despite these important insights, the rate and processes governing the earliest stages of diploidization remain poorly understood, and the relative importance of genetic drift, positive selection and relaxed purifying selection in the process of gene degeneration and loss is unclear. Here, we conduct whole genome resequencing in Capsella bursa-pastoris, a recently formed tetraploid with one of the most widespread species distributions of any angiosperm. Whole genome data provide strong support for recent hybrid origins of the tetraploid species within the last 100-300,000 years from two diploid progenitors in the Capsella genus. Major-effect inactivating mutations are frequent, but many were inherited from the parental species and show no evidence of being fixed by positive selection. Despite a lack of large-scale gene loss, we observe a decrease in the efficacy of natural selection genome-wide, due to the combined effects of demography, selfing and genome redundancy from whole genome duplication. Our results suggest that the earliest stages of diploidization are associated with quantitative genome-wide decreases in the strength and efficacy of selection rather than rapid gene loss, and that non-functionalization can receive a 'head start' through a legacy of deleterious variants and differential expression originating in parental diploid populations.

title: Evidence for widespread positive and negative selection in coding and conserved noncoding regions of Capsella grandiflora
abstract: The extent that both positive and negative selection vary across different portions of plant genomes remains poorly understood. Here, we sequence whole genomes of 13 Capsella grandiflora individuals and quantify the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we show that selection is strong in coding regions, but weak in most noncoding regions, with the exception of 5’ and 3’ untranslated regions (UTRs). However, estimates of selection in noncoding regions conserved across the Brassicaceae family show strong signals of selection. Additionally, we see reductions in neutral diversity around functional substitutions in both coding and conserved noncoding regions, indicating recent selective sweeps at these sites. Finally, using expression data from leaf tissue we show that genes that are more highly expressed experience stronger negative selection but comparable levels of positive selection to lowly expressed genes. Overall, we observe widespread positive and negative selection in coding and regulatory regions, but our results also suggest that both positive and negative selection in plant noncoding sequence are considerably rarer than in animal genomes.  Selection affects patterns of genomic variation, but it is unclear how much the genomic effects of selection vary across plant genomes, particularly in noncoding regions. To determine the strength and extent of selective signatures across the genome, we sequenced and analyzed genomes from 13 Capsella grandiflora individuals. C. grandiflora is well-suited for this task because it has experienced a large, stable effective population size, so we expect that selection signatures will not be overly distorted by demographic effects. Our analysis shows that positive and negative selection acting on new mutations have broadly shaped patterns of genomic diversity in coding regions but not in most noncoding regions. However, when we focus only on noncoding regions that show evidence of constraint across species, we see evidence for strong positive and negative selection. In addition, we find that genes with high expression experience stronger negative selection than genes with low expression, but the extent of positive selection appears to be equivalent across expression categories.

title: Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora
abstract: Despite long-standing interest in the contribution of cis-regulatory changes to adaptation, we still have a limited understanding of the selective importance and genomic determinants of cis-regulatory variation in natural populations. We use a combination of analyses of allele-specific expression and population genomic analyses to investigate the selective forces and genomic determinants of cis-regulatory variation in the outcrossing plant species Capsella grandiflora. We conclude that gene-specific functional constraints shape cis-regulatory variation and that genes with cis-regulatory variation are under relaxed purifying selection compared with other genes. Finally, we identify a link between gene body methylation and the extent of cis-regulatory constraint in natural populations. Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora. We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

title: Association mapping reveals the role of purifying selection in the maintenance of genomic variation in gene expression
abstract: Biologists have long sought to explain why we see genetic variation for traits in populations despite the expectation that selection will remove most variation. We address this question by using gene expression as a model trait and identifying the genetic loci that affect gene expression in a single, large population of the plant Capsella grandiflora. Alleles at loci that affect expression were rarer than expected under neutral expectations, and there was a negative correlation between phenotypic effect size and frequency of these alleles. These observations are consistent with the hypothesis that purifying selection acts on the genetic variation for expression. The evolutionary forces that maintain genetic variation in quantitative traits within populations remain poorly understood. One hypothesis suggests that variation is under purifying selection, resulting in an excess of low-frequency variants and a negative correlation between minor allele frequency and selection coefficients. Here, we test these predictions using the genetic loci associated with total expression variation (eQTLs) and allele-specific expression variation (aseQTLs) mapped within a single population of the plant Capsella grandiflora. In addition to finding eQTLs and aseQTLs for a large fraction of genes, we show that alleles at these loci are rarer than expected and exhibit a negative correlation between phenotypic effect size and frequency. Overall, our results show that the distribution of frequencies and effect sizes of the loci responsible for local expression variation within a single outcrossing population are consistent with the effects of purifying selection.

title: Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella
abstract: The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.

title: Trimmomatic: a flexible trimmer for Illumina sequence data
abstract: Although many next-generation sequencing (NGS) read preprocessing tools already existed, we could not find any tool or combination of tools that met our requirements in terms of flexibility, correct handling of paired-end data and high performance. We have developed Trimmomatic as a more flexible and efficient preprocessing tool, which could correctly handle paired-end data.  The value of NGS read preprocessing is demonstrated for both reference-based and reference-free tasks. Trimmomatic is shown to produce output that is at least competitive with, and in many cases superior to, that produced by other tools, in all scenarios tested. Availability and implementation: Trimmomatic is licensed under GPL V3. It is cross-platform (Java 1.5+ required) and available at http://www.usadellab.org/cms/index.php?page=trimmomatic Contact: usadel@bio1.rwth-aachen.de Supplementary information: Supplementary data are available at Bioinformatics online.

title: Fast and accurate short read alignment with Burrows–Wheeler transform
abstract: The enormous amount of short reads generated by the new DNA sequencing technologies call for the development of fast and accurate read alignment programs. A first generation of hash table-based methods has been developed, including MAQ, which is accurate, feature rich and fast enough to align short reads from a single individual. However, MAQ does not support gapped alignment for single-end reads, which makes it unsuitable for alignment of longer reads where indels may occur frequently. The speed of MAQ is also a concern when the alignment is scaled up to the resequencing of hundreds of individuals.  We implemented Burrows-Wheeler Alignment tool (BWA), a new read alignment package that is based on backward search with Burrows–Wheeler Transform (BWT), to efficiently align short sequencing reads against a large reference sequence such as the human genome, allowing mismatches and gaps. BWA supports both base space reads, e.g. from Illumina sequencing machines, and color space reads from AB SOLiD machines. Evaluations on both simulated and real data suggest that BWA is ∼10–20× faster than MAQ, while achieving similar accuracy. In addition, BWA outputs alignment in the new standard SAM (Sequence Alignment/Map) format. Variant calling and other downstream analyses after the alignment can be achieved with the open source SAMtools software package.

title: PoPoolationTE2: comparative population genomics of transposable elements using Pool-Seq
abstract: The evolutionary dynamics of transposable elements (TEs) are still poorly understood. One reason is that TE abundance needs to be studied at the population level, and despite recent advances in sequencing technologies, characterizing TE abundance in multiple populations by sequencing individuals separately is still too expensive. While sequencing pools of individuals (Pool-Seq) dramatically reduces sequencing costs, a comparison of TE abundance between pooled samples has been difficult, if not impossible, due to various biases. Here, we introduce a novel bioinformatic tool, PoPoolationTE2, which is specifically tailored for the comparison of TE abundance among pooled population samples or different tissues. Using computer simulations we demonstrate that PoPoolationTE2 not only faithfully recovers TE insertion frequencies and positions but, by homogenizing the power to identify TEs acrosss samples, it provides an unbiased comparison of TE abundance between pooled population samples. We anticipate that PoPoolationTE2 will greatly facilitate the analysis of TE insertion patterns in a broad range of applications.

title: The Capsella rubella genome and the genomic consequences of rapid mating system evolution
abstract: The shift from outcrossing to selfing is common in flowering plants, but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible <200,000 years ago. We report a C. rubella reference genome sequence and compare RNA expression and polymorphism patterns between C. rubella and its outcrossing progenitor Capsella grandiflora. We found a clear shift in the expression of genes associated with flowering phenotypes, similar to that seen in Arabidopsis, in which self fertilization evolved about 1 million years ago. Comparisons of the two Capsella species showed evidence of rapid genome-wide relaxation of purifying selection in C. rubella without a concomitant change in transposable element abundance. Overall we document that the transition to selfing may be typified by parallel shifts in gene expression, along with a measurable reduction of purifying selection.

title: A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data
abstract: Most existing methods for DNA sequence analysis rely on accurate sequences or genotypes. However, in applications of the next-generation sequencing (NGS), accurate genotypes may not be easily obtained (e.g. multi-sample low-coverage sequencing or somatic mutation discovery). These applications press for the development of new methods for analyzing sequence data with uncertainty.



We present a statistical framework for calling SNPs, discovering somatic mutations, inferring population genetical parameters and performing association tests directly based on sequencing data without explicit genotyping or linkage-based imputation. On real data, we demonstrate that our method achieves comparable accuracy to alternative methods for estimating site allele count, for inferring allele frequency spectrum and for association mapping. We also highlight the necessity of using symmetric datasets for finding somatic mutations and confirm that for discovering rare events, mismapping is frequently the leading source of errors.


AVAILABILITY
http://samtools.sourceforge.net.



hengli@broadinstitute.org.

title: Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM
abstract: BWA-MEM is a new alignment algorithm for aligning sequence reads or long query sequences against a large reference genome such as human. It automatically chooses between local and end-to-end alignments, supports paired-end reads and performs chimeric alignment. The algorithm is robust to sequencing errors and applicable to a wide range of sequence lengths from 70bp to a few megabases. For mapping 100bp sequences, BWA-MEM shows better performance than several state-of-art read aligners to date. 
Availability and implementation: BWA-MEM is implemented as a component of BWA, which is available at this http URL. 
Contact: hengli@broadinstitute.org

title: The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data.
abstract: Next-generation DNA sequencing (NGS) projects, such as the 1000 Genomes Project, are already revolutionizing our understanding of genetic variation among individuals. However, the massive data sets generated by NGS--the 1000 Genome pilot alone includes nearly five terabases--make writing feature-rich, efficient, and robust analysis tools difficult for even computationally sophisticated individuals. Indeed, many professionals are limited in the scope and the ease with which they can answer scientific questions by the complexity of accessing and manipulating the data produced by these machines. Here, we discuss our Genome Analysis Toolkit (GATK), a structured programming framework designed to ease the development of efficient and robust analysis tools for next-generation DNA sequencers using the functional programming philosophy of MapReduce. The GATK provides a small but rich set of data access patterns that encompass the majority of analysis tool needs. Separating specific analysis calculations from common data management infrastructure enables us to optimize the GATK framework for correctness, stability, and CPU and memory efficiency and to enable distributed and shared memory parallelization. We highlight the capabilities of the GATK by describing the implementation and application of robust, scale-tolerant tools like coverage calculators and single nucleotide polymorphism (SNP) calling. We conclude that the GATK programming framework enables developers and analysts to quickly and easily write efficient and robust NGS tools, many of which have already been incorporated into large-scale sequencing projects like the 1000 Genomes Project and The Cancer Genome Atlas.

title: From FastQ Data to High‐Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline
abstract: This unit describes how to use BWA and the Genome Analysis Toolkit (GATK) to map genome sequencing data to a reference and produce high‐quality variant calls that can be used in downstream analyses. The complete workflow includes the core NGS data‐processing steps that are necessary to make the raw data suitable for analysis by the GATK, as well as the key methods involved in variant discovery using the GATK.

title: The variant call format and VCFtools
abstract: The variant call format (VCF) is a generic format for storing DNA polymorphism data such as SNPs, insertions, deletions and structural variants, together with rich annotations. VCF is usually stored in a compressed manner and can be indexed for fast data retrieval of variants from a range of positions on the reference genome. The format was developed for the 1000 Genomes Project, and has also been adopted by other projects such as UK10K, dbSNP and the NHLBI Exome Project. VCFtools is a software suite that implements various utilities for processing VCF files, including validation, merging, comparing and also provides a general Perl API.

title: Demography and mating system shape the genome-wide impact of purifying selection in Arabis alpina
abstract: Plant mating systems have profound effects on levels and structuring of genetic variation, and can affect the impact of natural selection. While theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine crucifer Arabis alpina. We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing Scandinavian A. alpina show a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, both mating system and demography shape the impact of purifying selection on genomic variation in A. alpina. These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.  Intermediate outcrossing rates are theoretically predicted to maintain effective selection against harmful alleles, but few studies have empirically tested this prediction using genomic data. We used whole-genome resequencing data from alpine rock-cress to study how genetic variation and purifying selection vary with mating system. We find that populations with intermediate outcrossing rates have similar levels of genetic diversity as outcrossing populations, and that purifying selection against harmful alleles is efficient in mixed-mating populations. In contrast, self-fertilizing populations from Scandinavia have strongly reduced genetic diversity, and accumulate harmful mutations, likely as a result of demographic effects of postglacial colonization. Our results suggest that mixed-mating populations can avoid the negative evolutionary consequences of high self-fertilization rates.

title: On the number of segregating sites in genetical models without recombination.
abstract: The distribution is obtained for the number of segregating sites observed in a sample from a population which is subject to recurring, new, mutations but not subject to recombination. After allowance is made for the different effective population sizes, the results apply approximately to three population models, due to Wright, Burrows and Cockerham, and Moran. Included as extreme special cases are the distributions of the number of segregating sites in the whole population and of the number of heterozygous sites in a diploid individual. Some results of Fisher, Haldane, Kimura, and Ewens concerning the means of the distributions for different models are confirmed, but the variances, and the distributions themselves, are new.

title: Recent speciation of Capsella rubella from Capsella grandiflora, associated with loss of self-incompatibility and an extreme bottleneck
abstract: Flowering plants often prevent selfing through mechanisms of self-incompatibility (S.I.). The loss of S.I. has occurred many times independently, because it provides short-term advantages in situations where pollinators or mates are rare. The genus Capsella, which is closely related to Arabidopsis, contains a pair of closely related diploid species, the self-incompatible Capsella grandiflora and the self-compatible Capsella rubella. To elucidate the transition to selfing and its relationship to speciation of C. rubella, we have made use of comparative sequence information. Our analyses indicate that C. rubella separated from C. grandiflora recently (≈30,000–50,000 years ago) and that breakdown of S.I. occurred at approximately the same time. Contrasting the nucleotide diversity patterns of the 2 species, we found that C. rubella has only 1 or 2 alleles at most loci, suggesting that it originated through an extreme population bottleneck. Our data are consistent with diploid speciation by a single, selfing individual, most likely living in Greece. The new species subsequently colonized the Mediterranean by Northern and Southern routes, at a time that also saw the spread of agriculture. The presence of phenotypic diversity within modern C. rubella suggests that this species will be an interesting model to understand divergence and adaptation, starting from very limited standing genetic variation.

title: Rapid divergence and high diversity of miRNAs and miRNA targets in the Camelineae.
abstract: MicroRNAs (miRNAs) are short RNAs involved in gene regulation through translational inhibition and transcript cleavage. After processing from imperfect fold-back structures, miRNAs are incorporated into RNA-induced silencing complexes (RISCs) before targeting transcripts with varying degrees of complementarity. Some miRNAs are evolutionarily deep-rooted, and sequence complementarity with their targets is maintained through purifying selection. Both Arabidopsis and Capsella belong to the tribe Camelineae in the Brassicaceae, with Capsella rubella serving as an outgroup to the genus Arabidopsis. The genome sequence of C. rubella has recently been released, which allows characterization of its miRNA complement in comparison with Arabidopsis thaliana and Arabidopsis lyrata. Through next-generation sequencing, we identify high-confidence miRNA candidates specific to the C. rubella lineage. Only a few lineage-specific miRNAs have been studied for evolutionary constraints, and there have been no systematic studies of miRNA target diversity within or divergence between closely related plant species. Therefore we contrast sequence variation in miRNAs and their targets within A. thaliana, and between A. thaliana, A. lyrata and C. rubella. We document a surprising amount of small-scale variation in miRNA-target pairs, where many miRNAs are predicted to have species-specific targets in addition to ones that are shared between species. Our results emphasize that the transitive nature of many miRNA-target pairs can be observed even on a relatively short evolutionary time-scale, with non-random occurrences of differences in miRNAs and their complements in the miRNA precursors, the miRNA* sequences.

title: STAR: ultrafast universal RNA-seq aligner
abstract: Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases.



To align our large (>80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of >50 in mapping speed, aligning to the human genome 550 million 2 × 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80-90% success rate, corroborating the high precision of the STAR mapping strategy.



STAR is implemented as a standalone C++ code.

title: The Rate and Molecular Spectrum of Spontaneous Mutations in Arabidopsis thaliana
abstract: Evolution in Action Rates of evolution in gene and genome sequences have been estimated, but these estimates are subject to error because many of the steps of evolution over the ages are not directly measurable or are hidden under subsequent changes. Ossowski et al. (p. 92) now provide a more accurate measurement of how often spontaneous mutations arise in a nuclear genome. Mutations arising over 30 generations were compared by sequencing DNA from individual Arabidopsis thaliana plants. UV- and deamination-induced mutagenesis appeared to bias the type of mutations found. Rapid sequencing technologies allow a more accurate calculation of the mutation rate for plants. To take complete advantage of information on within-species polymorphism and divergence from close relatives, one needs to know the rate and the molecular spectrum of spontaneous mutations. To this end, we have searched for de novo spontaneous mutations in the complete nuclear genomes of five Arabidopsis thaliana mutation accumulation lines that had been maintained by single-seed descent for 30 generations. We identified and validated 99 base substitutions and 17 small and large insertions and deletions. Our results imply a spontaneous mutation rate of 7 × 10−9 base substitutions per site per generation, the majority of which are G:C→A:T transitions. We explain this very biased spectrum of base substitution mutations as a result of two main processes: deamination of methylated cytosines and ultraviolet light–induced mutagenesis.

title: The First Rule of Plant Transposable Element Silencing: Location, Location, Location
abstract: Transposable elements (TEs) are mobile units of DNA that comprise large portions of plant genomes. Besides creating mutations via transposition and contributing to genome size, TEs play key roles in chromosome architecture and gene regulation. TE activity is repressed by overlapping mechanisms of chromatin condensation, epigenetic transcriptional silencing, and targeting by small interfering RNAs. The specific regulation of different TEs, as well as their different roles in chromosome architecture and gene regulation, is specified by where on the chromosome the TE is located: near a gene, within a gene, in a pericentromere/TE island, or at the centromere core. In this Review, we investigate the silencing mechanisms responsible for inhibiting TE activity for each of these chromosomal contexts, emphasizing that chromosomal location is the first rule dictating the specific regulation of each TE.