The Arabidopsis lyrata genome sequence and the basis of rapid genome size change

On the Coevolution of Transposable Elements and Plant Genomes

Plant genomes are unique in an intriguing feature: the range of their size variation is unprecedented among living organisms. Although polyploidization contributes to this variability, transposable elements (TEs) seem to play the pivotal role. TEs, often considered intragenomic parasites, not only affect the genome size of the host, but also interact with other genes, disrupting and creating new functions and regulatory networks. Coevolution of plant genomes and TEs has led to tight regulation of TE activity, and growing evidence suggests their relationship became mutualistic. Although the expansions of TEs represent certain costs for the host genomes, they may also bring profits for populations, helping to overcome challenging environmental (biotic/abiotic stress) or genomic (hybridization and allopolyploidization) conditions. In this paper, we discuss the possibility that the possession of inducible TEs may provide a selective advantage for various plant populations.

Phytozome: a comparative platform for green plant genomics

The number of sequenced plant genomes and associated genomic resources is growing rapidly with the advent of both an increased focus on plant genomics from funding agencies, and the application of inexpensive next generation sequencing. To interact with this increasing body of data, we have developed Phytozome (http://www.phytozome.net), a comparative hub for plant genome and gene family data and analysis. Phytozome provides a view of the evolutionary history of every plant gene at the level of sequence, gene structure, gene family and genome organization, while at the same time providing access to the sequences and functional annotations of a growing number (currently 25) of complete plant genomes, including all the land plants and selected algae sequenced at the Joint Genome Institute, as well as selected species sequenced elsewhere. Through a comprehensive plant genome database and web portal, these data and analyses are available to the broader plant science research community, providing powerful comparative genomics tools that help to link model systems with other plants of economic and ecological importance.

i-ADHoRe 3.0--fast and sensitive detection of genomic homology in extremely large data sets

Comparative genomics is a powerful means to gain insight into the evolutionary processes that shape the genomes of related species. As the number of sequenced genomes increases, the development of software to perform accurate cross-species analyses becomes indispensable. However, many implementations that have the ability to compare multiple genomes exhibit unfavorable computational and memory requirements, limiting the number of genomes that can be analyzed in one run. Here, we present a software package to unveil genomic homology based on the identification of conservation of gene content and gene order (collinearity), i-ADHoRe 3.0, and its application to eukaryotic genomes. The use of efficient algorithms and support for parallel computing enable the analysis of large-scale data sets. Unlike other tools, i-ADHoRe can process the Ensembl data set, containing 49 species, in 1 h. Furthermore, the profile search is more sensitive to detect degenerate genomic homology than chaining pairwise collinearity information based on transitive homology. From ultra-conserved collinear regions between mammals and birds, by integrating coexpression information and protein–protein interactions, we identified more than 400 regions in the human genome showing significant functional coherence. The different algorithmical improvements ensure that i-ADHoRe 3.0 will remain a powerful tool to study genome evolution.

Widespread interspecific divergence in cis-regulation of transposable elements in the Arabidopsis genus

Transposable elements (TEs) are so abundant and variable that they count among the most important mutational sources in genomes. Nonetheless, little is known about the genetics of their variation in activity or silencing across closely related species. Here, we demonstrate that regulation of TE genes can differ dramatically between the two closely related Arabidopsis species A. thaliana and A. lyrata. In leaf and floral tissues of F1 interspecific hybrids, about 47% of TEs show allele-specific expression, with the A. lyrata copy being generally expressed at higher level. We confirm that TEs are generally expressed in A. lyrata but not in A. thaliana. Allele-specific differences in TE expression are associated with divergence in epigenetic modifications like DNA and histone methylation between species as well as with sequence divergence. Our data demonstrate that A. thaliana silences TEs much better than A. lyrata. For long terminal repeat retrotransposons, these differences are more pronounced for younger insertions. Interspecific differences in TE silencing may have a great impact on genome size changes.

Boechera, a model system for ecological genomics

The selection and development of a study system for evolutionary and ecological functional genomics (EEFG) depend on a variety of factors. Here, we present the genus Boechera as an exemplary system with which to address ecological and evolutionary questions. Our focus on Boechera is based on several characteristics as follows: (i) native populations in undisturbed habitats where current environments reflect historical conditions over several thousand years; (ii) functional genomics benefitting from its close relationship to Arabidopsis thaliana; (iii) inbreeding tolerance enabling development of recombinant inbred lines, near-isogenic lines and positional cloning; (iv) interspecific crosses permitting mapping for genetic analysis of speciation; (v) apomixis (asexual reproduction by seeds) in a genetically tractable diploid; and (vi) broad geographic distribution in North America, permitting ecological genetics for a large research community. These characteristics, along with the current sequencing of three Boechera species by the Joint Genome Institute, position Boechera as a rapidly advancing system for EEFG studies.

Higher intron loss rate in Arabidopsis thaliana than A. lyrata is consistent with stronger selection for a smaller genome

The number of introns varies considerably among different organisms. This can be explained by the differences in the rates of intron gain and loss. Two factors that are likely to influence these rates are selection for or against introns and the mutation rate that generates the novel intron or the intronless copy. Although it has been speculated that stronger selection for a compact genome might result in a higher rate of intron loss and a lower rate of intron gain, clear evidence is lacking, and the role of selection in determining these rates has not been established. Here, we studied the gain and loss of introns in the two closely related species Arabidopsis thaliana and A. lyrata as it was recently shown that A. thaliana has been undergoing a faster genome reduction driven by selection. We found that A. thaliana has lost six times more introns than A. lyrata since the divergence of the two species but gained very few introns. We suggest that stronger selection for genome reduction probably resulted in the much higher intron loss rate in A. thaliana, although further analysis is required as we could not find evidence that the loss rate increased in A. thaliana as opposed to having decreased in A. lyrata compared with the rate in the common ancestor. We also examined the pattern of the intron gains and losses to better understand the mechanisms by which they occur. Microsimilarity was detected between the splice sites of several gained and lost introns, suggesting that nonhomologous end joining repair of double-strand breaks might be a common pathway not only for intron gain but also for intron loss.

Developmental genetics and new sequencing technologies: the rise of nonmodel organisms

Much of developmental biology in the past decades has been driven by forward genetic studies in a few model organisms. We review recent work with relatives of these species, motivated by a desire to understand the evolutionary and ecological context for morphological innovation. Unfortunately, despite a number of shining examples, progress in nonmodel systems has often been slow. The current revolution in DNA sequencing has, however, enormous potential in extending the reach of genetics. We discuss how developmental biology will benefit from these advances, particularly by increasing the universe of study species.

Genome-wide comparison of nucleotide-binding site-leucine-rich repeat-encoding genes in Arabidopsis

Plants, like animals, use several lines of defense against pathogen attack. Prominent among genes that confer disease resistance are those encoding nucleotide-binding site-leucine-rich repeat (NB-LRR) proteins. Likely due to selection pressures caused by pathogens, NB-LRR genes are the most variable gene family in plants, but there appear to be species-specific limits to the number of NB-LRR genes in a genome. Allelic diversity within an individual is also increased by obligatory outcrossing, which leads to genome-wide heterozygosity. In this study, we compared the NB-LRR gene complement of the selfer Arabidopsis thaliana and its outcrossing close relative Arabidopsis lyrata. We then complemented and contrasted the interspecific patterns with studies of NB-LRR diversity within A. thaliana. Three important insights are as follows: (1) that both species have similar numbers of NB-LRR genes; (2) that loci with single NB-LRR genes are less variable than tandem arrays; and (3) that presence-absence polymorphisms within A. thaliana are not strongly correlated with the presence or absence of orthologs in A. lyrata. Although A. thaliana individuals are mostly homozygous and thus potentially less likely to suffer from aberrant interaction of NB-LRR proteins with newly introduced alleles, the number of NB-LRR genes is similar to that in A. lyrata. In intraspecific and interspecific comparisons, NB-LRR genes are also more variable than receptor-like protein genes. Finally, in contrast to Drosophila, there is a clearly positive relationship between interspecific divergence and intraspecific polymorphisms.

Evolution of the S-locus region in Arabidopsis relatives

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high-quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haplotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self-compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora.

Genomic determinants of protein evolution and polymorphism in Arabidopsis

Recent results from Drosophila suggest that positive selection has a substantial impact on genomic patterns of polymorphism and divergence. However, species with smaller population sizes and/or stronger population structure may not be expected to exhibit Drosophila-like patterns of sequence variation. We test this prediction and identify determinants of levels of polymorphism and rates of protein evolution using genomic data from Arabidopsis thaliana and the recently sequenced Arabidopsis lyrata genome. We find that, in contrast to Drosophila, there is no negative relationship between nonsynonymous divergence and silent polymorphism at any spatial scale examined. Instead, synonymous divergence is a major predictor of silent polymorphism, which suggests variation in mutation rate as the main determinant of silent variation. Variation in rates of protein divergence is mainly correlated with gene expression level and breadth, consistent with results for a broad range of taxa, and map-based estimates of recombination rate are only weakly correlated with nonsynonymous divergence. Variation in mutation rates and the strength of purifying selection seem to be major drivers of patterns of polymorphism and divergence in Arabidopsis. Nevertheless, a model allowing for varying negative and positive selection by functional gene category explains the data better than a homogeneous model, implying the action of positive selection on a subset of genes. Genes involved in disease resistance and abiotic stress display high proportions of adaptive substitution. Our results are important for a general understanding of the determinants of rates of protein evolution and the impact of selection on patterns of polymorphism and divergence.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana

Genetic differences between Arabidopsis thaliana accessions underlie the plant's extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

Whole-genome sequencing of multiple Arabidopsis thaliana populations

The plant Arabidopsis thaliana occurs naturally in many different habitats throughout Eurasia. As a foundation for identifying genetic variation contributing to adaptation to diverse environments, a 1001 Genomes Project to sequence geographically diverse A. thaliana strains has been initiated. Here we present the first phase of this project, based on population-scale sequencing of 80 strains drawn from eight regions throughout the species' native range. We describe the majority of common small-scale polymorphisms as well as many larger insertions and deletions in the A. thaliana pan-genome, their effects on gene function, and the patterns of local and global linkage among these variants. The action of processes other than spontaneous mutation is identified by comparing the spectrum of mutations that have accumulated since A. thaliana diverged from its closest relative 10 million years ago with the spectrum observed in the laboratory. Recent species-wide selective sweeps are rare, and potentially deleterious mutations are more common in marginal populations.

Complex evolutionary events at a tandem cluster of Arabidopsis thaliana genes resulting in a single-locus genetic incompatibility

Non-additive interactions between genomes have important implications, not only for practical applications such as breeding, but also for understanding evolution. In extreme cases, genes from different genomic backgrounds may be incompatible and compromise normal development or physiology. Of particular interest are non-additive interactions of alleles at the same locus. For example, overdominant behavior of alleles, with respect to plant fitness, has been proposed as an important component of hybrid vigor, while underdominance may lead to reproductive isolation. Despite their importance, only a few cases of genetic over- or underdominance affecting plant growth or fitness are understood at the level of individual genes. Moreover, the relationship between biochemical and fitness effects may be complex: genetic overdominance, that is, increased or novel activity of a gene may lead to evolutionary underdominance expressed as hybrid weakness. Here, we describe a non-additive interaction between alleles at the Arabidopsis thaliana OAK (OUTGROWTH-ASSOCIATED PROTEIN KINASE) gene. OAK alleles from two different accessions interact in F₁ hybrids to cause a variety of aberrant growth phenotypes that depend on a recently acquired promoter with a novel expression pattern. The OAK gene, which is located in a highly variable tandem array encoding closely related receptor-like kinases, is found in one third of A. thaliana accessions, but not in the reference accession Col-0. Besides recruitment of exons from nearby genes as promoter sequences, key events in OAK evolution include gene duplication and divergence of a potential ligand-binding domain. OAK kinase activity is required for the aberrant phenotypes, indicating it is not recognition of an aberrant protein, but rather a true gain of function, or overdominance for gene activity, that leads to this underdominance for fitness. Our work provides insights into how tandem arrays, which are particularly prone to frequent, complex rearrangements, can produce genetic novelty.

While intraspecific hybrids are vitally important in modern agriculture because they often perform better than their inbred parents, certain hybrid combinations fail to develop normally and are inferior to their parents. We have identified an Arabidopsis thaliana hybrid with several aberrant growth phenotypes that are caused by divergence at a single locus encoding the receptor-like kinase OUTGROWTH-ASSOCIATED PROTEIN KINASE (OAK). OAK belongs to a group of similar genes arranged in a tandem cluster that varies substantially between A. thaliana strains. OAK originated through duplication within the cluster with concurrent recruitment of coding sequences from nearby genes to form a new promoter with a novel expression pattern. Kinase activity of OAK is required for its effects, indicating that it is not recognition of an aberrant protein but rather a true gain of function that leads to the incompatibility. Most of the incompatibility seems to come from divergence within the extracellular ligand-binding domain of the OAK protein, indicating that heterodimers of OAK may have higher affinity for a natural substrate compared to either homodimer. Finally, mis-expression of the incompatible OAK alleles from the promoter present in the reference strain of A. thaliana also leads to genetic incompatibility, but with different phenotypic outcomes.

Improved haplotype-based detection of ongoing selective sweeps towards an application in Arabidopsis thaliana

Background

The increasing amount of genome information allows us to address various questions regarding the molecular evolution and population genetics of different species. Such genome-wide data sets including thousands of individuals genotyped at hundreds of thousands of markers require time-efficient and powerful analysis methods. Demography and sampling introduce a bias into present population genetic tests of natural selection, which may confound results. Thus, a modification of test statistics is necessary to introduce time-efficient and unbiased analysis methods.

Results

We present an improved haplotype-based test of selective sweeps in samples of unequally related individuals. For this purpose, we modified existing tests by weighting the contribution of each individual based on its uniqueness in the entire sample. In contrast to previous tests, this modified test is feasible even for large genome-wide data sets of multiple individuals. We utilize coalescent simulations to estimate the sensitivity of such haplotype-based test statistics to complex demographic scenarios, such as population structure, population growth and bottlenecks. The analysis of empirical data from humans reveals different results compared to previous tests. Additionally, we show that our statistic is applicable to empirical data from Arabidopsis thaliana. Overall, the modified test leads to a slight but significant increase of power to detect selective sweeps among all demographic scenarios.

Conclusions

The concept of this modification might be applied to other statistics in population genetics to reduce the intrinsic bias of demography and sampling. Additionally, the combination of different test statistics may further improve the performance of tests for natural selection.

Evolutionary genetics of plant adaptation

Plants provide unique opportunities to study the mechanistic basis and evolutionary processes of adaptation to diverse environmental conditions. Complementary laboratory and field experiments are important for testing hypotheses reflecting long-term ecological and evolutionary history. For example, these approaches can infer whether local adaptation results from genetic tradeoffs (antagonistic pleiotropy), where native alleles are best adapted to local conditions, or if local adaptation is caused by conditional neutrality at many loci, where alleles show fitness differences in one environment, but not in a contrasting environment. Ecological genetics in natural populations of perennial or outcrossing plants can also differ substantially from model systems. In this review of the evolutionary genetics of plant adaptation, we emphasize the importance of field studies for understanding the evolutionary dynamics of model and nonmodel systems, highlight a key life history trait (flowering time) and discuss emerging conservation issues.

Reference-guided assembly of four diverse Arabidopsis thaliana genomes

We present whole-genome assemblies of four divergent Arabidopsis thaliana strains that complement the 125-Mb reference genome sequence released a decade ago. Using a newly developed reference-guided approach, we assembled large contigs from 9 to 42 Gb of Illumina short-read data from the Landsberg erecta (Ler-1), C24, Bur-0, and Kro-0 strains, which have been sequenced as part of the 1,001 Genomes Project for this species. Using alignments against the reference sequence, we first reduced the complexity of the de novo assembly and later integrated reads without similarity to the reference sequence. As an example, half of the noncentromeric C24 genome was covered by scaffolds that are longer than 260 kb, with a maximum of 2.2 Mb. Moreover, over 96% of the reference genome was covered by the reference-guided assembly, compared with only 87% with a complete de novo assembly. Comparisons with 2 Mb of dideoxy sequence reveal that the per-base error rate of the reference-guided assemblies was below 1 in 10,000. Our assemblies provide a detailed, genomewide picture of large-scale differences between A. thaliana individuals, most of which are difficult to access with alignment-consensus methods only. We demonstrate their practical relevance in studying the expression differences of polymorphic genes and show how the analysis of sRNA sequencing data can lead to erroneous conclusions if aligned against the reference genome alone. Genome assemblies, raw reads, and further information are accessible through http://1001genomes.org/projects/assemblies.html.

Journey through the past: 150 million years of plant genome evolution

The genome sequence of the plant model organism Arabidopsis thaliana was presented in December of the year 2000. Since then, the 125 Mb sequence has revealed many of its evolutionary secrets. Through comparative analyses with other plant genomes, we know that the genome of A. thaliana, or better that of its ancestors, has undergone at least three whole genome duplications during the last 120 or so million years. The first duplication seems to have occurred at the dawn of dicot evolution, while the later duplications probably occurred <70 million years ago (Ma). One of those younger genome-wide duplications might be linked to the K-T extinction. Following these duplication events, the ancestral A. thaliana genome was hugely rearranged and gene copies have been massively lost. During the last 10 million years of its evolution, almost half of its genome was lost due to hundreds of thousands of small deletions. Here, we reconstruct plant genome evolution from the early angiosperm ancestor to the current A. thaliana genome, covering about 150 million years of evolution characterized by gene and genome duplications, genome rearrangements and genome reduction.

Strategies for silencing and escape: the ancient struggle between transposable elements and their hosts

Over the past several years, there has been an explosion in our understanding of the mechanisms by which plant transposable elements (TEs) are epigenetically silenced and maintained in an inactive state over long periods of time. This highly efficient process results in vast numbers of inactive TEs; indeed, the majority of many plant genomes are composed of these quiescent elements. This observation has led to the rather static view that TEs represent an essentially inert portion of plant genomes. However, recent work has demonstrated that TE silencing is a highly dynamic process that often involves transcription of TEs at particular times and places during plant development. Plants appear to use transcripts from silenced TEs as an ongoing source of information concerning the mobile portion of the genome. In contrast to our understanding of silencing pathways, we know relatively little about the ways in which TEs evade silencing. However, vast differences in TE content between even closely related plant species suggest that they are often wildly successful at doing so. Here, we discuss TE activity in plants as the result of a constantly shifting balance between host strategies for TE silencing and TE strategies for escape and amplification.